Align sequencing reads to long reference prokaryotic genome sequences using Bowtie2.
This App aligns the sequencing reads for a set of two or more samples to relatively long reference sequences, usually genomes, using Bowtie2 and outputs a set of alignments in BAM format. Alternatively, if the user has only a single sample of reads or a set of reads, this App will take single-end or paired-end library object instead of a sample set.
In addition, this App outputs the alignment statistics such as total reads, total reads mapped and unmapped, number of singletons, multiple alignments, and overall successful alignment rate in the table format.
Bowtie2 supports end-to-end read alignment mode by default. In this mode, it searches for alignments involving all of the read characters. This is also called an "untrimmed" or "unclipped" alignment. The App also supports "local" alignment mode, which produces alignments that might be "trimmed" (or "soft clipped") at one or both extremes in a way that optimizes the alignment score.
NOTE: This App is one of the steps in the Transcriptomics and Expression Analysis workflow in KBase, however it can also be run as a standalone for one or more RNA-seq samples.
Team members who developed & deployed algorithm in KBase: Srividya Ramakrishnan, Sunita Kumari, Shinjae Yoo, Priya Ranjan, Jim Thomason, and Vivek Kumar. For questions, please contact us.
- Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9: 357 359. doi:10.1038/nmeth.1923 , https://www.nature.com/articles/nmeth.1923
- Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10: R25. doi:10.1186/gb-2009-10-3-r25 , https://genomebiology.biomedcentral.com/articles/10.1186/gb-2009-10-3-r25
Module Commit: 6a2636c6f2a9723f50855e4c693a79a5ce2c464f