ERROR: ref genome sequence seems to aligned! remove and restart
Opened this issue · 3 comments
Hello,
When running parsnp v1.2 with the following command:
parsnp -r 1339277.3.fna -d genomes -p 1 -o out
I get the following strange error:
-->Reading Genbank file(s) for reference (.gbk) ..
|->[WARNING]: no genbank file provided for reference annotations, skipping..
ERROR: ref genome sequence 1339277.3.fna seems to aligned! remove and restart
However, if I specify a different genome with -r
, the program runs to completion. All genomes are in the genomes
directory. All are draft genomes from isolates from the same species.
I'm trying to incorporate parsnp into a pipeline, but this unpredictable error is problematic. Any help would be appreciated.
Thanks,
Stephen
Looking over the other issues more closely, I noticed this is a duplicate of the open issue #19
A solution to this would be much appreciated
Several programs have problems with the contig (I do not know if it happens with parsnp) but most of the time it is better to form pseudo-chromosomes. Also, open the .fna file and check it (looking for symbols or something wrong) or simplify it in the tag name
Has it been a solution or going-around to this issue?