Go Back   SEQanswers > General

Similar Threads
Thread Thread Starter Forum Replies Last Post
Qiagen SeqTarget system - Don't go there friend !! lek2k Sample Prep / Library Generation 2 08-18-2011 06:32 PM
Samtools don't read the whole data mikael Bioinformatics 6 08-09-2011 12:37 PM
Don't know if I'm right anna_ Introductions 4 12-13-2010 12:32 PM
Bam and Sam don't like my fasta file mindlessbrain Bioinformatics 2 12-09-2010 11:47 PM
something wrong with tophat-1.0.14 version, don't use it frankyue50 Bioinformatics 4 10-04-2010 09:33 AM

Thread Tools
Old 12-29-2011, 11:35 AM   #1
Location: USA

Join Date: Nov 2011
Posts: 28
Default BWA, I don't know what I'm doing


I'm not really sure I understand how to setup parameters in BWA.

I want to map 25bp non-paired reads to the mouse genome. I try to do that on galaxy...

I made an alignment using the following parameters:
- max edit distance (aln-n) : 0
- fraction on missing alignments given 2% uniform base error rate (aln-n) : 0.04
- max gap open (aln -o) : 1
(I wanted to allow 1 gap)

- max gap extension (aln -e) : -1
- disallow long del within (val)bp towards 3' end (aln -d) : 30
- disallow in/del within (val)bp towards the end (aln -i) : 30
(I guess that screws up my will to allow 1 gap since the reads are 25bp long and gaps are insertions, right?)

- # of 1st subsequences to take as seed (aln -l) : 30
- max edit distance in the seed (aln -k) : 2
- missmatch penalty (aln -M) : 3
- gap open penalty (aln -O) : 11
- gap extension penalty (aln -E) : 4
- max # of alignments to output in the XA tag (samse -n) : 3

I guess this settings should allow 0 gaps and 0 missmatches (MM), right?
However, I end up with an alignment where many reads have 1 or 2 MM.
Did I miss something?

what setting should I use to allow for 1 gap and 3 MM?

Thank you very much in advance for your help.
asheenlevrai is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 03:04 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO