Samtools faidx can't be run in succession on two reference files?
0
0
Entering edit mode
18 months ago
DdogBoss ▴ 20

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 • 796 views
ADD COMMENT
1
Entering edit mode

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

ADD REPLY
0
Entering edit mode

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.

ADD REPLY
1
Entering edit mode

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?

ADD REPLY

Login before adding your answer.

Traffic: 3000 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6