Skip to content

Invalid data for modkit dmr pair #604

@caballero

Description

@caballero

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

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions