Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,831 T→C 9.2% S347S (TCT→TCC cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,8310TC 9.2% 179.1 / 9.3 65S347S (TCT→TCCcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base T (30/29);  new base C (3/3);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

ACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATG  >  minE/385761‑385900
                                                                      |                                                                     
aCCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCtcttct                                                                       >  1:3546339/1‑71 (MQ=255)
   tGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTgg                                                                     >  1:1109884/1‑70 (MQ=255)
   tGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTgg                                                                     >  1:1779544/1‑70 (MQ=255)
    gTTCCCGTTCCTGTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCggta                                                                   >  1:3889583/1‑69 (MQ=255)
       cccGTTCCTGTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTAcga                                                                >  1:3531833/1‑69 (MQ=255)
             ccTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCg                                                           >  1:3569403/1‑70 (MQ=255)
             ccTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCg                                                           >  1:1577627/1‑70 (MQ=255)
                 tttATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAg                                                        <  1:3996159/69‑1 (MQ=255)
                    aTCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCt                                                   <  1:483455/71‑1 (MQ=255)
                    aTCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGc                                                    <  1:439657/70‑1 (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)
                               cTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAg                                        >  1:77299/1‑71 (MQ=255)
                                tgtgGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacga                                          <  1:2566111/68‑1 (MQ=255)
                                tgtgGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAgg                                       >  1:2097676/1‑71 (MQ=255)
                                 gtgGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGAcc                                                                 >  1:742290/1‑44 (MQ=255)
                                 gtgGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGAcc                                                                 >  1:3757692/1‑44 (MQ=255)
                                   ggCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAggg                                      >  1:2921285/1‑69 (MQ=255)
                                   ggCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAggg                                      >  1:3465857/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:54322/1‑71 (MQ=255)
                                            tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCtt                           <  1:2943414/71‑1 (MQ=255)
                                              tGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGc                                                >  1:3900882/1‑48 (MQ=255)
                                              tGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTa                         >  1:1125927/1‑71 (MQ=255)
                                                gCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCttt                          <  1:1556019/68‑1 (MQ=255)
                                                gCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCttt                          <  1:1100498/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:3036159/71‑1 (MQ=255)
                                                     cATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTa                  >  1:2520274/1‑71 (MQ=255)
                                                     cATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTa                  <  1:1995040/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)
                                                           cTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTgg            <  1:3569800/71‑1 (MQ=255)
                                                           cTGATCTCTTCTGGGACCACGCCGAAGATGCTCGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTgg            <  1:1793739/71‑1 (MQ=255)
                                                            tGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTAc                 <  1:2580754/65‑1 (MQ=255)
                                                            tGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTAc                 <  1:2372572/65‑1 (MQ=255)
                                                             gATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGa          >  1:1037604/1‑71 (MQ=255)
                                                             gATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGa          >  1:1857763/1‑71 (MQ=255)
                                                             gATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGa          >  1:3324890/1‑71 (MQ=255)
                                                             gATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGa          >  1:582615/1‑71 (MQ=255)
                                                              aTCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                                           <  1:98470/37‑1 (MQ=255)
                                                              aTCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATc                       <  1:686651/57‑1 (MQ=255)
                                                              aTCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGa          <  1:1050874/70‑1 (MQ=255)
                                                              aTCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGa          <  1:2580709/70‑1 (MQ=255)
                                                                ctctTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGt              >  1:1168887/1‑64 (MQ=255)
                                                                ctctTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTgg            <  1:332704/66‑1 (MQ=255)
                                                                ctctTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATg        >  1:493619/1‑70 (MQ=255)
                                                                ctctTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGt       >  1:1090924/1‑71 (MQ=255)
                                                                 tcttctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGc                             <  1:2250624/48‑1 (MQ=255)
                                                                 tcttctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGc                             <  1:2158458/48‑1 (MQ=255)
                                                                  cttctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAgg                                  >  1:1055739/1‑42 (MQ=255)
                                                                  cttctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGtt      >  1:2907151/1‑70 (MQ=255)
                                                                   ttctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCtt                           <  1:3226865/48‑1 (MQ=255)
                                                                   ttctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCtt                           <  1:3231774/48‑1 (MQ=255)
                                                                   ttctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTATGGTGGGATGtt      <  1:1360947/69‑1 (MQ=255)
                                                                   ttctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGAt         <  1:794879/66‑1 (MQ=255)
                                                                   ttctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATg        <  1:1417964/67‑1 (MQ=255)
                                                                   ttctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGtt      <  1:502095/69‑1 (MQ=255)
                                                                   ttctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTaa    >  1:1959911/1‑71 (MQ=255)
                                                                     ctGGGACCACGCCGAAGATGCTGGCGCACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGttta    >  1:236578/1‑67 (MQ=255)
                                                                     ctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAAt   >  1:3802626/1‑70 (MQ=255)
                                                                     ctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATg  >  1:2662586/1‑71 (MQ=255)
                                                                      |                                                                     
ACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATG  >  minE/385761‑385900

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