Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • 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 <==
    @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:

    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?

  • #2
    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.

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Techniques and Challenges in Conservation Genomics
      by seqadmin



      The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

      Avian Conservation
      Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
      03-08-2024, 10:41 AM
    • seqadmin
      The Impact of AI in Genomic Medicine
      by seqadmin



      Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
      02-26-2024, 02:07 PM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, 03-14-2024, 06:13 AM
    0 responses
    34 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-08-2024, 08:03 AM
    0 responses
    72 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-07-2024, 08:13 AM
    0 responses
    81 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-06-2024, 09:51 AM
    0 responses
    68 views
    0 likes
    Last Post seqadmin  
    Working...
    X