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%	183.1 / NA	57	I172I (ATC→ATT)	galK	galactokinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/26);  total (31/26)

GTACGTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACT  >  NC_000913/789330‑789598                                                                                                                                       |                                                                                                                                       gTACGTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATcc                                                                                                                                    <  1:381242/139‑1 (MQ=255) gTACGTTTGAAGTTACTGTTGATGATGACGACAACCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATcc                                                                                                                                    <  2:590479/139‑1 (MQ=255)  tACGTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATccc                                                                                                                                   >  1:433676/1‑139 (MQ=255)  tACGTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATcc                                                                                                                                    >  1:1319899/1‑138 (MQ=255)  tACGTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATcc                                                                                                                                    <  2:1319899/138‑1 (MQ=255)       ttGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGt                                                                                                                              >  1:62878/1‑139 (MQ=255)                cTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCACAGTGAGCGGCAATCGATCAGAAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTa                                                                                                                     >  1:1430649/1‑139 (MQ=255)                      atgatgaCGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACa                                                                                                                          >  1:239337/1‑128 (MQ=255)                      atgatgaCGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACa                                                                                                                          <  2:239337/128‑1 (MQ=255)                         atgaCGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGt                                                                                                            <  1:302774/139‑1 (MQ=255)                                           ttGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGt                                                                                          >  2:773096/1‑139 (MQ=255)                                              ggCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTaa                                                                                       >  1:1332449/1‑139 (MQ=255)                                                 aTGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAgcg                                                                                    <  1:222354/139‑1 (MQ=255)                                                  tGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAgcgc                                                                                   <  2:433676/139‑1 (MQ=255)                                                   ggAAACTGCTTTGGTCCCCAGTGAGCGGCTATCGATCAGCTAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTAACCGTTAAGcgcg                                                                                  >  2:760871/1‑139 (MQ=255)                                                   ggAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCg                                                                                                                                  <  2:516294/91‑1 (MQ=255)                                                   ggAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCg                                                                                                                                  >  1:516294/1‑91 (MQ=255)                                                              tGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGcc                                                                       >  1:663036/1‑139 (MQ=255)                                                              tGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGAAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGcc                                                                       >  2:785281/1‑139 (MQ=255)                                                               ggTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAAt                                                                                                                                      <  1:349026/75‑1 (MQ=255)                                                               ggTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCg                                                                      >  2:626829/1‑139 (MQ=255)                                                                 tCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTc                                                                    <  2:1332449/139‑1 (MQ=255)                                                                     cAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTgcgc                                                                        >  2:189830/1‑131 (MQ=255)                                                                     cAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTgcgc                                                                        <  1:189830/131‑1 (MQ=255)                                                                            cGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTAc                                                                                                                           <  1:854255/73‑1 (MQ=255)                                                                            cGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTAc                                                                                                                           >  2:854255/1‑73 (MQ=255)                                                                                  tCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGt                                                                                                            <  1:1383896/82‑1 (MQ=255)                                                                                  tCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGt                                                                                                            >  2:1383896/1‑82 (MQ=255)                                                                                  tCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATaaa                                                   >  2:54665/1‑139 (MQ=255)                                                                                         gCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGctgctg                                            <  1:54665/139‑1 (MQ=255)                                                                                          cAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCtgctgc                                           >  2:59746/1‑139 (MQ=255)                                                                                             ggCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAAt                                        >  2:1169505/1‑139 (MQ=255)                                                                                              gCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAAt                                        >  1:1275833/1‑138 (MQ=255)                                                                                              gCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAAt                                        <  2:1275833/138‑1 (MQ=255)                                                                                                aTGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACg                                     >  1:1408614/1‑139 (MQ=255)                                                                                                    tctttcttGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTgcgc                                                                        <  1:1262349/100‑1 (MQ=255)                                                                                                    tctttcttGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTgcgc                                                                        >  2:1262349/1‑100 (MQ=255)                                                                                                    tctttcttGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCCGCGGCAGATGATAAAGCTGCTGCAATACGGTTc                                 >  2:1150974/1‑139 (MQ=255)                                                                                                    tctttcttGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTc                                 >  1:1169987/1‑139 (MQ=255)                                                                                                        tcttGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTAc                                                                                                                           <  1:475035/45‑1 (MQ=255)                                                                                                        tcttGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTAc                                                                                                                           >  2:475035/1‑45 (MQ=255)                                                                                                        tcttGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGt                                                                                                            >  1:465693/1‑60 (MQ=255)                                                                                                        tcttGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGt                                                                                                            <  2:465693/60‑1 (MQ=255)                                                                                                         cttGCCGAGCGCGGAAATTAGCTGATCCATAATCCCg                                                                                                                                  >  1:334474/1‑37 (MQ=255)                                                                                                         cttGCCGAGCGCGGAAATTAGCTGATCCATAATCCCg                                                                                                                                  <  2:334474/37‑1 (MQ=255)                                                                                                         cttGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAAc                                                                                                                <  2:58491/55‑1 (MQ=255)                                                                                                         cttGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAAc                                                                                                                >  1:58491/1‑55 (MQ=255)                                                                                                            gCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACcgcg                         <  1:785281/139‑1 (MQ=255)                                                                                                               gAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACt                      <  2:1408614/139‑1 (MQ=255)                                                                                                                  cgcgGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTcc                   >  1:849727/1‑139 (MQ=255)                                                                                                                    cgGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCTCGATTTGTGCGCCGTCCATCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTACCAg                 >  2:433080/1‑139 (MQ=255)                                                                                                                        aaTTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACcgcg                         <  2:95083/127‑1 (MQ=255)                                                                                                                        aaTTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACcgcg                         >  1:95083/1‑127 (MQ=255)                                                                                                                                   ccATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGAttt                                                                              <  2:181116/63‑1 (MQ=255)                                                                                                                                   ccATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGAttt                                                                              >  1:181116/1‑63 (MQ=255)                                                                                                                                   ccATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACt  <  2:1169987/139‑1 (MQ=255)                                                                                                                                   ccATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACt  <  1:869822/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                       GTACGTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACT  >  NC_000913/789330‑789598

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