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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,8340GT17.6% 172.2 / 15.0 62G348G (GGG→GGTcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base G (27/24);  new base T (3/8);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 1.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGA  >  minE/385767‑385902
                                                                   |                                                                    
tCCCGTTCATGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGAc                                                                    >  1:2241936/1‑70 (MQ=255)
  ccGTTCCTGTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTAcga                                                                  <  1:2624561/70‑3 (MQ=255)
  ccGTTCCTGTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTAcga                                                                  <  1:198973/70‑3 (MQ=255)
   cGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCAc                                                                 >  1:3610705/1‑70 (MQ=255)
   cGTTCCTGTTGATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTacgacg                                                                >  1:3111843/1‑71 (MQ=255)
        cTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCAGGCCGaa                                                           <  1:513094/71‑1 (MQ=255)
          gTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCg                                                             >  1:871837/1‑67 (MQ=255)
          gTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTAcga                                                                  <  1:345865/62‑3 (MQ=255)
              aTCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCt                                                     <  1:2748866/71‑1 (MQ=255)
              aTCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGc                                                      <  1:3476404/70‑1 (MQ=255)
              aTCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCt                                                     >  1:793095/1‑71 (MQ=255)
               tCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTg                                                    <  1:1647770/71‑1 (MQ=255)
               tCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTg                                                    <  1:652499/71‑1 (MQ=255)
                  ccATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGc                                                      <  1:3153357/66‑1 (MQ=255)
                   cATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGGTGATCTCTTCTGGGACCACGc                                                               >  1:973459/1‑56 (MQ=255)
                    aTCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaa                                               <  1:609374/71‑1 (MQ=255)
                    aTCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaa                                               <  1:983120/71‑1 (MQ=255)
                     tCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaac                                              >  1:3689463/1‑71 (MQ=255)
                       gCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaatga                                            <  1:2135801/71‑1 (MQ=255)
                       gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacga                                            <  1:1849348/71‑1 (MQ=255)
                        ccTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaatgaa                                           <  1:2278416/71‑1 (MQ=255)
                        ccTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaatgaa                                           <  1:2030088/71‑1 (MQ=255)
                         cTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAg                                          >  1:2297925/1‑71 (MQ=255)
                         cTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAg                                          <  1:1613947/71‑1 (MQ=255)
                            tgGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaatgaaacc                                        >  1:3047714/1‑67 (MQ=255)
                            tgGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaatgaaac                                         <  1:3812158/69‑3 (MQ=255)
                            tgGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGc                                                               <  1:277114/47‑1 (MQ=255)
                            tgGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGc                                       <  1:2970180/71‑1 (MQ=255)
                              gcgcGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcga                                                <  1:910831/60‑1 (MQ=255)
                              gcgcGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAg                                     <  1:3158300/71‑1 (MQ=255)
                                   gtgtCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGt                                 >  1:2231017/1‑70 (MQ=255)
                                      tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCtt                             >  1:2903705/1‑71 (MQ=255)
                                      tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCtt                             >  1:2899356/1‑71 (MQ=255)
                                      tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCtt                             >  1:969233/1‑71 (MQ=255)
                                       cTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCtt                             >  1:3010349/1‑70 (MQ=255)
                                        tGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTg                                <  1:1400676/66‑1 (MQ=255)
                                        tGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTa                           >  1:611994/1‑71 (MQ=255)
                                          gCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTa                           <  1:3108896/69‑1 (MQ=255)
                                          gCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATc                         >  1:2382780/1‑71 (MQ=255)
                                           cTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTg                                >  1:3858316/1‑63 (MQ=255)
                                           cTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCg                        >  1:805865/1‑71 (MQ=255)
                                            ttCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcg                                                 >  1:1392380/1‑45 (MQ=255)
                                            ttCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCgg                       >  1:1162794/1‑71 (MQ=255)
                                              ccATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCTAACGAAGGGCAGGCGTGCTTTAt                          >  1:1428709/1‑66 (MQ=255)
                                              ccATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAgg                                         <  1:230622/51‑1 (MQ=255)
                                              ccATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGc                      >  1:819680/1‑70 (MQ=255)
                                                  gcgcTGATCTCTTCCGGTACGACGCCAAAACTGCTGGCTAATGAAACCGACGCGCGTTTCATCGGCTACgg                 >  1:3022534/1‑71 (MQ=255)
                                                    gcTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTg               <  1:3509131/71‑1 (MQ=255)
                                                     cTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTgg              >  1:3016323/1‑71 (MQ=255)
                                                     cTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTgg              <  1:1662811/71‑1 (MQ=255)
                                                     cTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTg               <  1:1583949/70‑1 (MQ=255)
                                                      tGATCTCTTCTGGGACCACTCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCttt                            <  1:2659821/56‑1 (MQ=255)
                                                      tGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCttt                            <  1:1582497/56‑1 (MQ=255)
                                                      tGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTgg              >  1:374929/1‑70 (MQ=255)
                                                       gATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacga                                            <  1:88656/39‑1 (MQ=255)
                                                       gATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacga                                            <  1:1591502/39‑1 (MQ=255)
                                                        aTCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTggg             <  1:3076262/69‑1 (MQ=255)
                                                         tctcTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTg               >  1:696242/1‑66 (MQ=255)
                                                             ttctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTaa      >  1:3139172/1‑71 (MQ=255)
                                                              tctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGt         >  1:515839/1‑67 (MQ=255)
                                                               ctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGc                               >  1:3580757/1‑44 (MQ=255)
                                                               ctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATg    >  1:2695806/1‑71 (MQ=255)
                                                                 gggACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATgg   >  1:1923530/1‑70 (MQ=255)
                                                                  ggACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATgg   <  1:2062311/69‑1 (MQ=255)
                                                                  ggACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGa  <  1:3520558/70‑1 (MQ=255)
                                                                   |                                                                    
TCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGA  >  minE/385767‑385902

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: