How can I count reads mapped to the ribosomal genes.
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?
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Differences between number of expressed genes obtained with RSEM and htseq-count
written by concetta •
-
Counting reads over genes with summarizeOverlaps
written by concetta •
-
very-sensitive-local bowtie2 or bwa_mem
written by Illusion •
-
How to calculate Mapping Quality in sequencing data?
written by --panda-- •
-
rRNA quantification: which is the best way?
written by A. Domingues
-
Bowtie2: missing forward reads after alignment
written by nikkihathi •
-
All reads with mapq value <= 1, htseq does not attribute counts to any feature
written by fhsantanna
-
count reads in ambiguous regions for both features
written by caner •
-
What is the acceptable % of reads that map to no feature (in RNASeq)?
written by noushin.farnoud
-
rRNA Removal In Rna-Seq Data
written by jockbanan
I think the option for
htseq-countto count multi-mapped reads is--nonunique all. Equivalent command for featureCounts is-M. Both programs will not count multi-mappers by default.Thank you so much for help.