GuoliangLi-HZAU / BatMeth2

BS-seq analysis pipeline
38 stars 12 forks source link

batmeth2 methgff #48

Open yhj-j opened 4 months ago

yhj-j commented 4 months ago

Chr01 332613000 . 22 610 Chr01.332613000.332613200 Chr01 332613201 . 22 610 Chr01.332613201.332613401 Chr01 332613401 . 22 610 Chr01.332613401.332613601 Chr01 332613601 . 22 610 Chr01.332613601.332613801 Chr01 332613801 . 22 610 Chr01.332613801.332614001 Chr01 332614001 . 22 610 Chr01.332614001.332614201 Chr01 332614201 . 22 610 Chr01.332614201.332614401 Chr01 332614401 . 22 610 Chr01.332614401.332614601 Chr01 332614601 . 22 610 Chr01.332614601.332614801 Chr01 332614801 . 22 610 Chr01.332614801.332615001 Chr01 332615001 . 22 610 Chr01.332615001.332615201 Chr01 332615201 . 22 610 Chr01.332615201.332615376 Chr02 1 . 22 610 Chr02.1.201 Chr02 201 . 22 610 Chr02.201.401 Chr02 401 . 22 610 Chr02.401.601 Chr02 601 . 22 610 Chr02.601.801 Chr02 801 . 22 610 Chr02.801.1001 Chr02 1001 . 22 610 Chr02.1001.1201 Chr02 1201 . 22 610 Chr02.1201.1401 Chr02 1401 . 22 610 Chr02.1401.1601 Chr02 1601 . 22 610 Chr02.1601.1801 Chr02 1801 . 22 610 Chr02.1801.2001 Chr02 2001 . 22 610 Chr02.2001.2201 Chr02 28000 . 47 47 Chr02.28000.28200 Chr02 57800 . 7 10 Chr02.57800.58000 Chr02 58000 . 2 10 Chr02.58000.58200 Chr02 59200 . 23 48 Chr02.59200.59400 Chr02 61800 . 12 24 Chr02.61800.62000 您好 老师,我在用batmeth2 methgff 计算区域甲基化水平时,对输出结果有点疑惑,软件输入的bed文件 0-based,输出结果中,不是很理解到底是1-based 还是0-based的,而且Chr01 最后的染色体数还超出了 我chr01界限1bp BatMeth2 methyGff --coverage 4 -nC 1 \ --genome ../../batmeth_idx/st_8_filter_Chr.genome.fa \ -b ../st_8_filter_Chr.genome.200bp_1based.bed3 \ --methratio ../../3.SiteAndSmallWindow/leaf_WGBS_rep1.methratio.txt \ --out leaf_WGBS_rep1.200bp.methy