SEQanswers

Go Back   SEQanswers > Applications Forums > Genomic Resequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
using bwa to map illumina paired end reads mikeworth Bioinformatics 6 08-13-2013 07:02 PM
paired-end reads mapped to genome.. gene with only one direction of paired-end reads? danwiththeplan Bioinformatics 2 09-22-2011 03:06 AM
Bowtie Illumina paired end reads alignment empyrean Bioinformatics 3 09-20-2011 10:51 AM
Sam flags for bwa-aligned paired end reads with identical + / - strand coordinates spark Bioinformatics 0 03-09-2011 05:00 AM
Problems with bwa on Q2 trimmed paired end reads curious_mapper Bioinformatics 2 05-06-2010 03:44 PM

Reply
 
Thread Tools
Old 06-07-2010, 04:39 AM   #1
AvinashP
Junior Member
 
Location: India

Join Date: Jun 2010
Posts: 7
Default BWA alignment for paired end reads

hi,

i am using BWA to align 115bp paired end solexa reads. I have the following doubts :-

1. Does BWA take the min and max insert sizes for paired end reads ?
2. What is maximum number of indels calculated by BWA ? Does keeping the indel as 3 will affect the performance ?
3. How do i get the unaligned reads seperately after using the 'sampe' command? The reads which are not paired and the reads which didn't map at all.
4. I plan to use samtools for SNP and indel detection. Can i get the reference as well as the consensus coordinate in the pile up format.

Any help for the above is appreciated.
Thanks
AvinashP is offline   Reply With Quote
Old 06-07-2010, 07:19 AM   #2
lh3
Senior Member
 
Location: Boston

Join Date: Feb 2008
Posts: 693
Default

1. bwa infers it
2. hard to say for paired-end reads. could be up to 20bp in some cases
3. see sam format
4. it depends on what you exactly want.
lh3 is offline   Reply With Quote
Old 06-08-2010, 04:11 AM   #3
AvinashP
Junior Member
 
Location: India

Join Date: Jun 2010
Posts: 7
Default

Thanks, for your immediate reply.

For (2.) can we restrict the indels with the parameters -o and -e in the 'aln' step? What is the ideal value for these parameters given 115bp reads ?

4. i want the consensus coordinate also in the pileup format with the other fields. How to get it using samtools ?
AvinashP is offline   Reply With Quote
Reply

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 11:05 PM.


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