Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • unrecognized .txt format

    Hi,

    I recently received data from illumina, and I am unable to determine the file format or how to convert it to a fastq for bwa. I found lots of documentation on a format called qseq, and several conversion scripts, but it doesn't seem to be the format I'm dealing with and the scripts don't work.

    The lines in the files look like this:

    Code:
    NAME:3:1:1138:1083#0/1:ATATATGAGTGGGGTATAAACGGAGATTGATTTCGT:abababbababbbb`bbaaa]b_ab`aa]^_a_``^

  • #2
    It looks like SCARF using FASTQ style encoding for the quality scores.

    Looks like you can just split the lines on the colon, something like this (in Python) should give you a FASTQ file:

    Code:
    #!/usr/bin/env python
    """
    Simple script to generate FASTQ from colon separated SCARF input 
    where the quality scores are ASCII encoded (FASTQ like).
    
    Use at the command line, piping the input and output.
    """
    import sys
    for line in sys.stdin:
        name, seq, qual = line.rstrip("\n").rsplit(":",2)
        assert len(seq) == len(qual), line
        sys.stdout.write("@%s\n%s\n+\n%s\n" % (name, seq, qual))
    The same idea would be equally trivial in Perl or your language of choice.
    Last edited by maubp; 11-17-2010, 05:51 AM.

    Comment


    • #3
      Thanks a lot! It seems to be working. Any idea why these files are in this format? They're supposed to be illumina files.

      Comment


      • #4
        SCARF is one of the many formats used by Solexa/Illumina, it gets confusing

        Comment

        Latest Articles

        Collapse

        • 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
        • seqadmin
          Current Approaches to Protein Sequencing
          by seqadmin


          Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
          04-04-2024, 04:25 PM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, Today, 08:47 AM
        0 responses
        10 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 04-11-2024, 12:08 PM
        0 responses
        60 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 04-10-2024, 10:19 PM
        0 responses
        57 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 04-10-2024, 09:21 AM
        0 responses
        53 views
        0 likes
        Last Post seqadmin  
        Working...
        X