Closed ghost closed 5 years ago
@drjsanger Has this been run over a BAM/CRAM file and confirmed to work as expected?
$ install/bin/mismatchQc -i /path/to/input.bam -p | samtools-1.3.1 view -H | grep -e '^@PG'
Results in:
...
@PG ID:PCAP-core-mismatchQC.1 PN:PCAP-core-mismatchQC PP:bammerge' CL:install/bin/mismatchQc -p -i /path/to/input.bam -o - DS:Marks a read as QCFAIL and adds aux tag 'mm' where the mismatch rate higher than the threshold VN:4.3.1
...
@keiranmraine Please comment