Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,948,123 T→C 33.8% T40A (ACC→GCC)  yjjB ← conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,948,1230TC33.8% 73.7 / 59.7 65T40A (ACC→GCC) yjjBconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (22/21);  new base C (6/16);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 1.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTACCGCCGAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGTACCGACCAGCA  >  minE/2948055‑2948185
                                                                    |                                                              
ttttACCGCCGAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTa                                                               >  1:3825701/1‑70 (MQ=255)
     ccgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACt                                                         <  1:159483/71‑1 (MQ=255)
      cgccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACt                                                         <  1:3044534/70‑1 (MQ=255)
       gccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACttt                                                       >  1:2101633/1‑71 (MQ=255)
       gccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACttt                                                       >  1:2432974/1‑71 (MQ=255)
       gccgAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACtt                                                        >  1:345474/1‑70 (MQ=255)
            aaTCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACTTTCGGAt                                                  <  1:1824204/71‑1 (MQ=255)
              tCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACTTTCGGATgc                                                <  1:11549/71‑1 (MQ=255)
                   gCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACTTTCGGATGCGCCa                                            <  1:2823278/70‑1 (MQ=255)
                     ggTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGAt                                         <  1:2269627/71‑1 (MQ=255)
                       tataAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATAc                                       <  1:3009252/71‑1 (MQ=255)
                           aGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACTTTCGGa                                                   >  1:2110812/1‑55 (MQ=255)
                           aGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAgc                                   >  1:2315141/1‑71 (MQ=255)
                             ccGATATGCCCGGGAACATAGGGATGACGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATACCAgcgc                                 >  1:187783/1‑71 (MQ=255)
                             ccGATATGCCCGGGAACATAGGGATAGCGGCCGCCACGGTAAACACttt                                                       >  1:60743/1‑49 (MQ=255)
                             ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACttt                                                       >  1:3898108/1‑49 (MQ=255)
                             ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACttt                                                       >  1:3699182/1‑49 (MQ=255)
                             ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACttt                                                       >  1:3360509/1‑49 (MQ=255)
                             ccGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATACCAgcgc                                 >  1:4011684/1‑71 (MQ=255)
                                atatGCCCGGGAACATAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCa                                            <  1:646339/57‑1 (MQ=255)
                                atatGCCCGGGAACATAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCa                                            <  1:598439/57‑1 (MQ=255)
                                atatGCCCGGGAACATAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCa                                            <  1:2955900/57‑1 (MQ=255)
                                  atGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCa                            >  1:953412/1‑71 (MQ=255)
                                   tGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCAt                           <  1:3213613/71‑1 (MQ=255)
                                   tGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCAt                           <  1:1400295/71‑1 (MQ=255)
                                   tGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCAt                           <  1:719794/71‑1 (MQ=255)
                                    gCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCAtt                          >  1:2665559/1‑71 (MQ=255)
                                    gCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGAGACCAtt                          >  1:2949392/1‑71 (MQ=255)
                                     cccGGGAACATAGGGATAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCAtt                          <  1:3380490/70‑1 (MQ=255)
                                     cccGGGAACATAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTg                         <  1:1171236/71‑1 (MQ=255)
                                      ccGGGAACATAGGGATAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGa                        <  1:3699965/71‑1 (MQ=255)
                                      ccGGGAACATAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTg                         <  1:1166379/70‑1 (MQ=255)
                                         ggAACATAGGGATAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATa                     >  1:3605087/1‑71 (MQ=255)
                                         ggAACATAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGcc                                             <  1:1154073/47‑1 (MQ=255)
                                         ggAACATAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGcc                                             <  1:911940/47‑1 (MQ=255)
                                          gAACATAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGcc                                             <  1:2164667/46‑1 (MQ=255)
                                           aaCATAGGGATAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATAcc                   <  1:1733726/71‑1 (MQ=255)
                                            aCATAGTGATAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCa                  <  1:2635669/71‑1 (MQ=255)
                                            aCATAGGGCTAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCa                  <  1:4145236/71‑1 (MQ=255)
                                            aCATAGGGATAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCa                  <  1:2701189/71‑1 (MQ=255)
                                            aCATAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCa                                            <  1:115949/45‑1 (MQ=255)
                                             cATAGGGATAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCaa                 >  1:283662/1‑71 (MQ=255)
                                              aTAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCATATACCAGCGCGACCATTGAATACCAAt                <  1:1637545/71‑1 (MQ=255)
                                              aTAGGGATAACGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAAt                <  1:2337674/71‑1 (MQ=255)
                                                    aTAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGTAcc          >  1:3042815/1‑71 (MQ=255)
                                                    aTAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGTAc           >  1:3175010/1‑70 (MQ=255)
                                                     tAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATAcc                   <  1:2118135/61‑1 (MQ=255)
                                                     tAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATAcc                   <  1:2426758/61‑1 (MQ=255)
                                                     tAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGTACCg         <  1:2526961/71‑1 (MQ=255)
                                                       aCGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAAt                <  1:4114021/62‑1 (MQ=255)
                                                       aCGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGt             >  1:3727008/1‑65 (MQ=255)
                                                       aCGGCCGCCACGGCAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGTa            >  1:340326/1‑66 (MQ=255)
                                                         ggCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGTACCGAc       >  1:3697277/1‑69 (MQ=255)
                                                         ggCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGTACCGAc       >  1:2237284/1‑69 (MQ=255)
                                                         ggCCGCCACGGCAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGTACCGAc       >  1:3042671/1‑69 (MQ=255)
                                                           ccgccACGGCAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCaa                 <  1:1740002/57‑1 (MQ=255)
                                                            cgccACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGTACCGACCAGCa  <  1:2740876/71‑1 (MQ=255)
                                                            cgccACGGCAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGTACCGACCAGCa  <  1:1862890/71‑1 (MQ=255)
                                                            cgccACGGCAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGTACCGACCAGCa  <  1:538108/71‑1 (MQ=255)
                                                            cgccACGGCAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGTACCGACCAGCa  <  1:114952/71‑1 (MQ=255)
                                                             gccACGGCAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATg               >  1:932892/1‑57 (MQ=255)
                                                               cACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGTACCg         <  1:3428277/61‑1 (MQ=255)
                                                                 cGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGTACCGAc       >  1:531790/1‑61 (MQ=255)
                                                                 cGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGTACCGAc       >  1:1258136/1‑61 (MQ=255)
                                                                 cGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGTACCGAc       >  1:933366/1‑61 (MQ=255)
                                                                    |                                                              
TTTTACCGCCGAAATCATTGCGGTATAAGCCGATATGCCCGGGAACATAGGGATAACGGCCGCCACGGTAAACACTTTCGGATGCGCCAGATACCAGCGCGACCATTGAATACCAATGGTACCGACCAGCA  >  minE/2948055‑2948185

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: