I am new to these types of quality control checks. I have not faced any error in my analyses but want to perform quality control. I have aligned my sncRNA reads with bowtie (-m 1 option) and after making the bam output from the alignment I performed HTSeq count. Now I want to see if most of the reads actually end up being represented in the HTSeq file as gene counts. How can I do that?
Should I just simply sum up the total of the counts reported by HTSeq and divide it by the total number of unique reads?
Any help will be appreciated, thanks
1 answer
https://htseq.readthedocs.io/en/master/count.html#usage at the bottom of the file you have some lines with some special counters you have to exclude for calculating the ratio. All of it has to sum up to the number of input reads. http://bioinf.wehi.edu.au/featureCounts/ is faster and gives you the same output, but the special counters are in a separate file.
Log in to answer this question.
-
Confirming intron length
written by Apex92 •
-
Bowtie error when indexing a reference file
written by Apex92 •
-
How to divide reads mapping to two different genes using featureCounts
written by Apex92 •
-
Samtools flagstat with sncRNA (single-end reads)
written by Apex92 •
-
Normalization of read counts from HTSeq by library size
written by Susmita Mandal
-
Conversion read counts to TPM for microRNAs
written by lessismore
-
fastQC for Oxford Nanopore reads
written by efratk •
-
What is the definition of "read depth" vs "coverage"? (again...)
written by ariel.balter
-
high aligned_not_unique number for HTSeq
written by bxia
-
How are expected error rates calculated?
written by Andrew •
