Entering edit mode
8.2 years ago
irritable_phd_syndrome
▴
130
I running a local blastx (v. 2.2.31) server with the most recent version of nr. I blastx some sequences and one of the entries returned is :
viral_c1 6e-16 48.91 79.5 85708;86385 7021350 AAF35296 Viruses Porcine circovirus 2;Porcine circovirus type 2-D replicase [Porcine circovirus-2]
where my query was formatted, i.e. " ... -outfmt '6 qseqid evalue pident bitscore staxids sgi sacc sskingdoms sscinames stitle'.
Why does my the 5th column, the taxids, have _two_ entries? I have a downstream analysis where I parse out the taxonomic IDs. Short of going through by hand and resolving each instance of these, I am not sure how I can programmatically choose / know the "correct" taxonomic ID. Could someone please explain why this is happening and the best way to resolve this?
