X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=bam_stat.c;h=a4522daa49257b920ab26586a6062d046a26d9e3;hb=02884d7b923597d15ccecd598e73675b5074a1eb;hp=142933666d687e261367664aaabcdbfcb4e3107f;hpb=091c8d2918d888c664dd8baaf9f2a84b2d2b6cd3;p=samtools.git diff --git a/bam_stat.c b/bam_stat.c index 1429336..a4522da 100644 --- a/bam_stat.c +++ b/bam_stat.c @@ -62,9 +62,9 @@ int bam_flagstat(int argc, char *argv[]) printf("%lld paired in sequencing\n", s->n_pair_all); printf("%lld read1\n", s->n_read1); printf("%lld read2\n", s->n_read2); - printf("%lld properly paired (%.2f%%)\n", s->n_pair_good, (float)s->n_pair_good / s->n_reads * 100.0); + printf("%lld properly paired (%.2f%%)\n", s->n_pair_good, (float)s->n_pair_good / s->n_pair_all * 100.0); printf("%lld with itself and mate mapped\n", s->n_pair_map); - printf("%lld singletons (%.2f%%)\n", s->n_sgltn, (float)s->n_sgltn / s->n_reads * 100.0); + printf("%lld singletons (%.2f%%)\n", s->n_sgltn, (float)s->n_sgltn / s->n_pair_all * 100.0); printf("%lld with mate mapped to a different chr\n", s->n_diffchr); printf("%lld with mate mapped to a different chr (mapQ>=5)\n", s->n_diffhigh); free(s);