Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 93,598 G→A 7.3% G145S (GGT→AGT)  murE → UDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate‑‑2, 6‑diaminopimelate ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091393,5980GA 7.3% 132.6 / 2.6 55G145S (GGT→AGT) murEUDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate‑‑2, 6‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base G (33/18);  new base A (2/2);  total (35/20)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.51e-01

GCGGGCCGCTTTTACCATGAACCCTCTGACAATTTACGTCTCGTGGGCGTAACGGGCACCAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAG  >  NC_000913/93457‑93738
                                                                                                                                             |                                                                                                                                            
gCGGGCCGCTTGGGCCATGGACCCTCTGACGATTGACGCCTCGTGGGCGTAACGGGCACCAACGGCAAAACCACGACTCCCCAGCGGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAGTCAGCGCGGTAATGGGCACCGTTgg                                                                                                                                             <  2:253869/143‑1 (MQ=255)
      cGCTTTTACCATGAACCCTCTGACAATTTACGTCTCGTGGGCGTAACGGGCACCAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACGGTTGGTAACgg                                                                                                                                       >  1:944117/1‑143 (MQ=255)
         ttttACCATGAACCCTCAGACAATTTACGGCTCGTGGGCGTAACGGGCACCAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGCACTCAGCGCGGGAATGGGCACCGCTGGTAACGGCct                                                                                                                                    >  2:969096/1‑143 (MQ=255)
                        gcttgcgaTTTACTGCGCGGGGGGGTAACGGGCGCCAACGGCAAAACCACGACGACCCAGCTGTTGGTGCAGTGAAACCAACTGGGTGGCGAAAGCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGAt                                                                                                                     <  2:963413/136‑1 (MQ=255)
                          tGACACTTTACGTCTCGTGGGCGTAACGGGCACCAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATcc                                                                                                                   <  1:609976/143‑1 (MQ=255)
                          tGACAATTTACGTCTCGTGGGCGTAACGGGCACCAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATcc                                                                                                                   >  1:914343/1‑143 (MQ=255)
                            aCAATTTACGTCTCGTGGGCGTAACGGGCACCAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTACTGGGCACCGTTGCTAGCGGCCTGCTGCTGGAAGTGATcgtg                                                                                                                 >  1:2799/1‑140 (MQ=255)
                               aTTTACGTCTCGTGGGCGTAACGGGCACCAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACa                                                                                                              >  2:363428/1‑143 (MQ=255)
                                 ttACGTCTCGTGGGCGTAACGGGCACCAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGACGGGGAAAGTGATCCCGACAgc                                                                                                            >  2:801625/1‑142 (MQ=255)
                                 ttACGTCTCGTGGGCGTAACGGGCACCAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGACGGCCACGGCCTGCGGGGGAACGTGCTCCCgcccgg                                                                                                            >  2:777587/1‑138 (MQ=255)
                                        tcGTGGGCGTAACGGGCACCAACGGCAAAACCGCGCCTACCCGGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAGATACa                                                                                                     <  2:918347/143‑1 (MQ=255)
                                         cGTGGGCGGACCGGGCACCAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACaa                                                                                                    <  1:359760/143‑1 (MQ=255)
                                             ggCGCAACGGGCACCAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCgg                                                                                                <  1:140686/143‑1 (MQ=255)
                                                 tAACGGGCACCAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCAGACAGAAAATACAACCGGTTCg                                                                                            >  2:980396/1‑143 (MQ=255)
                                                  aaCGGGCACCAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCgg                                                                                           <  1:598186/143‑1 (MQ=255)
                                                  aaCGGGCACCAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCgg                                                                                           >  1:777673/1‑143 (MQ=255)
                                                     gggCACCAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACGACCGGTTCGGccg                                                                                        >  1:1069639/1‑141 (MQ=255)
                                                     gggCACCAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCAGACAGAACATACGACCGGTCCGGCAg                                                                                        >  2:647056/1‑143 (MQ=255)
                                                       gCACCAACGGCAAAACCATGACTACCCAGCTGTTGGCGCAGGGGAGCAAACTGCTGGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTc                                                                                      <  1:71484/143‑1 (MQ=255)
                                                       gCACCAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCCACTGCTTGGCGAAATCAGCGCGGTAATGAGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTc                                                                                      >  1:1387/1‑143 (MQ=255)
                                                         aCCAACGGCAAAACCACGACTACCCAGCTGTTGCCGCAGCCGCGCAACCTAATTAACGAAATCAGCGCGGTGACGGGGACCGGGGCCAACGGCCAGGAGGGGAAAGTGCTCCCGCCAGCAAATCCAACCGGACCGGCAGTCGa                                                                                    >  2:439319/1‑143 (MQ=255)
                                                           cAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATg                                                                                  <  1:129450/143‑1 (MQ=255)
                                                           cAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTAGGCAGTCGATg                                                                                  >  1:595813/1‑143 (MQ=255)
                                                           cAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGAAATGCGCACCGTTGGTAACGACCAGCTGGGGAAAGTGATCCCGACAGGAAATACAGCCGGTGCGGCAGCCAATg                                                                                  >  2:72579/1‑143 (MQ=255)
                                                           cAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGAGCGGTAATGGGCACCGTTAGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATg                                                                                  >  2:901806/1‑143 (MQ=255)
                                                           cAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGAGCGGTAATGGGCACCGTTAGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATg                                                                                  >  2:899624/1‑143 (MQ=255)
                                                            aaCGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGAGGt                                                                                 >  1:301533/1‑143 (MQ=255)
                                                                 cAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTACCGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGACGAGGTTCGGc                                                                            >  2:285074/1‑143 (MQ=255)
                                                                       cACGACTGCCCAGCTGGTGGCGCAGTGGAGCCAACGGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGAGc                                                                      <  1:409066/143‑1 (MQ=255)
                                                                            ctgcccAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCg                                                                 <  2:1062326/140‑1 (MQ=255)
                                                                             tACCGAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGGAGTGGGCGCGCTAAAAAGCACCGTTGGGAACGGCCTGCGGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCgg                                                                <  1:207372/143‑1 (MQ=255)
                                                                             tACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGCTGTTCAGCATGAGCTGGCgg                                                                >  1:355923/1‑143 (MQ=255)
                                                                                ccAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGACAATACAACCGGTTCGGCAGTCGATGTCCAGCATGAGCTGACAGGGc                                                             >  2:17753/1‑143 (MQ=255)
                                                                                  aGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCGGGGCtg                                                           >  1:214972/1‑143 (MQ=255)
                                                                                          gcgcAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAATTACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCa                                                   <  1:578987/143‑1 (MQ=255)
                                                                                            gcAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAga                                                 >  2:117807/1‑142 (MQ=255)
                                                                                             cAGTGGAGCCAACTGCTTGGCGAAATCAGAGCGGTAATGGGCACCGTTAGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAggg                                                <  1:901806/143‑1 (MQ=255)
                                                                                             cAGTGGAGCCAACTGCTTGGCGAAATCAGAGCGGTAATGGGCACCGTTAGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAggg                                                <  1:899624/143‑1 (MQ=255)
                                                                                                       aaCTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGCTCAGCATGAGCTGGCGGGGCTGCCGGATCACGCCGCGACGttt                                      >  2:168373/1‑143 (MQ=255)
                                                                                                           gCTGGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTgcg                                  <  1:86373/143‑1 (MQ=255)
                                                                                                            cgcggCGAACCCAGCGCGGTAATGGGCACCGGTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAACTACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTgcgc                                 <  2:312837/140‑1 (MQ=255)
                                                                                                            cTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGCCAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTgcgc                                 >  1:901770/1‑143 (MQ=255)
                                                                                                             ttGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTccgcg                                >  1:772070/1‑142 (MQ=255)
                                                                                                               ggCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCact                              >  2:86473/1‑141 (MQ=255)
                                                                                                               ggCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAAt                              <  1:363428/143‑1 (MQ=255)
                                                                                                                gCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCGGGGCTGGTGGAGCAGGGCGCGACGTGTTGCGCCATg                             >  2:671817/1‑143 (MQ=255)
                                                                                                                  gAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGTGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGTGACGTTTTGCGCAGTGGa                           >  1:326236/1‑143 (MQ=255)
                                                                                                                  gAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAATGATCCCGACAGCAAATACAACCGGTTCGGCAGTCGATGTTCAGCATACGCTGCAGGCGCTAGTGAAACAGGCCGCGACGTTTTACGCCATgaa                           >  2:185387/1‑141 (MQ=255)
                                                                                                                   aaaTCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGaa                          >  1:854294/1‑143 (MQ=255)
                                                                                                                      tccgcgcGGTAATGCGCACCGTTGGTCACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGtt                       <  1:613768/140‑1 (MQ=255)
                                                                                                                         gcgcGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAGTGGAAGTTtcc                    >  2:423503/1‑143 (MQ=255)
                                                                                                                           gcgGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTtcctc                  <  2:767115/143‑1 (MQ=255)
                                                                                                                                aaTGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACg             <  1:950189/143‑1 (MQ=255)
                                                                                                                                     gCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCGGGGCTGGGGGATCAGCGCGCGACGCTTTGCGCCATGGGAGTTTGCTCCCACGGGcgg        >  2:22868/1‑141 (MQ=255)
                                                                                                                                        ccGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTa     >  1:701857/1‑143 (MQ=255)
                                                                                                                                           ttGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAg  >  1:177866/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                            
GCGGGCCGCTTTTACCATGAACCCTCTGACAATTTACGTCTCGTGGGCGTAACGGGCACCAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCCTCCCACGGGCTGGTACAG  >  NC_000913/93457‑93738

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: ‑