Number of aligned reads in exon, intron ...
3
3
Entering edit mode
8.0 years ago
samuelrivero ▴ 50

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 • 7.2k views
ADD COMMENT
0
Entering edit mode

Welcome to Biostars!

Please check previous questions: Extracting Read Count For Each Gene/Exon From Rna-Seq Bam Files

ADD REPLY
2
Entering edit mode
8.0 years ago
Hi! I'm usually doing this with Picard suite, see http://picard.sourceforge.net/picard-metric-definitions.shtml
ADD COMMENT
1
Entering edit mode
8.0 years ago
Ido Tamir 5.2k

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

ADD COMMENT
1
Entering edit mode
8.0 years ago
samuelrivero ▴ 50

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

ADD COMMENT
1
Entering edit mode

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

ADD REPLY

Login before adding your answer.

Traffic: 985 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