SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
1000 Genomes VCF Files ashkot Bioinformatics 8 02-25-2013 02:12 PM
confused by all the 1000 Genomes fastq files efoss Bioinformatics 3 10-11-2011 06:43 AM
bam files of exomes from 1000 genomes jorjial Bioinformatics 3 06-06-2011 03:16 AM
fastq files comparison to 1000 gp data Masta Bioinformatics 0 02-22-2011 10:27 PM
1000 genomes fastq format uam Bioinformatics 1 12-10-2008 07:51 AM

Reply
 
Thread Tools
Old 10-15-2009, 01:37 PM   #1
tgenahmet
Member
 
Location: Phoenix

Join Date: Apr 2009
Posts: 11
Default using BWA to align SOLiD fastq files from 1000 Genomes

Hi All,

I've been playing with BWA to align some color space data from 1000 Genomes. So far all my alignments are TNNNNNNNNNNNNNNNN's.

I'm trying to figure out the correct procedure for using BWA with (color spaced) fastq files.

In this thread lh3 says:

Quote:
Alignment is done in color space. After alignment, color reads are translated to nucleotides essentially with the code in maq's csmap2nt.
But in this thread it sounds like the input needs to be converted to nucleotide space.

Quote:
bwa uses fastQ files with colors represented as ACGT, perhaps the 1000 genomes fastq files represents it as 0123? Also, bwa does not use the first color so you may have to strip it or use the solid2fastq script.
Here is what my input looks like (straight from 1000 Genomes FTP):
Quote:
==> SRR008434_2.fastq <==
@VAB_SOLiD0097_20080822_2_Pilot_1_Ceph_11994_A_Lib_1_2Kb_MP1561_6_54
T3232322122230111000000220
+
!*(/&$'0++%/##3'%#(#&$%)#(
@VAB_SOLiD0097_20080822_2_Pilot_1_Ceph_11994_A_Lib_1_2Kb_MP1561_6_60
T0233233322023332023032223
+
!$##((#%$#%$,$#('$'#.#,+%#
@VAB_SOLiD0097_20080822_2_Pilot_1_Ceph_11994_A_Lib_1_2Kb_MP1561_6_93
T0102012101023023220003000
My output in sam format:

Quote:
VAB_SOLiD0097_20080822_2_Pilot_1_Ceph_11994_A_Lib_1_2Kb_MP1561_6_54 77 * 0 0 * * 0 0 GNNNNNNNNNNNNNNNNNNNNNNNNN !##$&&&$)#.*(#%#.%$&'*%#$#
VAB_SOLiD0097_20080822_2_Pilot_1_Ceph_11994_A_Lib_1_2Kb_MP1561_6_54 141 * 0 0 * * 0 0 TNNNNNNNNNNNNNNNNNNNNNNNNN !*(/&$'0++%/##3'%#(#&$%)#(
VAB_SOLiD0097_20080822_2_Pilot_1_Ceph_11994_A_Lib_1_2Kb_MP1561_6_60 77 * 0 0 * * 0 0 GNNNNNNNNNNNNNNNNNNNNNNNNN !&%%%%&5)#,'%#$&)$%'$&,%(%
VAB_SOLiD0097_20080822_2_Pilot_1_Ceph_11994_A_Lib_1_2Kb_MP1561_6_60 141 * 0 0 * * 0 0 TNNNNNNNNNNNNNNNNNNNNNNNNN !$##((#%$#%$,$#('$'#.#,+%#
VAB_SOLiD0097_20080822_2_Pilot_1_Ceph_11994_A_Lib_1_2Kb_MP1561_6_93 77 * 0 0 * * 0 0 GNNNNNNNNNNNNNNNNNNNNNNNNN !&8';#'$.'&(?>:&',+6<<,56&
I indexed the reference genome using "bwa index -c", aligned using "bwa aln -c", then a "bwa sampe". No errors were reported in any of the steps.

Any ideas?
tgenahmet is offline   Reply With Quote
Old 10-15-2009, 05:19 PM   #2
pliang
Junior Member
 
Location: Canada

Join Date: Aug 2009
Posts: 9
Default

I had to agree with what chipper said that "bwa uses fastQ files with colors represented as ACGT". I had to give it up using the color space sequence directly. Instead of doing your own alignment, you may try use the alignment data they provided. Sorry I couldn't be more helpful.
pliang 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 03:59 PM.


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