I want to map reads to AMGs sequences to get AMGs abundance information. But I do not know how to get AMGs sequences file. I try to use identified protein ID to separate it from prodigal.ffn file. But some protein ID have "fragment" . Would you have some ways to solve this questions.
Rather than pulling from the prodigal.ffn file you can pull from the phages_combined.ffn file. The names should match up with "fragment" prophages and whole phages.
I want to map reads to AMGs sequences to get AMGs abundance information. But I do not know how to get AMGs sequences file. I try to use identified protein ID to separate it from prodigal.ffn file. But some protein ID have "fragment" . Would you have some ways to solve this questions.