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_000913788,0310TC97.8% 152.0 / ‑3.7 45K269K (AAA→AAGgalMaldose 1‑epimerase; type‑1 mutarotase; galactose mutarotase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (15/29);  total (16/29)
Fisher's exact test for biased strand distribution p-value = 3.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

CGGGCTGTCCGGTAGAAACTCGCTTTCCAGAGCCAGCCCTTGCCAGTCGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGG  >  NC_000913/787887‑788172
                                                                                                                                                |                                                                                                                                             
cGGGCTGTCCGGTAGAAACTCGCTTTCCAGAGCCAGCCCTTGCCAGTCGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCa                                                                                                                                         <  2:126520/151‑1 (MQ=255)
    cTGTCCGGTAGAAACTCGCTTTCCAGAGCCAGCCCTTGCCAGTCGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATctg                                                                                                                                     >  1:426466/1‑151 (MQ=255)
         cGGTAGAAACTCGCTTTCCAGAGCCAGCCCTGGCCAGTCGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCAGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGAc                                                                                                                                <  2:653535/151‑1 (MQ=255)
                    cGCTTTCCAGAGCCAGCCCTTGCCAGTCGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGCGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATgcgc                                                                                                                     >  1:565383/1‑151 (MQ=255)
                         tCCAGAGCCAGCCCTTGCCAGTCGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCgccgcc                                                                                                                 <  2:283546/150‑1 (MQ=255)
                          ccAGAGCCAGCCCTTGCCAGTCGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCa                                                                                                                <  2:322576/150‑1 (MQ=255)
                           cAGAGCCAGCCCTTGCCAGTCGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGCCCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACt                                                                                                              <  1:774440/151‑1 (MQ=255)
                             gagCCAGCCCTTGCCAGTCGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACttt                                                                                                            <  1:591704/151‑1 (MQ=255)
                                ccagccCTTGCCAGTCGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCtt                                                                                                         <  2:103634/151‑1 (MQ=255)
                                    cccTTGCCAGTCGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCa                                                                                                     <  1:220368/151‑1 (MQ=255)
                                    cccTTGCCAGTCGGCGTAAGGGTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGCCCAGACATGCGCCGCCACTTTTTTGCCa                                                                                                     <  1:18727/151‑1 (MQ=255)
                                     ccTTGCCAGTCGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCAt                                                                                                    >  2:133429/1‑151 (MQ=255)
                                        tGCCAGTCGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCg                                                                                                  <  1:141400/150‑1 (MQ=255)
                                            aGTCGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCttt                                                                                             >  2:716181/1‑151 (MQ=255)
                                               cGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTgg                                                                                           <  2:155366/150‑1 (MQ=255)
                                               cGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGc                                                                                          <  2:424889/151‑1 (MQ=255)
                                               cGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGc                                                                                          <  2:16596/151‑1 (MQ=255)
                                                           cggttcCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGaa                                                                              <  2:277945/151‑1 (MQ=255)
                                                               tcCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGc                                                                           <  2:447502/150‑1 (MQ=255)
                                                                 cccGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTg                                                                        >  2:144871/1‑151 (MQ=255)
                                                                  ccGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGa                                                                       <  1:321519/151‑1 (MQ=255)
                                                                    gcgACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATc                                                                     <  1:240798/151‑1 (MQ=255)
                                                                            gtgCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCt                                                              <  1:465413/150‑1 (MQ=255)
                                                                                 gccgAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACt                                                        >  2:137283/1‑151 (MQ=255)
                                                                                   cgAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACttt                                                      >  1:719963/1‑151 (MQ=255)
                                                                                      ggAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTgcg                                                   >  1:465967/1‑151 (MQ=255)
                                                                                           ttGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGAt                                              >  1:139334/1‑151 (MQ=255)
                                                                                             gCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGAtcg                                            <  2:47301/151‑1 (MQ=255)
                                                                                             gCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGAtcg                                            <  2:220211/151‑1 (MQ=255)
                                                                                             gCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGAtc                                             <  1:59504/150‑1 (MQ=255)
                                                                                                       aaTTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGc                                    <  2:394315/149‑1 (MQ=255)
                                                                                                          tGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGa                                >  2:206847/1‑150 (MQ=255)
                                                                                                           gCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGaaa                              <  2:752591/151‑1 (MQ=255)
                                                                                                             agagCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACt                            <  1:769176/151‑1 (MQ=255)
                                                                                                             agagCCGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACt                            >  2:499932/1‑151 (MQ=255)
                                                                                                                 ccGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACt                        >  2:647792/1‑151 (MQ=255)
                                                                                                                  cGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTg                       >  2:292191/1‑151 (MQ=255)
                                                                                                                  cGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTg                       <  2:342342/151‑1 (MQ=255)
                                                                                                                  cGGAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTg                       <  2:654018/151‑1 (MQ=255)
                                                                                                                    gAGCGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGc                     <  1:661572/151‑1 (MQ=255)
                                                                                                                       cGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCgat                  <  1:716107/151‑1 (MQ=255)
                                                                                                                       cGGTGGTGTAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTCCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGACGTCGGCAAGAAACTCACTgg                      >  2:317939/1‑147 (MQ=255)
                                                                                                                               tAGACCTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTgg          >  1:430803/1‑151 (MQ=255)
                                                                                                                                    cTTCAGCTGCAACTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTg     >  1:552989/1‑151 (MQ=255)
                                                                                                                                        aGCTGCAACTTTGCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCgg  <  1:348762/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                             
CGGGCTGTCCGGTAGAAACTCGCTTTCCAGAGCCAGCCCTTGCCAGTCGGCGTAAGGTTCGGTTCCCCGCGACGGTGTGCCGCCGAGGAAGTTGCCGGAGTAGAATTGCAGAGCCGGAGCGGTGGTGTAGACCTTCAGCTGCAATTTTTCATCTGCTGACCAGACATGCGCCGCCACTTTCTTGCCATCGCCTTTGGCCTGTAACAAGAATGCGTGATCGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGG  >  NC_000913/787887‑788172

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: