cerebis / qc3C

Reference-free quality assessment for Hi-C sequencing data
GNU Affero General Public License v3.0
12 stars 1 forks source link

Provide greater detail when BAM sort order is incorrect. #60

Closed cerebis closed 2 years ago

cerebis commented 2 years ago

Hi Matthew,

thanks for your reply. Indeed I can confirm that my bam file does not contain the expected header:

$> samtools view -H L79851_Track-123454.aligned.sorted.bam| grep '^@HD'
@HD     VN:1.6  SO:unknown

The bam files were generated as follows:

samtools sort -n [aligned.bam] -o [aligned.sorted.bam]

However, the logs for that particular step are missing so I will regenerate just to be sure.

Originally posted by @agilly in https://github.com/cerebis/qc3C/issues/58#issuecomment-1034540210