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

GTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACT  >  NC_000913/789334‑789598                                                                                                                                   |                                                                                                                                        gTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGTGAGGCAGTCGATCAGCAAGGCATGATCTTTCTTGCCGATCGCGGAAATTATCTGATCCATAATCCCGCa                                                                                                                                 >  2:587989/1‑139 (MQ=255) gTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCa                                                                                                                                 >  1:177705/1‑139 (MQ=255)  tttGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGTGGAAATTAGCTGATCCATAATCCCGCAg                                                                                                                                <  2:552719/139‑1 (MQ=255)    tGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGtt                                                                                                                              >  1:804189/1‑139 (MQ=255)           aCTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCt                                                                                                                       >  2:187558/1‑139 (MQ=255)                 gatgatgaCGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAAc                                                                                                                 <  2:804189/139‑1 (MQ=255)                   tgatgaCGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTg                                                                                                               >  1:475965/1‑139 (MQ=255)                          gacAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCt                                                                                                        <  2:429883/139‑1 (MQ=255)                            cAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGc                                                                                                      >  1:1198657/1‑139 (MQ=255)                                  cacCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTg                                                                                                <  1:1136269/139‑1 (MQ=255)                                       ttGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTACTTCCTGACCGt                                                                                           >  1:537978/1‑139 (MQ=255)                                               ggAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTcc                                                                                                  >  1:538600/1‑124 (MQ=255)                                               ggAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTcc                                                                                                  <  2:538600/124‑1 (MQ=255)                                                      gCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTtgtg                                                                           <  1:289100/139‑1 (MQ=255)                                                        tttGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTgcg                                                                         <  2:475965/139‑1 (MQ=255)                                                        tttGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTACTTCCTGACCGTTAAGCGCGA‑TTTGTgcg                                                                         <  2:537978/139‑1 (MQ=255)                                                        tttGGTCCCCAGTGAGCGGCAATCGATCAGCAAGACATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTgcg                                                                         <  2:1536397/139‑1 (MQ=255)                                                         ttGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACaaa                                                                                                                  <  2:1070605/98‑1 (MQ=255)                                                         ttGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACaaa                                                                                                                  >  1:1070605/1‑98 (MQ=255)                                                         ttGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTgcgc                                                                        <  2:177705/139‑1 (MQ=255)                                                           ggTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTGCGCCg                                                                      >  1:1606887/1‑139 (MQ=255)                                                                 cAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTGCGCCGTCCAGc                                                                >  2:597539/1‑139 (MQ=255)                                                                    tGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTGCGCCGTCCAGCGGc                                                             >  1:817153/1‑139 (MQ=255)                                                                     gAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTGCGCCGTCCAGCGGCa                                                            <  2:817153/139‑1 (MQ=255)                                                                        cGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCATTAAGCGCGATTTTGTGCGCCGTCCAGCGGCAga                                                          >  1:1229710/1‑139 (MQ=255)                                                                         ggCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTGCGCCGTCCAGCGGCa                                                            <  2:1085195/135‑1 (MQ=255)                                                                         ggCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTGCGCCGTCCAGCGGCa                                                            >  1:1085195/1‑135 (MQ=255)                                                                             aTCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTGCGCCGTCCAGCGGCAGATGATaa                                                    >  1:1298631/1‑139 (MQ=255)                                                                                    aGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTgcg                                                                         <  1:1368208/111‑1 (MQ=255)                                                                                    aGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTgcg                                                                         >  2:1368208/1‑111 (MQ=255)                                                                                            aTGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTgcgc                                                                        >  2:83772/1‑104 (MQ=255)                                                                                            aTGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTgcgc                                                                        <  1:83772/104‑1 (MQ=255)                                                                                              gATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATa                                       >  2:1112905/1‑135 (MQ=255)                                                                                              gATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATa                                       <  1:1112905/135‑1 (MQ=255)                                                                                              gATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGt                                   >  1:839613/1‑139 (MQ=255)                                                                                                tctttcttGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTGCGCCGTCCAGCGGCAGATGATAAAgctg                                               <  1:1230858/125‑1 (MQ=255)                                                                                                tctttcttGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTGCGCCGTCCAGCGGCAGATGATAAAgctg                                               >  2:1230858/1‑125 (MQ=255)                                                                                                    tcttGCCGAGCGCGGAAATTAGCTGATCCATAATccc                                                                                                                                    >  1:835453/1‑37 (MQ=255)                                                                                                    tcttGCCGAGCGCGGAAATTAGCTGATCCATAATccc                                                                                                                                    <  2:835453/37‑1 (MQ=255)                                                                                                     cttGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTGCGCCGTc                                                                    <  2:1447434/99‑1 (MQ=255)                                                                                                     cttGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTGCGCCGTc                                                                    >  1:1447434/1‑99 (MQ=255)                                                                                                        gCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACcgcg                         <  2:1606887/139‑1 (MQ=255)                                                                                                        gCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACcgcg                         <  2:1298631/139‑1 (MQ=255)                                                                                                        gCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACcgcg                         >  1:1447327/1‑139 (MQ=255)                                                                                                        gCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACcgcg                         <  1:597539/139‑1 (MQ=255)                                                                                                          cGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGAc                       >  2:742724/1‑139 (MQ=255)                                                                                                           gAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACt                      >  1:743649/1‑139 (MQ=255)                                                                                                            aGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTgcg                                                                         >  2:558895/1‑87 (MQ=255)                                                                                                            aGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTgcg                                                                         <  1:558895/87‑1 (MQ=255)                                                                                                                     aTTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAg            <  1:609303/139‑1 (MQ=255)                                                                                                                               ccATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACt  <  2:743649/139‑1 (MQ=255)                                                                                                                               ccATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACt  <  2:1198657/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                                        GTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGA‑TTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACT  >  NC_000913/789334‑789598

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