Question: Mapping reads to nuclear, chloroplast, and mitochondria genome separately but all resulting bam files are similar sizes
0
gravatar for kristina.mahan
11 months ago by
kristina.mahan110 wrote:

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?

ADD COMMENTlink written 11 months ago by kristina.mahan110
3
gravatar for ATpoint
11 months ago by
ATpoint41k
Germany
ATpoint41k wrote:

Because a BAM contains all reads (mapped and unmapped) that were in the original fastq file, so this is normal and expected.

ADD COMMENTlink written 11 months ago by ATpoint41k

so for further process just remove unmapped reads from those bam files?

ADD REPLYlink written 11 months ago by kristina.mahan110

Technically you can do that to reduce file size. If you want that is on you.

ADD REPLYlink modified 11 months ago • written 11 months ago by ATpoint41k

This would be dependent on aligner options. Some aligners will allow you to drop unmapped reads.

ADD REPLYlink written 11 months ago by genomax92k
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: 1696 users visited in the last hour