Discrepancy in calls when running individual modified base pileup and combined modbase pileup

Valutazione

Direzione di ricerca

Descrizione

Dear @ArtRand,

I ran modkit pileup on the modbam files for analyzing 5mC and 5hmC individually as well as using the --combine-mods commands. The total global methylation reported is:

I calculated this % using: % = (Total modified reads) / (Total modified + Total canonical reads) × 100

It is odd that the methylation % for 5mC and 5hmC are higher than those reported in the modkit output using the --combine-mods option. The complete command I ran for combined modification calls was:
$ modkitv0.6.1/dist_modkit_v0.6.1_481e3c9/modkit pileup --cpg --combine-strands --sampling-frac 1.0 --ref /PlasmoDB64_Pfalciparum3D7_Genome.fasta -t 8 /P55_sorted.bam --combine-mods ./P55_combinemh_CpG.bed --log ./P55.log --modified-bases 5mC 5hmC

The initial run generating individual 5mC and 5hmC values were run with default parameters (--filter-percentile 10).

If the case is such that the combined run used a stricter probability threshold, and so, reads that were weakly called as "m" or "h" individually got rejected when competing against each other... then how do we explain this phenomenon biologically?
Ideally, the combined run should give a methylation % either similar to or higher than that of the individual runs gave.

I request your assistance on the same.

Thank you.

Preparare l'ambiente

Come iniziare

1. Leggi tutta la issue e poi la guida ai contributi del progetto.
2. Commenta sulla issue per dire che te ne occupi tu — evita che due persone facciano lo stesso lavoro.
3. Fai un fork del repository e lavora su un branch.
4. Apri una pull request che faccia riferimento al numero della issue.

Altre issue di nanoporetech/modkit

• Test failure due to concurrency issue

  Aperta

  ⋯Segnala un problema con questa issue

  • Qualcosa non va?

  • Già presa in carico
  • Non adatta ai principianti
  • Obsoleta o non valida
  • Altro…
  • Invia Annulla
  bug

  Difficoltà 2/5 1-3 ore Idoneità per principianti 68/100

  nanoporetech/modkit#520 · 2 commenti · Aperta il 14/10/2025

• Samtools code for modkit repair

  Aperta

  ⋯Segnala un problema con questa issue

  • Qualcosa non va?

  • Già presa in carico
  • Non adatta ai principianti
  • Obsoleta o non valida
  • Altro…
  • Invia Annulla
  documentation

  Difficoltà 1/5 Meno di un'ora Idoneità per principianti 65/100

  nanoporetech/modkit#336 · 1 commento · Aperta il 9/1/2025

• modkit probability threshold

  Aperta

  ⋯Segnala un problema con questa issue

  • Qualcosa non va?

  • Già presa in carico
  • Non adatta ai principianti
  • Obsoleta o non valida
  • Altro…
  • Invia Annulla

  Difficoltà 5/5 Più di una settimana Idoneità per principianti 25/100

  nanoporetech/modkit#723 · Aperta il 18/9/2026

• `dmr pair` region output has `effect_size` sign inverted relative to the documentation and to single-site output

  Aperta

  ⋯Segnala un problema con questa issue

  • Qualcosa non va?

  • Già presa in carico
  • Non adatta ai principianti
  • Obsoleta o non valida
  • Altro…
  • Invia Annulla

  Difficoltà 3/5 1-2 giorni Idoneità per principianti 66/100

  nanoporetech/modkit#721 · 1 commento · Aperta il 9/9/2026

• pileup: worker panic (base_mods_adapter.rs:80 "entered unreachable code") is silently swallowed in whole-genome mode — presents as an infinite zero-CPU hang

  Aperta

  ⋯Segnala un problema con questa issue

  • Qualcosa non va?

  • Già presa in carico
  • Non adatta ai principianti
  • Obsoleta o non valida
  • Altro…
  • Invia Annulla

  Difficoltà 4/5 3-5 giorni Idoneità per principianti 55/100

  nanoporetech/modkit#719 · Aperta il 25/8/2026

Tutte le issue di nanoporetech/modkit

Issue simili

• Investigate remaining allow(dead_code) in thermal_model_iterative.rs

  Aperta

  ⋯Segnala un problema con questa issue

  • Qualcosa non va?

  • Già presa in carico
  • Non adatta ai principianti
  • Obsoleta o non valida
  • Altro…
  • Invia Annulla

  Difficoltà 2/5 1-3 ore Idoneità per principianti 72/100

  anchapin/fluxion#4107 · Aperta il 27/9/2026

  I maintainer di solito rispondono entro 1 giorno

• cranelift/README.md references removed test-all.sh

  Aperta

  ⋯Segnala un problema con questa issue

  • Qualcosa non va?

  • Già presa in carico
  • Non adatta ai principianti
  • Obsoleta o non valida
  • Altro…
  • Invia Annulla

  Difficoltà 1/5 Meno di un'ora Idoneità per principianti 92/100

  bytecodealliance/wasmtime#14433 · Aperta il 27/9/2026

  I maintainer di solito rispondono entro 1 giorno

• Specify Chinese fonts for "用户指南" in https://just.systems/

  Aperta

  ⋯Segnala un problema con questa issue

  • Qualcosa non va?

  • Già presa in carico
  • Non adatta ai principianti
  • Obsoleta o non valida
  • Altro…
  • Invia Annulla

  Difficoltà 2/5 1-3 ore Idoneità per principianti 68/100

  casey/just#3721 · Aperta il 27/9/2026

• Formatting a file that only contains comments deletes the contents of the file

  Aperta

  ⋯Segnala un problema con questa issue

  • Qualcosa non va?

  • Già presa in carico
  • Non adatta ai principianti
  • Obsoleta o non valida
  • Altro…
  • Invia Annulla
  bug

  Difficoltà 2/5 1-3 ore Idoneità per principianti 76/100

  nushell/nufmt#231 · 1 commento · Aperta il 27/9/2026

• bug(python): all-score rerankers lose the input Arrow schema

  Aperta

  ⋯Segnala un problema con questa issue

  • Qualcosa non va?

  • Già presa in carico
  • Non adatta ai principianti
  • Obsoleta o non valida
  • Altro…
  • Invia Annulla

  Difficoltà 2/5 1-3 ore Idoneità per principianti 88/100

  lancedb/lancedb#4331 · Aperta il 27/9/2026

  I maintainer di solito rispondono entro 1 giorno

Altre issue su Rust