Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,414,689 C→G 12.9% intergenic (+17/‑58) ackA → / → pta acetate kinase/phosphate acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,414,6890CG12.9% 105.6 / 13.2 47intergenic (+17/‑58)ackA/ptaacetate kinase/phosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (23/18);  new base G (3/3);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.49e-02

CGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGG  >  NC_000913/2414547‑2414828
                                                                                                                                              |                                                                                                                                           
cGCAACCGGGCTTCACGTTCCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCcagctc                                                                                                                                             <  2:574798/143‑1 (MQ=255)
 gCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGGAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCcagctca                                                                                                                                            <  1:164697/143‑1 (MQ=255)
         gCTGCAGGTTTCGTCAAATCTGGTTTCATCAACAAAGAGGGTGCCCGTCCTGCGGTGGTTATCCCAACCAACGAATAACTGGTTATCGCGCAAGGCGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTGAGCTGGCgg                                                                                                                                    <  2:378870/143‑1 (MQ=255)
         gCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCGCACCGCCAGCTCAGCTCGCgg                                                                                                                                    >  2:156103/1‑143 (MQ=255)
         gCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAACTGGCgg                                                                                                                                    <  1:473487/143‑1 (MQ=255)
         ccTGCACGTTTCGGCAAATCTGGTTTAATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCgg                                                                                                                                    <  1:323220/142‑1 (MQ=255)
          cTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTGAGCTGGCGGt                                                                                                                                   <  1:759956/143‑1 (MQ=255)
            gCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCGCACCGCCAGCTGAGCTGGCGGTgt                                                                                                                                 <  1:363048/143‑2 (MQ=255)
            gCACGTTTCGGCAAATCTGGTTTCATCAACAAAAAAGGTATCCGTCTTTCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGc                                                                                                                                 >  1:378880/1‑143 (MQ=255)
              aCGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTg                                                                                                                               <  1:195251/143‑1 (MQ=255)
                 tttCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACAAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGgtt                                                                                                                            >  1:68488/1‑143 (MQ=255)
                        aaaTCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACCCCGCCCGCTCAGCTGGCGCTGCAGCTTATTAAcc                                                                                                                     >  1:536321/1‑143 (MQ=255)
                                tttCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAAt                                                                                                             <  1:382877/143‑1 (MQ=255)
                                   caccgacaaGGACGGTACCCGTCCTGCGGTGGTCATCCCAACCGACGAAGAGCTGGTTACCGCGCAAGACGCGCGCCGCCTGATTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATcgg                                                                                                          <  2:801808/138‑1 (MQ=255)
                                    ctcaACAAAGAAGGTACCCGTCGTGCGGTGGTTATCCCAACCAACGAAGATCTGGTTGTCGCGCCCGTGGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATcggc                                                                                                         <  2:220063/142‑1 (MQ=255)
                                    atcaACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCGGTTTTTTAACCCGCCAAATccgc                                                                                                         >  2:556015/1‑143 (MQ=255)
                                    atcaACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCCGATGTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAGCCCGCCAACCcggc                                                                                                         >  2:508360/1‑143 (MQ=255)
                                     tcaACAAAGAGGGGACCCGTCCTGCGGGGGTTATCCCAACCAAGGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATcggcg                                                                                                        <  2:651346/143‑1 (MQ=255)
                                       aacaaAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGGAGAACTGGCTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACGCGCCGAATcggcggg                                                                                                      >  2:297322/1‑142 (MQ=255)
                                       aacaaAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCGAAGTcggcggc                                                                                                      >  1:390928/1‑142 (MQ=255)
                                         caaCGAAGGTGCCCGTCCTGCGGGGGTTGTCCCGACCAACGAAGAACGGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCTGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTaa                                                                                                    <  2:668546/143‑1 (MQ=255)
                                         caaAGAGGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTaa                                                                                                    <  2:649879/143‑1 (MQ=255)
                                                gggaCCCGTCCTGCGGTGGTTAGCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAaga                                                                                             <  1:22535/140‑1 (MQ=255)
                                                 gTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACAGCCGGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAACTCGGCGCGAACgcagggg                                                                                            >  2:118750/1‑137 (MQ=255)
                                                  tACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGCTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGCACGAgg                                                                                           >  2:123624/1‑143 (MQ=255)
                                                      cGTCCTGCGGTGGTTATCCCAACCAACGAAGACCTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATaa                                                                                       <  2:545390/143‑1 (MQ=255)
                                                        tCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAAc                                                                                     >  1:758424/1‑143 (MQ=255)
                                                               gtggtTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTcc                                                                              <  1:356728/143‑1 (MQ=255)
                                                                  gtTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAACCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGt                                                                           >  1:665238/1‑143 (MQ=255)
                                                                         caaccaacGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGtattatta                                                                    >  2:320633/1‑143 (MQ=255)
                                                                          aaccaacGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGGGGATGAACCGTGTCCCGtattattat                                                                   >  1:47348/1‑143 (MQ=255)
                                                                               acGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATCATGATGCTGa                                                              >  2:698777/1‑143 (MQ=255)
                                                                                 gaagaaCTGGTTATCGCGCAAGACGCGAGCCGCCTGGCTGCCTGATCTCACACCGCCAGCTCCGCTGGCGGTGCTGTTTCGTAACCCGCCAACGCGGCGCTAGCGAAAGAGGCTAAACCGGATCCCGTCCTCTAATGCTCATc                                                            >  2:489651/1‑143 (MQ=255)
                                                                                        tGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGGATTATTATGCTTATCCCTACCg                                                     >  1:678270/1‑143 (MQ=255)
                                                                                               cgcgCAAGCCGCGAGCGGCCTGACGGCCGGATTTCACACCGCCAGCTCAGCTGGCGGTCCTGTTTTGTAACCCGCCAAATCGGCGGTAACG‑AAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGc                                             <  1:690792/143‑1 (MQ=255)
                                                                                                    aaGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCAAAATCGGCGGTCACGAAAGAGGATAAACCGTGTCCCGTATGATGATGCTGATCCCTACCGGGACCAGCGTCg                                         >  1:454472/1‑143 (MQ=255)
                                                                                                          cgggccgccTGACTGCCTGATTGCCCACCGCCAGCTCAGCTGGCGGTGCTGTTTTGCAACCCGACAGATCGGCGGTAACGAGAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGCACCAGCGTGGGTCTGa                                   <  2:115681/140‑1 (MQ=255)
                                                                                                             gccgccTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGCACCAGCGTCGGTCTGACCa                                >  1:755679/1‑143 (MQ=255)
                                                                                                             gccgccTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGGAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCa                                <  1:595074/143‑1 (MQ=255)
                                                                                                             gccgccTGACTGCCTGATTTAACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAACTCGACGGTAACGAAAGAGGATAAACCGCGTCCCGTATTGTTATGCTGATACCTACCGGACCAAGCGTCGGTCTgcaca                                >  2:233384/1‑139 (MQ=255)
                                                                                                                       tGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAAGCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCGGACCAGCGGCAGCCt                      >  2:808417/1‑143 (MQ=255)
                                                                                                                               tccaccCCGCCAGCGCAGCGGGCGGTGCTGTTTTGTAGCCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGa              <  1:354959/141‑1 (MQ=255)
                                                                                                                                tCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGGACCAGCGTCGGTCTGACCAGCGTCAGCCTTGTCGTgct             >  1:334996/1‑141 (MQ=255)
                                                                                                                                   cacCGCCAGCTGAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGCTCCTACCAGCGTCAGCCTTGGCCTCATCCg          >  2:865001/1‑143 (MQ=255)
                                                                                                                                    cccGCCATCTCAGCTGGCGTTGCTGGTTTAGAACGTGCCAAATCGGCCGGAACGAAAGAGGATAAACCGTGGCCCGTACTATTATGCTCATCCCTCCCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGt         <  1:466718/142‑1 (MQ=255)
                                                                                                                                     ccGCCAGCTGAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTg        >  2:551774/1‑143 (MQ=255)
                                                                                                                                     ccGCCAGCTGAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGAGCCGTg        >  2:18977/1‑143 (MQ=255)
                                                                                                                                     ccGCCAGCCCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTg        <  1:436952/143‑1 (MQ=255)
                                                                                                                                        tcagctcagctGGCGGCGCTGTCTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCGACCGGAACCAGCGTCGGTCTGTCTAGTGTCAGCCTTGGCGTGATCCGTGCaa     <  1:636290/142‑1 (MQ=255)
                                                                                                                                        ccAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCaa     >  2:773561/1‑143 (MQ=255)
                                                                                                                                         cagctcagctGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAAt    <  1:724493/143‑1 (MQ=255)
                                                                                                                                           gcgcagctGGCGGTGCTGTTTTGTATCCCGCCAAATCGGCGGTAACGAAAGAGGATAAAGCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATgg  <  1:57407/140‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                           
CGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGG  >  NC_000913/2414547‑2414828

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑