This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Samtools faidx can't be run in succession on two reference files?

I am running a snakemake pipeline, and I had two reference genomes. One with fasta extension, and the other without.

I also had two rules to index these files using samtools faidx. The first rule had the .fai file created, and the second rule did not.

Why is this?

samtools

what were the command lines , what were the error messages ?

I figured out the error by only using one index instead of two, but I was just wondering why that was the case that two index rules won't necessarily work.

In response to Istvan, I will post the actual rules later.

Some bioinformaticians, when confronted with a problem, think

"I know, I'll use snakemake."

Now they have two problems.


That being said, you have to show the rules otherwise how could people figure out what is wrong?

0 answers

No answers yet.

Log in to answer this question.