Retrieving Gene Names From 3'-End Aligned Reads
Hello, I am planing an experiment, where I will collect RNA-seq data on a 3'-polyA capture library. My question is, after alignment with Bowtie, what is the best way to retrieve the gene names of the loci to which the data aligned? These are short reads that map to the 3'-end of transcripts and that's it. Any idea, will be appreciated. G.
• 2,564 views
•
link
updated
by
Damian Kao
•
written
by
GPR
•
written
by
GPR
1 answer
You need a genome annotation file preferably in gtf or gff format for the genes you want to look at.
I used HTSeq-count a lot for counting reads that correspond to an annotation. http://www-huber.embl.de/users/anders/HTSeq/doc/count.html#count
• 1 views
•
link
written
by
Damian Kao
Log in to answer this question.
More posts like this
-
Mouse miRNA mapping
written by AellSidious •
-
Extract primer-containing reads from RNA-Seq
written by Ruixuan •
-
RNA-seq read alignment best practices for low-input Smart-seq2 samples
written by goldberg.jm
-
Strand Specific Tophat Output
written by GPR
-
Varscan2 Strandness Filter
written by GPR
-
Tophat2, Strand-Specific And Genes.Gtf
written by GPR
-
Bowtie Sam Output Not Parsed
written by GPR
-
Bowtie Output "Not Enough Fields"
written by GPR
-
Microrna Reads Alignment
written by GPR
-
Bwa And Splicing Junctions
written by GPR
Do you have a gene model file for the organism? Then after mapping, its basically the same procedure. You could just count the number of reads mapped to each gene, albeit it maps to the 3' region of your gene. If necessary, to filter out certain mapping errors, you could also choose a window size from the 3' region of your gene and obtain counts over that region. Those with counts are the ones you are interested in (typically with a certain threshold over the number of read counts per million, cpm).