...
- Prepare the vibCho reference index for bowtie2 from GenBank records
- Align reads using bowtie2, producing a SAM file
- Convert the SAM file to a BAM file (samtools view)
- Sort the BAM file by genomic location (samtools sort)
- Index the BAM file (samtools index)
- Gather simple alignment statistics (samtools flagstat and samtools idxstatidxstats)
Obtaining the GenBank records
...