Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

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

Open
timeresistance1996 opened this issue Mar 27, 2024 · 0 comments
Open

Comments

@timeresistance1996
Copy link

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?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant