Header Leaderboard Ad

Collapse

sff 454 files into fasta

Collapse

Announcement

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

  • sff 454 files into fasta

    Dear users, I have question:

    I am trying to convert a .sff 454 file into fasta, using sff_extract in python2.3.4 that includes Biophyton-1.49. I have used the following command:

    $ python sff_extract my_file.sff

    and this results in a syntax error:


    File "sff_extract", line 405

    qual_line = ' '.join((str(q) for q in qual))

    ^

    SyntaxError: invalid syntax


    Does anyone could suggest what's wrong?

    Thank you in advance

    P.R.

  • #2
    That looks to me like some Python syntax requiring Python 2.4 or later. You could change this:

    qual_line = ' '.join((str(q) for q in qual))

    to:

    qual_line = ' '.join([str(q) for q in qual])

    but there are probably going to be more issues. I would urge you to update to Python 2.6, Python 2.5 or Python 2.4 (in that order of preference assuming you have no other constraints).

    You mentioned Biopython 1.49. You should update (e.g. for FASTQ support). However, the more recent releases also required Python 2.4 or later.

    Also, although there is SFF in the pipeline for Biopython is isn't included as yet (i.e. it is not even in the brand new Biopython 1.53 released yesterday). Its on a branch on github, but if you (or anyone else) would like to help test it I am confident we can have this in the next release of Biopython.

    Comment


    • #3
      Originally posted by maubp View Post
      You mentioned Biopython 1.49. You should update (e.g. for FASTQ support). However, the more recent releases also required Python 2.4 or later.

      Also, although there is SFF in the pipeline for Biopython is isn't included as yet (i.e. it is not even in the brand new Biopython 1.53 released yesterday). Its on a branch on github, but if you (or anyone else) would like to help test it I am confident we can have this in the next release of Biopython.
      In case anyone is interested, the SFF support has been merged into the Biopython main branch, and should be included with our next release (Biopython 1.54). Anyone interested in testing it before then can grab the code from our repository.
      http://www.biopython.org/wiki/SourceCode

      Comment


      • #4
        I know it wont help if you are setting up something very customized, but for just extracting the sequences sffinfo -s should do the trick.

        And its included in the recent 2.3 software.

        Comment


        • #5
          Yeah, if you have the Roche Linux only tools, sffinfo does offer SFF to FASTA (or QUAL, with or without trimming). It was there in v2.0 as well, and I expect in v1 too. See also this thread for getting the tools:
          http://seqanswers.com/forums/showthread.php?t=114

          Comment

          Latest Articles

          Collapse

          • seqadmin
            A Brief Overview and Common Challenges in Single-cell Sequencing Analysis
            by seqadmin


            ​​​​​​The introduction of single-cell sequencing has advanced the ability to study cell-to-cell heterogeneity. Its use has improved our understanding of somatic mutations1, cell lineages2, cellular diversity and regulation3, and development in multicellular organisms4. Single-cell sequencing encompasses hundreds of techniques with different approaches to studying the genomes, transcriptomes, epigenomes, and other omics of individual cells. The analysis of single-cell sequencing data i...

            01-24-2023, 01:19 PM
          • seqadmin
            Introduction to Single-Cell Sequencing
            by seqadmin
            Single-cell sequencing is a technique used to investigate the genome, transcriptome, epigenome, and other omics of individual cells using high-throughput sequencing. This technology has provided many scientific breakthroughs and continues to be applied across many fields, including microbiology, oncology, immunology, neurobiology, precision medicine, and stem cell research.

            The advancement of single-cell sequencing began in 2009 when Tang et al. investigated the single-cell transcriptomes
            ...
            01-09-2023, 03:10 PM

          ad_right_rmr

          Collapse
          Working...
          X