Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Reformatting BLAST+ pairwise alignment?

    Greetings!

    The default for BLAST+ output is a pairwise alignment broken up into 60bp increments:

    Code:
    Query  1     TTATTTTTCTTGATAGAATTAATTCTAATATTAGCACATTAAAAAATACTAAGTC-AAAT  59
                 ||||||| ||||||||  |||||| |||||||| |||||||||||||| |||||| ||||
    Sbjct  542   TTATTTTGCTTGATAGCGTTAATTATAATATTACCACATTAAAAAATATTAAGTCAAAAT  601
    
    Query  60    TGTAATATTTCATTTGATATTGTATGCATACAATTTT-ATC---T--TACCTTTAAGAAT  113
                 ||||||||||||||||||  ||||||||||||||||| |||   |  ||||||||| |||
    Sbjct  602   TGTAATATTTCATTTGAT-CTGTATGCATACAATTTTAATCTTATCATACCTTTAATAAT  660
    
    Query  114   ATGTATAA-TAATTTATACCCATGTAGTT-TATTAAAT--A-AACAT-ATTTGTGCGTAT  167
                 |||||||| |||| |||| | | ||| |   || ||||  |    || |||| | |||||
    Sbjct  661   ATGTATAATTAATATATATCAACGTATTCGAATCAAATGGATTTAATGATTTTTTCGTAT  720
    
    Query  168   CCAATGCGTTTATCTATTGTTATTACTA-T--TA-T--TATTGTTACTGCTATATCACTG  221
                   ||| |||||||||||||||||||||| |  || |  |||| ||| ||| || ||||||
    Sbjct  721   TTAATACGTTTATCTATTGTTATTACTATTACTATTACTATTATTATTGCCATCTCACTG  780
    I inquired with NCBI about a command for BLAST+ (neƩ StandALone BLAST) to have the aligment break at different increments or not at all. Like this:

    Code:
    Query  1     TTATTTTTCTTGATAGAATTAATTCTAATATTAGCACATTAAAAAATACTAAGTC-AAATTGTAATATTTCATTTGATATTGTATGCATACAATTTT-ATC---T--TACCTTTAAGAAT 113
                 ||||||| ||||||||  |||||| |||||||| |||||||||||||| |||||| ||||||||||||||||||||||  ||||||||||||||||| |||   |  ||||||||| |||
    Sbjct  542   TTATTTTGCTTGATAGCGTTAATTATAATATTACCACATTAAAAAATATTAAGTCAAAATTGTAATATTTCATTTGAT-CTGTATGCATACAATTTTAATCTTATCATACCTTTAATAAT  660
    They forwarded the request on to the development team, but don't have that functionality at the time in BLAST+. I was hoping someone might have had the same thought and wrote a little script to reformat these alignments?

    Thanks!
    Last edited by Genomics101; 12-04-2013, 04:42 PM.

  • #2
    Use a different output format (ASN, CSV, or XML). The default text-based format is only meant to be for viewing, not data manipulation, and can change without notice.

    Comment


    • #3
      Save your blast result as XML and process it with a XSLT stylesheet.

      Examples:

      * http://www.biostars.org/p/6635/

      * http://www.biostars.org/p/85258/

      Comment

      Latest Articles

      Collapse

      • seqadmin
        Quality Control Essentials for Next-Generation Sequencing Workflows
        by seqadmin




        Like all molecular biology applications, next-generation sequencing (NGS) workflows require diligent quality control (QC) measures to ensure accurate and reproducible results. Proper QC begins at nucleic acid extraction and continues all the way through to data analysis. This article outlines the key QC steps in an NGS workflow, along with the commonly used tools and techniques.

        Nucleic Acid Quality Control
        Preparing for NGS starts with isolating the...
        02-10-2025, 01:58 PM
      • seqadmin
        An Introduction to the Technologies Transforming Precision Medicine
        by seqadmin


        In recent years, precision medicine has become a major focus for researchers and healthcare professionals. This approach offers personalized treatment and wellness plans by utilizing insights from each person's unique biology and lifestyle to deliver more effective care. Its advancement relies on innovative technologies that enable a deeper understanding of individual variability. In a joint documentary with our colleagues at Biocompare, we examined the foundational principles of precision...
        01-27-2025, 07:46 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by seqadmin, 02-07-2025, 09:30 AM
      0 responses
      67 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 02-05-2025, 10:34 AM
      0 responses
      104 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 02-03-2025, 09:07 AM
      0 responses
      83 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 01-31-2025, 08:31 AM
      0 responses
      45 views
      0 likes
      Last Post seqadmin  
      Working...
      X