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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,4110GT13.7% 171.2 / 10.6 58G207G (GGG→GGTcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base G (24/26);  new base T (6/2);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 2.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

AAGCCCTGACTCATAGCCCGTGGGGAACATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTT  >  minE/385341‑385481
                                                                      |                                                                      
aaGCCCTGACTCATAGCCCGTGGGGAACATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATgggg                                                                        <  1:242675/71‑1 (MQ=255)
aaGCCCTGACTCATAGCCCGTGGGGAACATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATgggg                                                                        <  1:3131026/71‑1 (MQ=255)
aaGCCCTGACTCATAGCCCGTGGGGAACATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATgggg                                                                        <  1:2831227/71‑1 (MQ=255)
            aTAGCCCGTGGGGAACATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCtcc                                                              >  1:2972575/1‑67 (MQ=255)
            aTAGCCCGTGGGGAACATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTgc                                                              >  1:1020216/1‑69 (MQ=255)
                      gggAACATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTc                                                  <  1:1058188/71‑1 (MQ=255)
                         aaCATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCgg                                               >  1:92002/1‑71 (MQ=255)
                         aaCATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCgg                                               <  1:2936589/71‑1 (MQ=255)
                          aCATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCggg                                              <  1:3534153/71‑1 (MQ=255)
                          aCATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCggg                                              >  1:3784579/1‑71 (MQ=255)
                            atatACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGc                                            >  1:1397155/1‑71 (MQ=255)
                               taCCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTAc                                                                  >  1:2171876/1‑46 (MQ=255)
                               taCCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGt                                          <  1:1355792/70‑1 (MQ=255)
                               taCCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGt                                          <  1:2984431/70‑1 (MQ=255)
                               taCCGTTGCGTTCACCATTCCGCTCGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGt                                          <  1:3746699/70‑1 (MQ=255)
                                aCCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGTGATCTACCTGCGCTATCTGCGTCCGGGGCGTAt                                        >  1:1814932/1‑71 (MQ=255)
                                aCCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCt                                                                    <  1:3576874/43‑1 (MQ=255)
                                aCCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCt                                                                    <  1:1222103/43‑1 (MQ=255)
                                 ccGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTAt                                                         >  1:2354023/1‑53 (MQ=255)
                                 ccGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTAtt                                       >  1:3842965/1‑71 (MQ=255)
                                      gCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCg                                   <  1:618196/70‑1 (MQ=255)
                                       cGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTa                                                          <  1:128096/46‑1 (MQ=255)
                                         ttCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATc                                                        >  1:1176646/1‑46 (MQ=255)
                                           cACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCg                                   >  1:66222/1‑65 (MQ=255)
                                            aCCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAgtgt                             <  1:1025092/70‑1 (MQ=255)
                                            aCCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAgtgt                             <  1:2085597/70‑1 (MQ=255)
                                            aCCATTCCGCTGGCGCTGTTTAGGGGGATCTACCTGCGCTATCTGCGTcc                                                 <  1:2694765/50‑1 (MQ=255)
                                             ccATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCg                           <  1:3562301/71‑1 (MQ=255)
                                              cATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCt                                                           <  1:3481994/38‑1 (MQ=255)
                                              cATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCt                                                           <  1:3478850/38‑1 (MQ=255)
                                              cATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCgg                          <  1:1839672/71‑1 (MQ=255)
                                                ttCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGt                         >  1:1935080/1‑70 (MQ=255)
                                                 tCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCa                       >  1:1795004/1‑71 (MQ=255)
                                                   cGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATc                     <  1:2778612/71‑1 (MQ=255)
                                                   cGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATc                     <  1:2284494/71‑1 (MQ=255)
                                                     attgcgcTGTTTATGGGTATCTACATGCGCTTTATTCGTCCGGGGCGTGTGGGTGAAGTCTCTGTCATTgg                   >  1:3002723/4‑71 (MQ=255)
                                                     attgcgcTGTTTATGGGTATCTACATGCGCTTTATTCGTCCGGGGCGTGTGGGTGAAGTCTCTGTCATTgg                   >  1:3567367/4‑71 (MQ=255)
                                                      ttgcgcTGTTTATGGGTATCTACATGCGCTTTATTCGTCCGGGGCGTGTGGGTGAAGTCTCTGTCa                       <  1:188346/64‑1 (MQ=255)
                                                      tGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCgg                   >  1:957080/1‑70 (MQ=255)
                                                        gcgcTGTTTATGGGTATCTACATGCGCTTTATTCGTCCGGGGCGTGTGGGTGAAGTCTCTGTCAtt                     >  1:1990209/1‑65 (MQ=255)
                                                        gcgcTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTc                 >  1:3244522/1‑70 (MQ=255)
                                                        gcgcTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTc                 >  1:2757177/1‑70 (MQ=255)
                                                        gcgcTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGAGTATTGGCGAAGTGTCGGTCATCGGTc                 >  1:1934121/1‑70 (MQ=255)
                                                        gcgcTGTTTATGGGGATCTACCTGCACTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGt                  <  1:3056243/69‑1 (MQ=255)
                                                         cgcTGTTTATGGGTATCTACATGCGCTTTATTCGTCCGGGGCGTGTGGGTGAAGTCTCTGTCATTGGtatc               >  1:75250/1‑68 (MQ=255)
                                                         cgcTGTTTATGGGGGTCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTg               >  1:2216504/1‑71 (MQ=255)
                                                         cgcTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTg                                      >  1:195986/1‑48 (MQ=255)
                                                         cgcTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTg               <  1:467657/71‑1 (MQ=255)
                                                         cgcTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTg               <  1:1875192/71‑1 (MQ=255)
                                                         cgcTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTg               <  1:37189/71‑1 (MQ=255)
                                                          gcTGTTTATGGGTATCTACATGCGCTTTATTCGTCCGGGGCGTGTGGGTGAAGtctc                            <  1:2799748/57‑1 (MQ=255)
                                                          gcTGTTTATGGGTATCTACATGCGCTTTATTCGTCCGGGGCGTGTGGGTGAAGTCTCTGTCAttg                    >  1:2751892/1‑63 (MQ=255)
                                                              tttATGGGTATCTTCATGCGCTTTATTCGTCCGGGGCGTGTGGGTGAAgt                               >  1:707868/1‑50 (MQ=255)
                                                              tttATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCAt                      <  1:3002841/59‑1 (MQ=255)
                                                              tttATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCAt                      <  1:2042298/59‑1 (MQ=255)
                                                              tttATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCAt                      <  1:2046690/59‑1 (MQ=255)
                                                              tttATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTAt           >  1:3533357/1‑70 (MQ=255)
                                                               ttATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTc                            >  1:1104422/1‑52 (MQ=255)
                                                                tATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAg                                >  1:3270449/1‑47 (MQ=255)
                                                                 aTGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCt       >  1:3800526/1‑71 (MQ=255)
                                                                  tGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATTGGTCTGGTATTCCt       >  1:3697298/1‑70 (MQ=255)
                                                                     ggATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGAtt   <  1:335181/71‑1 (MQ=255)
                                                                      gATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGAttt  <  1:331608/71‑1 (MQ=255)
                                                                      |                                                                      
AAGCCCTGACTCATAGCCCGTGGGGAACATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTT  >  minE/385341‑385481

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: