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

Imitervirales assignment although not well supported by ViPhOGs #110

Open
hoelzer opened this issue Jul 14, 2023 · 0 comments
Open

Imitervirales assignment although not well supported by ViPhOGs #110

hoelzer opened this issue Jul 14, 2023 · 0 comments

Comments

@hoelzer
Copy link
Collaborator

hoelzer commented Jul 14, 2023

This issue was discovered during the introduction of the new voting scheme and adjustment of the virus taxonomy: #103

We ran the assembly.fasta test file. In the 08-final/gff/assembly_virify.gff file, we got for pos.phage.5 the following CDS and ViPhOG hits:

pos.phage.5	Prodigal	CDS	1934	3370	.	+	0	ID=pos.phage.5_3;gbkey=CDS;viphog=ViPhOG1981;viphog_taxonomy=Arquatrovirinae
pos.phage.5	Prodigal	CDS	13526	15382	.	+	0	ID=pos.phage.5_19;gbkey=CDS;viphog=ViPhOG9461;viphog_taxonomy=Twortvirinae
pos.phage.5	Prodigal	CDS	15397	18063	.	+	0	ID=pos.phage.5_20;gbkey=CDS;viphog=ViPhOG9903;viphog_taxonomy=Picovirinae
pos.phage.5	Prodigal	CDS	18063	19595	.	+	0	ID=pos.phage.5_21;gbkey=CDS;viphog=ViPhOG9462;viphog_taxonomy=Twortvirinae
pos.phage.5	Prodigal	CDS	22518	23501	.	+	0	ID=pos.phage.5_26;gbkey=CDS;viphog=ViPhOG27646;viphog_taxonomy=Picovirinae
pos.phage.5	Prodigal	CDS	37511	37756	.	+	0	ID=pos.phage.5_54;gbkey=CDS;viphog=ViPhOG9436;viphog_taxonomy=Klosneuvirinae
pos.phage.5	Prodigal	CDS	37765	38172	.	+	0	ID=pos.phage.5_55;gbkey=CDS;viphog=ViPhOG1970;viphog_taxonomy=Sepvirinae
pos.phage.5	Prodigal	CDS	40637	40816	.	+	0	ID=pos.phage.5_62;gbkey=CDS;viphog=ViPhOG4695;viphog_taxonomy=Twortvirinae

but in the final taxonomy assignment in 08-final/taxonomy/high_confidence_viral_contigs_prodigal_annotation_taxonomy.tsv we got

contig_ID	genus	subfamily	family	order	class
pos.phage.5			1	Imitervirales	Megaviricetes

Why do we get based on these ViPhOGs Imitervirales as annotation?

Actually, Klosneuvirinae is a subfamily within Imitervirales but it doesn't make sense that the pipeline classifies it with only one ViPhOG hit.

@guille0387 was looking into that and probably found the reason.

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