Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Help!, strange fastq format

    Dear All

    I download some fastq reads from SRA, but the format is wired. The reads are all number and dots.
    Please give me some pointers if you have seen this before? Thank you in advance.

    example:
    Code:
    @SRR1175538.1 1_11_419_F5-RNA length=35
    G.12.1.21.2.20.1.3013..321..20233033
    +SRR1175538.1 1_11_419_F5-RNA length=35
    !!@@!@!@@!@!@@!@!@@@@!!@@@!!@@@?4<@/
    @SRR1175538.2 1_12_62_F5-RNA length=35
    G.02.0.01.2301.1.1222.1333.111102223
    +SRR1175538.2 1_12_62_F5-RNA length=35
    !!@@!@!@@!@@@=!@!@@@@!@@2@!@?>?@@@@@
    @SRR1175538.3 1_12_580_F5-RNA length=35
    G.10.0.13.0032.2.3123.1121.022330003
    +SRR1175538.3 1_12_580_F5-RNA length=35
    !!<@!@!@?!@@@8!@!@2@@!@;@@!@;@@@2@@@
    @SRR1175538.4 1_12_1917_F5-RNA length=35
    G.00.2.12.3.10.1.1133.2333.100133010
    +SRR1175538.4 1_12_1917_F5-RNA length=35
    !!@@!@!@@!@!@@!@!@@@@!@@?@!@@@@@@@@@
    @SRR1175538.5 1_13_1110_F5-RNA length=35
    G132.2120.22320023122.3200.320102100
    +SRR1175538.5 1_13_1110_F5-RNA length=35
    !@@@!@?@A!5=833052/3-!@?.2!@-55/@6.?
    @SRR1175538.6 1_13_1767_F5-RNA length=35
    G002.3231.03320230131.1020.131301132
    +SRR1175538.6 1_13_1767_F5-RNA length=35
    !@@@!@@8@!@@//76*?<0@!@20@!@62/@@=-?
    @SRR1175538.7 1_14_646_F5-RNA length=35
    G232.1032.21300133222.2022.233022033
    +SRR1175538.7 1_14_646_F5-RNA length=35
    !@@@!@@@@!@@@@<@@@@@@!@@@/!@@;@@@/=@
    @SRR1175538.8 1_15_1325_F5-RNA length=35
    G331.2201.30111132130.0011.202011230
    +SRR1175538.8 1_15_1325_F5-RNA length=35
    !@@=!@*@@!@@@@@@@@@@>!@@@@!@@@?@@@?@
    @SRR1175538.9 1_15_1850_F5-RNA length=35
    G110.3000.03303201101.1013.001113220
    +SRR1175538.9 1_15_1850_F5-RNA length=35
    !:@0!9:@*!=2;.6/[email protected]!@*0:!@*.@;@.3:
    Last edited by GenoMax; 05-19-2017, 09:23 AM.

  • #2
    It's data from the SOLiD platform.

    They use colorspace rather than basespace,
    so you get a sequence of 0,1,2,3 rather than A,T,G,C.

    The . indicate missing base calls.

    Comment


    • #3
      Thank you.

      Originally posted by mastal View Post
      It's data from the SOLiD platform.

      They use colorspace rather than basespace,
      so you get a sequence of 0,1,2,3 rather than A,T,G,C.

      The . indicate missing base calls.

      Comment


      • #4
        If you want to map those to a reference:
        (1) Don't. Delete the file and pretend like it never existed.
        (2) No, really, I mean it. Delete the file.
        (3) Okay, the reason there are so many "." in the sequence is that the SOLiD gave you reads from every bead it found, no matter how bad the data was. Also the .fastq's started at the very top of the slide, near the edge, where the data is worst.
        (4) But you don't need to know this, right? Please tell me you took my advice and deleted the file?
        (5) I mean, you can't just take the numbers and convert them to sequence. Any string of numbers could encode any of 4 sequences. It is the base at the beginning of each sequence that tells you the right conversion path. But if any "base" in the read was incorrect, it corrupts the conversion path. So if you did map these, you would want to use an old-time mapper in "colorspace" mode to do the mapping. BWA and Bowtie probably can still do this type of mapping.
        (6) But you don't need to know that because you deleted the file, right? If not, well, I did warn you...

        --
        Phillip

        Comment


        • #5
          Haha

          But, yes, pmiguel is completely correct and I suggest you follow his advice.

          Well, actually, I wrote really good software for mapping and variant-calling of Solid reads. But then I deleted them because, well, the platform was obsolete. And unfortunately, UTSW is very strict about preventing employees from providing knowledge to the rest of the world... their legal department told me very bluntly that all of the taxpayer- and grant-funded development that I did was their property, which they choose to keep secret.

          So! The Solid platform is terrible. I wrote software that actually makes it useful. But I'm not allowed to share it with you.

          Still, no matter how good the software is, Illumina is vastly better than Solid, so it's best to discard the Solid data and use an alternative platform.
          Last edited by Brian Bushnell; 05-19-2017, 10:55 AM.

          Comment


          • #6
            I completely agree with pmiguel:



            Discussion of next-gen sequencing related bioinformatics: resources, algorithms, open source efforts, etc

            Comment

            Latest Articles

            Collapse

            • seqadmin
              Strategies for Sequencing Challenging Samples
              by seqadmin


              Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
              03-22-2024, 06:39 AM
            • 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

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by seqadmin, Yesterday, 06:37 PM
            0 responses
            7 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, Yesterday, 06:07 PM
            0 responses
            7 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-22-2024, 10:03 AM
            0 responses
            49 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-21-2024, 07:32 AM
            0 responses
            66 views
            0 likes
            Last Post seqadmin  
            Working...
            X