Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Split long fasta seq into smaller seqs by removing 'n'

    Hello,

    My fasta file has a long consensus sequence (gigabases long) that is padded with 'n' between the actual sequences. Like this:

    actgggacnnnnnnnnnnnnnnnnnnnnnnnnnnactgacgtggattgc
    aatnnnnnnnnnnnnnnnaccaattggatagagaccnnn

    I've searched rather intensively to see if the internet can solve my problem. Not really, though they were close. For example, people report splitting the long sequence into smaller files but I do not want that. I want the n removed and have all sequences in the same file. Better yet, with the start and end of the sequences.

    Desired output (>seqname_start_end) :
    >seq1_1_8
    actgggac

    >seq2_35_52
    actgacgtggattgcaat

    >seq3_68_85
    accaattggatagagacc

    If anyone could point me towards a right tool (bioperl, etc) or give me a pseudo-code in perl, I would appreciate it.

    Thanks.

  • #2
    IDBA contains a script called split_scaffold

    Contribute to loneknightpy/idba development by creating an account on GitHub.


    Here's what it does:

    IN =>

    >test
    acgacgacgacgacgacgacgacagcannnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnacgacgacgacagcagtagatgatgatagtag


    OUT =>

    >teste_0
    ACGACGACGACGACGACGACGACAGCA
    >teste_1
    ACGACGACGACAGCAGTAGATGATGATAGTAG

    Comment


    • #3
      Originally posted by azneto View Post
      IDBA contains a script called split_scaffold

      Contribute to loneknightpy/idba development by creating an account on GitHub.

      Great! Thanks a lot! Worked flawlessly.

      It doesn't get me the start/end of the base positions but that I can live with.

      Comment

      Latest Articles

      Collapse

      • seqadmin
        Understanding Genetic Influence on Infectious Disease
        by seqadmin




        During the COVID-19 pandemic, scientists observed that while some individuals experienced severe illness when infected with SARS-CoV-2, others were barely affected. These disparities left researchers and clinicians wondering what causes the wide variations in response to viral infections and what role genetics plays.

        Jean-Laurent Casanova, M.D., Ph.D., Professor at Rockefeller University, is a leading expert in this crossover between genetics and infectious...
        09-09-2024, 10:59 AM
      • seqadmin
        Addressing Off-Target Effects in CRISPR Technologies
        by seqadmin






        The first FDA-approved CRISPR-based therapy marked the transition of therapeutic gene editing from a dream to reality1. CRISPR technologies have streamlined gene editing, and CRISPR screens have become an important approach for identifying genes involved in disease processes2. This technique introduces targeted mutations across numerous genes, enabling large-scale identification of gene functions, interactions, and pathways3. Identifying the full range...
        08-27-2024, 04:44 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by seqadmin, 09-11-2024, 02:44 PM
      0 responses
      12 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 09-06-2024, 08:02 AM
      0 responses
      145 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 09-03-2024, 08:30 AM
      0 responses
      152 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 08-27-2024, 04:40 AM
      0 responses
      163 views
      0 likes
      Last Post seqadmin  
      Working...
      X