Question: Base by base alignment from BLAST to reference
0
gravatar for kspata
2.4 years ago by
kspata70
Chicago
kspata70 wrote:

Hi All,

I have aligned a 5k sequence to a reference sequence longer than 5k. The Blast alignment report shows alignment not from the start position of the reference but from bases from where alignment starts. I want an alignment report where reference which does not match to the consensus is also mentioned in the report.

For example (Desired Output);

Query      ~~~~~~~~~~  ~~~~~~~~~~ ~~~~~~~~~~  AGCTTGGCCC  ATTGCATACG  TTGTATCCAT 
Reference  AAACGAGCTC  TGCTTTTATA GGCGCCCACC  AGCTTGGCCC  ATTGCATACG  TTGTATCCAT

Query      ~~~~~~~~~~ ~~~~~~~~~ 
Reference  AATTAAGCGC GCAACGATG

Blast only reports bases which are aligned. Example:

Query     AGCTTGGCCC  ATTGCATACG  TTGTATCCAT 
Reference AGCTTGGCCC  ATTGCATACG  TTGTATCCAT

I can only use BLAST to do this.

How can I obtain such alignment report using Blast?

Thanks in advance!!!

blast alignment • 738 views
ADD COMMENTlink modified 2.4 years ago by h.mon31k • written 2.4 years ago by kspata70
3
gravatar for h.mon
2.4 years ago by
h.mon31k
Brazil
h.mon31k wrote:

Blast is a local alignment tool, which means it won't do what you want. You can play with gap and indel opening and extension parameters and try your luck. References:

Unable to blast entire coding sequence of gene

Force blast to align the whole query read

ADD COMMENTlink written 2.4 years ago by h.mon31k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1897 users visited in the last hour