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]]

