Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,026,302 C→A 10.8% G190G (GGG→GGT malI ← DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,026,3020CA10.8% 168.4 / 4.0 54G190G (GGG→GGTmalIDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (30/18);  new base A (4/2);  total (34/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTG  >  minE/1026233‑1026370
                                                                     |                                                                    
cAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGcc                                                                     <  1:2229401/71‑1 (MQ=255)
cAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGcc                                                                     <  1:1309172/71‑1 (MQ=255)
  ggCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCt                                                                   >  1:1663103/1‑71 (MQ=255)
   gCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTc                                                                  >  1:1503653/1‑71 (MQ=255)
       aaTTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCccag                                                              >  1:789849/1‑71 (MQ=255)
       aaTTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCccag                                                              >  1:1789455/1‑71 (MQ=255)
        aTTTTAGTAGAGTTGCACAATAGCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCccagc                                                             >  1:2010586/1‑70 (MQ=255)
              gtagAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGACCTCCCAGCCAgg                                                         <  1:2500748/69‑1 (MQ=255)
                   gTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCg                                                  >  1:2083544/1‑71 (MQ=255)
                   gTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGACCTCCCAGCCAGGCGATCCg                                                  >  1:1285583/1‑71 (MQ=255)
                      gCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTc                                                                  >  1:1283890/1‑52 (MQ=255)
                      gCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCt                                                >  1:2054008/1‑70 (MQ=255)
                      gCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCt                                                <  1:1412239/70‑1 (MQ=255)
                       cacaATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGACCTCCCAGCCAGGCGATCCGCTg                                               <  1:2742208/70‑1 (MQ=255)
                        acaATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGAt                                             <  1:1922430/71‑1 (MQ=255)
                          aaTAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCccag                                                              <  1:1919419/52‑1 (MQ=255)
                          aaTAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGACCTCCCAGCCAGGCGATCCGCTGATGc                                           >  1:2603204/1‑71 (MQ=255)
                              gCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGc                                                 >  1:3491560/1‑61 (MQ=255)
                                  ccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCg                                   <  1:1238507/71‑1 (MQ=255)
                                   cAACCCGTTCTGCACGGGTTAATGAGGAACTTTGc                                                                      >  1:1221726/1‑35 (MQ=255)
                                   cAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCg                                   >  1:424506/1‑70 (MQ=255)
                                   cAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCg                                   >  1:3583973/1‑70 (MQ=255)
                                   cAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCg                                   >  1:3365254/1‑70 (MQ=255)
                                   cAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGa                                  <  1:2172825/71‑1 (MQ=255)
                                    aaCCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGaa                                 >  1:1184902/1‑71 (MQ=255)
                                       ccGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTg                                               <  1:836067/54‑1 (MQ=255)
                                        cGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATgag                             <  1:1937077/71‑1 (MQ=255)
                                        cGTTCTGCACGGGTTAATGAGGAACTTTGACCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAAt                                >  1:2395582/1‑68 (MQ=255)
                                         gTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAAt                                >  1:3501025/1‑67 (MQ=255)
                                          ttCTGCACGGGTTAATGAGGAACTTTGCCCTCccagc                                                             <  1:238771/37‑1 (MQ=255)
                                            cTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGAGATCCGCTGATGCCCATTGCGAATGAGATGc                         >  1:3865367/1‑71 (MQ=255)
                                             tGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGa                                              >  1:3628231/1‑49 (MQ=255)
                                              gCACGGGTTAATGAGGAACTTTGCCCTCccagccag                                                          >  1:3129795/1‑36 (MQ=255)
                                               cACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTcc                      >  1:3382371/1‑71 (MQ=255)
                                               cACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTcc                      >  1:3386974/1‑71 (MQ=255)
                                                  gggTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAAt                                <  1:507955/58‑1 (MQ=255)
                                                      tAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCt                                                >  1:381423/1‑38 (MQ=255)
                                                      tAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTcc                      >  1:1305072/1‑64 (MQ=255)
                                                      tAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTcaac                >  1:3522542/1‑70 (MQ=255)
                                                      tAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTcaac                >  1:2251922/1‑70 (MQ=255)
                                                         tGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTcaac                <  1:3455317/67‑1 (MQ=255)
                                                          gAGGAACTTTGCCCTCGCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACt            <  1:1583639/70‑1 (MQ=255)
                                                          gAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTcaac                <  1:3694202/66‑1 (MQ=255)
                                                           aGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACtg           <  1:1305170/70‑1 (MQ=255)
                                                            ggAACTTTGCCCTCCCATCCAGGCGATCCGCTGATGCCCAtt                                      >  1:3575082/1‑42 (MQ=255)
                                                            ggAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCAtt                                      >  1:2021687/1‑42 (MQ=255)
                                                            ggAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACtgtg         <  1:252787/71‑1 (MQ=255)
                                                            ggAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACtgtg         <  1:1867891/71‑1 (MQ=255)
                                                             taaCTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTcaacaa              <  1:2130884/64‑1 (MQ=255)
                                                             gAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTc                   <  1:3535180/60‑1 (MQ=255)
                                                              aaCTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACtgtg         >  1:1839693/1‑69 (MQ=255)
                                                               aCTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGAt                           >  1:2788995/1‑50 (MQ=255)
                                                                 tttGACCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGcc    >  1:3877911/1‑71 (MQ=255)
                                                                  ttGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTcaa                 >  1:1371500/1‑57 (MQ=255)
                                                                  ttGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTcaa                 >  1:3582647/1‑57 (MQ=255)
                                                                   tGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTg  >  1:1416591/1‑71 (MQ=255)
                                                                    gCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTg  <  1:1203645/70‑1 (MQ=255)
                                                                     |                                                                    
CAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTG  >  minE/1026233‑1026370

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: