GEvo Blastn Bug

Bug Description

When using blastn in GEvo, a large HSP appear in the middle of the region when the query sequence length is changed by 1 nucleotide.

Visualization

Difference between analyses is that addition of 1 nucleotide to the top panel of the top analysis

Report of "disappearing" HSP in blast file

Top portion of blast report which contains the large HSP (problem HSP starts at "Query 129"):

BLASTN 2.2.24+

Query=  AT1G75520 
Length=478

Subject=  Bra008203 
Length=17361


 Score = 68.0 bits (35),  Expect = 3e-14
 Identities = 76/94 (80%), Gaps = 4/94 (4%)
 Strand=Plus/Plus

Query  1      CTAGGTTTCGTGTTCCACTGATCAAAGATTTGAAAAAAAACATATACTTAGTAAACTTCA  60
              |||||||| ||||||||||||||||||| ||    || |||||||    |  |||||| |
Sbjct  9943   CTAGGTTTTGTGTTCCACTGATCAAAGAGTT----AAGAACATATTTCAATAAAACTTTA  9998

Query  61     AGCAATTTTTATATTACCCAATTGAATTTCTCCA  94
              |  |||||||||| |||||| |||||||||||||
Sbjct  9999   AATAATTTTTATACTACCCAGTTGAATTTCTCCA  10032


 Score = 54.5 bits (28),  Expect = 3e-10
 Identities = 32/34 (94%), Gaps = 0/34 (0%)
 Strand=Plus/Plus

Query  442    TTGGCGTAGATGAATGTAAACGGATGGTAATATA  475
              |||||||| |||||||||||||||| ||||||||
Sbjct  10342  TTGGCGTATATGAATGTAAACGGATAGTAATATA  10375


 Score = 50.7 bits (26),  Expect = 4e-09
 Identities = 77/95 (81%), Gaps = 5/95 (5%)
 Strand=Plus/Plus

Query  129    ATATATATATACCCAACAACTGAGAAAAGATGGAAAAAGTTTAGTTAAAAACTGGTCCTG  188
              ||||||||||| |||||| ||  |||||||||||| || || |||||||||     | | 
Sbjct  10072  ATATATATATAACCAACATCTCGGAAAAGATGGAACAAATT-AGTTAAAAAAA---CATT  10127

Query  189    GGCGGCTTTAAATTATATTTATGCACTTAAATTTA  223
              ||||||||||||| ||||| ||| | |||||||||
Sbjct  10128  GGCGGCTTTAAATCATATTCATGTA-TTAAATTTA  10161


 Score = 48.8 bits (25),  Expect = 2e-08
 Identities = 47/58 (81%), Gaps = 0/58 (0%)
 Strand=Plus/Plus

Query  316    GGTTTTTACTTAGATAATATCGTGTCATTCCATCTAGATTCAACCCCTGTCTACAATA  373
              |||||| ||| |||| |||| ||  || | |  ||||||||||||||| |||||||||
Sbjct  10207  GGTTTTCACTGAGATGATATTGTTCCAGTTCCACTAGATTCAACCCCTCTCTACAATA  10264


 Score = 25.7 bits (13),  Expect = 0.15
 Identities = 13/13 (100%), Gaps = 0/13 (0%)
 Strand=Plus/Plus

Query  127   TAATATATATATA  139
             |||||||||||||
Sbjct  7327  TAATATATATATA  7339

Top portion of blast report which DOES NOT contains the large HSP :

BLASTN 2.2.24+


Query=  AT1G75520 
Length=479

Subject=  Bra008203 
Length=17361


 Score = 69.9 bits (36),  Expect = 7e-15
 Identities = 77/95 (81%), Gaps = 4/95 (4%)
 Strand=Plus/Plus

Query  1      ACTAGGTTTCGTGTTCCACTGATCAAAGATTTGAAAAAAAACATATACTTAGTAAACTTC  60
              ||||||||| ||||||||||||||||||| ||    || |||||||    |  |||||| 
Sbjct  9942   ACTAGGTTTTGTGTTCCACTGATCAAAGAGTT----AAGAACATATTTCAATAAAACTTT  9997

Query  61     AAGCAATTTTTATATTACCCAATTGAATTTCTCCA  95
              ||  |||||||||| |||||| |||||||||||||
Sbjct  9998   AAATAATTTTTATACTACCCAGTTGAATTTCTCCA  10032


 Score = 54.5 bits (28),  Expect = 3e-10
 Identities = 32/34 (94%), Gaps = 0/34 (0%)
 Strand=Plus/Plus

Query  443    TTGGCGTAGATGAATGTAAACGGATGGTAATATA  476
              |||||||| |||||||||||||||| ||||||||
Sbjct  10342  TTGGCGTATATGAATGTAAACGGATAGTAATATA  10375


 Score = 48.8 bits (25),  Expect = 2e-08
 Identities = 47/58 (81%), Gaps = 0/58 (0%)
 Strand=Plus/Plus

Query  317    GGTTTTTACTTAGATAATATCGTGTCATTCCATCTAGATTCAACCCCTGTCTACAATA  374
              |||||| ||| |||| |||| ||  || | |  ||||||||||||||| |||||||||
Sbjct  10207  GGTTTTCACTGAGATGATATTGTTCCAGTTCCACTAGATTCAACCCCTCTCTACAATA  10264


 Score = 25.7 bits (13),  Expect = 0.15
 Identities = 13/13 (100%), Gaps = 0/13 (0%)
 Strand=Plus/Plus

Query  128   TAATATATATATA  140
             |||||||||||||
Sbjct  7327  TAATATATATATA  7339