This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Number of aligned reads in exon, intron ...

Hello,

I am new in RNA-seq. I would like to know how calculate the number (%) of aligned reads in exons, introns in a bam file.

Thanks

exon-coverage rna-seq
• 10,240 views • link
updated by Ram Ram • written by samuelrivero • samuelrivero

3 answers

Hi! I'm usually doing this with Picard suite, see http://picard.sourceforge.net/picard-metric-definitions.shtml
• 0 views • link
written by mikhail.shugay mikhail.shugay

You could try http://rseqc.sourceforge.net/#read-distribution-py also have a look at the other programs from this suite.

• 0 views • link
updated by Ram Ram • written by Ido Tamir Ido Tamir

Thank you so much for the help. I will check those suites. Do you think using samtools in this way will work?

samtools view -b -F exons.bed aligned.bam > aligned.exons.bam

Thank you again.

Sam

• 0 views • link
updated by Ram Ram • written by samuelrivero • samuelrivero

You don't need to explicitly use samtools for CollectRNASeqMetrics from Picard (http://picard.sourceforge.net/command-line-overview.shtml#CollectRnaSeqMetrics)

the command line should be following (after downloading and unpacking Picard tools):

java -Xmx4G -jar picard/CollectRnaSeqMetrics.jar REF_FLAT=exome.refFlat INPUT=aligned.bam CHART_OUTPUT=rnaSeqCoverage.pdf OUTPUT=rnaSeqMetrics.txt

the exome.refFlat file for H. Sapiens could be downloaded from here.

for other species check http://hgdownload.cse.ucsc.edu/downloads.html

PS. To reply to a specific answer please further use "comment" option

• 0 views • link
updated by Ram Ram • written by mikhail.shugay mikhail.shugay

Log in to answer this question.

More posts like this