Question: Counting reads aligned to forward or reverse strand
0
gravatar for caggtaagtat
2.0 years ago by
caggtaagtat1.1k
caggtaagtat1.1k wrote:

Hello,

I mapped reads (stranded library prep) with STAR to my reference genome and learned, that you can count reads, which aligned to the forward or reverse strand of the reference seuquence, using the FLAG option.

Reads mapped to Forward strand

samtools view -c -F 20 FILE.bam >> Output_file.txt

Reads mapped to Reverse strand

samtools view -c -f 16 FILE.bam >> Output_file.txt

This showed me in all files a strong inbalance of read distribution between the strands, which I did not expect, with a few hundred thousend reads on the plus strand and a few million reads on the reverese strand. The library preparation should not make any difference in this, that's why I am wondering, if I may have something wrong in the code?

There are no multimappers in the data and the data only comes from mitochondrial RNA.

Edit: a word

rna-seq • 1.7k views
ADD COMMENTlink modified 2.0 years ago by michael.ante3.6k • written 2.0 years ago by caggtaagtat1.1k
2
gravatar for michael.ante
2.0 years ago by
michael.ante3.6k
Austria/Vienna
michael.ante3.6k wrote:

Hi caggtaagtat,

I'd check the gene expression with htseq-count or featureCounts setting the right stranded parameter for your library prep. In human for instance, more than 75% of the mt-genes are located on the + strand. Having little expression on the minus-strand genes would explain your result.

Cheers,

Michael

ADD COMMENTlink written 2.0 years ago by michael.ante3.6k

Thank you,

I will try to use htseq-count and the RSeQC now. I was just wondering, if I can really access the stranded coverage with these Flags in samtools view or if I got something wrong.

ADD REPLYlink modified 2.0 years ago • written 2.0 years ago by caggtaagtat1.1k

Ok, so using the HTseq-count script, I got a much more believable result, with most reads on the heavy strands.

ADD REPLYlink written 2.0 years ago by caggtaagtat1.1k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 828 users visited in the last hour