Hi,
I have used a reference genome for assembling. For the Draft assembly fasta file is to big to mapping to the reference genome. So for some fasta sequence I only choose the first 1,000,000bp for mapping. I have used the parameter "blastn -outfmt 7". And the names of the blast result file are in right format.
When I run lachesis, it shows the following err massege:
Lachesis: TextFileParsers.cc:488: void ParseBlastAlignmentFiles(const std::vector<std::basic_string >&, const std::vector&, const std::vector<std::basic_string >&, const string&): Assertion `stop_on_Q <= query_lengths[query_ID]' failed.
/opt/gridview//pbs/dispatcher/mom_priv/jobs/714734.admin.SC: line 1: 30684
I don't know why it shows this. And how can I fix it.
looking forward for help
Hi, I have used a reference genome for assembling. For the Draft assembly fasta file is to big to mapping to the reference genome. So for some fasta sequence I only choose the first 1,000,000bp for mapping. I have used the parameter "blastn -outfmt 7". And the names of the blast result file are in right format. When I run lachesis, it shows the following err massege: Lachesis: TextFileParsers.cc:488: void ParseBlastAlignmentFiles(const std::vector<std::basic_string >&, const std::vector&, const std::vector<std::basic_string >&, const string&): Assertion `stop_on_Q <= query_lengths[query_ID]' failed.
/opt/gridview//pbs/dispatcher/mom_priv/jobs/714734.admin.SC: line 1: 30684
I don't know why it shows this. And how can I fix it. looking forward for help