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

Evaluación

Línea de trabajo

Descripción

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.

Preparar el entorno

Primeros pasos

1. Lee el issue completo y luego la guía de contribución del proyecto.
2. Comenta en el issue que vas a ocuparte — evita que dos personas hagan lo mismo.
3. Haz un fork del repositorio y trabaja en una rama.
4. Abre un pull request que haga referencia al número del issue.

Más de nanoporetech/modkit

• Test failure due to concurrency issue

  Abierto

  ⋯Informar de un problema con este issue

  • ¿Algo no está bien?

  • Ya está asignado
  • No es para principiantes
  • Obsoleto o no válido
  • Otra cosa…
  • Enviar Cancelar
  bug

  Dificultad 2/5 1-3 horas Aptitud para principiantes 68/100

  nanoporetech/modkit#520 · 2 comentarios · Abierto el 14/10/2025

• Samtools code for modkit repair

  Abierto

  ⋯Informar de un problema con este issue

  • ¿Algo no está bien?

  • Ya está asignado
  • No es para principiantes
  • Obsoleto o no válido
  • Otra cosa…
  • Enviar Cancelar
  documentation

  Dificultad 1/5 Menos de una hora Aptitud para principiantes 65/100

  nanoporetech/modkit#336 · 1 comentario · Abierto el 9/1/2025

• modkit probability threshold

  Abierto

  ⋯Informar de un problema con este issue

  • ¿Algo no está bien?

  • Ya está asignado
  • No es para principiantes
  • Obsoleto o no válido
  • Otra cosa…
  • Enviar Cancelar

  Dificultad 5/5 Más de una semana Aptitud para principiantes 25/100

  nanoporetech/modkit#723 · Abierto el 18/9/2026

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

  Abierto

  ⋯Informar de un problema con este issue

  • ¿Algo no está bien?

  • Ya está asignado
  • No es para principiantes
  • Obsoleto o no válido
  • Otra cosa…
  • Enviar Cancelar

  Dificultad 3/5 1-2 días Aptitud para principiantes 66/100

  nanoporetech/modkit#721 · 1 comentario · Abierto el 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

  Abierto

  ⋯Informar de un problema con este issue

  • ¿Algo no está bien?

  • Ya está asignado
  • No es para principiantes
  • Obsoleto o no válido
  • Otra cosa…
  • Enviar Cancelar

  Dificultad 4/5 3-5 días Aptitud para principiantes 55/100

  nanoporetech/modkit#719 · Abierto el 25/8/2026

Todos los issues de nanoporetech/modkit

Issues similares

• A report-only `scan -g` tells you to run `socket-patch scan --mode agent [PATHS]` without `-g`, so following the hint scans the cwd project instead of the global install

  Abierto

  ⋯Informar de un problema con este issue

  • ¿Algo no está bien?

  • Ya está asignado
  • No es para principiantes
  • Obsoleto o no válido
  • Otra cosa…
  • Enviar Cancelar
  agent:triaged bug bughunt pm:npm priority:p1

  Dificultad 2/5 1-3 horas Aptitud para principiantes 82/100

  SocketDev/socket-patch#464 · 1 comentario · Abierto el 1/10/2026

  Los mantenedores suelen responder en 1 día

• MD005 compares byte columns, so a tab-indented sibling is reported

  Abierto

  ⋯Informar de un problema con este issue

  • ¿Algo no está bien?

  • Ya está asignado
  • No es para principiantes
  • Obsoleto o no válido
  • Otra cosa…
  • Enviar Cancelar
  bug

  Dificultad 2/5 1-3 horas Aptitud para principiantes 76/100

  akiomik/mado#482 · Abierto el 1/10/2026

  Los mantenedores suelen responder en 1 día

• Incorrect operator in `Deno.lint.BinaryExpression` declaration

  Abierto

  ⋯Informar de un problema con este issue

  • ¿Algo no está bien?

  • Ya está asignado
  • No es para principiantes
  • Obsoleto o no válido
  • Otra cosa…
  • Enviar Cancelar

  Dificultad 1/5 Menos de una hora Aptitud para principiantes 90/100

  denoland/deno#36934 · Abierto el 1/10/2026

  Los mantenedores suelen responder en 3 días

• `tx edit` fails on transaction XDR that is wrapped across multiple lines

  Abierto

  ⋯Informar de un problema con este issue

  • ¿Algo no está bien?

  • Ya está asignado
  • No es para principiantes
  • Obsoleto o no válido
  • Otra cosa…
  • Enviar Cancelar
  bug

  Dificultad 2/5 1-3 horas Aptitud para principiantes 78/100

  stellar/stellar-cli#2773 · Abierto el 1/10/2026

  Los mantenedores suelen responder en 2 días

• RVF SDK: deleteByFilter omits required native filter fields

  Abierto

  ⋯Informar de un problema con este issue

  • ¿Algo no está bien?

  • Ya está asignado
  • No es para principiantes
  • Obsoleto o no válido
  • Otra cosa…
  • Enviar Cancelar

  Dificultad 2/5 1-3 horas Aptitud para principiantes 78/100

  ruvnet/RuVector#1105 · Abierto el 1/10/2026

  Los mantenedores suelen responder en 1 día

Más issues de Rust