Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	788,031	T→C	K269K (AAA→AAG)	galM ←	aldose 1‑epimerase; type‑1 mutarotase; galactose mutarotase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	788,031	0	T	C	100.0%	163.2 / NA	44	K269K (AAA→AAG)	galM	aldose 1‑epimerase; type‑1 mutarotase; galactose mutarotase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/24);  total (20/24)

GTCCGGTAGAAACTCGCTTTCCAGAGCCAGCCCTTGCCAGTCGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATC  >  NC_000913/787893‑788177                                                                                                                                           |                                                                                                                                                   gTCCGGTAGAAACTCGCTTTCCAGAGCCAGCCCTTGCCAGTCGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATctgc                                                                                                                                         <  1:193891/150‑1 (MQ=255) gTCCGGTAGAAACTCGCTTTCCAGAGCCAGCCCTTGCCAGTCGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATctg                                                                                                                                          <  2:838825/149‑1 (MQ=255) gTCCGGTAGAAACTCGCTTTCCAGAGCCAGCCCTTGCCAGTCGGCGTAAGGTTCGGTTCCACGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATctgc                                                                                                                                         <  1:285786/150‑1 (MQ=255)      gTAGAAACTCGCTTTCCAGAGCCAGCCCTTGCCAGTCGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGAcc                                                                                                                                    >  1:722306/1‑150 (MQ=255)      gTAGAAACTCGCTTTCCAGAGCCAGCCCTTGCCAGTCGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCa                                                                                                                                   <  2:314541/151‑1 (MQ=255)       tAGAAACTCGCTTTCCAGAGCCAGCCCTTGCCAGTCGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGAcc                                                                                                                                    >  1:679076/1‑149 (MQ=255)                gCTTTCCAGAGCCAGCCCTTGCCAGTCGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGcgcc                                                                                                                         >  1:6722/1‑151 (MQ=255)                     ccAGAGCCAGCCCTTGCCAGTCGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCAc                                                                                                                    <  1:404671/151‑1 (MQ=255)                     ccAGAGCCAGCCCTTGCCAGTCGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCAc                                                                                                                    >  1:815707/1‑151 (MQ=255)                         agccagccCTTGCCAGTCGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACttt                                                                                                                 >  1:201008/1‑150 (MQ=255)                                ccTTGCCAGTCGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCAt                                                                                                         <  1:63057/151‑1 (MQ=255)                                    gCCAGTCGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGCCCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGc                                                                                                      <  2:120723/150‑1 (MQ=255)                                    gCCAGTCGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGcc                                                                                                     >  1:541812/1‑151 (MQ=255)                                            gCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCt                                                                                             <  2:242668/151‑1 (MQ=255)                                             cGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTg                                                                                            >  2:355783/1‑151 (MQ=255)                                                 aGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAAc                                                                                        <  1:569626/151‑1 (MQ=255)                                                  ggttcggttcCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACa                                                                                       <  1:577164/151‑1 (MQ=255)                                                      cggttcCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGa                                                                                    <  2:138031/150‑1 (MQ=255)                                                           ccccGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCg                                                                               <  1:576652/150‑1 (MQ=255)                                                              cgcgACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGAt                                                                           >  2:41075/1‑151 (MQ=255)                                                               gcgACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATc                                                                          <  2:51494/151‑1 (MQ=255)                                                                 gACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGAAATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCATGTAACAAGAATGCGTGATCGt                                                                        >  1:721508/1‑151 (MQ=255)                                                                    ggTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAAc                                                                     >  1:38450/1‑151 (MQ=255)                                                                    ggTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAAc                                                                     >  1:243386/1‑151 (MQ=255)                                                                         gccgccGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTc                                                                >  2:721640/1‑151 (MQ=255)                                                                          ccgccaAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTATAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGGCATGCGCCGCCGCTTTCTTGCCATCGCCTTTGGCCTGTAACGAGAATGCGTGATCGTAAccg                                                                   <  2:621957/147‑2 (MQ=255)                                                                               gAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTg                                                          <  2:878433/151‑1 (MQ=255)                                                                                  gAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTgcgc                                                       <  1:167159/151‑1 (MQ=255)                                                                                         ccGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGAtcgt                                                <  2:842448/151‑1 (MQ=255)                                                                                          cGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGAtcgtc                                               >  1:579189/1‑151 (MQ=255)                                                                                             agtagtATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGAtcgtcgtc                                            >  1:714430/1‑151 (MQ=255)                                                                                             aGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGAtcgtcgtc                                            <  1:288775/151‑1 (MQ=255)                                                                                                  aaTTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCa                                        >  1:671826/1‑150 (MQ=255)                                                                                                       cAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAAc                                  >  1:514056/1‑151 (MQ=255)                                                                                                       cAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTATCACTTTGCGCTGATCGTCGTCGGCAAGAAAc                                  <  1:116011/151‑1 (MQ=255)                                                                                                        agagCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAAc                                  <  2:638202/150‑1 (MQ=255)                                                                                                            ccGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACt                             <  2:180688/151‑1 (MQ=255)                                                                                                            ccGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACt                             >  1:738866/1‑151 (MQ=255)                                                                                                             cGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTg                            <  1:875897/151‑1 (MQ=255)                                                                                                                         gtAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGAtttt                 <  2:61270/150‑1 (MQ=255)                                                                                                                            gACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGgcg             <  2:150347/151‑1 (MQ=255)                                                                                                                               cTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGACACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGgcgc            >  2:354195/1‑149 (MQ=255)                                                                                                                                   aGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCgg       >  1:808436/1‑150 (MQ=255)                                                                                                                                        cAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATc  >  2:469405/1‑150 (MQ=255)                                                                                                                                           |                                                                                                                                                   GTCCGGTAGAAACTCGCTTTCCAGAGCCAGCCCTTGCCAGTCGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATC  >  NC_000913/787893‑788177

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