Can't locate Bowtie2 index files
1
3
Entering edit mode
3.8 years ago
kdc15 ▴ 40

Hi,

I am trying to align my sequences to a reference genome but Bowtie2 is having trouble locating the index files that I downloaded from Illumina. I am working on a server and have read instructions about moving the files to the bowtie2 index subdirectory, however, this does not exist. It also says that bowtie2 is not a directory.

Does anyone know how I can solve this problem?

bowtie2 ChIP-Seq alignment • 4.9k views
ADD COMMENT
0
Entering edit mode

Can you provide a link where you downloaded the sequence/annotation/index bundle from? Bundles from iGenomes contain indexes for multiple aligners. You don't need to move any files into directories etc as long as you provide right paths.

ADD REPLY
0
Entering edit mode

I got the files from this page under pre-built indexes: http://bowtie-bio.sourceforge.net/tutorial.shtml#preb

ADD REPLY
0
Entering edit mode

I am trying to align my sequences to a reference genome but Bowtie2 is having trouble locating the index files that I downloaded from Illumina

@Istvan has a simple guide on how you would use bowtie2 here.

ADD REPLY
0
Entering edit mode

Thank you, this helped!!

ADD REPLY
3
Entering edit mode
3.8 years ago
flogin ▴ 280

To use the basic mapping function with bowtie2 you can do this:

$ bowtie2-build reference_sequence.fasta reference_sequence.fasta.index
$ bowtie2 -x reference_sequence.fasta.index -U dataset_foward.fastq,dataset_reverse.fastq -S output.sam

Parameters: -x = input index, -U=fastq datasets separated by coma (e.g.case of paired-end datasets); -S= output in sam format.

If you want to visualize your mapping analysis on IGV or Artemis (e.g.), you need to use Samtools to convert sam in bam and create a sorted and index files:

$ samtools view -Sb output.sam -o output.bam
$ samtools sort output.bam -o output.sorted.bam
$ samtools index output.sorted.bam

I hope that it helps u.

ADD COMMENT
1
Entering edit mode

Thank you, I have managed to do this now. Looks like I was using Bowtie1 instead of Bowtie2

ADD REPLY

Login before adding your answer.

Traffic: 956 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