-
Notifications
You must be signed in to change notification settings - Fork 5
Description
I use metaDMG with "./metaDMG-cpp getdamage -r 1" and "/metaDMG-cpp dfit .bdamage.gz --showfits 1".
.bdamage.gz.dfit.gz seems no error exit. But the reference sequence names has been renamed to "taxid"(the first column of dfit result).
May I ask the renaming rules so that I can map the sequence names back to .bdamage.gz.dfit.gz.
Header of .bdamage.gz.dfit.gz:
taxid A q c phi llh ncall sigmaD Zfit fwdx0 fwdxConf0 fwdx1 fwdxConf1 fwdx2 fwdxConf2 fwdx3 fwdxConf3 fwdx4 fwdxConf4 bwdx0 bwdxConf0 bwdx1 bwdxConf1 bwdx2 bwdxConf2 bwdx3 bwdxConf3 bwdx4 bwdxConf4
0 0.053594 0.538223 0.006153 52610.912454 -9.006416 3983.000000 0.039217 1.366616 0.059748 0.041272 0.030902 0.048002 0.017582 0.028685 0.011431 0.019745 0.008590 0.017445 0.059748 0.063354 0.030902 0.036091 0.017582 0.022222 0.011431 0.020852 0.008590 0.018461
1 0.019535 0.477313 0.000000 16.948712 -4.725200 817.000000 0.044401 0.439969 0.019535 0.044401 0.010211 0.032640 0.005337 0.021941 0.002790 0.015574 0.001458 0.011491 0.019535 0.044401 0.010211 0.027679 0.005337 0.021790 0.002790 0.014911 0.001458 0.011576
Header of my .sam:
@hd VN:1.0 SO:unsorted
@sq SN:k127_375571 LN:1077
@sq SN:k127_150230 LN:2244
By the way, which one is the recomand damage valu: A or q or any other value?