Entering edit mode
5.0 years ago
kristina.mahan
▴
170
I mapped reads onto a nuclear, chloroplast, and mitochondrial genome separately using bwa mem and then sorted using samtools. Entire genome: 179 Mb, chloroplast: 169Kb, mitochondria: 43Kb. The resulting bam files are all similar in sizes (23G, 27G, and 27G). Is their an error? Why would this be the case?
so for further process just remove unmapped reads from those bam files?
Technically you can do that to reduce file size. If you want that is on you.
This would be dependent on aligner options. Some aligners will allow you to drop unmapped reads.