Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • #16
    I have questions regarding to reads from Solid 5500 mapping.

    Dear all,
    I have recently been working on dataset (75bps single end reads/fragment) output from solid 5500. I used bowtie to map to mouse genome, but the result was very unreasonable. There were only about 10 percent of total reads could be mapped to mouse genome.

    Then, I used blat from UCSC to blat the unmappable reads back to mouse genome. I came across that unmappable reads from bowtie output had many output after I "blatted". Then, I took a look at the result from blat, I noticed that most of them could map/align to the mouse genome with shorter length. For instance, the original 75bps read that could not map to mouse genome with bowtie, I got some output with shorter reads (46 bps) from blat. I wonder where the rest of the reads went?

    p.s: almost all the unmappable reads from bowtie got output from blat with shorter reads length aligned to mouse genome.

    Mike

    Comment


    • #17
      Are you mapping in color space? If not, were these ECC reads, or just normal chemistry? If at all possible, you want to do your mapping in color space. Bowtie allows you to do this.

      Have you run the data set through a QC checking program? Eg, FastQC? Might be worth it to see if the quality dropped off dramatically at some point in the run.

      --
      Phillip

      Comment


      • #18
        The csfasta and qual file obtained from the solid XSQ barcode separation process, then color space was converted to nucleotide (convert the file to fastsanger format in galaxy then used "manipulate FASTQ" tool in galaxy to convert to NT.)

        They were not ECC reads, it was just normal chemistry which followed ABI protocol. By the way, is it normal to amplify to 2^20 cycle for a library ? thanks in advance.


        Mike

        Comment


        • #19
          Originally posted by acadMT View Post
          For instance, the original 75bps read that could not map to mouse genome with bowtie, I got some output with shorter reads (46 bps) from blat. I wonder where the rest of the reads went?
          They went off into colorspace. If you absolutely need to map your reads in basespace you could try bowtie2 with the --local option, but mapping in colorspace is the only way to get reasonable results.

          Comment


          • #20
            Hi Mike,

            as others have said map in colour space using the csfasta and qual files. Alternatively to bowtie you could try bwa, though both don't seem great for 60 bp data from the 5500xl (33% and 40% mapping rate for human SE data). It seems the bwa developers are considering stopping their

            Bioscope and Clc Genomics workbench seem a lot better (~85 % mapping rate), in part because they (or at least Bioscope) iteratively trim the sequences.

            Novoalign-CS is also supposed to be very promising, but we are just in the process of getting a license.

            Comment

            Latest Articles

            Collapse

            • seqadmin
              Recent Advances in Sequencing Analysis Tools
              by seqadmin


              The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...
              05-06-2024, 07:48 AM
            • seqadmin
              Essential Discoveries and Tools in Epitranscriptomics
              by seqadmin




              The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...
              04-22-2024, 07:01 AM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by seqadmin, Yesterday, 06:35 AM
            0 responses
            14 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 05-09-2024, 02:46 PM
            0 responses
            19 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 05-07-2024, 06:57 AM
            0 responses
            18 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 05-06-2024, 07:17 AM
            0 responses
            19 views
            0 likes
            Last Post seqadmin  
            Working...
            X