SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Nextera Expanded Exome on MiSeq? JosieSeq Illumina/Solexa 3 10-01-2013 07:29 AM
bwa aln fails mike2vandy Bioinformatics 4 12-08-2012 08:23 AM
too small output .sai file created by bwa aln? ..leads to bwa sampe hanging? Stina Bioinformatics 12 12-05-2012 07:32 AM
Questions about bwa aln -n and -o komais Bioinformatics 0 10-19-2011 07:56 PM
Bwa aln aleferna Bioinformatics 1 07-25-2010 10:12 PM

Reply
 
Thread Tools
Old 10-03-2013, 11:32 PM   #1
gmarco
Member
 
Location: Spain

Join Date: Oct 2012
Posts: 36
Unhappy Nextera Expanded Exome BWA aln

Hello,

I'm trying to align Nextera Expanded Exome (62Mb) samples from HiSeq 2000 (Paired end ~99bp insert size) with BWA aln.

I'm having some troubles. Calculate SA coordinate seems to take a lot more time than for TruSeq 71Mb exome. More than 24 hours at this moment to align R1 and still has not finished. All my R1 and R2 are split in 4M reads (output from Illumina CASSAVA software).

Some aligns are stuck here:
Code:
[bwa_aln_core] calculate SA coordinate...
Others have advanced a bit more bit still unusually slow for more than 24 hours:

Code:
[bwa_aln] 17bp reads: max_diff = 2
[bwa_aln] 38bp reads: max_diff = 3
[bwa_aln] 64bp reads: max_diff = 4
[bwa_aln] 93bp reads: max_diff = 5
[bwa_aln] 124bp reads: max_diff = 6
[bwa_aln] 157bp reads: max_diff = 7
[bwa_aln] 190bp reads: max_diff = 8
[bwa_aln] 225bp reads: max_diff = 9
[bwa_read_seq] 6.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 61080.65 sec
[bwa_aln_core] write to the disk... 0.08 sec
[bwa_aln_core] 262144 sequences have been processed.
[bwa_read_seq] 6.7% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 60561.20 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 524288 sequences have been processed.
[bwa_read_seq] 7.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 59024.93 sec
[bwa_aln_core] write to the disk... 0.08 sec
[bwa_aln_core] 786432 sequences have been processed.
[bwa_read_seq] 7.8% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 58316.68 sec
[bwa_aln_core] write to the disk... 0.07 sec
[bwa_aln_core] 1048576 sequences have been processed.
[bwa_read_seq] 8.0% bases are trimmed.
[bwa_aln_core] calculate SA coordinate... 57773.30 sec
[bwa_aln_core] write to the disk... 0.08 sec
[bwa_aln_core] 1310720 sequences have been processed.
[bwa_read_seq] 5.5% bases are trimmed.
[bwa_aln_core] calculate SA coordinate...
Thanks

Thank you

Last edited by gmarco; 10-03-2013 at 11:34 PM. Reason: missing data
gmarco is offline   Reply With Quote
Reply

Tags
align, bwa, exome, expanded, nextera

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:05 AM.


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