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

Plasmids classified as viruses #6

Open
asl opened this issue Nov 16, 2020 · 4 comments
Open

Plasmids classified as viruses #6

asl opened this issue Nov 16, 2020 · 4 comments
Assignees

Comments

@asl
Copy link
Member

asl commented Nov 16, 2020

Consider the attached FASTA file. It contains two sequences that match a plasmid from RefSeq at 99% IDY. However, there are classified as virus by viralVerify. Looks strange.
ic5-metaplasmidspades.plasmid.fasta.zip

@mikeraiko
Copy link
Collaborator

Interesting. They both predicted as viruses by other tools, such as VirSorter, and match phages P1 and P7 with IDY 96-98% and span 86%.
And this is fine because these phages are somewhat special - they can exist as a plasmids in the cell, and have both plasmid- and phage-specific genetic features.

@asl
Copy link
Member Author

asl commented Nov 17, 2020

FWIW Platon (https://github.com/oschwengers/platon) predicts it as plasmid.

@Dmitry-Antipov
Copy link
Contributor

I'm pretty sure that plasmidVerify would also classify it as plasmid.
For the sequences that similar to both P and V and not similar to C a plasmid prediction tool would say they are plasmid, and a viral prediction tool would say they are viral

Current logic for viralVerify is to compare V vs P+C and then only if P+C compare P vs C. It is vulnerable to such cases, but I suppose it is not a common situation..

@Dmitry-Antipov
Copy link
Contributor

Dmitry-Antipov commented Nov 17, 2020

And these contigs are plasmids and phages in the same time.

So, the right answer for such contigs should be virus_or_plasmid that is not currently supported by viralverify.

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

3 participants