Skip to content

The rule for renaming reference sequence names to "taxid".  #23

@timeresistance1996

Description

@timeresistance1996
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?

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions