Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,816 T→C 12.1% H342H (CAT→CAC cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,8160TC12.1% 177.7 / 15.2 66H342H (CAT→CACcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base T (32/26);  new base C (4/4);  total (36/30)
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.83e-01

GTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACG  >  minE/385748‑385886
                                                                    |                                                                      
gTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATgc                                                                      >  1:3956390/1‑71 (MQ=255)
  aTGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATgcgc                                                                    >  1:1966939/1‑71 (MQ=255)
     gACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTg                                                                  <  1:3286345/70‑1 (MQ=255)
     gACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTg                                                                  <  1:3891542/70‑1 (MQ=255)
     gACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGa                                                                 >  1:3595653/1‑71 (MQ=255)
      aCCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGAt                                                                <  1:1313880/71‑1 (MQ=255)
          gTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGAtctc                                                             >  1:1286010/1‑70 (MQ=255)
          gTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGAt                                                                >  1:777853/1‑67 (MQ=255)
          gTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATctct                                                            <  1:3176107/71‑1 (MQ=255)
           tAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCtctt                                                           <  1:2134460/71‑1 (MQ=255)
             aCCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCtcttct                                                         >  1:3546339/1‑71 (MQ=255)
             aCCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCtcttc                                                          <  1:1236067/70‑1 (MQ=255)
                tGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTgg                                                       >  1:1109884/1‑70 (MQ=255)
                tGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTgg                                                       >  1:1779544/1‑70 (MQ=255)
                tGTTCCCGTTCCTGTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCtcttcc                                                         <  1:2596108/68‑2 (MQ=255)
                 gTTCCCGTTCCTGTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCggta                                                     >  1:3889583/1‑69 (MQ=255)
                  ttCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCtcttc                                                          >  1:3576218/1‑65 (MQ=255)
                    cccGTTCCTGTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTAcga                                                  >  1:3531833/1‑69 (MQ=255)
                      cGTTCCTGTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATctct                                                            >  1:3208517/1‑59 (MQ=255)
                       gTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTg                                                                  <  1:1885310/52‑1 (MQ=255)
                          ccTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCg                                             >  1:3569403/1‑70 (MQ=255)
                          ccTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCg                                             >  1:1577627/1‑70 (MQ=255)
                              tttATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATgc                                                                      <  1:3088359/41‑1 (MQ=255)
                              tttATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAg                                          <  1:3996159/69‑1 (MQ=255)
                                 aTCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCt                                     <  1:483455/71‑1 (MQ=255)
                                 aTCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGc                                      <  1:439657/70‑1 (MQ=255)
                                    aCCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGa                                                                 >  1:132746/1‑40 (MQ=255)
                                      cATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCt                                     <  1:164617/66‑1 (MQ=255)
                                        tCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaat                              <  1:1062723/71‑2 (MQ=255)
                                        tCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaat                              <  1:2842094/71‑2 (MQ=255)
                                        tCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaat                              >  1:2736575/1‑70 (MQ=255)
                                          gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                             >  1:3856479/1‑70 (MQ=255)
                                          gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                             <  1:2907269/70‑1 (MQ=255)
                                          gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                             <  1:158737/70‑1 (MQ=255)
                                          gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacaa                            <  1:2581713/71‑3 (MQ=255)
                                            cTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCtcttc                                                          >  1:4100515/1‑39 (MQ=255)
                                            cTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAg                          >  1:77299/1‑71 (MQ=255)
                                             tgtgGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacga                            <  1:2566111/68‑1 (MQ=255)
                                             tgtgGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAgg                         >  1:2097676/1‑71 (MQ=255)
                                              gtgGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGAcc                                                   >  1:3757692/1‑44 (MQ=255)
                                              gtgGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGAcc                                                   >  1:742290/1‑44 (MQ=255)
                                                ggCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAggg                        >  1:3465857/1‑69 (MQ=255)
                                                ggCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAggg                        >  1:2921285/1‑69 (MQ=255)
                                                   gcgGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGaa                                           <  1:743418/47‑1 (MQ=255)
                                                      gtgtCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCg                  >  1:414362/1‑69 (MQ=255)
                                                      gtgtCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGt                 >  1:4022824/1‑70 (MQ=255)
                                                      gtgtCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTg                >  1:829044/1‑71 (MQ=255)
                                                         tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGaa                                           <  1:2734030/41‑1 (MQ=255)
                                                         tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTgg                                   >  1:345771/1‑49 (MQ=255)
                                                         tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacga                            <  1:3645840/56‑1 (MQ=255)
                                                         tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCtt             <  1:2943414/71‑1 (MQ=255)
                                                         tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCtt             >  1:54322/1‑71 (MQ=255)
                                                           tGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGc                                  >  1:3900882/1‑48 (MQ=255)
                                                           tGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTa           >  1:1125927/1‑71 (MQ=255)
                                                             gCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCttt            <  1:1100498/68‑1 (MQ=255)
                                                             gCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCttt            <  1:1556019/68‑1 (MQ=255)
                                                             gCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTa           <  1:547433/69‑1 (MQ=255)
                                                             gCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTAt          >  1:2635685/1‑70 (MQ=255)
                                                             gCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATc         >  1:3676676/1‑71 (MQ=255)
                                                              cTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCttt            <  1:3661634/67‑1 (MQ=255)
                                                              cTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCg        >  1:647106/1‑71 (MQ=255)
                                                                 ccATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGc      >  1:4145850/1‑70 (MQ=255)
                                                                  cATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTa    >  1:2520274/1‑71 (MQ=255)
                                                                  cATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTa    <  1:1995040/71‑1 (MQ=255)
                                                                  cATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTa    <  1:3036159/71‑1 (MQ=255)
                                                                  cATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTa    <  1:3602948/71‑1 (MQ=255)
                                                                    tGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACg  <  1:2023850/71‑1 (MQ=255)
                                                                    tGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACg  <  1:2020927/71‑1 (MQ=255)
                                                                    |                                                                      
GTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACG  >  minE/385748‑385886

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