Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,843 G→A 11.1% P351P (CCG→CCA cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,8430GA11.1% 173.7 / 10.3 63P351P (CCG→CCAcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base G (28/28);  new base A (2/5);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

CTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGT  >  minE/385775‑385911
                                                                    |                                                                    
cTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCAGGCCGaa                                                                    <  1:513094/71‑1 (MQ=255)
  gTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCg                                                                      >  1:871837/1‑67 (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:652499/71‑1 (MQ=255)
       tCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTg                                                             <  1:1647770/71‑1 (MQ=255)
          ccATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGc                                                               <  1:3153357/66‑1 (MQ=255)
            aTCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaa                                                        <  1:983120/71‑1 (MQ=255)
            aTCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaa                                                        <  1:609374/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)
                    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:969233/1‑71 (MQ=255)
                              tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCtt                                      >  1:2899356/1‑71 (MQ=255)
                              tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCtt                                      >  1:2903705/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:1662811/71‑1 (MQ=255)
                                             cTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTgg                       >  1:3016323/1‑71 (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:1591502/39‑1 (MQ=255)
                                               gATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacga                                                     <  1:88656/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)
                                                              cACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCt      <  1:3416980/71‑1 (MQ=255)
                                                               aCGCCGAAGATGGTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGa                     >  1:1521606/1‑55 (MQ=255)
                                                                cGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTc    >  1:2964384/1‑71 (MQ=255)
                                                                cGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTc    >  1:424096/1‑71 (MQ=255)
                                                                  ccGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGa                     <  1:46536/52‑1 (MQ=255)
                                                                  ccGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCg   <  1:3240639/70‑1 (MQ=255)
                                                                  ccGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGt  >  1:3810750/1‑71 (MQ=255)
                                                                    gAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATgg            <  1:1053513/59‑1 (MQ=255)
                                                                    |                                                                    
CTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGT  >  minE/385775‑385911

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: