SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
BWA: Is alignment summary available? Protaeus Bioinformatics 2 07-05-2011 06:12 PM
BWA alignment followed by TopHat slny Bioinformatics 9 06-06-2011 06:19 AM
Spliced alignment with BWA telos SOLiD 7 10-06-2010 06:32 AM
Alignment statistics from bwa? menenuh Bioinformatics 1 01-29-2010 01:14 PM
BWA Alignment Errors Seq1234 Bioinformatics 0 12-15-2009 08:59 PM

Reply
 
Thread Tools
Old 10-21-2011, 11:43 AM   #1
CTGF
Junior Member
 
Location: US

Join Date: Oct 2011
Posts: 7
Default bwa alignment failed

Hello, I am trying to align my Illumina hiseq data using bwa:

./bwa index -c -a bwtsw -p hg19 hg19.fa

./bwa aln -t 4 -c hg19 quatrim.fq > quatrim.sai

./bwa samse -f quatrim.sam hg19 quatrim.sai quatrim.fq

Then .sam was convered to .bam

I used samstool to check the alignment:

./samtools idxstats quatrim.bam

chr1 249250621 1 0
chr2 243199373 0 0
chr3 198022430 2 0
chr4 191154276 1 0
chr5 180915260 1 0
chr6 171115067 0 0
chr7 159138663 0 0
chr8 146364022 1 0
chr9 141213431 0 0
chr10 135534747 0 0
chr11 135006516 1 0
chr12 133851895 0 0
chr13 115169878 0 0
chr14 107349540 1 0
chr15 102531392 0 0
chr16 90354753 0 0
chr17 81195210 0 0
chr18 78077248 0 0
chr19 59128983 3 0
chr20 63025520 1 0
chr21 48129895 0 0
chr22 51304566 0 0
chrX 155270560 0 0
chrY 59373566 0 0
chrM 16571 0 0
* 0 0 3081338

When I ran the above steps, everything looked fine. The quality scores of the reads are also very good. I blasted some, and they mapped perfectly to the genome. I am wondering which step(s) might have problems. Thanks!
CTGF is offline   Reply With Quote
Old 10-24-2011, 12:35 AM   #2
Bukowski
Senior Member
 
Location: UK

Join Date: Jan 2010
Posts: 390
Default

bwa aln -c is for colorspace data (i.e. from SOLiD platforms). As is bwa index -c

You state you have Illumina data.

I suggest removing the -c options from both commands.
Bukowski 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 01:12 AM.


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