Skip to content

Empty results of mfg #77

@zhangaicen

Description

@zhangaicen

Dear Guo,
Thanks for developing this great tool! Recently I am studying every detail function of CGmaptools, but when I used the mfg parameter (cgmaptools mfg -i CGmap.gz -r testbin.bed -c 2 -x CG > test.mfg),I got a blank file while the log file reported that it was successfully completed. the CGmap.gz file:
LG1 G 6782 CHH CA 0.33 1 3
LG1 G 6784 CHH CT 0.33 1 3
LG1 G 6806 CHH CT 0.67 2 3
LG1 G 6807 CHH CC 0.00 0 3
LG1 G 6815 CHH CT 0.67 2 3
LG1 G 6819 CHH CA 0.67 2 3
LG1 G 6824 CHH CT 0.33 1 3
LG1 G 6826 CHH CA 0.33 1 3
LG1 G 6843 CHH CT 0.00 0 5
LG1 C 6845 CHH CA 0.33 1 3
LG1 G 6853 CHH CT 0.00 0 4
LG1 G 6856 CHH CT 0.60 3 5
LG1 C 6861 CHH CC 0.00 0 4
LG1 C 6862 CHH CT 0.25 1 4
LG1 G 6867 CHH CT 0.40 2 5
LG1 C 6868 CG CG 1.00 4 4
LG1 G 6869 CG CG 1.00 5 5
LG1 G 6873 CHH CT 0.57 4 7
LG1 G 6874 CHH CC 0.38 3 8
LG1 C 6875 CHH CA 0.75 3 4
and the testbin.bed:
#chr chr U1 U2 U3 U4 U5 R1 R2 R3 R4 R5 D1 D2 D3 D4 End
LG1 - 23681 23631 23581 23531 23481 23431 23144 22857 22570 22283 21996 21946 21896 21846 21796
LG1 - 24569 24519 24469 24419 24369 24319 24128 23938 23748 23558 23368 23318 23268 23218 23168
LG1 + 31325 31375 31425 31475 31525 31575 31876 32177 32478 32779 33081 33131 33181 33231 33281
LG1 - 45248 45198 45148 45098 45048 44998 44696 44395 44094 43793 43492 43442 43392 43342 43292
LG1 - 47405 47355 47305 47255 47205 47155 46940 46726 46512 46298 46084 46034 45984 45934 45884
LG1 - 50415 50365 50315 50265 50215 50165 49837 49510 49183 48856 48529 48479 48429 48379 48329
LG1 - 53401 53351 53301 53251 53201 53151 52635 52120 51605 51090 50575 50525 50475 50425 50375
LG1 + 55274 55324 55374 55424 55474 55524 55825 56126 56427 56728 57030 57080 57130 57180 57230
LG1 + 57084 57134 57184 57234 57284 57334 58157 58981 59804 60628 61452 61502 61552 61602 61652
LG1 + 62957 63007 63057 63107 63157 63207 63502 63798 64093 64389 64685 64735 64785 64835 64885
LG1 - 67050 67000 66950 66900 66850 66800 66389 65978 65568 65157 64747 64697 64647 64597 64547
LG1 + 68688 68738 68788 68838 68888 68938 69554 70170 70786 71402 72019 72069 72119 72169 72219

I checked the bed file, and there were no negative values.
Could you please offer me some suggestions for solutions? Hope your response, thank you very much.

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions