Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • lindenb
    replied
    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/

    Leave a comment:


  • gringer
    replied
    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.

    Leave a comment:


  • Genomics101
    started a topic Reformatting BLAST+ pairwise alignment?

    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.

Latest Articles

Collapse

  • seqadmin
    Exploring the Dynamics of the Tumor Microenvironment
    by seqadmin




    The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
    07-08-2024, 03:19 PM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 07-25-2024, 06:46 AM
0 responses
9 views
0 likes
Last Post seqadmin  
Started by seqadmin, 07-24-2024, 11:09 AM
0 responses
26 views
0 likes
Last Post seqadmin  
Started by seqadmin, 07-19-2024, 07:20 AM
0 responses
160 views
0 likes
Last Post seqadmin  
Started by seqadmin, 07-16-2024, 05:49 AM
0 responses
127 views
0 likes
Last Post seqadmin  
Working...
X