Go Back   SEQanswers > Bioinformatics > Bioinformatics

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

Thread Tools
Old 10-15-2009, 02:37 PM   #1
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:

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.

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):
==> SRR008434_2.fastq <==
My output in sam format:

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, 06:19 PM   #2
Junior Member
Location: Canada

Join Date: Aug 2009
Posts: 9

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

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 04:08 PM.

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