Question: How can I count reads mapped to the ribosomal genes.
0
gravatar for biplab
20 months ago by
biplab100
University of California, Davis
biplab100 wrote:

There are two ribosomal genes(RDN37-1, RDN37-2) in my genomic sequence file, both have same sequence. Therefore when I did the alignment by bowtie2, reads that mapped to ribosomal genes mapped both regions thus low mapping quality for those reads. After mapping when I used htseq-count in order to count read mapped to features RDN37-1 and RDN37-2 give 0 count because low quality mapping. In this case, what will be best way to count reads for RDN37-1?

rna-seq alignment next-gen • 610 views
ADD COMMENTlink written 20 months ago by biplab100
2

I think the option for htseq-count to count multi-mapped reads is --nonunique all. Equivalent command for featureCounts is -M. Both programs will not count multi-mappers by default.

ADD REPLYlink written 20 months ago by genomax70k

Thank you so much for help.

ADD REPLYlink modified 19 months ago • written 19 months ago by biplab100
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: 1724 users visited in the last hour