I'm using samtool/1.9.0 and I'm trying to merge 2 files, but what I observed is when I give the inputs in different order I get very different output.
I have f1.bam 3G and f2.bam 200M. When I do samtools merge output.bam f1.bam f2.bam, I get a that's slightly larger than 3G. When I do samtools merge output.bam f2.bam f1.bam, I get a that's slightly larger than 200M. I'm sure they are not just f1.bam or f2.bam, but I wonder what could have been wrong or is it a issue with samtools/1.9.0? I also observed that the total counts I get from the output bam files of samtools merge output.bam f1.bam f2.bam are 10 time larger than that from samtools merge output.bam f2.bam f1.bam This is the flastat output for samtools merge output.bam f1.bam f2.bam 164862366 + 0 in total (QC-passed reads + QC-failed reads) 0 + 0 secondary 0 + 0 supplementary 0 + 0 duplicates 164862366 + 0 mapped (100.00% : N/A) 0 + 0 paired in sequencing 0 + 0 read1 0 + 0 read2 0 + 0 properly paired (N/A : N/A) 0 + 0 with itself and mate mapped 0 + 0 singletons (N/A : N/A) 0 + 0 with mate mapped to a different chr 0 + 0 with mate mapped to a different chr (mapQ>=5) 164862366 (edited) [white_check_mark][eyes][raised_hands] 1:23<https://uchicagomoskolab.slack.com/archives/GQJPM6K7B/p1642015420010500> This is the flastat output for samtools merge output.bam f2.bam f1.bam 11009638 + 0 in total (QC-passed reads + QC-failed reads) 0 + 0 secondary 0 + 0 supplementary 0 + 0 duplicates 11009638 + 0 mapped (100.00% : N/A) 0 + 0 paired in sequencing 0 + 0 read1 0 + 0 read2 0 + 0 properly paired (N/A : N/A) 0 + 0 with itself and mate mapped 0 + 0 singletons (N/A : N/A) 0 + 0 with mate mapped to a different chr 0 + 0 with mate mapped to a different chr (mapQ>=5) 11009638
_______________________________________________ Samtools-help mailing list Samtools-help@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/samtools-help
