Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,228,724 A→G T216A (ACT→GCT)  ycgM → putative isomerase/hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,228,7240AG96.2% 177.9 / ‑0.3 52T216A (ACT→GCT) ycgMputative isomerase/hydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (25/25);  minor base T (1/1);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.40e-01

CTGATCGCTTATATGAGCAAGTTTTTTACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAACTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCTGATGAGCGATGTACCTTTCTGG  >  NC_000913/1228592‑1228850
                                                                                                                                    |                                                                                                                              
cTGATCGCTTATATGAGCAAGTTTTTTACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTcgcg                                                                                                                          <  2:1268685/139‑1 (MQ=255)
cTGATCGCTTATATGAGCAAGTTTTTTACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTcgcg                                                                                                                          <  2:1045348/139‑1 (MQ=255)
 tGATCGCTTATATGAGCAAGTTTTTTACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGt                                                                                                                         <  2:780549/139‑1 (MQ=255)
   aTCGCTTATATGAGCAAGTTTTTTACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGttt                                                                                                                       >  2:961935/1‑139 (MQ=255)
          atatGAGCAAGTTTTTTACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTcgcg                                                                                                                          >  2:442961/1‑129 (MQ=255)
          atatGAGCAAGTTTTTTACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTcgcg                                                                                                                          <  1:442961/129‑1 (MQ=255)
               aGCAAGTTTTTTACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACtttt                                                                                                           >  1:667198/1‑139 (MQ=255)
                  aaGTTTTTCACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTCATACTTTTTgc                                                                                                        >  2:1190621/1‑139 (MQ=255)
                     ttttttACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTgccgc                                                                                                     >  1:288006/1‑139 (MQ=255)
                              cTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCTGCCCGTTGCAAAGCGGTTATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAgcg                                                                                            <  1:1190621/139‑1 (MQ=255)
                                 aaGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAgcggcg                                                                                         >  1:120730/1‑139 (MQ=255)
                                   ggCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCggcggc                                                                                       <  2:120730/139‑1 (MQ=255)
                                      cGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAgcg                                                                                            <  2:1634092/131‑1 (MQ=255)
                                      cGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAgcg                                                                                            >  1:1634092/1‑131 (MQ=255)
                                           aCGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGt                                                                                                                         >  2:1753119/1‑97 (MQ=255)
                                           aCGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGt                                                                                                                         <  1:1753119/97‑1 (MQ=255)
                                               tgtgTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATc                                                                           >  1:241425/1‑139 (MQ=255)
                                                gtgtTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCg                                                                          <  1:1354253/139‑1 (MQ=255)
                                                   ttGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCtg                                                                       >  2:1454494/1‑139 (MQ=255)
                                                    tGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCtgt                                                                      >  2:1291082/1‑139 (MQ=255)
                                                       cAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGcc                                                                   <  2:288006/139‑1 (MQ=255)
                                                        aGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCa                                                                  >  2:176163/1‑139 (MQ=255)
                                                           cACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGt                                                                                                                    <  2:139017/86‑1 (MQ=255)
                                                           cACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGt                                                                                                                    >  1:139017/1‑86 (MQ=255)
                                                           cACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGAc                                                               <  1:1291082/139‑1 (MQ=255)
                                                             cGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTg                                                             >  2:1032434/1‑139 (MQ=255)
                                                                             cccGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGttt                                             <  2:241425/139‑1 (MQ=255)
                                                                                 ttGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGAc                                                               <  1:882661/117‑1 (MQ=255)
                                                                                 ttGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGAc                                                               >  2:882661/1‑117 (MQ=255)
                                                                                           gTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTaa                                          <  1:453358/128‑1 (MQ=255)
                                                                                           gTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTaa                                          >  2:453358/1‑128 (MQ=255)
                                                                                            tgatgaGCTGACAGTCACTTTCGATGGGCATTCTTTGACATCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCtgtg                                                                     <  1:802392/100‑1 (MQ=255)
                                                                                            tgatgaGCTGACAGTCACTTTCGATGGGCATTCTTTGACATCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCtgtg                                                                     >  2:802392/1‑100 (MQ=255)
                                                                                            tgatgaGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGtt                                                                                                                        >  2:835595/1‑49 (MQ=255)
                                                                                            tgatgaGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGtt                                                                                                                        <  1:835595/49‑1 (MQ=255)
                                                                                            tgatgaGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGttt                                             >  1:364923/1‑124 (MQ=255)
                                                                                            tgatgaGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGttt                                             <  2:364923/124‑1 (MQ=255)
                                                                                            tgatgaGCTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGa                              >  2:990246/1‑139 (MQ=255)
                                                                                                   cTGACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCtg                       >  1:1015125/1‑139 (MQ=255)
                                                                                                     gACAGTCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCtgat                     <  1:990246/139‑1 (MQ=255)
                                                                                                          tCACTTTCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCTGATGAGCg                <  2:122833/139‑1 (MQ=255)
                                                                                                                tCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCtgtg                                                                     >  2:1644059/1‑80 (MQ=255)
                                                                                                                tCGATGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCtgtg                                                                     <  1:1644059/80‑1 (MQ=255)
                                                                                                                    tGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAacac                          <  2:1384647/119‑1 (MQ=255)
                                                                                                                    tGGGCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAacac                          >  1:1384647/1‑119 (MQ=255)
                                                                                                                      ggCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGAc                                                               >  2:410514/1‑80 (MQ=255)
                                                                                                                      ggCATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGAc                                                               <  1:410514/80‑1 (MQ=255)
                                                                                                                        cATTCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCTGATGAGCGATGTACCTTTCTgg  <  2:1015125/139‑1 (MQ=255)
                                                                                                                          ttCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACg                                        <  2:1748831/99‑1 (MQ=255)
                                                                                                                          ttCTTTGACAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACg                                        >  1:1748831/1‑99 (MQ=255)
                                                                                                                                 aCAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCgg                               <  2:253561/101‑1 (MQ=255)
                                                                                                                                 aCAGCTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCgg                               >  1:253561/1‑101 (MQ=255)
                                                                                                                                    |                                                                                                                              
CTGATCGCTTATATGAGCAAGTTTTTTACCCTCAAGGCCGGTGACGTTGTGTTGACAGGCACGCCTGATGGCGTCGGCCCGTTGCAAAGCGGTGATGAGCTGACAGTCACTTTCGATGGGCATTCTTTGACAACTCGCGTTTTGTAATACTTTTTGCCGCCTGAAAGCGGCGGCAAAACTTGCATCGCTGTGCCAGACTGGTTATAAGGTGCGTTTTAACGTAATGGCGGAACACCTGATGAGCGATGTACCTTTCTGG  >  NC_000913/1228592‑1228850

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: ‑