Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • SOLiD paired end data and BWA - cleaning/mapping

    Hello everyone,

    I notice there has been the odd post about this general topic before but I am still confused!

    I have SOLiD single and paired end data which I would like to take through the following steps:
    1) Quality assessment - Fastqc
    2) Data cleaning - Quality trimming/removing adapters if necessary
    3) Mapping using BWA

    At the moment I am converting the csfasta and qual files in to double-encoded (fakespace) fastq files for BWA. I reverse the F5 reads so BWA treats them correctly. However, can the same fastq files be used with Fastqc? Also, can CutAdapt and fastq_quality_trimmer handle this format? If anyone else out there still uses SOLiD data and can give me some details of their workflows I would be very grateful.

    Best wishes,
    Kath

  • #2
    I used to work with SOLiD data, and I have to ask why you are going to all of this trouble. In my experience, the best way to map SOLiD data is to stay in colorspace, and use Bfast, or mapreads (part of lifescope/bioscope).


    Brian

    Comment


    • #3
      Thanks for your reply. I have Lifescope mapped BAM files but I'm trying to dedup using Picard and am getting some very high duplication rates. I've been told that Picard seems to overestimate duplication rates in Lifescope mapped data and that BWA provides better compatibility with Picard. I did look into using Bfast but previous benchmarking studies suggested BWA would be better.

      Anyway, it looks like there is a script - solid-trimmer.py - that will do what I want.

      Comment


      • #4
        OK. Are you working with DNA-seq, or RNA-seq? If it is DNA seq, I would be worried about not working in colorspace, if it is RNA-seq, removing duplicates from SOLiD short read data may not be a good idea. Just my experience.

        Brian

        Comment


        • #5
          We tested all SOLiD mappers very thoroughly, and NovoalignCS (commercial but not too expensive) won hands-down. I would really recommend it.

          BWA dropped color space support a long time ago and did not perform very well in terms of # mapped reads. Bfast wasn't much fun either.

          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
          82 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