Here is a complete list of some important parameters used in both aligners to control the alignment of FastQ data.
BWA
|
Bowtie
|
||
Parameter
|
Description
|
Parameter
|
Description
|
Bwa
alignment options
|
There are different parameters for
alignment which are listed below
|
-u
20
|
This option causes Bowtie to exit
after processing first 20 reads from file
|
-n NUM
|
Maximum edit distance if the value is INT, or
the fraction of missing alignments given 2% uniform base error rate if FLOAT
|
-S
|
To obtain alignment in SAM format
|
-o INT
|
Maximum number of gap opens
|
-l
|
To confirm that output files are
created
|
-e INT
|
Maximum number of gap extensions
|
--suppress
|
Suppresses all outputs if columns
number is not mentioned otherwise suppresses the mentioned columns
|
-a STR
|
Algorithm for constructing BWT index
|
-a
|
Instructs bowtie to report all valid alignments
subject to the alignment policy
|
-d INT
|
Disallow a long deletion within INT bp towards the 3’-end |
-k
|
Reports all the valid alignments
|
-i INT
|
Disallow an indel within INT bp towards the
ends
|
--best
|
Guarantees that bowtie will find and
print the alignments for a given read in best-to-worst order
|
-p STR
|
Prefix of the output database
|
-p
|
Value greater than 1 allows bowtie to
use more than computer processor
|
-l INT
|
Take the first INT subsequence as seed. If INT
is larger than the query sequence, seeding will be disabled.
|
--offrate
|
Adjusts the error message about memory
having been exhausted
|
-k INT
|
Maximum edit distance in the seed
|
-C
|
Enables alignment of color space reads
produced by ABI SOLiD
|
-t INT
|
Number of threads (multi-threading mode)
|
-v <int>
|
Only alignment up to <int>
mismatches in the entire alignment are considered “valid”
|
-O INT
|
Gap open penalty
|
-n <int>
|
Only alignment with up to <int>
mismatches in the seed portion of alignment are considered “valid”
|
-I
|
The input is in the Illumina 1.3+ read format
(quality equals ASCII-64)
|
-l <int>
|
When -n is specified, -l sets the
length of the seed portion of the alignment to <int>
|
-E INT
|
Gap extension penalty
|
-e <int>
|
When -n is specified, -e sets the
limit
|
-b
|
Specify the input read sequence file is the BAM
format. For paired-end data, two ends in a pair must be grouped together and
options -1 or -2 are usually applied to specify which
end should be mapped.
|
-m <int>
|
If bowtie finds more than <int>
valid arguments for a read, suppress all alignments for a read
|
-M INT
|
Mismatch penalty. BWA will not search for
suboptimal hits with a score lower than (bestScore-misMsc)
|
-M <int>
|
If bowtie finds more than <int>
valid alignments for a read, suppress the alignments but report one of them
at random
|
samse
|
Generate alignments in the SAM format given
single-end reads. Repetitive hits will be randomly chosen.
|
-1 is specified for single end
alignment
|
|
sampe
|
Generate alignments in the SAM format given
paired-end reads. Repetitive read pairs will be placed randomly.
|
-2
|
Specified for paired end alignment
|
No comments:
Post a Comment