Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,311,576	A→G	Y163C (TAT→TGT)	tonB →	Ton complex subunit TonB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,311,576	0	A	G	100.0%	93.5 / NA	28	Y163C (TAT→TGT)	tonB	Ton complex subunit TonB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/12);  total (16/12)

GAGTCGCGTCCGGCATCACCGTTTGAAAATACGGCACCGGCACGCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAA  >  NC_000913/1311440‑1311708                                                                                                                                         |                                                                                                                                     gagTCGCGTCCGGCATCACCGTTTGAAAATACGGCACCGGCACGCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTGTcc                                                                                                                                   <  1:851268/140‑1 (MQ=255)   gTCGCGTCCGGCATCACCGTTTGAAAATACGGCACCGGCACGCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTGTCCgg                                                                                                                                 <  2:14066/140‑1 (MQ=255)    tCGCGTCCGGCATCACCGTTTGAAAATACGGCACCGGCACGCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTGTCCGGc                                                                                                                                >  1:227479/1‑140 (MQ=255)                        tGAAAATACGGCACCGGCACGCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTGTCCGGCACGAGCACAGGCATTGCGCa                                                                                                            <  1:403309/140‑1 (MQ=255)                        tGAAAATACGGCACCGGCACGCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTGTCCGGCACGAGCACAGGCATTGCGCa                                                                                                            <  2:492148/140‑1 (MQ=255)                                cggcaccggcacGCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTGTCCGGCACGAGCACAGGCATTGCGCATTGAAggg                                                                                                    <  2:227479/140‑1 (MQ=255)                                cggcaccggcacGCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTGTCCGGCACGAGCACAGGCATTGCGCATTGAAggg                                                                                                    <  2:834745/140‑1 (MQ=255)                                 ggcaccggcccGCCTGACATCACGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTGTCCGGCACGAGCACAGGCATTGCGCATTGAAGGGc                                                                                                   >  2:758075/1‑140 (MQ=255)                                      cggcacGCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTGTCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGgtt                                                                                              >  1:537185/1‑140 (MQ=255)                                           cGCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTGTCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGgttaaagt                                                                                         >  2:392011/1‑140 (MQ=255)                                             ccTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTGTCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGgttaaagtta                                                                                       >  2:476991/1‑140 (MQ=255)                                                       aGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTGTCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGt                                                                             <  2:656836/140‑1 (MQ=255)                                                        gTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTGTCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGtt                                                                            >  2:656407/1‑140 (MQ=255)                                                           cAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTGTCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACg                                                                         <  1:656407/140‑1 (MQ=255)                                                             gCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTGTCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGcc                                                                       >  2:682951/1‑140 (MQ=255)                                                              cAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTGTCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCg                                                                      >  2:408901/1‑140 (MQ=255)                                                                 cGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGGGTCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGAt                                                                   >  2:708873/1‑140 (MQ=255)                                                                  ggCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGGAATCAGCCGCAGTGTCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATg                                                                  >  1:94046/1‑140 (MQ=255)                                                                       cAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTGTCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTcgc                                                             <  2:537185/140‑1 (MQ=255)                                                                         aCCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTGTCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGt                                                           <  1:476991/140‑1 (MQ=255)                                                                              cAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTGTCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATa                                                      <  1:408901/140‑1 (MQ=255)                                                                                   cGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTGTCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGa                                                        >  2:213688/1‑133 (MQ=255)                                                                                   cGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTGTCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGa                                                        <  1:213688/133‑1 (MQ=255)                                                                                   cGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTGTCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTa                                                 >  2:213047/1‑140 (MQ=255)                                                                                             tgtgGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTGTCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCctct                                       >  1:479070/1‑140 (MQ=255)                                                                                                  cTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTGTCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGcc                                  >  2:670743/1‑140 (MQ=255)                                                                                                             cgcgCATTAAGCCGTAATCAGCCGCAGTGTCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGaa                       >  2:985941/1‑140 (MQ=255)                                                                                                                           tAATCAGCCGCAGTGTCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTg         >  1:973219/1‑140 (MQ=255)                                                                                                                                  ccGCAGTGTCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGaa  >  2:405677/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                     GAGTCGCGTCCGGCATCACCGTTTGAAAATACGGCACCGGCACGCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAA  >  NC_000913/1311440‑1311708

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