You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Hi ! You're right, something's wrong here.. :-)
It is however relatively innocuous and an easy fix: by default, VirSorter extends the prophage sequence by 50 nucleotide beyond the last gene in 5' and 3' (to include potential att sites and not end a contig right on a start / stop codon). I just forgot to include a check to make sure we don't extend past the contig, i.e. if the last gene of the prophage is at the end of the contig, the coordinate from VirSorter will be 50 nucleotides beyond the contig (hence 63493 vs 63443).
I'll fix it asap, but in the meantime you can safely use these results anyway, as the prediction is overall correct (also, the genbank file you get from this prophage is accurate too, since it automatically adjust to the actual sequence length).
Hi! I just observed a possible issue where a predicted prophage sequence's stop is larger than the actual contig size:
Test assembly:
kleiner_2015.fasta.gz
Command used:
Observation
It seems that VirSorter predicts a prophage in a range that is actually larger than the contig size. Example:
>NODE_51_length_63443_cov_50.479870
So contig
NODE_51
has a length of 63443 nt.Now VirSorter predicts a prophage for this contig from position 19922-63493:
So the predicted prophage's stop position is larger than the actual contig size when I understand the output correctly?
The text was updated successfully, but these errors were encountered: