Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,834 G→T 6.5% G348G (GGG→GGT cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,8340GT 6.5% 253.2 / 2.5 76G348G (GGG→GGTcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base G (34/37);  new base T (2/3);  total (36/40)
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.35e-01

CCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAA  >  minE/385768‑385903
                                                                  |                                                                     
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:2736575/1‑70 (MQ=255)
                    tCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaat                                               <  1:2842094/71‑2 (MQ=255)
                    tCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaat                                               <  1:1062723/71‑2 (MQ=255)
                      gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                                              <  1:2907269/70‑1 (MQ=255)
                      gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                                              <  1:158737/70‑1 (MQ=255)
                      gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                                              >  1:3856479/1‑70 (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: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: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:2520274/1‑71 (MQ=255)
                                              cATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTa                     <  1:3036159/71‑1 (MQ=255)
                                              cATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTa                     <  1:3602948/71‑1 (MQ=255)
                                              cATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTa                     <  1:1995040/71‑1 (MQ=255)
                                                tGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACg                   <  1:2020927/71‑1 (MQ=255)
                                                tGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACg                   <  1:2023850/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:3324890/1‑71 (MQ=255)
                                                      gATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGa             >  1:1037604/1‑71 (MQ=255)
                                                      gATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGa             >  1:582615/1‑71 (MQ=255)
                                                      gATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGa             >  1:1857763/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)
                                                                 ggACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGaa  <  1:184914/71‑1 (MQ=255)
                                                                  |                                                                     
CCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAA  >  minE/385768‑385903

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: