Hello, is there a tool for processing the blast output and perform a statistical analysis of the alignment? For example a plot of the position of the mismatches? I cant believe that there is no existing tool. Here is a example of my blast output:
>ref|NC_018914.1| Homo sapiens chromosome 3, alternate assembly CHM1_1.0, whole
genome shotgun sequence
gb|CM001611.1| Homo sapiens cell-line CHM1htert chromosome 3, whole genome shotgun
sequence
Length=197815086
Score = 54.0 bits (27), Expect = 3e-05
Identities = 33/35 (94%), Gaps = 0/35 (0%)
Strand=Plus/Minus
Query 1 CTAGAATGCACACTCCTACCTCCTTTACCAAACGT 35
|||||||||||| |||||| |||||||||||||||
Sbjct 169676496 CTAGAATGCACATTCCTACTTCCTTTACCAAACGT 169676462