Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 635,841 A→T 11.2% E22V (GAA→GTA) ‡ mpl → UDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947635,8410AT11.2% 142.5 / 7.0 53E22V (GAA→GTA) ‡mplUDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base A (23/24);  new base T (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

TTACCAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCAT  >  NC_002947/635694‑635988
                                                                                                                                                    |                                                                                                                                                    
ttACCAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGaa                                                                                                                                                     <  6:314499/149‑1 (MQ=255)
 tACCAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTCCCCTGCCTCGAGGT‑TTGGCGAGGTCGACGCTGCATATTCACGTTGGGGGTATTTGCGGGAGTTTGATGGGTTCGCTGGCGGTGCTGGCCAAGGacc                                                                                                                                                    >  4:183476/1‑147 (MQ=255)
    cAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTg                                                                                                                                                <  3:32418/149‑1 (MQ=255)
      ggggCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGc                                                                                                                                              >  5:271242/1‑149 (MQ=255)
        ggCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGTTCC‑TTGGCCa                                                                                                                                            <  5:168568/149‑1 (MQ=255)
           gACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGTATTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCAAAGGAAC‑TTGGCCAcc                                                                                                                                          >  1:88689/1‑149 (MQ=255)
                aGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGTTCC‑TTGGCCag                                                                                                                                           <  5:62615/142‑2 (MQ=255)
                aGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACcgcg                                                                                                                                       <  6:112858/146‑1 (MQ=255)
                aGCGAAACAGCGGGTGCGATTTGTCGCGACCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTATCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGTTCC‑TTGGCCag                                                                                                                                           >  6:62615/1‑141 (MQ=255)
                 gCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCAc                                                                                                                                   >  6:310909/1‑149 (MQ=255)
                   gAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCg                                                                                                                                 >  6:283425/1‑149 (MQ=255)
                       cAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGctc                                                                                                                             <  6:163266/149‑1 (MQ=255)
                        aGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGctct                                                                                                                            >  1:69895/1‑149 (MQ=255)
                          cGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGa                                                                                                                          <  6:76040/149‑1 (MQ=255)
                           gggTGCGATTTGTCGCGCCCTGGGCTACCCTGTCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGAc                                                                                                                         >  1:263302/1‑149 (MQ=255)
                           gggTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGAc                                                                                                                         <  6:159197/149‑1 (MQ=255)
                           gggTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGAc                                                                                                                         <  5:191829/149‑1 (MQ=255)
                              tGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGcg                                                                                                                      <  4:108704/149‑2 (MQ=255)
                                       tCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTAt                                                                                                             >  5:109718/1‑149 (MQ=255)
                                              cTGGGCTACCCTGTCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGa                                                                                                      <  2:263302/149‑1 (MQ=255)
                                              cTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGa                                                                                                      <  3:5404/149‑1 (MQ=255)
                                                 ggCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGa                                                                                                   <  6:109718/149‑1 (MQ=255)
                                                 ggCTACCCTGCCTCGATGT‑TTTGCGAGGTCGAACATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGa                                                                                                   <  5:177716/149‑1 (MQ=255)
                                                     aCCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCTGTATTTGCTGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGTCTCTGACGCCAACGTCTATCCCCCGATGAGCa                                                                                                <  5:310909/148‑1 (MQ=255)
                                                         tGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAg                                                                                           >  2:308302/1‑149 (MQ=255)
                                                           ccTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGc                                                                                          <  2:69895/148‑1 (MQ=255)
                                                               gATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCg                                                                                       >  5:78403/1‑147 (MQ=255)
                                                                         cGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAggg                                                                             <  8:141063/148‑1 (MQ=255)
                                                                          gAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCa                                                                                <  2:276896/144‑1 (MQ=255)
                                                                          gAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCa                                                                                >  1:276896/1‑144 (MQ=255)
                                                                            ggTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATc                                                                         <  6:271242/149‑1 (MQ=255)
                                                                                            tCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGaa                                                                                     <  3:261159/121‑1 (MQ=255)
                                                                                            tCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGaa                                                                                     >  4:261159/1‑121 (MQ=255)
                                                                                               caTTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTAt                                                       >  7:84967/1‑148 (MQ=255)
                                                                                                  tctcGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATACCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGAcc                                                   >  7:124818/1‑149 (MQ=255)
                                                                                                    tcGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAACGTTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGAccc                                                  >  1:244205/1‑149 (MQ=255)
                                                                                                    tcGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCg                                                 >  5:144440/1‑149 (MQ=255)
                                                                                                      ggTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGc                                               <  6:144440/149‑1 (MQ=255)
                                                                                                        tATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGcc                                              >  4:41741/1‑148 (MQ=255)
                                                                                                                   cTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCg                                  <  8:284038/149‑1 (MQ=255)
                                                                                                                        aTGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGGACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGgcgc                              <  6:78403/148‑1 (MQ=255)
                                                                                                                               cGCTGGCGGTGCCGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCtg                      >  5:285066/1‑149 (MQ=255)
                                                                                                                                    gCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGGCGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTggtggt                 >  5:291189/1‑149 (MQ=255)
                                                                                                                                    gCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCtggtg                   <  2:88689/147‑1 (MQ=255)
                                                                                                                                       gcggtGGCCAAGTTCC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATc                                    <  2:57015/123‑1 (MQ=255)
                                                                                                                                       gcggtGGCCAAGTTCC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATc                                    >  1:57015/5‑127 (MQ=255)
                                                                                                                                           tGGCCAAGTTCC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCgg           >  4:297791/1‑148 (MQ=255)
                                                                                                                                            ggCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGc          >  7:43030/1‑148 (MQ=255)
                                                                                                                                             gCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGc          <  1:308302/147‑1 (MQ=255)
                                                                                                                                               cAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACACAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACg      >  2:260433/1‑149 (MQ=255)
                                                                                                                                                 aGGAAC‑TTGGCCCCCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGCACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGc     >  6:280691/1‑148 (MQ=255)
                                                                                                                                                  ggAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCa   >  5:121348/1‑149 (MQ=255)
                                                                                                                                                   gAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCt                       <  6:251637/128‑1 (MQ=255)
                                                                                                                                                   gAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCt                       >  5:251637/1‑128 (MQ=255)
                                                                                                                                                   gAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCAt  <  3:183476/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTACCAGGGGCGACTGAGCGAAACAGCGGGTGCGATTTGTCGCGCCCTGGGCTACCCTGCCTCGATGT‑TTTGCGAGGTCGACCATGCATATTCACATTCTCGGTATTTGCGGCACTTTCATGGGTTCGCTGGCGGTGCTGGCCAAGGAAC‑TTGGCCACCGCGTCACCGGCTCTGACGCCAACGTCTATCCCCCGATGAGCACCCAGCTCGAAGCCCAGGGTATCGAACTGACCCAAGGCTATGACCCGGCCCAGCTGGATCCGGCGCCAGACCTGGTGGTCATCGGCAACGCCAT  >  NC_002947/635694‑635988

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