I have a question regarding the read count. After the sequencing, i found that i have 513k reads that were passed. But when i align (minimap2) and check how much reads i have in my bam file using Samtools, i found myself with 1.7M reads. After that i was looking for the mapped reads only and i have 280k. Than when i use featureCount to quantify and count my reads, it tells me that i have 1.7M reads, but assign only 268k reads. Can you help me please ?
Another thing : I was also wondering how to quantify rRNAs in my sequencing using featureCount. I tried to retrieve all rRNAs feature in the count file generated by featureCount, but the number of reads aligning to the rRNAs is greater than the reads assigned by featureCount.
Thank you very much for your help!!