Entering edit mode
mauede@alice.it
▴
870
@mauedealiceit-3511
Last seen 9.7 years ago
I used function gene.seq to get the 3utr sequence relative to a vector
containing ""ensembl_transcript_id" items
> gene.seq <- getSequence (id=gene.map[,"ensembl_transcript_id"],
type="ensembl_transcript_id",
seqType="3utr", mart=hmart)
> length(gene.map[,"ensembl_transcript_id"])
[1] 246
Surprisingly, I get fewer 3utr sequences than the input number of
ensembl_transcript_id.
> length(gene.seq[,"3utr"])
[1] 216
How come ?
Thank you in advance,
Maura
tutti i telefonini TIM!
[[alternative HTML version deleted]]
