Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 387,864 T→G 16.0% A43A (GCT→GCG)  U712_01810 → Surfactin synthase subunit 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP006881387,8640TG16.0% 97.2 / 3.5 57A43A (GCT→GCG) U712_01810Surfactin synthase subunit 2
Reads supporting (aligned to +/- strand):  ref base T (30/17);  new base G (8/2);  total (38/19)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATATGAGCAAAAAATCGATTCAAAAGGTGTACGCACTGACACCAATGCAGGAGGGAATGCTGTATCATGCGATGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCGGTCATACGATATTCTCCGTACGGTATTTGTGCATCAGCAGCTGCAAAAACCGCGTCAGGTCGTGTTAGCGGAACGCAAAACAAAAGTGCATTAT  >  CP006881/387733‑388002
                                                                                                                                   |                                                                                                                                          
cATATGAGCAAAAAATCGATTCAAAAGGTGTACGCACTGACACCAATGCAGGAGGGAATGCTGTATCATGCGATGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCTTTTGATc                                                                                                                                     >  7:41842/1‑139 (MQ=255)
                        aaGGTGTACGCACTGACACCAATGCAGGAGGGAATGCTGTATCATGCGATGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCa                                                                                                             <  3:14522/139‑1 (MQ=255)
                             gtccGCACTGAAACCAATGCAGGAGGGAATGCTGTATCATGCGCTGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTCGAGCTTGGGATTCACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGaa                                                                                                        <  5:160175/136‑1 (MQ=255)
                             gtACGCACTGACACCAATGCAGGAGGGAATGCTGTATCATGCGATGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCCCGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGaa                                                                                                        >  6:213433/1‑139 (MQ=255)
                             gtACGCACTGACACCAATGCAGGAGGGAATGCTGTATCATGCGATGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGaa                                                                                                        >  4:266928/1‑139 (MQ=255)
                                cGCACTGACACCAATGCAGGAGGGAATGCTGTATCATGCGATGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACt                                                                                                      >  5:107810/1‑138 (MQ=255)
                                   aCTGACACCAATGCTGGAGTGAATTCTGTATCATGCGATGCTTGATCGGTATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGTTTCACGGCGCTTTTGATCTTGAAATATTTGAGAAAAGCGTCAATGAACTGAtt                                                                                                  <  6:14546/139‑1 (MQ=255)
                                     tGACACCAATGCAGGAGGGAATGCTGTATCATGCGATGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCTTTTGTTCTTGAAATATTTTAGAAAAGCGGCAATGAACTGATTCg                                                                                                >  2:277195/1‑139 (MQ=255)
                                     tGACAACAATGCAGGAGGGAATGCTGTATCATGCGATGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCTTTTGATCTTGAAATCTTTTAGAAAAGAGTCAATGAACTGATTct                                                                                                >  6:135011/1‑138 (MQ=255)
                                      gACACCAATGCAGGAGGGAATGCTGTATCATGCGATGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCgg                                                                                               >  6:62596/1‑139 (MQ=255)
                                      gACACCAATGCAGGAGGGAATGCTGTATCATGCGATGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCgg                                                                                               >  2:199619/1‑139 (MQ=255)
                                      gACACCAATGCAGGAGGGAATGCTGTATCATGCGATGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCgg                                                                                               >  4:116129/1‑139 (MQ=255)
                                      gACACCAATGCAGGAGGGAATGCTGTATCATGCGATGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCGTTTTATATTTAAATATTTGTGAAAAGCGTCAATGAACTGATTCgg                                                                                               >  6:157855/1‑139 (MQ=255)
                                        cacCAAAGCAGGACGGAATGCTGTATCATGCGACGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACCGCGCGATTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCGGTc                                                                                             <  3:269365/139‑1 (MQ=255)
                                        caaaaaTGCAGTTGTGATTGATGTATCATGCGATGCTTGATCCGCATTCTTCCTCGAACTTCACACAATTCGTGCTTGGGATTCACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCGGTc                                                                                             <  2:238302/135‑1 (MQ=255)
                                           cAATGCAGGAGGGAATGCTGTATCATGCGATGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCGGTCATa                                                                                          <  5:62596/139‑1 (MQ=255)
                                            aaTTCAGGAGGGAATGCTGTATCCTGCTATGCTTGATCCGCTTTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCGGTCATAc                                                                                         <  3:6691/139‑1 (MQ=255)
                                             aTGCAGGAGGGAATGCTGTATCATGCGATGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCGGTCATAc                                                                                         >  1:156367/1‑138 (MQ=255)
                                             aTGCAGGAGGGAATGCTGTATCATGCGATGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCGGTCATAc                                                                                         >  1:35477/1‑138 (MQ=255)
                                                cAGGAGTGAATGCCGTATCATGCGATGCTTGATGCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGAATGACGGCTCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCGGTCATACGat                                                                                      <  5:44296/138‑1 (MQ=255)
                                                cAGGAGGGAATGCTGTATCATGCGATGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCGTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCTGTCATACGat                                                                                      <  5:294534/138‑1 (MQ=255)
                                                   gaggGAATGCTGTATCATGCGATGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCGTTTGATCTTGAAATCTTTGATAAAAGCGTCAATTAACTGATTCGGTCATACGaaaat                                                                                   >  3:42327/1‑134 (MQ=255)
                                                         aTGCTGTATCATGCGATGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCGGTCATACGATATTCTCCGTa                                                                            <  2:244522/139‑1 (MQ=255)
                                                            cTGTATCATGCGATGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCGGTCATACGATATTCTCCGTACg                                                                          <  6:107810/138‑1 (MQ=255)
                                                               tATCATGCGATGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCTTTTTATATTGAAATCTTTTATAAAAGCGTCAATGAACTGATTAGGTCATACGAAATTATACGTACGGTAt                                                                      >  3:255206/1‑139 (MQ=255)
                                                                aTCATGCGATGCTTGATCCGCATTCTTAATCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCTTTTGATCTTTAAATATTTCTGAAAATCATCAATGAAATTATTCGGTCATACCTTATTATCCGTACGGGAtt                                                                     >  2:74002/1‑139 (MQ=255)
                                                                 tCATGCGATGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCTTTTGATGTTGAAATCTTTGAGAAAAGCGTCAATGATCTGCTTCGGTGATACGATATTCTCCGTTCGGTAttt                                                                    >  4:31617/1‑139 (MQ=255)
                                                                    tGCGATGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGTGATTCACGGCGCGTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGGACAGAGTCGGGCATAGCATATTCTACGTACGTGATTtttg                                                                 >  2:250220/1‑139 (MQ=255)
                                                                    tGCGATGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCGGTCATACGATATTGTGCGTTCGGTATTtgtg                                                                 >  4:232956/1‑139 (MQ=255)
                                                                    tGCGATGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCGGTCATACGATATTCTCCGTACGGTATTtgtg                                                                 >  4:114606/1‑139 (MQ=255)
                                                                    tGCGATGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCGGTCATACGATATTCTCCGTACGGTATTtgtg                                                                 >  1:72721/1‑139 (MQ=255)
                                                                     gCGATGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCGTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCGGTCATACGATATTCTCCGTACGGTATTtgtg                                                                 >  7:86253/1‑138 (MQ=255)
                                                                      cGATGCTTGATCCGCATTCTTACTCGTAGTTCACACAATTAGAGCTGGGGATTCACGGCGCTTTTGATATGGAAAAAGTTGAGAAAACGGTCATTTACCTGATTCGGTGATACGAAATTCTCCTTAAGGTATTTGTGc                                                                >  8:227464/1‑138 (MQ=255)
                                                                         tGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCTTTTGAGCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCGGTCATACGATATTCTCCGTACGGTATTTGTGCAt                                                              >  3:250288/1‑137 (MQ=255)
                                                                         tGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCGTTTTATATTGAAATATTTTAGAAAAGCGTCAATGAACATCTTCGGTGATACGATATTATACGTACGGTATTTTTGCAt                                                              >  5:26470/1‑137 (MQ=255)
                                                                           cTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCTTTTGATGGTGAAATCTTTGGGGAAAGCGGCAAGGCACTGATGGGGGCATTCGATATTCTGCGTTCGGGATGTGTGCATccgc                                                          >  2:20416/1‑139 (MQ=255)
                                                                                    cATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCGGTCATACGATATTCTCCGTACGGTATTTGTGCATCAGCAGCTGCaa                                                  >  7:5117/1‑138 (MQ=255)
                                                                                     aTTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCGTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCTGTCATACGATATTCTCCGTACGGTATTTGTGCATCAGCAGCTGCaaa                                                 >  6:160771/1‑138 (MQ=255)
                                                                                      ttcttcCTCGTACATCACACAATTACCGCTTGGGATTCACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCGGTCATACGATATTCTCCGTACGGTATTTGTGCATCAGCAGCTGCaaaaa                                               <  2:147570/139‑1 (MQ=255)
                                                                                                  cTTCACACAATTAGAGCTTGGGATTCACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCGGTCATACGATATTCTCCGTACGGTATTTGTGCATCAGCAGCTGCAAAAACCGCGTCAGGt                                    <  7:23754/138‑1 (MQ=255)
                                                                                                  cTTCACACAATTAGAGCTTGGGATTCACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCGGTCATACGATATTCTCCGTACGGTATTTGTGCATCAGCAGCTGCAAAAACCGCGTCAGGt                                    <  2:35477/138‑1 (MQ=255)
                                                                                                  cTTCACACAATTAGAGCTTGGGATTCACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCGGTCATACGATATTCTCCGTACGGTATTTGTGCATCAGCAGCTGCAAAAACCGCGTCAGGt                                    <  2:72721/138‑1 (MQ=255)
                                                                                                   ttCACACAATTAGAGCTTGGGATTCACGGCGCGTTTGGTCTTGAAATCTTTGAGAAAAGCGGCCATGAACGGATTCGGTGATACGATATTCTCCGTACGGTATTTGTGCATCAGCCGCTGCCAAAACCGCGTCAGGGCg                                  >  6:255121/1‑139 (MQ=255)
                                                                                                                agCTTGGGACTCACGGCTCTTTTTATCTTGAACTCTTTGAGAAAAGCGTCAATGAACTGAATCGGTCATACGATATTCTACGTACGGTATTTGTGCATCAGCAGCTGCAAAAACCGCGTCATGTCGTGTTTGCGGAACg                     <  7:227464/139‑1 (MQ=255)
                                                                                                                  cTTGGGATTCACGGCGCTTTTTATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTAGGTCATACGATATTATCCGTACGGTATTTGTGCATCAGCAGCTGCAAAAACCGCGTCAGGTCGTGTTAGCGGAACGc                    >  5:53155/1‑138 (MQ=255)
                                                                                                                    tGGGATTCACGGCGCTTTTGTTCTTGAAATCTTTTAGAAAAGCGTCAATGAAATGATTCGGTCATACGATATTCTCCGTACGGTATTTTTGCATCAGCAGCTGCAAAAACCGCGTCAGGTAGTGTTAGCAGAACGcaa                  >  7:239075/1‑138 (MQ=255)
                                                                                                                    tGGGATTCACGGCGCTTTTGTTCTTGAAATCTTTGAGAAAAGCGGCAATGAACTGATTCGGGGATACGATATTCTCCGGACGGTATTTGTGCATCCGCAGCTGCAAAAACCGCGTCATGTCGTGGTTGGGGAACGcaa                  >  6:219283/1‑138 (MQ=255)
                                                                                                                    tGGGATTCACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCGGTCATa                                                                                          >  7:95225/1‑66 (MQ=255)
                                                                                                                    tGGGATTCACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCGGTCATa                                                                                          <  8:95225/66‑1 (MQ=255)
                                                                                                                    tGGGATTCACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCGGTCATACGATATTGTCCGTACGGTATTTGTGCATCAGCAGCTGCAAAAACCGCGTCAGGTCGTGTTTGCGGACCGcaa                  >  8:34271/1‑138 (MQ=255)
                                                                                                                    tGGGATTCACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCGGTCATACGATATTATCCGTACGGTATTTGTTCATCAGCAGCTGCAAAAACCGCGTCAGGTCGTGTTAGCGGAACGcaa                  >  7:34660/1‑138 (MQ=255)
                                                                                                                    tGGGATTCACGGCGCTTTTGATCATGAAATCATTGAGAAAAGCGTCAATGAACTGATTAGGTGATACGATATTCTCCGTACGGTAATTGTGCATCAGCAGCTGCAAAAACCGCGTCACGTAGTGTTAGCCGCACgaac                  >  8:24904/1‑135 (MQ=255)
                                                                                                                    tGGGATTCACGGCGCGTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCGGTCATACGATATTCTCCGTACGGTATTTGTGCATCAGCAGCTGCAAAAACCGCGGCCGGTCGTGGTAGCGGACCGcaa                  >  2:211488/1‑138 (MQ=255)
                                                                                                                     gggATTCACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCGGTCATACGATATTCTCCGTACGGTATTTGTGCATCAGCAGCTGCAAAAACCGCGTCAGGTCGTGTTAGCGGAACGcaa                  <  8:5117/137‑1 (MQ=255)
                                                                                                                      ggATTCACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCGGTCATACGATATTATCCGTACGGTATTTGTGCATCAGCAGCTGCAAAAACCGCGTCAGGTCGTGTTAGCGGAACGcaaaa                >  7:213242/1‑138 (MQ=255)
                                                                                                                           cACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCGGTCATACGATATTCTCCGTACGGTATTTGTGCATCAGCAGCTGCAAAAACCGCGTCAGGTCGTGTTAGCGGAACGCAAAACAAAAg          <  1:250220/139‑1 (MQ=255)
                                                                                                                                cgcTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCGGTCATACGATATTCTCCGTACGGTATTTGTGCATCAGCAGCTGCAAAAACCGCGTCAGGTCGTGTTAGCGGAACGCAAAACAAAAGTGCAt     >  4:205586/1‑139 (MQ=255)
                                                                                                                                   ttttGATCTTGTACTCTTTGAGAATAGCGTCAATGAACTGATTCGGTCATACGATATTCTCCGTACGGTATTTGTGCATCAGCAGCTGCAAAAACCGCGTCAGGTCGTGTTAGCGGAACGCAAAACAAAAGTGCATTAt  <  1:93023/139‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                          
CATATGAGCAAAAAATCGATTCAAAAGGTGTACGCACTGACACCAATGCAGGAGGGAATGCTGTATCATGCGATGCTTGATCCGCATTCTTCCTCGTACTTCACACAATTAGAGCTTGGGATTCACGGCGCTTTTGATCTTGAAATCTTTGAGAAAAGCGTCAATGAACTGATTCGGTCATACGATATTCTCCGTACGGTATTTGTGCATCAGCAGCTGCAAAAACCGCGTCAGGTCGTGTTAGCGGAACGCAAAACAAAAGTGCATTAT  >  CP006881/387733‑388002

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