Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,953,303 T→C M201T (ATG→ACG)  cmoA → carboxy‑SAM synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,953,3030TC98.0% 171.9 / ‑4.0 51M201T (ATG→ACG) cmoAcarboxy‑SAM synthase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (26/24);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.42e-01

GGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGAC  >  NC_000913/1953158‑1953450
                                                                                                                                                 |                                                                                                                                                   
gggCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTgg                                                                                                                                                <  2:471288/151‑1 (MQ=255)
  gCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGCGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTg                                                                                                                                                 >  1:111021/1‑148 (MQ=255)
    tGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGaaa                                                                                                                                             >  1:295773/1‑150 (MQ=255)
     gAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAAc                                                                                                                                           <  2:766932/151‑1 (MQ=255)
      aaCCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACg                                                                                                                                          >  2:587625/1‑151 (MQ=255)
       aCCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGt                                                                                                                                         >  1:614283/1‑151 (MQ=255)
          ccGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTCAt                                                                                                                                      <  1:796374/151‑1 (MQ=255)
               ggcgCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACAACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTg                                                                                                                                  >  2:117287/1‑150 (MQ=255)
                gcgcgcTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGAc                                                                                                                                >  1:127986/1‑151 (MQ=255)
                  gcgcTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCg                                                                                                                              <  1:542363/151‑1 (MQ=255)
                        gTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCg                                                                                                                        <  2:303250/151‑1 (MQ=255)
                         tGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGt                                                                                                                       >  2:238356/1‑151 (MQ=255)
                          gCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTg                                                                                                                      >  2:267413/1‑151 (MQ=255)
                           cTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTg                                                                                                                      <  1:415424/150‑1 (MQ=255)
                              tCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGaa                                                                                                                   >  2:369486/1‑150 (MQ=255)
                                ggAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAAcc                                                                                                                <  2:147611/151‑1 (MQ=255)
                                   aaaaTTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCAt                                                                                                             >  2:428405/1‑151 (MQ=255)
                                   aaaaTTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCAt                                                                                                             >  1:830125/1‑151 (MQ=255)
                                    aaaTTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATa                                                                                                            <  2:121043/151‑1 (MQ=255)
                                           gTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACg                                                                                                     >  2:825913/1‑151 (MQ=255)
                                               cGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTg                                                                                                 <  1:369462/151‑1 (MQ=255)
                                                gAAGATGCCAAAGTTGGTGAACTGGTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGc                                                                                                >  1:518370/1‑151 (MQ=255)
                                                      gCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACg                                                                                          <  2:748394/151‑1 (MQ=255)
                                                        cAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGcc                                                                                        <  2:613619/151‑1 (MQ=255)
                                                        cAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGcc                                                                                        >  2:258301/1‑151 (MQ=255)
                                                           aGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGAGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGAACGCTGGAAAACGTGATGCTGACCGAATCCGCGGAAACCCATAAAGCACGCCTGCATAACGCCgg                                                                                      >  1:14818/1‑150 (MQ=255)
                                                                   aaCTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGc                                                                              >  2:209121/1‑150 (MQ=255)
                                                                     cTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCAt                                                                            >  1:316961/1‑150 (MQ=255)
                                                                         tgttCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTg                                                                        <  1:801029/150‑1 (MQ=255)
                                                                            tCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCt                                                                    >  1:743529/1‑151 (MQ=255)
                                                                                 aTGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGt                                                               <  2:108470/151‑1 (MQ=255)
                                                                                          gACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCt                                                      >  2:390511/1‑151 (MQ=255)
                                                                                           aCTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCtt                                                     <  1:840709/151‑1 (MQ=255)
                                                                                               tAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAATCCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTaa                                                  <  2:803121/150‑1 (MQ=255)
                                                                                                  aCGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTGAACttt                                              <  1:684051/151‑1 (MQ=255)
                                                                                                      gCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGt                                           <  1:141604/150‑1 (MQ=255)
                                                                                                        cAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCa                                        <  2:718562/151‑1 (MQ=255)
                                                                                                           cGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACtg                                     >  2:131964/1‑151 (MQ=255)
                                                                                                                 cAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGc                                >  1:261012/1‑150 (MQ=255)
                                                                                                                      aaCTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTaaa                          <  2:255386/151‑1 (MQ=255)
                                                                                                                          ggAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAg                        <  2:651418/149‑1 (MQ=255)
                                                                                                                           gagaTCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAATCag                     <  1:254599/151‑1 (MQ=255)
                                                                                                                           gagaTCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCag                     <  1:203702/151‑1 (MQ=255)
                                                                                                                                 aGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACg               <  2:110750/151‑1 (MQ=255)
                                                                                                                                   ccAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCt             >  1:835698/1‑151 (MQ=255)
                                                                                                                                    cAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTg            >  2:586056/1‑151 (MQ=255)
                                                                                                                                     aGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGc           >  2:195478/1‑151 (MQ=255)
                                                                                                                                      gAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGc           >  1:433761/1‑150 (MQ=255)
                                                                                                                                        aGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATg        <  1:392830/151‑1 (MQ=255)
                                                                                                                                        aGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGCGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATg        >  1:268295/1‑151 (MQ=255)
                                                                                                                                              gcaCGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGAc  >  2:590531/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGAC  >  NC_000913/1953158‑1953450

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