User Tools

Site Tools


mkatari-bioinformatics-august-2013-bowtienotes

This is an old revision of the document!


Back to Manny's Bioinformatics Workshop Home

Creating and Running Bowtie using Cassava

Make you are in "interactive mode" 1) Download cassava genome from phytozome

wget ftp://ftp.jgi-psf.org/pub/compgen/phytozome/v9.0/Mesculenta/assembly/Mesculenta_147.fa.gz

2) created a directory called cassava

mkdir cassava

3) move fasta sequence in the directory

mv Mesculenta_147.fa.gz cassava/

4) uncompress using gunzip

cd cassava
gunzip Meculenta_147.fa.gz

5) Count number of scaffold in the file

grep ">" Mesculenta_147.fa | wc

6) load bowtie module

module load bowtie2

7) Create the bowtie index.

bowtie2-build Mesculenta_147.fa cassava

8) Run Bowtie using single end fastq as input. See bowtie2 for instructions on how to run it on pair-end sequences.

bowtie2 -x cassava/cassava -U test.fastq -S test.sam
mkatari-bioinformatics-august-2013-bowtienotes.1398860966.txt.gz · Last modified: 2014/04/30 12:29 by mkatari