Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,962,329 A→C 11.2% D390A (GAC→GCC)  lysA‑II → diaminopimelate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,962,3290AC11.2% 131.8 / 7.8 51D390A (GAC→GCC) lysA‑IIdiaminopimelate decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (23/22);  new base C (3/3);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCTGCGAGACTGGCGACTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTGGACCTGGTCAGC  >  NC_002947/5962181‑5962478
                                                                                                                                                    |                                                                                                                                                     
aTCTGCGAGACTGGCGACTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGAc                                                                                                                                                      >  2:1361494/1‑150 (MQ=255)
 tCTGCGAGACTGGCGACTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACg                                                                                                                                                     <  2:1173363/150‑1 (MQ=255)
  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                                                                            >  1:369729/1‑138 (MQ=255)
                                                                                      cTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGa                                                                            <  2:369729/138‑1 (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                                                                                                                                    >  2:350719/1‑57 (MQ=255)
                                                                                                               tACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGaa                                                                                                                                    <  1:350719/57‑1 (MQ=255)
                                                                                                                    cacCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCAtgggctgg                                  <  2:416079/150‑1 (MQ=255)
                                                                                                                       ccGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGATTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCa                               >  2:1028045/1‑150 (MQ=255)
                                                                                                                        cgcgGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCAtgggc                                     <  1:546107/143‑1 (MQ=255)
                                                                                                                        cgcgGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCAtgggc                                     >  2:546107/1‑143 (MQ=255)
                                                                                                                            gCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGAt                                             <  1:472271/131‑1 (MQ=255)
                                                                                                                            gCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGAt                                             >  2:472271/1‑131 (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                     <  1:205245/138‑1 (MQ=255)
                                                                                                                                             cTGGTCGCCGTCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCAt                     >  2:205245/1‑138 (MQ=255)
                                                                                                                                              tGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTGGACCTg        <  2:1297362/150‑1 (MQ=255)
                                                                                                                                                    aCGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTGGAGCTGGTcagc  >  1:200296/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATCTGCGAGACTGGCGACTTCCTCGGCAAGGACCGCGTCCTCAACCTGGCCGAAGGTGACCTGCTGGCCGTAGAGTCCGCGGGCGCCTATGGTTTTGTCATGAGCTCCAACTACAACACCCGCGGCCGTTGCGCTGAAATCCTGGTCGACGGCGACCAGGCGTTCGAAGTACGCCGCCGCGAGACCGTCGCCGAACTGTATGCTGGCGAAAGCCTGCTGCCGGAGTAAGCCCATGCTGCTGCGTTTTACCAAGATGCATGGGCTGGGCAACGACTTCATGGTCCTGGACCTGGTCAGC  >  NC_002947/5962181‑5962478

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