Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,618 G→C 11.3% L276L (CTG→CTC cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,6180GC11.3% 138.3 / 8.1 53L276L (CTG→CTCcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base G (25/22);  new base C (2/4);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

GCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTT  >  minE/385548‑385686
                                                                      |                                                                    
gCGTGCAGTTGACCTGGATGCTGGTGGGTTACTGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTActgctg                                                                      >  1:3156345/1‑71 (MQ=255)
gCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTActgctg                                                                      >  1:3949890/1‑71 (MQ=255)
     cAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGccc                                                                  <  1:1764602/70‑1 (MQ=255)
      aGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCcg                                                                <  1:3206180/71‑1 (MQ=255)
        ttGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGATGGCCCcccg                                                              >  1:852744/1‑71 (MQ=255)
          gACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGt                                                             >  1:5671959/1‑70 (MQ=255)
          gACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGt                                                             >  1:1315308/1‑70 (MQ=255)
            ccTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGAc                                                          >  1:5090822/1‑71 (MQ=255)
               ggATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTAc                                                       >  1:4665770/1‑71 (MQ=255)
               ggATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCACCGCGTg                                                            >  1:5331024/1‑66 (MQ=255)
                gATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGcccc                                                                 >  1:4819517/1‑60 (MQ=255)
                gATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTAcc                                                      >  1:1069388/1‑71 (MQ=255)
                gATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTAcc                                                      >  1:4752566/1‑71 (MQ=255)
                gATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTAc                                                       <  1:4483431/70‑1 (MQ=255)
                 aTGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACct                                                     >  1:4723507/1‑71 (MQ=255)
                      ggtggGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTa                                                >  1:4706999/1‑71 (MQ=255)
                      ggtggGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTa                                                >  1:144673/1‑71 (MQ=255)
                       gtggGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTAc                                               >  1:2173621/1‑71 (MQ=255)
                        tggGTTACGGTTTTGTGGCAGCGGTGCTGCCGGTGTGGTTActgctg                                                                      <  1:2753046/47‑1 (MQ=255)
                           gTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTc                                           <  1:3594466/71‑1 (MQ=255)
                             tACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCt                                         >  1:4062551/1‑71 (MQ=255)
                             tACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCt                                         >  1:556715/1‑71 (MQ=255)
                              aCGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTAcc                                                      >  1:4165621/1‑57 (MQ=255)
                                ggTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGaa                                      <  1:833257/71‑1 (MQ=255)
                                 gTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGaaa                                     <  1:1130967/71‑1 (MQ=255)
                                    ttGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTcc                                          <  1:2931034/63‑1 (MQ=255)
                                    ttGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATc                                  >  1:3968050/1‑71 (MQ=255)
                                            gcggTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATc                                  <  1:5542898/63‑1 (MQ=255)
                                            gcggTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATc                                  <  1:2092252/63‑1 (MQ=255)
                                            gcggTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATc                                  <  1:3265248/63‑1 (MQ=255)
                                            gcggTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATc                                  <  1:5671157/63‑1 (MQ=255)
                                             cactgctgcCAGTGTGGCTGATCCTCGCACCGCGCGACTATCTGGCAACCTTCCTGaaaa                                    >  1:907676/4‑60 (MQ=255)
                                                tgctgcCAGTGTGGCTGATCCTCGCACCGCGCGACTATCTGGCAACCTTCCTGa                                       <  1:3358470/54‑1 (MQ=255)
                                                 gctgcCAGTGTGGCTGATCCTCGCACCGCGCGACTATCTGGCAACCTTCCTGAAAATCGGCGTTAtcgtcg                     >  1:3116366/1‑71 (MQ=255)
                                                  ctgcCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTa                           >  1:1879722/1‑64 (MQ=255)
                                                  ctgcCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCgg                    >  1:1053986/1‑71 (MQ=255)
                                                  ctgcCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCgg                    >  1:3163480/1‑71 (MQ=255)
                                                  ctgcCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCgg                    <  1:4568831/71‑1 (MQ=255)
                                                  ctgcCAGTGTGGCTGATCCTCGCACCGCGCGACTATCTGGCAACCTTCCTGAAAATCGGCGTTAtcgtcg                     <  1:5292467/70‑1 (MQ=255)
                                                  ctgcCAGTGTGGCTGATCCTCGCACCGCGCGACTATCTGGCAACCTTCCTGAAAATCGGCGTTAtcgtc                      <  1:1211863/69‑1 (MQ=255)
                                                  ctgcCAGTGTGGCTGATCCTCGCACCGCGCGACTATCTGGCAACCTTCCTGAAAATCGGCGTTAtcgtc                      <  1:5149104/69‑1 (MQ=255)
                                                     ccGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGa                                       <  1:1612282/49‑1 (MQ=255)
                                                     ccGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCgg                    <  1:5232363/68‑1 (MQ=255)
                                                     ccGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCgg                    <  1:1553514/68‑1 (MQ=255)
                                                       ggTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTg                <  1:1446727/70‑1 (MQ=255)
                                                         tgtgGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGc              <  1:5483196/70‑1 (MQ=255)
                                                         tgtgGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCg             <  1:457091/71‑1 (MQ=255)
                                                          gtgGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCgg            >  1:2327305/1‑71 (MQ=255)
                                                           tgGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGaa                                      <  1:2177273/44‑1 (MQ=255)
                                                             gTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATGGTCGGt                   <  1:5998/61‑1 (MQ=255)
                                                             gTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCTGTCTGGCGGTAg         <  1:1606371/71‑1 (MQ=255)
                                                             gTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAg         >  1:3797468/1‑71 (MQ=255)
                                                              ttACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAgg        >  1:2483854/1‑71 (MQ=255)
                                                                aCTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATc                                  >  1:5690135/1‑43 (MQ=255)
                                                                aCTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGAc                             >  1:2137469/1‑48 (MQ=255)
                                                                aCTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTg                <  1:5039293/61‑1 (MQ=255)
                                                                  tgctgGCCCCGCGTGACTACCTCTGTACCTTCCTGAAAATCGGGACTATCGTCGGt                   <  1:3512216/56‑1 (MQ=255)
                                                                  tgctgGCCCCGCGTGACTACCTCTCTACCTTCCTGaaa                                     <  1:1899124/38‑1 (MQ=255)
                                                                  tgctgGCCCCGCGTGACTACCTCTCTACCTTCCTGaaa                                     <  1:1940238/38‑1 (MQ=255)
                                                                  tgctgGCCCCGCGTGACTACCTCTCTACCTTCCTGaaa                                     <  1:4448642/38‑1 (MQ=255)
                                                                  tgctgGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGTGACTATCGTCGGt                   <  1:1516846/56‑1 (MQ=255)
                                                                  tgctgGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGt                   <  1:2109879/56‑1 (MQ=255)
                                                                    ctgGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCAtttt  >  1:2013277/1‑71 (MQ=255)
                                                                    ctgGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCAtt    <  1:4871700/69‑1 (MQ=255)
                                                                      |                                                                    
GCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTT  >  minE/385548‑385686

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: