Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,827,811 G→C 11.5% G223A (GGG→GCG) ‡ mtgA → monofunctional biosynthetic peptidoglycan transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,827,8110GC11.5% 142.0 / 10.8 52G223A (GGG→GCG) ‡mtgAmonofunctional biosynthetic peptidoglycan transglycosylase
Reads supporting (aligned to +/- strand):  ref base G (26/20);  new base C (3/3);  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 = 7.06e-01

GTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATC  >  NC_002947/5827666‑5827951
                                                                                                                                                 |                                                                                                                                            
gTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCCATATTCCCCGGCAGATGAGCCAGCTGGGCGGGccc                                                                                                                                          >  2:63107/1‑150 (MQ=255)
     cGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCCCGCCCAGCtg                                                                                                                                     <  2:615088/150‑2 (MQ=255)
      gCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCCCGCCCAGCtgg                                                                                                                                    <  1:465571/150‑3 (MQ=255)
                       ccAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGActc                                                                                                                   >  1:149690/1‑150 (MQ=255)
                            aggcagCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCa                                                                                                              <  1:359247/150‑1 (MQ=255)
                                  cgcagcTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGc                                                                                                        <  1:368685/150‑1 (MQ=255)
                                     cgcTGCCAGCTGCGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCttt                                                                                                     <  1:630336/149‑1 (MQ=255)
                                     agcTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTTATGCAGCTGGACTCCTCCAGAAAGCtt                                                                                                      >  1:449890/1‑149 (MQ=255)
                                     agcTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCtt                                                                                                      <  2:449890/149‑1 (MQ=255)
                                      cttGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTg                                                                                                    >  2:735109/3‑150 (MQ=255)
                                                        aaGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGAc                                                                                                                     <  1:130581/115‑1 (MQ=255)
                                                        aaGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGAc                                                                                                                     >  2:130581/1‑115 (MQ=255)
                                                        aaGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCCGAAAGCTTTGAATGCAGGCTCGCGCCtt                                                                                  >  2:442801/1‑150 (MQ=255)
                                                          gCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTc                                                                                >  2:359471/1‑150 (MQ=255)
                                                              aaTCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCtt                                                                                  <  2:319048/144‑1 (MQ=255)
                                                              aaTCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCtt                                                                                  >  1:319048/1‑144 (MQ=255)
                                                              aaTCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCgg                                                                            >  2:502184/1‑150 (MQ=255)
                                                              aaTCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCgg                                                                            >  2:306583/1‑150 (MQ=255)
                                                              aaTCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCgg                                                                            <  2:202222/150‑1 (MQ=255)
                                                                    gTGGAGTGCCAGCCGACCGAGTTCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAAc                                                                      <  1:75049/150‑1 (MQ=255)
                                                                    gTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCttt                                                                                                     <  2:715622/119‑1 (MQ=255)
                                                                    gTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCttt                                                                                                     >  1:715622/1‑119 (MQ=255)
                                                                        aGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGc                                                                  <  1:502184/150‑1 (MQ=255)
                                                                        aGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGc                                                                  <  1:735109/150‑1 (MQ=255)
                                                                          tGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTc                                                                >  1:191077/1‑150 (MQ=255)
                                                                            ccAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCCTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCt                                                              <  2:626337/150‑1 (MQ=255)
                                                                              aGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCCCGCCCAGCtg                                                                                                                                     >  2:471561/1‑76 (MQ=255)
                                                                              aGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCCCGCCCAGCtg                                                                                                                                     <  1:471561/77‑2 (MQ=255)
                                                                              aGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCGTGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGGAAACCCGCTCCTAc                                                            >  1:63167/1‑150 (MQ=255)
                                                                                ccgaccgaGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGAAAGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACgg                                                          <  1:63107/150‑1 (MQ=255)
                                                                                ccgaccgaGTGCCTACGTGGCCAGCCGGGCCGGCTCGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTCGAATTCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACgg                                                          <  2:426067/150‑1 (MQ=255)
                                                                                  gaccgaGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGc                                                        >  2:286666/1‑150 (MQ=255)
                                                                                        gTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGc                                                  >  2:723643/1‑150 (MQ=255)
                                                                                            cTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATcc                                              >  1:723611/1‑150 (MQ=255)
                                                                                                 tGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGc                                         <  1:256712/150‑1 (MQ=255)
                                                                                                  ggCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGcc                                        <  2:723611/150‑1 (MQ=255)
                                                                                                     cAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCTCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAATCGCATCCCCGGgctcc                                     >  1:112723/1‑150 (MQ=255)
                                                                                                      aGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGcctcct                                    >  2:5024/1‑150 (MQ=255)
                                                                                                          gggccggCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTcgcg                                >  2:740143/1‑150 (MQ=255)
                                                                                                                 cTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAAc                         <  2:149690/150‑1 (MQ=255)
                                                                                                                  tGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAAcc                        >  1:497161/1‑150 (MQ=255)
                                                                                                                  tGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGAGGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCACCTCGCGGGCAAAcc                        >  2:554085/1‑150 (MQ=255)
                                                                                                                            ggCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAACACTGCTCATAGTTTGTCGTGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCAGGCCGACCCGCGCCTACa              >  2:229176/1‑150 (MQ=255)
                                                                                                                             gCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCCAACCCGCGCCTACAg             <  1:626020/150‑1 (MQ=255)
                                                                                                                             gCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAg             <  1:442801/150‑1 (MQ=255)
                                                                                                                               aGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGc           >  1:467082/1‑150 (MQ=255)
                                                                                                                                      gtaGCTGGGCCCGCCCAGCTGGCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCa    >  2:667742/3‑150 (MQ=38)
                                                                                                                                      ccAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCa    >  2:343268/1‑150 (MQ=255)
                                                                                                                                      ccAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCACATCCCCGTCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCa    <  1:359471/150‑1 (MQ=255)
                                                                                                                                       taGCTGGGCCCGCCCAGCTACCTGATGCAGCTGGGCTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACTGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCAt   >  1:530310/2‑150 (MQ=255)
                                                                                                                                        aGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATc  >  2:55749/1‑150 (MQ=255)
                                                                                                                                        aGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATc  >  1:160576/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
GTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATC  >  NC_002947/5827666‑5827951

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