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
dies with Invalid reference index -1 #46
Comments
Hi,
I need to get this genome polished as soon as possible, is there a way to fix my issue? |
Hi Dario, The last log output line "Processing 000237F-010-01:1-12490" and the stack trace indicate it was trying to fetch reads mapping to that contig (000237F-010-01). Does that contig exist in the BAM header (@sq record)? If not, then that's the problem, though certainly Pilon should be more graceful about it (e.g., ignore that input contig). |
Indeed, it was missing, together with 130 (!) other contigs. I took them out and had a look: they are still related to my species. So I started aligning the reads with bwa mem (default settings) instead of bowtie2: do you think it will help? |
Hello,
I am having issues running Pilon on a genome where I could finish a run previously.
The command is
java -Xmx158G -jar /opt/pilon-1.18.jar --genome ../path/pdac_fu_124_pa_ctg.fa --frags aln_for_pilon_s.bam --output ../path/pdac_fu_124_pa_ctg_pilon.fa --diploid --threads 12 --changes
the reads were aligned with bowtie2, the error looks like:
I can't understand how to fix it.
thanks!
The text was updated successfully, but these errors were encountered: