Open cgjosephlee opened 5 years ago
Hi,
The output files from call-methylation are zero-based, not 1-based. This is for compatibility with bisulfite sequencing data, which is usually in the 0-based bed format. The example bisulfite data we provide in the tutorial (bisulfite.ENCFF835NTC.example.tsv) should align with our calls.
Jared
Thanks for the reply.
The first base should be [0,1]
in zero-based format, right? (ref)
As the second example, the same sites are
bisulfite.ENCFF835NTC.example.tsv
chr20 5005939 5005940 . 12 + 5005939 5005940 255,0,0 12 100
methylation_frequency.tsv
chr20 5005939 5005939 1 5 5 1.000 GAGATCGCACC
which is confusing.
Indeed, the end coordinates are encoded slightly differently. I will consider changing our end coordinate to match the BED output.
Thank you, this will be really helpful.
I ran the example in https://nanopolish.readthedocs.io/en/latest/quickstart_call_methylation.html. However, I noticed that some reported coordinates were not pointed to C/G base and probably shift -1 base.
For example, in
methylation_frequency.tsv
but the mC in at 4988400bp
Another example,
I roughly counted that there were ~16k sites had this sort of behavior. Did I interpret wrong?