A question about the output file

PengNi / ccsmeth

Detecting DNA methylation from PacBio CCS reads

BSD 3-Clause Clear License

72 stars 11 forks source link

A question about the output file #23

Closed hunkui1 closed 2 years ago

hunkui1 commented 2 years ago

hi，I used this command： ccsmeth call_freqb --input_bam /home/lyx/wangaokai/data/pacbio_analyse/tden-mas/03.call_mods/r64114_20201106_074159.hifi.pbmm2.call_mods.modbam.bam --ref /home/lyx/wangaokai/data/genome/tden/Tden.genome.fasta --output ./r64114_20201106_074159.hifi.pbmm2.call_mods.modbam.feaq --threads 20 --sort I got： a1f6b147fc12b0996429a47f9c2b7f9 21389c1987c446652b1add0198f19df This is different from the reference format in README. How do I understand it？

PengNi commented 2 years ago

@hunkui1 , the output has the first 10 columns (chrom to modification_frequency), but with prob_0_sum/prob_1_sum replaced with "-".

hunkui1 commented 2 years ago

I see. Thank you for your answer.