Hi folks, Please correct me if I am wrong. I am dealing with how to get the cordinates of different genome regions such as UTR/intergenic/intragenic etc... and from the genePred format (http://genome.ucsc.edu/FAQ/FAQformat.html#format9), I think I can get them like follow:
If Strand = '+': 3UTR = txEnd-cdsEnd 5UTR = cdsStart-txStart Intragenic(i) = exonEnds(i)-exonStarts(i) Intergenic = all regions that do not overlap with gene cordinates (between txStart and txEnd) For Strand = '-', everything should be reversed, such as 3UTR = cdStart-txStart etc... Thank you very much in advance, D. _______________________________________________ Genome maillist - [email protected] https://lists.soe.ucsc.edu/mailman/listinfo/genome
