Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,569,883 C→T 10.9% S51S (AGC→AGT yibI → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,569,8830CT10.9% 212.4 / 14.0 73S51S (AGC→AGTyibIpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (34/31);  new base T (3/5);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGC  >  minE/2569816‑2569947
                                                                   |                                                                
cTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTg                                                               <  1:2881446/71‑1 (MQ=255)
 ttATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGt                                                              >  1:197004/1‑71 (MQ=255)
   aTCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGtt                                                             >  1:990417/1‑70 (MQ=255)
    tCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGttt                                                            <  1:1208724/70‑1 (MQ=255)
     cTGATTGCCAAAAAACGTAATCATCCTCATGCCGACGCCATTCATGTTGCTGGTTGGGTGAGTCTTTTTAc                                                          <  1:1610422/71‑1 (MQ=255)
      tGATTGCCAAAAAACGTAATCATCCTCATGCCGACGCCATTCATGTTGCTGGTTGGGTGAGTCTTTTTAc                                                          >  1:1228502/1‑70 (MQ=255)
      tGATTGCCAAAAAACGTAATCATCCTCATGCCGACGCCATTCATGTTGCTGGTTGGGTGAGTCTTTTTAc                                                          >  1:1917940/1‑70 (MQ=255)
        aTTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTa                                                           <  1:3277672/67‑1 (MQ=255)
         ttGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTg                                                      <  1:1446602/71‑1 (MQ=255)
         ttGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTg                                                      <  1:514100/71‑1 (MQ=255)
         ttGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTg                                                      >  1:1638252/1‑71 (MQ=255)
           gCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCa                                                    >  1:2509197/1‑71 (MQ=255)
             cAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCAt                                                   <  1:2097289/70‑1 (MQ=255)
               aaaaGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATg                                                  <  1:2705871/69‑1 (MQ=255)
               aaaaGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGtt                                                <  1:1744519/71‑1 (MQ=255)
                 aaGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTa                                               >  1:2790158/1‑70 (MQ=255)
                  cgcgcAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATc                                             <  1:2449483/70‑1 (MQ=255)
                     gcAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGtt                                                             <  1:521538/52‑1 (MQ=255)
                        atcATCCTCATGCCGACGCCATTCATGTTGCTGGTTGGGTGAGTCTTTTTACGTTACATGTTATCTGGCCg                                       >  1:2212166/3‑71 (MQ=255)
                          cATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGt                                      <  1:3314130/70‑1 (MQ=255)
                            tCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTc                                   <  1:1419141/71‑1 (MQ=255)
                              cccATGCCGACGCTATTCATACGGCGGTCTGGGTGAGCCTGTTTACTCTGCATGtt                                                >  1:3788282/1‑56 (MQ=255)
                              cccATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCt                                  <  1:708545/70‑1 (MQ=255)
                               tcATGCCGACGCCATTCATGTTGCTGGTTGGGTGAGTCTTTTTACGTTACATGt                                                 <  1:3151205/53‑1 (MQ=255)
                               ccATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTc                                   <  1:939965/68‑1 (MQ=255)
                                cATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGttt                                    >  1:4115960/1‑66 (MQ=255)
                                cATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCtgtg                               >  1:3910992/1‑71 (MQ=255)
                                cATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCtgtg                               >  1:1727044/1‑71 (MQ=255)
                                cATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCtgtg                               >  1:1575045/1‑71 (MQ=255)
                                cATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCtgt                                >  1:2191511/1‑70 (MQ=255)
                                aaTGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCtgtg                               >  1:502231/2‑71 (MQ=255)
                                 aTGCCGACGCTATTCATACGGCGGGCTGGGTTAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCtgtg                               >  1:194766/1‑70 (MQ=255)
                                 aTGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGt                                                 <  1:1738963/52‑1 (MQ=255)
                                 aTGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTgg                              <  1:427964/71‑1 (MQ=255)
                                 aTGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTgg                              <  1:633390/71‑1 (MQ=255)
                                  tGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCtgtg                               >  1:1332663/1‑69 (MQ=255)
                                    ccGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGttt                                    >  1:1139567/1‑62 (MQ=255)
                                    ccGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATc                           >  1:210277/1‑71 (MQ=255)
                                    ccGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATc                           >  1:2339690/1‑71 (MQ=255)
                                     cGACTCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCt                          <  1:2711162/71‑1 (MQ=255)
                                      gACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCa                                                    >  1:2551776/1‑44 (MQ=255)
                                      gACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTg                         <  1:225772/71‑1 (MQ=255)
                                      gACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTg                         <  1:2407866/71‑1 (MQ=255)
                                       aCGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATc                           >  1:1413131/1‑68 (MQ=255)
                                       aCGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATc                           >  1:3326168/1‑68 (MQ=255)
                                       aCGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATc                           >  1:2201061/1‑68 (MQ=255)
                                       aCGCCATTCATGTTGCTGGTTGGGTGAGTCTTTTTACGTTACATGTTATCTGGCCGTTTTTGTGGATTTg                         <  1:664420/70‑1 (MQ=255)
                                       aCGCCATTCATGTTGCTGGTTGGGTGAGTCTTTTTACGTTACATGTTATCTGGCCGTTTTTGTGGATTTg                         <  1:2047322/70‑1 (MQ=255)
                                         gCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCtgtgag                             <  1:2844112/64‑3 (MQ=255)
                                         gCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGc                      >  1:1555806/1‑71 (MQ=255)
                                         gCTATTCATACGGCGGGCTGGCGGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCtgtgag                             <  1:811604/64‑3 (MQ=255)
                                           tATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGt                                                 <  1:599802/42‑1 (MQ=255)
                                            aTTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGa                             <  1:4035978/61‑1 (MQ=255)
                                            aTTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGAc                   >  1:4000188/1‑71 (MQ=255)
                                            aTTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGAc                   >  1:3525466/1‑71 (MQ=255)
                                             ttCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTg                         <  1:2566614/64‑1 (MQ=255)
                                               cATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCa                                                    >  1:307500/1‑35 (MQ=255)
                                               cATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCa                                                    >  1:1275227/1‑35 (MQ=255)
                                               cATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGct                >  1:2754663/1‑71 (MQ=255)
                                               cATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGc                 >  1:670692/1‑70 (MQ=255)
                                                aTACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGctc               >  1:650214/1‑71 (MQ=255)
                                                aTACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGctc               <  1:1827599/71‑1 (MQ=255)
                                                 tACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTgg                              <  1:1251511/55‑1 (MQ=255)
                                                   cggcggGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCtgt                                >  1:3263133/1‑51 (MQ=255)
                                                   cggcggGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCtgt                                >  1:3424403/1‑51 (MQ=255)
                                                   cggcggGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCtgt                                >  1:2853821/1‑51 (MQ=255)
                                                      cggGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGct                <  1:1708824/64‑1 (MQ=255)
                                                      cggGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGct                <  1:1857656/64‑1 (MQ=255)
                                                       gggCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAAc        <  1:222475/71‑1 (MQ=255)
                                                         gttGGGTGAGTCTTTTTACGTTACATGTTATCTGGCCGTTTTTGTGGATTTGGGCCACGctt               <  1:1649038/60‑2 (MQ=255)
                                                          cTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATc                           >  1:42408/1‑49 (MQ=255)
                                                           tGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGa                    <  1:345542/55‑1 (MQ=255)
                                                             ggTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGc  >  1:1226190/1‑71 (MQ=255)
                                                                   |                                                                
CTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGC  >  minE/2569816‑2569947

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: