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