Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,523,411 T→A 11.0% C25S (TGC→AGC)  murG → UDP‑N‑acetylglucosamine‑‑N‑acetylmuramyl‑ (pentapeptide) pyrophosphoryl‑undecaprenol N‑acetylglucosamine transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,523,4110TA11.0% 153.2 / 8.6 53C25S (TGC→AGC) murGUDP‑N‑acetylglucosamine‑‑N‑acetylmuramyl‑ (pentapeptide) pyrophosphoryl‑undecaprenol N‑acetylglucosamine transferase
Reads supporting (aligned to +/- strand):  ref base T (23/24);  new base A (3/3);  total (26/27)
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.53e-01

GGAAAGCCGCACGCACCTGGGCAGTGAGGAGCACGAATTCAGCGAAAGCGATTTTGCCGAGGAGACCAGTCATGGCCGCTGAAGGCAAAAATGTATTGATCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGAACGAGCTGGTGCCGCAGGCGGGTTTGCCATTGCATCTGATTCAGGTCAGCGGCTTGCGTGGCAAGGGCAAGCTGTCGTTGCTC  >  NC_002947/1523268‑1523560
                                                                                                                                               |                                                                                                                                                     
ggAAAGCCGCACGCATCTGGGCAGGGAGGAGCACGAATTCAGCGAAAGCGATTTTGCCGAGGAGACCAGTCATGGCCGCTGAAGGCAAAAATGTATTGATCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGccc                                                                                                                                                 >  2:59658/1‑150 (MQ=255)
 gAAAGCCGCACGCACCTGGGCAGTGAGGAGCACGAATTCAGCGAAAGCGATTTTGCCGAGGAGACCAGTCATGGCCGCTGAAGGCAAAAATGTATTGATCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCg                                                                                                                                                >  1:362273/1‑150 (MQ=255)
  aaaGCCGCACGCACCTGGGCAGTGAGGAGCACGAATTCAGCGAAAGCGATTTTGCCGAGGAGACCAGTCATGGCCGCTGAAGGCAAAAATGTATTGATCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGt                                                                                                                                               <  1:775124/150‑1 (MQ=255)
          acgcacCTGGGCAGTGAGGAGCACGAATTCAGCGAAAGCGATTTTGCCGAGGAGACCAGTCATGGCCGCTGAAGGCAAAAATGTATTGATCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCa                                                                                                                                       >  2:664686/1‑150 (MQ=255)
           cgcacCTGGGCAGTGAGGAGCACGAATTCAGCGAAAGCGATTTTGCCGAGGAGACCAGTCATGGCCGCTGAAGGCAAAAATGTATTGATCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAg                                                                                                                                      <  2:296043/150‑1 (MQ=255)
              acCTGGGCAGTGAGGAGCACGAATTCAGCGAAAGCGATTTGGCCGAGGAGACCAGTCATGGCCGCTGAAGGCAAAAATGTATTGATCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCGCGTGAGTTCCAGgcg                                                                                                                                   >  1:1359471/1‑150 (MQ=255)
                            gagCACGAATTCAGCGAAAGCGATTTTGCCGAGGAGACCAGTCATGGCCGCTGAAGGCAAAAATGTATTGATCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGt                                                                                                                     <  1:664686/150‑1 (MQ=255)
                                   aaTTCAGCGAAAGCGATTTTGCCGAGGAGACCAGTCATGGCCGCTGAAGGCAAAAATGTATTGATCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCActgg                                                                                                              <  1:425171/150‑1 (MQ=255)
                                             aaGCGATTTTGCCGAGGAGACCAGTCATGGCCGCTGAAGGCAAAAATGTATTGATCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTc                                                                                                    >  2:225081/1‑150 (MQ=255)
                                              aGCGATTTTGCCGAGGAGACCAGTCATGGCCGCTGAAGGCAAAAATGTATTGATCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGtt                                                                                                                            >  2:1360211/1‑125 (MQ=255)
                                              aGCGATTTTGCCGAGGAGACCAGTCATGGCCGCTGAAGGCAAAAATGTATTGATCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGtt                                                                                                                            <  1:1360211/125‑1 (MQ=255)
                                                 gATTTTGCCGAGGAGACCAGTCATGGCCGCTGAAGGCAAAAATGTATTGATCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCcgcg                                                                                                >  2:1143794/1‑150 (MQ=255)
                                                  aTTTTGCCGAGGAGACCAGTCATGGCCGCTGAAGGCAAAAATGTATTGATCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGt                                                                                               <  1:1246928/150‑1 (MQ=255)
                                                   ttttGCCGAGGAGACCAGTCATGGCCGCTGAAGGCAAAAATGTATTGATCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTg                                                                                              >  1:1323573/1‑150 (MQ=255)
                                                        ccGAGGAGACCAGTCATGGCCGCTGAAGGCAAAAATGTATTGATCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATc                                                                                         <  1:225081/150‑1 (MQ=255)
                                                        ccGAGGAGACCAGTCATGGCCGCTGAAGGCAAAAATGTATTGATCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATc                                                                                         <  1:267788/150‑1 (MQ=255)
                                                             gagaCCAGTCATGGCCGCTGAAGGCAAAAATGTATTGATCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGaa                                                                                    <  1:827670/150‑1 (MQ=255)
                                                               gaCCAGTCATGGCCGCTGAAGGCAAAAATGTATTGATCATGGCCGGCGGCCCCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGt                                                                                                                             <  2:924526/107‑1 (MQ=255)
                                                               gaCCAGTCATGGCCGCTGAAGGCAAAAATGTATTGATCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGt                                                                                                                             >  1:924526/1‑107 (MQ=255)
                                                                aCCAGTCATGGCCGCTGAAGGCAAAAATGTATTGATCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGAACGa                                                                                 <  2:362273/150‑1 (MQ=255)
                                                                  cAGTCATGGCCGCTGAAGGCAAAAATGTATTGATCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGAACGAg                                                                                >  1:109950/1‑149 (MQ=255)
                                                                  cAGTCATGGCCGCTGAAGGCAAAAATGTATTGATCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGAACGAg                                                                                <  2:109950/149‑1 (MQ=255)
                                                                    gTCATGGCCGCTGAAGGCAAAAATGTATTGATCATGGCCGGCGGCACAGGGGGGCACGTATTCTCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGGGTGCAGTGGCTGGGCACTCCGCGGGGTTTCGAGACGGAGCTg                                                                             >  2:123106/1‑150 (MQ=255)
                                                                     tCATGGCCGCTGAAGGCAAAAATGTATTGATCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGAACGAGCTgg                                                                            >  1:70626/1‑150 (MQ=255)
                                                                         ggCCGCTGAAGGCAAAAATGTATTGATCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGAACGAGCTGGTGcc                                                                        >  1:1043940/1‑150 (MQ=255)
                                                                          gCCGCTGAAGGCAAAAATGTATTGATCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCActggctgg                                                                                                          >  1:659077/1‑115 (MQ=255)
                                                                          gCCGCTGAAGGCAAAAATGTATTGATCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCActggctgg                                                                                                          <  2:659077/115‑1 (MQ=255)
                                                                                 aaGGCAAAAATGTATTGATCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGAACGAGCTGGTGCCGCAGGCgg                                                                <  1:59658/150‑1 (MQ=255)
                                                                                 aaGGCAAAAATGTATTGATCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGAACGAGCTGGTGCCGCAGGCgg                                                                <  1:1143794/150‑1 (MQ=255)
                                                                                        aaaTGTATTGATCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGAACGAGCTGGTGCCGCAGGCGGGTTTGcc                                                         >  1:662886/1‑150 (MQ=255)
                                                                                                  aTCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGAACGAGCTGGTGCCGCAGGCGGGTTTGCCATTGCATCTg                                               >  1:471469/1‑150 (MQ=255)
                                                                                                   tCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGAACGAGCTGGTGCCGCAGGCGGGTTTGCCATTGCATCTGa                                              <  2:1323573/150‑1 (MQ=255)
                                                                                                     aTGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGAACGAGCTGGTGCCGCAGGCGGGTTTGCCATTGCATCTGAtt                                            >  2:1208269/1‑150 (MQ=255)
                                                                                                         ccGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGGTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGAACGAGCTGGTGCCGCAGGCGGGTTTGCCATTGCATCTGATTCAgg                                        >  1:1038670/1‑150 (MQ=255)
                                                                                                                 aCCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGAACGAGCTGGTGCCGCAGGCGGGTTTGCCATTGCATCTGATTCAGGTCAGCGGc                                >  2:140056/1‑150 (MQ=255)
                                                                                                                       gggCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGAACGAGCTGGTGCCGCAGGCGGGTTTGCCATTGCATCTGATTCAGGTCAGCGGCTTGCGt                          <  2:232818/150‑1 (MQ=255)
                                                                                                                              tATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGAACGAGCTGGTGCCGCAGGCGGGTTTGCCATTGCATCTGATTCAGGTCAGCGGCTTGCGTggcaagg                   <  2:662886/150‑1 (MQ=255)
                                                                                                                              actcacgggCGCAGGCCAGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGAACGAGCTGGTGCCGCAGGCGGGTTTGCCATTGCATCTGATTCAGGTCAGCGGCTTGCGTggcaagg                   <  1:1165745/143‑1 (MQ=255)
                                                                                                                              actcacgggCGCAGGCCAGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGAACGAGCTGGTGCCGCAGGCGGGTTTGCCATTGCATCTGATTCAGGTCAGCGGCTTGCGTggcaagg                   >  2:176016/8‑150 (MQ=255)
                                                                                                                                ttCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGAACGAGCTGGTGCCGCAGGCGGGTTTGCCATTGCATCTGATTCAGGTCAGCGGCTTGCGTggcaagggc                 <  2:563053/150‑1 (MQ=255)
                                                                                                                                tcacgggCGCAGGCCAGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGAACGAGCTGGTGCCGCAGGCGGGTTTGCCATTGCa                                                   >  2:1034879/6‑116 (MQ=255)
                                                                                                                                tcacgggCGCAGGCCAGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGAACGAGCTGGTGCCGCAGGCGGGTTTGCCATTGCa                                                   <  1:1034879/111‑1 (MQ=255)
                                                                                                                                  cccGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGAACGAGCTGGTGCCGCAGGCGGGTTTGCCATTGCATCTGATTCAGGTCAGCGGCTTGCGTggcaagggcaa               <  2:346037/150‑1 (MQ=255)
                                                                                                                                   ccGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCATTCCGCGTGGTATCGAGAACGAGCTGGTGCCGCAGGCGGGTTTGCCATTGCATCTGATTCAGGTCAGCgg                                 >  2:600263/1‑131 (MQ=255)
                                                                                                                                   ccGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCATTCCGCGTGGTATCGAGAACGAGCTGGTGCCGCAGGCGGGTTTGCCATTGCATCTGATTCAGGTCAGCgg                                 <  1:600263/131‑1 (MQ=255)
                                                                                                                                   ccGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGAACGAGCTGGTGCCGCAGGCGGGTTTGCCATTGCATCTGATTCAGGTCAGCGGCTTGc                            <  1:991652/136‑1 (MQ=255)
                                                                                                                                   ccGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGAACGAGCTGGTGCCGCAGGCGGGTTTGCCATTGCATCTGATTCAGGTCAGCGGCTTGc                            >  2:991652/1‑136 (MQ=255)
                                                                                                                                   ccGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGACACTCCGCGTGGTATCGAGAACGAGCTGGTGCCGCAGGCGGGTTTGCCATTGCATCTGATTCAGGTCAGCGGCTTGCGTggcaagggcaag              >  1:1328899/1‑150 (MQ=255)
                                                                                                                                      gcgcTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGACACTCCGCGTGGTATCGAGAACGAGCTGGTGCCGCAGGCGGGTTTGCCATTGCATCTGATTCAGGTCAGCGGCTTGCGTGGCAAGGGCAAGCTg           <  2:1328899/150‑1 (MQ=255)
                                                                                                                                      gcgcAGGCCAGCTCCCGTTAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGAACGAGCTGGTGCCGCAGGCGGGTTTGCCATTGCATCTGATTCAGGTCAGCGGCTTGCGTGGCAAGGGCAAGCTg           >  2:380464/1‑150 (MQ=255)
                                                                                                                                          tGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGAACGAGCTGGTGCCGCAGGCGGGTTTGCCATTGCATCTGATTCAGGTCAGCGGCTTGCGTGGCAAGGGCAAGCTGTCGt       >  2:316517/1‑150 (MQ=255)
                                                                                                                                          aggCCAGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGAACGAGCTGGTGCCGCAGGCGGGTTTGCCATTGCATCTGATTCAGGTCAGCGGCTTGCGTGGCAAGGGCAAGCTGTCGt       <  1:176016/149‑1 (MQ=255)
                                                                                                                                           ggCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGAACGAGCTGGTGCCGCAGGCGGGTTTGCCATTGCATCTGATTCAGGTCAGCGGCTTGCGTGGCAAGGGCAAGCTGTCGtt      <  2:1043940/150‑1 (MQ=255)
                                                                                                                                               tGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGAACGAGCTGGTGCCGCAGGCGGGTTTGCCATTGCATCTGATTCAGGTCAGCGGCTTGCGTGGCAAGGGCAAGCTGTCGTTGCTc  >  2:1184943/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
GGAAAGCCGCACGCACCTGGGCAGTGAGGAGCACGAATTCAGCGAAAGCGATTTTGCCGAGGAGACCAGTCATGGCCGCTGAAGGCAAAAATGTATTGATCATGGCCGGCGGCACCGGGGGGCACGTATTCCCGGCGCTGGCCTGCGCCCGTGAGTTCCAGGCGCGCGGTTACAGCGTGCACTGGCTGGGCACTCCGCGTGGTATCGAGAACGAGCTGGTGCCGCAGGCGGGTTTGCCATTGCATCTGATTCAGGTCAGCGGCTTGCGTGGCAAGGGCAAGCTGTCGTTGCTC  >  NC_002947/1523268‑1523560

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