Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 93,668 A→G E168G (GAG→GGG)  murE → UDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate:meso‑ diaminopimelate ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091393,6680AG98.0% 175.4 / NA 50E168G (GAG→GGG) murEUDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate:meso‑ diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (19/30);  minor base C (1/0);  total (20/30)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGA  >  NC_000913/93520‑93815
                                                                                                                                                    |                                                                                                                                                   
ggCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGc                                                                                                                                                   <  1:825789/151‑1 (MQ=255)
   aaaaCCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTgg                                                                                                                                                <  2:72789/151‑1 (MQ=255)
      aCCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCgg                                                                                                                                             <  2:710105/151‑1 (MQ=255)
      aCCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCgg                                                                                                                                             >  1:235288/1‑151 (MQ=255)
              tACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCtgggg                                                                                                                                     >  2:852672/1‑149 (MQ=255)
               aCCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCtggtgg                                                                                                                                    >  2:265463/1‑151 (MQ=255)
                   aGCGGTTGGCTCAGTAGGGCCGAGTGCTTGACGAAATAAGCGAGGTAATGGGCACCGTTGGTAACGGCCAGTTGGGGAATTTGATGCGGACAGGAAATACAAACGGTTCTGCAGTCGATGTTCAGCATGGGCTGGGGGGACTGATGGATc                                                                                                                                 <  2:428651/150‑1 (MQ=255)
                    gCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAg                                                                                                                               <  2:514588/151‑1 (MQ=255)
                      tGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAggg                                                                                                                             >  1:542570/1‑151 (MQ=255)
                       gTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGTGGCTGGTGGATCAgg                                                                                                                              <  2:667184/149‑1 (MQ=255)
                          ggCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGgcgc                                                                                                                          <  1:498300/150‑1 (MQ=255)
                           gcgcAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGa                                                                                                                        <  1:326230/151‑1 (MQ=255)
                                   gAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTgc                                                                                                                <  2:236655/151‑1 (MQ=255)
                                          cTGCTTGGCGAAATCAGAGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGTAAATACAACTGGTTCGGCAGCCGGTGTTCAGCATGCGCTGGCGGGGCGGGTGGATCAGGGcgcgtaggtatcagcact                                                                                                           >  2:196062/1‑135 (MQ=255)
                                                  cGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTtcc                                                                                                 <  2:894599/151‑1 (MQ=255)
                                                     aaTCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGGGGATCAGGGCGCGACGTTTTGCGCAATGGAAGGTtcctcc                                                                                              >  2:132945/1‑151 (MQ=255)
                                                                 aaTGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTa                                                                                  <  2:765242/151‑1 (MQ=255)
                                                                  aTGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTAc                                                                                 <  1:566420/151‑1 (MQ=255)
                                                                        aCCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCAcc                                                                           >  1:195862/1‑151 (MQ=255)
                                                                          cGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCgt                                                                         <  1:319254/151‑1 (MQ=255)
                                                                           gTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGGTCAGGGCGCGACGGTTTGCGCAATGGAAGGTTCCTCCCCCGGGCTGGGACAGCACCgtg                                                                        >  2:375847/1‑151 (MQ=255)
                                                                           gTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATTAGGGCGCGACGGTTTGGGCAATGGAAGGTTCCTCCCACCGGCTGGGACAGCAcc                                                                           >  2:121427/1‑148 (MQ=255)
                                                                            ttGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGGTTTGCGCAATGGGAGTTTCCTCCCACGGGCTGGGACAGCACCgtgt                                                                       >  2:395608/1‑151 (MQ=255)
                                                                             tGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCgtgt                                                                       <  2:420648/150‑1 (MQ=255)
                                                                               gTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTggc                                                                    <  1:560305/151‑1 (MQ=255)
                                                                                 aaCGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTggcg                                                                   <  2:542520/150‑1 (MQ=255)
                                                                                 aaCGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTggcg                                                                   <  1:132967/150‑1 (MQ=255)
                                                                                   cGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCa                                                                <  2:169306/151‑1 (MQ=255)
                                                                                       ctgctgGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGGGGATCAGGGCGCGACGTTTTGCGCAATGGGAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCAtt                                                              >  2:285178/1‑149 (MQ=255)
                                                                                         gctgGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGaaa                                                          <  2:434455/151‑1 (MQ=255)
                                                                                          ctgGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGaaa                                                          >  1:203880/1‑150 (MQ=255)
                                                                                               gAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTgcg                                                    <  1:150747/151‑1 (MQ=255)
                                                                                                   gTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGCTGTTCAGCATGGGCTGGCGGTGGTGGTGGATCAGGGCGCGACGGTTTGCGCAATGGAAGTTTCCTCCCCCGGGCTGGTACAGCACCGCGTGGCGGCATTTAAATTTgc                                                     >  1:443354/1‑146 (MQ=255)
                                                                                                         ccGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCt                                          <  2:239023/151‑1 (MQ=255)
                                                                                                          cGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCtt                                         >  1:261358/1‑151 (MQ=255)
                                                                                                            aCAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTa                                       <  2:373545/151‑1 (MQ=255)
                                                                                                             cAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTAc                                      <  2:625483/151‑1 (MQ=255)
                                                                                                              aGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGGTCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGGACAGCACCGTGTGGCGGCATTGAAAATTGCGGCGTCGGGCTTTAcc                                     >  2:347191/1‑151 (MQ=255)
                                                                                                                   aTACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACt                                 <  1:245242/150‑1 (MQ=255)
                                                                                                                   aTACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCAGCGTCGGTCTTTACCAACtt                                <  1:229004/151‑1 (MQ=255)
                                                                                                                   aTACAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCCGGTGGGTCGGGGGGCGCGGTTTTGGGGAGTGGGAGTTTCCTCCCCCCGGCTGGTACAGGCCCGGGTGGGGGGCATGAAATTTGAGGGGTCGGGCTTTAc                                      >  2:511974/1‑145 (MQ=255)
                                                                                                                      cAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTTGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTGTCCTCCCTCGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTa                               <  2:533778/149‑1 (MQ=255)
                                                                                                                      cAACCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAg                             <  1:100379/151‑1 (MQ=255)
                                                                                                                       aaCCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGc                            <  2:389972/151‑1 (MQ=255)
                                                                                                                        aCCGGTTCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGcc                           >  1:120270/1‑151 (MQ=255)
                                                                                                                              tCGGCAGTCGATGTTCAGCATGGGCTGGCGGGGCTGGGGGATCAGGGCGCGACGTTTTTCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCCGGTGGCGGCATTGAAAATTGCGGCGTCGGTCTTTTCCAACTTAAGCCGCGATc                     >  2:206235/1‑151 (MQ=255)
                                                                                                                                    gTCGATGTTCAGCATGGGCTGGCGGGGCTGGGGGATCAGGGCGCGACGGTTTGCGCCATGGAAGTTTCCTCCCACGGGCTGGGACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGGCTTTTCCCACTTAAGCCGCGATCACCtt                >  2:142091/1‑150 (MQ=255)
                                                                                                                                        aTGTTCAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGAtt            >  1:902677/1‑150 (MQ=255)
                                                                                                                                             cAGCATGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATg      <  2:347275/151‑1 (MQ=255)
                                                                                                                                                 aTGGGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGa  <  1:765493/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAGCACCGTGTGGCGGCATTGAAATTTGCGGCGTCGGTCTTTACCAACTTAAGCCGCGATCACCTTGATTATCATGGTGA  >  NC_000913/93520‑93815

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