Strange results appeared when extracting modification information from DRS data

[rain-zjg]

Hello,
Recently, I've been analyzing RNA modification data from DRS (Direct RNA Sequencing) and would like to extract information for a specific type of RNA modification individually. I used the following code to extract modification information of inosine , but the results do not include inosine modifications.

${modkit} pileup --ref ${ref} --preload-references --modified-bases A:17596 --log-filepath ${output}/modkit_${name}.log --header --threads ${thread} ${workdir}/${bam} ${output}/test.pileup.bed

Here is the output:

And here is the information of modification tags:

Then I tried extracting all available modification types, and the results were normal.

Why is this happening? Looking forward to the developer's response. Thank you very much!

[ArtRand]

Hello @rain-zjg,

Thanks for noticing this, there is a bug in the routine that picks which sub-routine to use. Another user had noticed a similar bug. I've attached a build with a fix. There is a hidden flag --use-dynamic which will also produce the correct output with v0.6.0. This change will make it into v0.6.1 shortly.

modkit_v0.6.1_4d52624.zip