This function runs the bwa program using the "aln" method for
aligning sequences. It is appropriate for matching most types of
short-read data.
The GUI is little more than a wrapper around command line tools,
which can essentially be repeatedly manually as follows.
-
Calculate and save the consensus for all contigs in the database in
fastq format.
-
Index the consensus sequence using "bwa index".
-
Map our input data against the bwa index using "bwa aln".
Repeat for reverse matches too.
-
Generate SAM format from the alignments using "bwa samse" or "bwa
sampe".
-
Convert to BAM and sort by position.
-
Import the BAM file, appending to the existing gap5 database
(equivalent to
tg_index -a
).
This function runs the bwa program using the "dbwtsw" method for
aligning sequences. This should be used when attempting to align
longer sequences or data with lots of indels.
The GUI is little more than a wrapper around command line tools,
which can essentially be repeatedly manually as follows.
-
Calculate and save the consensus for all contigs in the database in
fastq format.
-
Index the consensus sequence using "bwa index".
-
Map our input data against the bwa index using "bwa dbwtsw".
-
Convert to BAM and sort by position.
-
Import the BAM file, appending to the existing gap5 database
(equivalent to
tg_index -a
).
Last generated on 25 November 2011.