Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,699,077 T→C 38.9% S319S (AGT→AGC csiE → stationary phase inducible protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,699,0770TC38.9% 25.7 / 37.3 59S319S (AGT→AGCcsiEstationary phase inducible protein
Reads supporting (aligned to +/- strand):  ref base T (26/10);  new base C (11/12);  total (37/22)
Fisher's exact test for biased strand distribution p-value = 9.68e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

CACCACCCGCGCGGCGCTGGCCGGATTTGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGAT  >  minE/1699011‑1699142
                                                                  |                                                                 
caccacCCGCGCGGCGCTGGCCGGATTTGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGtg                                                                  <  1:1188665/68‑1 (MQ=255)
caccacCCGCGCGGCGCTGGCCGGATTTGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGtg                                                                  <  1:1109687/68‑1 (MQ=255)
caccacCCGCGCGGCGCTGGCCGGATTTGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGtg                                                                  <  1:2843942/68‑1 (MQ=255)
caccacCCGCGCGGCGCTGGCCGGATTTGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGtg                                                                  <  1:2359789/68‑1 (MQ=255)
caccacCCGCGCGGCGCTGGCCGGATTTGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGtg                                                                  <  1:2199800/68‑1 (MQ=255)
caccacCCGCGCGGCGCTGGCCGGATTTGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGtg                                                                  <  1:1351751/68‑1 (MQ=255)
caccacCCGCGCGGCGCTGGCCGGATTTGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGtg                                                                  <  1:3111178/68‑1 (MQ=255)
caccacCCGCGCGGCGCTGGCCGGATTTGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGtg                                                                  <  1:2052624/68‑1 (MQ=255)
caccacCCGCGCGGCGCTGGCCGGATTTGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGtg                                                                  <  1:3270521/68‑1 (MQ=255)
caccacCCGCGCGGCGCTGGCCGGATTTGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGtg                                                                  <  1:1694819/68‑1 (MQ=255)
caccacCCGCGCGGCGCTGGCCGGATTTGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGtg                                                                  <  1:490956/68‑1 (MQ=255)
  ccaccCGCGCGGCGCTGGCCGGATTTGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGtggtc                                                               >  1:2509248/1‑69 (MQ=255)
    accCGCGCGGCGCTGGCCGGATTTGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGtggtctg                                                             >  1:298753/1‑69 (MQ=255)
    accCGCGCGGCGCTGGCCGGATTTGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGtggtctg                                                             >  1:1932447/1‑69 (MQ=255)
            ggcgCTGGCCGGATTTGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGt                                                     <  1:497981/69‑1 (MQ=255)
             gcgcTGGCCGGATTTGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTg                                                    >  1:358808/1‑69 (MQ=255)
             gcgcTGGCCGGATTTGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTg                                                    >  1:3358450/1‑69 (MQ=255)
                      ggATTTGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGt                                           <  1:516141/69‑1 (MQ=255)
                       gATTTGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGtg                                                                  <  1:394747/45‑1 (MQ=255)
                           tGAAAGTTAATACGGCGTCCATCTTTCTGATGAGGAAAGtg                                                                  <  1:2266225/41‑1 (MQ=255)
                           tGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGtg                                                                  <  1:3255499/41‑1 (MQ=255)
                           tGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTg                                      <  1:930618/69‑1 (MQ=255)
                           tGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTg                                      <  1:3224477/69‑1 (MQ=255)
                           tGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTg                                      <  1:573378/69‑1 (MQ=255)
                            gAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTg                                      <  1:3106256/68‑1 (MQ=255)
                            gAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTg                                      <  1:3086634/68‑1 (MQ=255)
                               aGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCt                                   <  1:1686863/68‑1 (MQ=255)
                                      aCGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAg                            >  1:3007300/1‑68 (MQ=255)
                                      aCGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAg                            >  1:1338400/1‑68 (MQ=255)
                                       cGGCGTCCATCTTTCTGATGAGGAAGGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAgg                           >  1:1474800/1‑68 (MQ=255)
                                       cGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATTCAgg                           >  1:122002/1‑68 (MQ=255)
                                       cGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAgg                           >  1:1995739/1‑68 (MQ=255)
                                       cGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAgg                           >  1:1236708/1‑68 (MQ=255)
                                         gCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGaa                         <  1:2726659/68‑1 (MQ=255)
                                            tCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAAc                      >  1:840798/1‑68 (MQ=255)
                                            tCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACg                     >  1:1343292/1‑69 (MQ=255)
                                              cATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGAt                   >  1:3333019/1‑69 (MQ=255)
                                              cATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGAt                   >  1:288099/1‑69 (MQ=255)
                                              cATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGAt                   >  1:445025/1‑69 (MQ=255)
                                              cATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGAt                   >  1:2185535/1‑69 (MQ=255)
                                                tCTTTCTGATGAGGAAAGTGGTCTGGTCGCGTTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATc                  >  1:2409344/1‑68 (MQ=255)
                                                tCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGAt                   >  1:802418/1‑67 (MQ=255)
                                                tCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATc                  >  1:1485925/1‑68 (MQ=255)
                                                tCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATc                  >  1:2772686/1‑68 (MQ=255)
                                                tCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATc                  >  1:1658932/1‑68 (MQ=255)
                                                tCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATc                  >  1:1888167/1‑68 (MQ=255)
                                                tCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATc                  >  1:307662/1‑68 (MQ=255)
                                                 cTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTg                >  1:1522807/1‑69 (MQ=255)
                                                 cTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTg                >  1:2365187/1‑69 (MQ=255)
                                                 cTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTg                >  1:2236664/1‑69 (MQ=255)
                                                 cTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTg                >  1:1893693/1‑69 (MQ=255)
                                                 cTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTg                >  1:2100424/1‑69 (MQ=255)
                                                  tttCTGATGAGGAAAGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGaa                         >  1:2547270/1‑59 (MQ=255)
                                                  tttCTGATGAGGAAAGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGaa                         >  1:774549/1‑59 (MQ=255)
                                                  tttCTGATGAGGAAAGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGaa                         >  1:1879509/1‑59 (MQ=255)
                                                  tttCTGATGAGGAAAGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGaa                         >  1:316880/1‑59 (MQ=255)
                                                  tttCTGATGAGGAAAGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGaa                         >  1:2519207/1‑59 (MQ=255)
                                                  tttCTGATGAGGAAAGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGaa                         >  1:3040980/1‑59 (MQ=255)
                                                  tttCTGATGAGGAAAGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGaa                         >  1:2849114/1‑59 (MQ=255)
                                                  tttCTGATGAGGAAAGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGaa                         >  1:2347926/1‑59 (MQ=255)
                                                  tttCTGATGAGGAAAGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGaa                         >  1:2903541/1‑59 (MQ=255)
                                                  tttCTGATGAGGAAAGCGGTCTGGTCGCGGTGAATTTTGGt                                           >  1:1559450/1‑41 (MQ=37)
                                                  tttCTGATGAGCAAAGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGaa                         >  1:2086869/1‑59 (MQ=255)
                                                         tgatgaAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTg                                      <  1:1637531/39‑1 (MQ=255)
                                                               aaGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGAAAACGACCTGCATGAGAAACAGat  <  1:1870099/69‑1 (MQ=255)
                                                               aaGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGAAAACGACCTGCATGAGAAACAGat  <  1:1718305/69‑1 (MQ=255)
                                                               aaGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGAAAACGACCTGCATGAGAAACAGat  <  1:1700660/69‑1 (MQ=255)
                                                               aaGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGAAAACGACCTGCATGAGAAACAGat  <  1:2966101/69‑1 (MQ=255)
                                                               aaGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGAAAACGACCTGCATGAGAAACAGat  <  1:1581362/69‑1 (MQ=255)
                                                               aaGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGAAAACGACCTGCATGAGAAACAGat  <  1:2158215/69‑1 (MQ=255)
                                                               aaGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGAAAACGACCTGCATGAGAAACAGat  <  1:1523948/69‑1 (MQ=255)
                                                               aaGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGAAAACGACCTGCATGAGAAACAGat  <  1:197341/69‑1 (MQ=255)
                                                               aaGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGAAAACGACCTGCATGAGAAACAGat  <  1:2472251/69‑1 (MQ=255)
                                                               aaGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGAAAACGACCTGCATGAGAAACAGat  <  1:2631683/69‑1 (MQ=255)
                                                               aaGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGAAAACGACCTGCATGAGAAACAGat  <  1:2140985/69‑1 (MQ=255)
                                                               aaGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGAAAACGACCTGCATGAGAAACAGat  <  1:948341/69‑1 (MQ=255)
                                                                  |                                                                 
CACCACCCGCGCGGCGCTGGCCGGATTTGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGAT  >  minE/1699011‑1699142

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