Error in Freebayes

Hi

I have been attempting to identify polymorphic DNA sites according to this workflow: https://usegalaxy.org/u/coursera/p/genomic-data-science-with-galaxyidentify-polymorphic-sites.

I am trying to use Freebayes to identify polymorphic sites from the SAM file that I cleaned up. However, I keep getting an “unable to find FASTA index” error irrespective of the version I use. Any feedback would be helpful. Thanks!

1 Like

Hi and welcome,

one suggestion first: it is way more helpful if, for reporting errors, you click on that little bug icon on the failed datset and copy/paste the complete error message you find there, instead of posting screenshots.

Now for your actual issue: the error is not about a missing FASTA index, but about specific entries in there. In other words, one or more of the chromosome names listed in your SAM/BAM inputs could not be found in the reference genome you’ve been trying to call variants against.
For your dataset 31 that makes lots of sense since Group1 is most definitely not a chromosome identifier anywhere. For dataset 32, the name is cut off in the screenshot, but it is up to you anyway to check whether this is a valid chromosome name in your analysis.

Best wishes,
Wolfgang

1 Like

That makes sense. Thanks!