Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	789,464	G→A	I172I (ATC→ATT)	galK ←	galactokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	789,464	0	G	A	100.0%	168.7 / NA	53	I172I (ATC→ATT)	galK	galactokinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/27);  total (26/27)

GGGTACGTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAA  >  NC_000913/789328‑789593                                                                                                                                         |                                                                                                                                  gggTACGTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAAt                                                                                                                                 <  2:1103239/139‑1 (MQ=255) gggTACGTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAAt                                                                                                                                 <  1:1304652/139‑1 (MQ=255)     aCGTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCg                                                                                                                             <  1:677623/139‑1 (MQ=255)     aCGTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCg                                                                                                                             <  1:366317/139‑1 (MQ=255)     aCGTTTGAAGTTACTGTTGAAGATTACGACAGCCAGACCGTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAAGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCg                                                                                                                             <  2:1140250/139‑1 (MQ=255)             aGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACa                                                                                                                     <  2:507291/139‑1 (MQ=255)                 aCTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCt                                                                                                                 >  2:247992/1‑139 (MQ=255)                         tgatgaCGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAAc                                                                                                           >  2:48409/1‑137 (MQ=255)                         tgatgaCGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAAc                                                                                                           <  1:48409/137‑1 (MQ=255)                            tgaCGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCg                                                                                                                             >  1:1004051/1‑116 (MQ=255)                            tgaCGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCg                                                                                                                             <  2:1004051/116‑1 (MQ=255)                              acgacAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGtttt                                                                                                    >  2:553412/1‑139 (MQ=255)                               cgacAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTc                                                                                                   >  2:978446/1‑139 (MQ=255)                                          ccTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGAc                                                                                        >  1:132505/1‑139 (MQ=255)                                             ttGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGt                                                                                     >  1:1049070/1‑139 (MQ=255)                                               gggCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAAc                                                                                                           <  2:638336/115‑1 (MQ=255)                                               gggCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAAc                                                                                                           >  1:638336/1‑115 (MQ=255)                                                     ggAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATcc                                                                                                                               <  1:938994/89‑1 (MQ=255)                                                     ggAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGcgcg                                                                             >  2:565735/1‑139 (MQ=255)                                                     ggAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGcgcg                                                                             >  2:820990/1‑139 (MQ=255)                                                     ggAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGcgcg                                                                             >  2:26103/1‑139 (MQ=255)                                                     ggAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCAGATCCATAATcc                                                                                                                               >  2:938994/1‑89 (MQ=255)                                                         aCTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGAttt                                                                         <  1:565735/139‑1 (MQ=255)                                                             cTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTgc                                                                     >  2:511566/1‑139 (MQ=255)                                                              tttGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTgcg                                                                    <  1:511566/139‑1 (MQ=255)                                                              tttGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTgcg                                                                    <  1:463820/139‑1 (MQ=255)                                                                tGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGcc                                                                  >  2:351833/1‑139 (MQ=255)                                                                 ggTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGttt                                                                                                     <  1:748909/103‑1 (MQ=255)                                                                 ggTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGttt                                                                                                     >  2:748909/1‑103 (MQ=255)                                                                          tGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCg                                                                                                                             <  2:1134448/70‑1 (MQ=255)                                                                          tGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCg                                                                                                                             >  1:1134448/1‑70 (MQ=255)                                                                          tGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGc                                                        <  2:132505/139‑1 (MQ=255)                                                                              cGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAgat                                                    >  2:1062320/1‑139 (MQ=255)                                                                                 cAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAgatgat                                                 <  1:1062320/139‑1 (MQ=255)                                                                                     cGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAg                                             <  1:688688/139‑1 (MQ=255)                                                                                        tCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAgctg                                          <  1:978446/139‑1 (MQ=255)                                                                                           gCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGt                                                                                                       <  1:1266107/75‑1 (MQ=255)                                                                                           gCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGt                                                                                                       >  2:1266107/1‑75 (MQ=255)                                                                                              aGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTcc                                                              <  2:999880/113‑1 (MQ=255)                                                                                              aGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTcc                                                              >  1:999880/1‑113 (MQ=255)                                                                                               ggCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAgc                                            <  2:255721/130‑1 (MQ=255)                                                                                               ggCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAgc                                            >  1:255721/1‑130 (MQ=255)                                                                                               ggCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAAt                                   >  2:420761/1‑139 (MQ=255)                                                                                                     aTCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGtt                             >  1:1117952/1‑139 (MQ=255)                                                                                                         ttcttGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGa                         <  2:589058/139‑1 (MQ=255)                                                                                                         ttcttGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGa                         >  1:1167487/1‑139 (MQ=255)                                                                                                              gCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAg                                                            <  2:1231156/99‑1 (MQ=255)                                                                                                              gCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAg                                                            >  1:1231156/1‑99 (MQ=255)                                                                                                              gCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACcgcg                    <  1:351833/139‑1 (MQ=255)                                                                                                               ccGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGa                   >  1:105237/1‑139 (MQ=255)                                                                                                                    cgcgGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTcc              >  2:370290/1‑139 (MQ=255)                                                                                                                      cgGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAg            <  1:1196041/139‑1 (MQ=255)                                                                                                                                cTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGgaa  <  1:247992/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                  GGGTACGTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAA  >  NC_000913/789328‑789593

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