Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,517,281 C→A 10.9% A157A (GCC→GCA)  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_0029471,517,2810CA10.9% 127.9 / 4.6 45A157A (GCC→GCA) murEUDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate‑‑2, 6‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base C (19/21);  new base A (2/3);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

CACCAACGGCAAGACAAGTGTTACCCAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACATCGCGGTGATGACCAACCTGTCCCGCGACCACCTC  >  NC_002947/1517137‑1517425
                                                                                                                                                   |                                                                                                                                                
cACCAACGGCAAGACAAGTGTTACCCAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTg                                                                                                                                               >  3:288285/1‑148 (MQ=255)
  ccAACGGCAAGACAAGTGTTACCCAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTg                                                                                                                                               <  8:218976/146‑1 (MQ=255)
  ccAACGGCAAGACAAGTGTTACCCAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTg                                                                                                                                               >  7:218976/1‑146 (MQ=255)
   cAACGGCAAGACAAGTGTTACCCAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAAt                                                                                                                                           >  7:58117/1‑149 (MQ=255)
    aaCGGCAAGACAAGTGTTACCCAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATc                                                                                                                                          <  2:224316/149‑1 (MQ=255)
      cgccAAGACAAGTGTTACCACACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGATTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCACGGCAATCgg                                                                                                                                        <  8:236694/146‑2 (MQ=255)
          aaGACAAGTGTTACCCAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCACGGCAATCgggtcg                                                                                                                                    <  8:152938/149‑6 (MQ=255)
                    ttACCCAACTGGTCGCCAAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTc                                                                                                                          <  8:45100/149‑1 (MQ=255)
                        ccAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCATAAGGCGAGTTGAGATAGGCCGCCGCCTGACCACCCCCGACCAGATTGGCGTGCAATCGACCCTCAACCACCTc                                                                                                                          >  1:256411/1‑148 (MQ=255)
                                              gACCTGCTTGGCCTGCGCTGCGGGCTGATTGGCACCCTGGGTACCGTTTTC‑TACTGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCGCGACCGGATTGCCGTGGTATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGc                                                                                                <  2:259998/149‑1 (MQ=255)
                                              gACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTgg                                                                                                 <  4:288285/148‑1 (MQ=255)
                                                ccTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCa                                                                                              >  4:54407/1‑149 (MQ=255)
                                                 cTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTg                                                                                                                                               <  1:145916/99‑1 (MQ=255)
                                                 cTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTg                                                                                                                                               >  2:145916/1‑99 (MQ=255)
                                                    cTTGGCCAGCGCTGCGGCCTGATTGGCAACCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATgg                                                                                           >  8:146203/1‑148 (MQ=255)
                                                          cAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTc                                                                                    >  1:232599/1‑149 (MQ=255)
                                                          cAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCACGGCAATc                                                                                                                                          <  8:258672/95‑1 (MQ=255)
                                                          cAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCACGGCAATc                                                                                                                                          >  7:258672/1‑95 (MQ=255)
                                                              gcTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACTACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTcc                                                                                <  8:58117/149‑1 (MQ=255)
                                                                    gCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGccac                                                                                                   <  2:299146/124‑3 (MQ=255)
                                                                    gCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGccac                                                                                                   >  1:299146/1‑122 (MQ=255)
                                                                         aTTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAgtgt                                                                                     >  3:122646/1‑133 (MQ=255)
                                                                         aTTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAgtgt                                                                                     <  4:122646/133‑1 (MQ=255)
                                                                             gCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTg                                                                                                                                               <  5:327354/71‑1 (MQ=255)
                                                                             gCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTg                                                                                                                                               >  6:327354/1‑71 (MQ=255)
                                                                              cACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAg                                                                <  2:256411/149‑1 (MQ=255)
                                                                                        aCCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCGAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGc                                                      >  3:177628/1‑149 (MQ=255)
                                                                                            gTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCa                                                  >  8:6985/1‑149 (MQ=255)
                                                                                               tC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACt                                                <  3:54407/148‑1 (MQ=255)
                                                                                                c‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACTg                                               <  3:216347/148‑1 (MQ=255)
                                                                                                c‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACTg                                               <  1:119567/148‑1 (MQ=255)
                                                                                                c‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACTg                                               <  2:232599/148‑1 (MQ=255)
                                                                                                   aCGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCa                                                                  >  5:223379/1‑127 (MQ=255)
                                                                                                   aCGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCa                                                                  <  6:223379/127‑1 (MQ=255)
                                                                                                    cGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACTgg                                              <  7:146203/146‑1 (MQ=255)
                                                                                                           ttG‑CAGA‑GCGGCGGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACAt                                    >  4:197429/1‑149 (MQ=255)
                                                                                                               cAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACATcgc                                 >  1:169108/1‑149 (MQ=255)
                                                                                                                  a‑gCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACATCGCGGt                              >  5:254166/1‑149 (MQ=255)
                                                                                                                  a‑gCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACATCGCGGt                              <  7:6985/149‑1 (MQ=255)
                                                                                                                               aCCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCGAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACATCGCGGTGATGACCAACCt                  <  4:177628/149‑1 (MQ=255)
                                                                                                                                      ccgacccgaTTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACATCGCGGAGATGACAAACCTGTCCcgc           >  8:236248/1‑149 (MQ=255)
                                                                                                                                           tcgaTTGCACGGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACCTCGCGGTGATGACCAACCTGTCCCGCGccca      >  2:230862/2‑149 (MQ=255)
                                                                                                                                           ccgaTTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACATCGCGGTGATGACCaa                     <  6:225257/134‑1 (MQ=255)
                                                                                                                                           ccgaTTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATCCCCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACATCGCGGTGATGACCaa                     >  5:225257/1‑134 (MQ=255)
                                                                                                                                               ttGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACATCGCGGTGATGACCAACCTGTCCCGCGACCACCTc  <  6:254166/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CACCAACGGCAAGACAAGTGTTACCCAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTC‑TACGGCGAGTTG‑CAGA‑GCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGCCCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACATCGCGGTGATGACCAACCTGTCCCGCGACCACCTC  >  NC_002947/1517137‑1517425

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