Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,050,350	C→A	S108S (TCC→TCA)	insB1 →	IS1 transposase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,050,350	0	C	A	100.0%	19.7 / NA	13	S108S (TCC→TCA)	insB1	IS1 transposase B
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (10/3);  total (10/3)

TTATGAGCCTGCTGTCACCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGTATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGC  >  NC_000913/1050283‑1050483                                                                    |                                                                                                                                      ttctGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGTATAACCTGAATCTg                                                               >  1:128047/4‑139 (MQ=11)           tgctgtggccTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGTATAACCTGAATCTGAGGCAGCAc                                                      >  2:1208941/9‑139 (MQ=11)           tgctgtggccTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGTATAACCTGAATCTGAGGCAGCAc                                                      <  2:1288779/131‑1 (MQ=11)            gctgtggccTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCCGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGTATAACCTGAATCTGAGGCAGCAcc                                                     >  2:188568/8‑139 (MQ=11)                 ggccTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCGAGCGATATACGCAGCGAATTGAGCGGTATAACCTGAATCTGAGGCAGCACCTGGCa                                                >  2:1003426/3‑139 (MQ=11)                 ggccTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGTATAACCTGAATCTGAGGCAGCACCTGGCa                                                >  1:399318/3‑139 (MQ=11)                             cgtggtataaaTGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGTATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACgg                                    >  2:405408/11‑139 (MQ=11)                                            ggATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGTATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTtc                     >  2:399222/1‑139 (MQ=18)                                                ggctggcCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGTATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCaa                 >  2:105973/1‑139 (MQ=17)                                                         cTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCGAGCGATATACGCAGCGAATTGAGCGGTATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTgg        <  1:1003426/139‑1 (MQ=18)                                                          tGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGTATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGa       >  1:45053/1‑139 (MQ=17)                                                           gTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGTATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAg      <  1:157428/139‑1 (MQ=18)                                                               gAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGTATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGc  >  2:731037/1‑139 (MQ=17)                                                                    |                                                                                                                                      TTATGAGCCTGCTGTCACCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGTATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGC  >  NC_000913/1050283‑1050483

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑