Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 93,579 C→A 7.8% S138R (AGC→AGA)  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,5790CA 7.8% 133.5 / 2.7 53S138R (AGC→AGA) murEUDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate‑‑2, 6‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base C (31/17);  new base A (2/2);  total (34/19)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.54e-01

GAGCGTTTATCTGCACTGGCGGGCCGCTTTTACCATGAACCCTCTGACAATTTACGTCTCGTGGGCGTAACGGGCACCAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCC  >  NC_000913/93439‑93720
                                                                                                                                            |                                                                                                                                             
gAGCGTTTATCTGCACTGGCGGGCCGCTTTTACCATGAACCCTCTGACGATTTACGTCTCGTGGGCGTACCGGGCACCAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAgcgc                                                                                                                                             <  1:302876/143‑1 (MQ=255)
     tttATCTGCACTGGCGGGCCGCTTTTACCATGAACCCTCTGACAATTTACGTCTCGTGGGCGTAACGGGCACCAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTaa                                                                                                                                        >  2:598186/1‑143 (MQ=255)
             cACTGGCGGGCCGCTTTTACCATGAACCCTCTGACAATTTACGTCTCGTGGGCGTAACGGGCACCAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCAcc                                                                                                                                >  1:934331/1‑143 (MQ=255)
                tGGCGGGCCGCTTTTACCATGAACCCTCTGACAATTTACGTCTCGTGGGCGTAACGGGCACCAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGtt                                                                                                                             <  1:725363/143‑1 (MQ=255)
                  gCGGGCCGCTTTTACCATGAACCCTCTGACAATTTACGTGTCGTGGGCCTCACGGGCACCGACGGCAACCACGCGACTACCCAGCTGTCGCCGCAGTCGAGCCAAACGCGCGGCGACATCGGGGCGGCAATCGGCAACGttcg                                                                                                                           >  2:239706/1‑141 (MQ=255)
                  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:899624/1‑143 (MQ=255)
                                                                             cAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGAGCGGTAATGGGCACCGTTAGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATg                                                                >  2:901806/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)
                                                                                                                                            |                                                                                                                                             
GAGCGTTTATCTGCACTGGCGGGCCGCTTTTACCATGAACCCTCTGACAATTTACGTCTCGTGGGCGTAACGGGCACCAACGGCAAAACCACGACTACCCAGCTGTTGGCGCAGTGGAGCCAACTGCTTGGCGAAATCAGCGCGGTAATGGGCACCGTTGGTAACGGCCTGCTGGGGAAAGTGATCCCGACAGAAAATACAACCGGTTCGGCAGTCGATGTTCAGCATGAGCTGGCGGGGCTGGTGGATCAGGGCGCGACGTTTTGCGCAATGGAAGTTTCC  >  NC_000913/93439‑93720

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