-
Notifications
You must be signed in to change notification settings - Fork 27
Invalid data for modkit dmr pair #604
Description
Dear devs,
Thanks for the great tool. I am doing some DMR analysis, but I found an error when modkit loads the data as it is comparing the valid_coverage only with count_methylated + count_cannonical, in the bedMehtyl definition, it should be Nvalid_cov = Nmod + Nother_mod + Ncanonical:
For example, this is my error:
batch failed: invalid data, valid coverage (8) is not equal to the sum of canonical and modified counts (6), [BedMethylLine { chrom: "1", interval: Interval { start: 3259525, stop: 3259526, val: () }, raw_mod_code: Code('m'), strand: Positive, count_methylated: 4, valid_coverage: 8, count_canonical: 2, count_other: 2, count_delete: 0, count_fail: 4, count_diff: 0, count_nocall: 0 }] chrom: 1 starting at 3259525, stopping
The line in the bedMethyl is:
1 3259525 3259526 m 8 + 3259525 3259526 255,0,0 8 50.00 4 2 2 04 0 0
valid_coverage = 8 (col 5 and 10)
count_methylated = 4 (col 12)
count_cannonical = 2 (col 13)
count_other_mod = 2 (col 14)
I am using modkit_v0.6.1_481e3c9