Counting reads on genome
1
0
Entering edit mode
8.7 years ago
EVR ▴ 610

Hi,

I aligned my reads to de novo transcriptome using bowtie. Now the unaligned reads from the previous step was mapped to genome using tophat2. Now I would like to count the reads mapped to every scaffold and along the coverage region. How can I achieve this as tools like tophat,feature counts requires gtf file which has location of transcripts in every scaffold.But here the unaligned reads are mapped outside of the transcripts and hence I dont know whether using GTF file will be appropriate.

Any guidance wil be highly useful.

Thanks in advance

RNA-Seq genome transcriptome • 1.8k views
ADD COMMENT
0
Entering edit mode
8.7 years ago
mastal511 ★ 2.1k

Samtools idxstats should give you the number of reads mapped to every reference sequence in your genome file.

See the samtools manual:

http://www.htslib.org/doc/samtools.html

ADD COMMENT
0
Entering edit mode

Qualimap also provides a nice visual summary for BAM files.

ADD REPLY
0
Entering edit mode

@mastal511 and @genomax2 Thanks for your guidance. I estimated the reads like you said. Now how can I get the coverage information i.e. in scaffold_XX the region between 23445 and 23670 contains 200 reads. Can I make use of genome coverage? Kindly guide me

ADD REPLY
0
Entering edit mode

What question are you trying to answer after this step (or in general)?

ADD REPLY

Login before adding your answer.

Traffic: 1639 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6