Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Splitting fasta AND quality on masked sequence

    Hello,
    does anybody have a script or any program that allows me to split fasta sequences and the according quality file on masked bases such as "N"? Everything I find works perfectly on fasta files but the quality files are never touched. I think even FastQ format would be ok because conversion is easy these days using Galaxy.

    Any help would be appreciated.

    Cheers,

    Markus

    Example:

    >fasta
    GACTAGCGATCGACGATCGNNNNNNGACGAGCGACGAGCGA
    >qual
    20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 10 10 10 10 10 10 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20

    >fasta_new1
    GACTAGCGATCGACGATCG
    >fasta_new2
    GACGAGCGACGAGCGA
    >qual_new1
    20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20
    >qual_new2
    20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20

  • #2
    I don't think you'll find any tool. Sounds like is time to put your scripting skills to work ?
    -drd

    Comment


    • #3
      You can convert to fastQ using a perl script helpfully provided in these forums - I've done this in the past, you'll have to search for it. Then as you say manipulation is easier.

      Comment


      • #4
        My scripting skills are the problem... I know some basic shell stuff like grep, sed, awk etc. I come more from the wetlab. Using pre-existing tools and some mild reformatting is ok for me. But I am not sure how to approach this...
        I am still trying to learn some biopython in my spare time...

        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
        10 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, Yesterday, 06:07 PM
        0 responses
        9 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 03-22-2024, 10:03 AM
        0 responses
        50 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 03-21-2024, 07:32 AM
        0 responses
        67 views
        0 likes
        Last Post seqadmin  
        Working...
        X