Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,399,666 G→A R105C (CGC→TGC)  mreC ← cell wall structural complex MreBCD transmembrane component MreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,399,6660GA98.1% 199.6 / ‑4.7 52R105C (CGC→TGC) mreCcell wall structural complex MreBCD transmembrane component MreC
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (28/23);  total (29/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GATGACCGGCTGGCCTTCATAAACGCCATTAACGCTACCTTTATCGATAACAACTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGA  >  NC_000913/3399523‑3399811
                                                                                                                                               |                                                                                                                                                 
gatgaCCGGCTGGCCTTCATAAACGCCATTAACGCTACCTTTATCGATAACAACTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGAcgc                                                                                                                                            <  1:292139/151‑1 (MQ=255)
gatgaCCGGCTGGCCTTCATAAACGCCATTAACGCTACCTTTATCGATAACAACTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGAc                                                                                                                                              <  1:684024/149‑1 (MQ=255)
gatgaCCGGCTGGCCTTCATAAACGCCATTAACGCTACCTTTATCGATAACAACTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGAc                                                                                                                                              <  2:287600/149‑1 (MQ=255)
   gaCCGGCTGGCCTTCATAAACGCCATTAACGCTACCTTTATCGATAACAACTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGAcgcgcg                                                                                                                                         <  2:740058/151‑1 (MQ=255)
       ggCTGGCCTTCATAAACGCCATTAACGCTACCTTTATCGATAACAACTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTtc                                                                                                                                      >  2:479341/1‑150 (MQ=255)
          tGGCCTTCATAAACGCCATTAACGCTACCTTTATCGATAACAACTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCt                                                                                                                                  >  2:927925/1‑151 (MQ=255)
          tGGCCTTCATAAACGCCATTAACGCTACCTTTATCGATAACAACTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCt                                                                                                                                  >  2:636772/1‑151 (MQ=255)
               ttCATAAACGCCATTAACGCTACCTTTATCGATAACAACTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTATCCTGTTTg                                                                                                                             >  1:951082/1‑151 (MQ=255)
                tCATAAACGCCATTAACGCTACCTTTATCGATAACAACTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTg                                                                                                                             <  2:819231/150‑1 (MQ=255)
                 cATAAACGCCATTAACGCTACCTTTATCGATAACAACTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGt                                                                                                                            <  1:963889/150‑1 (MQ=255)
                    aaaCGCCATTAACGCTACCTTTATCGATAACAACTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTg                                                                                                                        >  2:7421/1‑151 (MQ=255)
                         ccATTAACGCTACCTTTATCGATAACAACTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCa                                                                                                                    >  2:906777/1‑150 (MQ=255)
                            ttAACGCTACCTTTATCGATAACAACTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGc                                                                                                                 >  1:479907/1‑150 (MQ=255)
                                cGCTACCTTTATCGATAACAACTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATc                                                                                                              >  2:518805/1‑149 (MQ=255)
                                      cTTTATCGATAACAACTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCATCAGTTc                                                                                                      >  1:215245/1‑151 (MQ=255)
                                        ttATCGATAACAACTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCAc                                                                                                    >  2:1045312/1‑151 (MQ=255)
                                          aTCGATAACAACTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACt                                                                                                   <  1:761805/150‑1 (MQ=255)
                                            cGATAACAACTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGt                                                                                                 <  2:241574/150‑1 (MQ=255)
                                              aTAACAACTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGtttt                                                                                              <  2:895694/151‑1 (MQ=255)
                                                aacaacTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGttttt                                                                                             <  2:311416/150‑1 (MQ=255)
                                                aacaacTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTc                                                                                            >  1:497206/1‑151 (MQ=255)
                                                 acaacTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCa                                                                                           <  1:959231/151‑1 (MQ=255)
                                                    acTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGc                                                                                         >  2:310321/1‑150 (MQ=255)
                                                       tGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCaa                                                                                       <  2:899998/149‑1 (MQ=255)
                                                              tataAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTg                                                                              >  1:412708/1‑151 (MQ=255)
                                                               ataAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCt                                                                               <  1:578903/149‑1 (MQ=255)
                                                               ataAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGa                                                                             <  1:831717/151‑1 (MQ=255)
                                                                     aTCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAAc                                                                       >  1:835975/1‑151 (MQ=255)
                                                                     aTCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAAc                                                                       >  1:174387/1‑151 (MQ=255)
                                                                       cGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACg                                                                      <  2:372529/150‑1 (MQ=255)
                                                                         ttAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGccc                                                                   >  1:573694/1‑151 (MQ=255)
                                                                         ttAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGccc                                                                   >  2:105474/1‑151 (MQ=255)
                                                                           aaCCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCgg                                                                 <  1:964172/151‑1 (MQ=255)
                                                                           aaCCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCgg                                                                 >  2:123522/1‑151 (MQ=255)
                                                                             ccGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGtt                                                               >  1:393034/1‑151 (MQ=255)
                                                                                        aCCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCGAGTTCTa                                                    <  2:1066674/151‑1 (MQ=255)
                                                                                          cTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAAt                                                  >  1:184653/1‑151 (MQ=255)
                                                                                            gAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAAt                                                  <  1:907925/149‑1 (MQ=255)
                                                                                             aGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTgg                                               <  2:1090305/151‑1 (MQ=255)
                                                                                                 aCCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCAc                                           >  1:235165/1‑151 (MQ=255)
                                                                                                  ccATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCAcg                                          >  1:763058/1‑151 (MQ=255)
                                                                                                        tCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAgg                                     <  2:596402/150‑1 (MQ=255)
                                                                                                        tCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGc                                    >  2:81861/1‑151 (MQ=255)
                                                                                                         cTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGcc                                   >  2:329934/1‑151 (MQ=255)
                                                                                                            cTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGc                                >  1:215610/1‑151 (MQ=255)
                                                                                                                tCCTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAGGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTc                             >  1:100405/1‑150 (MQ=255)
                                                                                                                 ccTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCt                            <  1:167218/150‑1 (MQ=255)
                                                                                                                 ccTGACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTg                           >  2:415695/1‑151 (MQ=255)
                                                                                                                    gACGCAGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGc                          <  1:680032/149‑1 (MQ=255)
                                                                                                                         aGCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGc                   <  1:1057532/151‑1 (MQ=255)
                                                                                                                          gCGGGGAACCCAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGcc                  >  2:1097664/1‑151 (MQ=255)
                                                                                                                                  cccAGCAGCTCGCACAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGc           >  2:422977/1‑150 (MQ=255)
                                                                                                                                            cgcaCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACga  <  1:24182/149‑1 (MQ=255)
                                                                                                                                            cgcaCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACga  <  1:514279/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
GATGACCGGCTGGCCTTCATAAACGCCATTAACGCTACCTTTATCGATAACAACTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGA  >  NC_000913/3399523‑3399811

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