This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Bowtie error: Could not locate a Bowtie index corresponding to basename

Hi everyone! I am an undergrad and pretty new in ChIP-seq. And I need to use bowtie with hg19. I indexed hg19 with bowtie-build and completed it. I would like to run my single-end fastq files with bowtie. I am using this to run bowtie.

bowtie --large-index </home/cilga/Desktop/VPC/ref/bowtie_indexes/bt1/hg19> -q /home/cilga/Desktop/VPC/chipseq/samples/dht_treatment/SRR3728821_pass.fastq  -S SRR3728821.sam

I am not sure whether I do wrong. And I got this error:

Could not locate a Bowtie index corresponding to basename "/home/cilga/Desktop/VPC/chipseq/samples/dht_treatment/SRR3728821_pass.fastq" Command: bowtie --wrapper basic-0 -S /home/cilga/Desktop/VPC/chipseq/samples/dht_treatment/SRR3728821_pass.fastq SRR3728821.sam

Is there anyone can help me? Thank you for your helps.

PS.. can someone write an example command for bowtie?

bowtie chip-seq single-end hg19 alignment
• 8,833 views • link
updated by h.mon h.mon • written by cilgaiscan • cilgaiscan

I believe the error message does not match with the command you pasted above. Can you check again that it is effectively the command you used?

• 1 views • link
written by WouterDeCoster WouterDeCoster

yes i double checked. i use same command and got same error..

• 1 views • link
written by cilgaiscan • cilgaiscan

The error says:

Command: bowtie --wrapper basic-0 -S /home/cilga/Desktop/VPC/chipseq/samples/dht_treatment/SRR3728821_pass.fastq SRR3728821.sam

Are you sure you didn't swap some arguments there?

• 1 views • link
written by WouterDeCoster WouterDeCoster

i am.i used the arguments i shared in my question :(

• 1 views • link
written by cilgaiscan • cilgaiscan

cilgaiscan : You are not actually including <> around the large index file location in your command, correct?

• 1 views • link
written by GenoMax GenoMax

1 answer

Try the following:

bowtie -S -q --large-index /home/cilga/Desktop/VPC/ref/bowtie_indexes/bt1/hg19 /home/cilga/Desktop/VPC/chipseq/samples/dht_treatment/SRR3728821_pass.fastq > SRR3728821.sam

The problem, I believe, is for bowtie -S is just a boolean flag, and you are using it as in bowtie2, where -S file.sam creates the output file, instead of using stdout.

• 1 views • link
written by h.mon h.mon

Log in to answer this question.

More posts like this