Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,962,332 G→T 10.1% G391V (GGC→GTC)  lysA‑II → diaminopimelate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,962,3320GT10.1% 140.0 / 3.8 49G391V (GGC→GTC) lysA‑IIdiaminopimelate decarboxylase
Reads supporting (aligned to +/- strand):  ref base G (23/21);  new base T (2/3);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

CTGCGAGACTGGCGACTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTGGACCTGGTCAGCCA  >  NC_002947/5962183‑5962480
                                                                                                                                                     |                                                                                                                                                    
cTGCGAGACTGGCGACTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACgg                                                                                                                                                      >  1:1400650/1‑150 (MQ=255)
   cGAGACTGGCGACTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGa                                                                                                                                                   >  2:1354434/1‑150 (MQ=255)
       aCTGGCGACTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGCCGTCGACCAg                                                                                                                                               <  1:550969/150‑1 (MQ=255)
         tGGCGACTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGc                                                                                                                                             <  1:868452/150‑1 (MQ=255)
            cGACTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGTCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGtt                                                                                                                                          <  1:1361494/150‑1 (MQ=255)
               cTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGa                                                                                                                                       >  1:1376492/1‑150 (MQ=255)
                ttCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGaa                                                                                                                                      >  2:470173/1‑150 (MQ=255)
                    tCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTAc                                                                                                                                  <  1:681781/150‑1 (MQ=255)
                          aGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGCCGTCGACCAg                                                                                                                                               <  1:1338956/131‑1 (MQ=255)
                          aGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGCCGTCGACCAg                                                                                                                                               >  2:1338956/1‑131 (MQ=255)
                              ccGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGaa                                                                                                                                      <  1:172208/136‑1 (MQ=255)
                              ccGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGaa                                                                                                                                      >  2:172208/1‑136 (MQ=255)
                                     cTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGATTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTc                                                                                                                 >  1:824748/1‑150 (MQ=255)
                                             ggCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACt                                                                                                         <  1:110449/150‑1 (MQ=255)
                                                       gACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGc                                                                                               <  1:577953/150‑1 (MQ=255)
                                                           tgctgGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGaaa                                                                                           <  2:1376492/150‑1 (MQ=255)
                                                           tgctgGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGGGCGAAGTACGCCGCCGCGGGACCGTCGCCGAACTGTATGCTGGCGaaa                                                                                           >  1:86725/1‑150 (MQ=255)
                                                                 ccGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCctgc                                                                                     >  2:206413/1‑150 (MQ=255)
                                                                 ccGTAGAGTCCGCGAGCGCCTATGGGTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCctgc                                                                                     >  1:392302/1‑150 (MQ=255)
                                                                  cGTAGAGTCCGCGGGCGCCGATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCctgct                                                                                    >  2:574779/1‑150 (MQ=255)
                                                                   gTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTg                                                                                                 <  1:921587/136‑1 (MQ=255)
                                                                   gTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTg                                                                                                 >  2:921587/1‑136 (MQ=255)
                                                                         tCCGCGGGCGCCTATGGTTTTGTCATGAGCCCAAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAg                                                                             <  2:824748/150‑1 (MQ=255)
                                                                          ccGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTTCCGGAGt                                                                            <  1:1354434/150‑1 (MQ=255)
                                                                                gcgcCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGccc                                                                      <  1:712639/150‑1 (MQ=255)
                                                                                    cTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGa                                                                              <  2:369729/138‑1 (MQ=255)
                                                                                    cTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGa                                                                              >  1:369729/1‑138 (MQ=255)
                                                                                     tATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCAtgct                                                                 >  1:1137157/1‑150 (MQ=255)
                                                                                     tATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCAtgct                                                                 >  1:1297362/1‑150 (MQ=255)
                                                                                               tCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCAGACATGAATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGtttt                                                       >  2:38374/1‑150 (MQ=255)
                                                                                                cATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTa                                                      <  1:206413/150‑1 (MQ=255)
                                                                                                cATGAGCTCCAACTACAACACCCGCGGCCGTTGAGCTGCACCCCCGGTCCCCGGGGACCAGGCGTTCGAAGGACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTa                                                      >  1:388023/1‑150 (MQ=255)
                                                                                                  tGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTAcc                                                    >  2:341664/1‑150 (MQ=255)
                                                                                                             tACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGaa                                                                                                                                      <  1:350719/57‑1 (MQ=255)
                                                                                                             tACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGaa                                                                                                                                      >  2:350719/1‑57 (MQ=255)
                                                                                                                  cacCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCAtgggctgg                                    <  2:416079/150‑1 (MQ=255)
                                                                                                                     ccGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGATTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCa                                 >  2:1028045/1‑150 (MQ=255)
                                                                                                                      cgcgGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCAtgggc                                       >  2:546107/1‑143 (MQ=255)
                                                                                                                      cgcgGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCAtgggc                                       <  1:546107/143‑1 (MQ=255)
                                                                                                                          gCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGAt                                               >  2:472271/1‑131 (MQ=255)
                                                                                                                          gCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGAt                                               <  1:472271/131‑1 (MQ=255)
                                                                                                                          gCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGAc                            <  2:392302/150‑1 (MQ=255)
                                                                                                                            cGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACtt                          >  1:767164/1‑150 (MQ=255)
                                                                                                                                gcgcTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATg                      <  1:133526/150‑1 (MQ=255)
                                                                                                                                    tGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTcc                  >  2:1237318/1‑150 (MQ=255)
                                                                                                                                      aaaTCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTg                <  2:388023/150‑1 (MQ=255)
                                                                                                                                           cTGGTCGCCGTCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCAt                       >  2:205245/1‑138 (MQ=255)
                                                                                                                                           cTGGTCGCCGTCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCAt                       <  1:205245/138‑1 (MQ=255)
                                                                                                                                            tGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTGGACCTg          <  2:1297362/150‑1 (MQ=255)
                                                                                                                                                  aCGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTGGAGCTGGTcagc    >  1:200296/1‑150 (MQ=255)
                                                                                                                                                    ggCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTGTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTGGACCTGGTcagcca  <  2:1133189/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTGCGAGACTGGCGACTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTGGACCTGGTCAGCCA  >  NC_002947/5962183‑5962480

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