Repository navigation
mlml "error: chrom or pos mismatch in both files" #135
Description
While using dnmtools (v 1.2.4 by command line) mlml with mock bisulfite and oxidated bisulfite data. Some, but not all, samples yielded this error:
error: chrom or pos mismatch in both files
The input files were generated from dnmtools counts from bam files from Bismark alignment.
I found that the reference had some small contig "chromosomes" one of which apparently had at least 1 read for one library time but not the other, causing this mismatch in the chromosomes represented in the two .meth files.
The error was resolved when I filtered out rows in the .meth file from that contig with awk, but it would be better if mlml solved this case.
Activity
andrewdavidsmith commentedon Sep 12, 2023
@christianwake I think we can make this more convenient. My preference would be to work harder to ensure that the sites included in the counts files depend on the reference genome, and not the mapped reads files. And in addition have an option for
mlmlthat could-forceit to ignore a discrepancy if the user were confident the two files are from the same reference and have a consistent sorted order. I think we've made some progress in the former, but it's not yet in a release.andrewdavidsmith commentedon Sep 12, 2023
Relevant here: PR #118 "counts: reporting chromosomes without mapped reads" So using
dnmtools countsfrom the source repo since about 2 weeks ago would have hopefully avoided the issue.andrewdavidsmith commentedon Oct 5, 2023
@christianwake This is probably solved for most users in dnmtools v1.4.0 because it will ensure that the
.countsinput files have all chromosomes from the reference genome, even if no read maps to them. I'm not closing this yet because I'm considering updatingmlmlto make it robust to missing chromosomes, but it's tricky sincemlmlcan take 3 inputs, so distinguishing inconsistent sorting vs. missing chroms is more involved. In case you are using Conda, we plan to release 1.4.1 on Conda this week.