Hello everybody.
i'm a new user of bwa to align solexa reads to the reference genome hg19.i have finished bwa index for the hg19 genome,but not very clear how to use bwa aln command,especially all the paremeters in bwa aln:
i want to know the following conception such as:
1. what's seed length? whether it is equal to reads length,my reads length is 75bp
2. how to set number of threads? How much is best?
3. how to set -q INT quality threshold for read trimming down to 35bp [0]
why 35bp?
how to set this paremeter for my 75bp reads.
4. which options are essential just like the following paremeters?
thank you all very much. yours sincerely Alex. your helps are highly appreciated.
my email is [email protected]
bwa aln ********************
Options: -n NUM max #diff (int) or missing prob under 0.02 err rate (float) [0.04]
-o INT maximum number or fraction of gap opens [1]
-e INT maximum number of gap extensions, -1 for disabling long gaps [-1]
-i INT do not put an indel within INT bp towards the ends [5]
-d INT maximum occurrences for extending a long deletion [10]
-l INT seed length [32]
-k INT maximum differences in the seed [2]
-m INT maximum entries in the queue [2000000]
-t INT number of threads [1]
-M INT mismatch penalty [3]
-O INT gap open penalty [11]
-E INT gap extension penalty [4]
-R INT stop searching when there are >INT equally best hits [30]
-q INT quality threshold for read trimming down to 35bp [0]
-c input sequences are in the color space
-L log-scaled gap penalty for long deletions
-N non-iterative mode: search for all n-difference hits (slooow)
i'm a new user of bwa to align solexa reads to the reference genome hg19.i have finished bwa index for the hg19 genome,but not very clear how to use bwa aln command,especially all the paremeters in bwa aln:
i want to know the following conception such as:
1. what's seed length? whether it is equal to reads length,my reads length is 75bp
2. how to set number of threads? How much is best?
3. how to set -q INT quality threshold for read trimming down to 35bp [0]
why 35bp?
how to set this paremeter for my 75bp reads.
4. which options are essential just like the following paremeters?
thank you all very much. yours sincerely Alex. your helps are highly appreciated.
my email is [email protected]
bwa aln ********************
Options: -n NUM max #diff (int) or missing prob under 0.02 err rate (float) [0.04]
-o INT maximum number or fraction of gap opens [1]
-e INT maximum number of gap extensions, -1 for disabling long gaps [-1]
-i INT do not put an indel within INT bp towards the ends [5]
-d INT maximum occurrences for extending a long deletion [10]
-l INT seed length [32]
-k INT maximum differences in the seed [2]
-m INT maximum entries in the queue [2000000]
-t INT number of threads [1]
-M INT mismatch penalty [3]
-O INT gap open penalty [11]
-E INT gap extension penalty [4]
-R INT stop searching when there are >INT equally best hits [30]
-q INT quality threshold for read trimming down to 35bp [0]
-c input sequences are in the color space
-L log-scaled gap penalty for long deletions
-N non-iterative mode: search for all n-difference hits (slooow)
Comment