Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,569,922 C→T 12.3% I64I (ATC→ATT yibI → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,569,9220CT12.3% 193.8 / 14.8 73I64I (ATC→ATTyibIpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (32/32);  new base T (3/6);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGAAAGCGA  >  minE/2569854‑2569992
                                                                    |                                                                      
gACGCCATTCATGTTGCTGGTTGGGTGAGTCTTTTTACGTTACATGTTATCTGGCCGTTTTTGTGGATTTg                                                                      <  1:3418814/71‑1 (MQ=255)
gACGCCATTCATGTTGCTGGTTGGGTGAGTCTTTTTACGTTACATGTTATCTGGCCGTTTTTGTGGATTTg                                                                      <  1:2600880/71‑1 (MQ=255)
gACGCCATTCATGTTGCTGGTTGGGTGAGTCTTTTTACGTTACATGTTATCTGGCCGTTTTTGTGGATTTg                                                                      <  1:611473/71‑1 (MQ=255)
 aCGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTgg                                                                     >  1:1384040/1‑71 (MQ=255)
      aTTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGa                                                                 >  1:3498529/1‑70 (MQ=255)
       ttCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACg                                                               >  1:3337287/1‑71 (MQ=255)
       ttCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACg                                                               >  1:1590569/1‑71 (MQ=255)
          aTACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGc                                                                   <  1:3270038/64‑1 (MQ=255)
          aTACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGctc                                                            >  1:3766557/1‑71 (MQ=255)
          aTACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGctc                                                            <  1:381442/71‑1 (MQ=255)
           tACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGctct                                                           <  1:3087948/71‑1 (MQ=255)
            aCGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGctct                                                           >  1:3399315/1‑70 (MQ=255)
            aCGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTa                                                          >  1:2907231/1‑71 (MQ=255)
             cggcggGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGctc                                                            >  1:3786883/1‑68 (MQ=255)
             cggcggGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGctc                                                            >  1:840233/1‑68 (MQ=255)
             cggcggGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGctc                                                            >  1:1410986/1‑68 (MQ=255)
               gcggGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCa                                                       <  1:2277035/71‑1 (MQ=255)
                  ggCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAAcc                                                    <  1:1521807/71‑1 (MQ=255)
                   gttGGGTGAGTCTTTTTACGTTACATGTTATCTGGCCGTTTTTGTGGATTTGGGCCACGCTTTa                                                          >  1:1399135/3‑64 (MQ=255)
                     tGGGTGAGTCTTTTTACGTTACATGTTATCTGGCCGTTTTTGTGGATTTGGGCCACGCTTTa                                                          >  1:935435/1‑62 (MQ=255)
                     tGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCg                                                   >  1:2309242/1‑69 (MQ=255)
                      gggTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTg                                                                      <  1:727625/49‑1 (MQ=255)
                       ggTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAg                                                <  1:1096122/70‑1 (MQ=255)
                       ggTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGc                                               >  1:619740/1‑71 (MQ=255)
                        gTGAGTCTTTTTACGTTACATGTTATCTGGCCGTTTTTGTGGATTTGGGCCACGCTTTACCGCCCGGAGCg                                              <  1:791423/71‑1 (MQ=255)
                         tGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTAt                                                         >  1:1754319/1‑59 (MQ=255)
                         tGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGt                                             >  1:767200/1‑71 (MQ=255)
                           aGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTgg                                           >  1:3481944/1‑71 (MQ=255)
                           aGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTgg                                           >  1:2586938/1‑71 (MQ=255)
                              cTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTAt                                                         >  1:2444902/1‑54 (MQ=255)
                                gTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTa                                                          >  1:2757702/1‑51 (MQ=255)
                                 tttACTCTGCATGTTATCTGGCCGTTTCTTTGGATCTGGGCGACGCTCTATCAACCGGa                                                 <  1:2788372/59‑1 (MQ=255)
                                 tttACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCaa                                                      <  1:3919277/54‑1 (MQ=255)
                                 tttACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTgg                                       <  1:2841122/69‑1 (MQ=255)
                                 tttACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGCCGACGCTCTATCAACCGGAGCGTGGCTgg                                       <  1:3581725/69‑1 (MQ=255)
                                    aCTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTa                                   <  1:3701735/70‑1 (MQ=255)
                                     ctctGCATGTTATCTGGCCTTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTAt                                  <  1:1233029/70‑1 (MQ=255)
                                      tctGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGctc                                                            <  1:618787/43‑1 (MQ=255)
                                      tctGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCa                                                       <  1:1926712/48‑1 (MQ=255)
                                      tctGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTggg                                      <  1:1138669/65‑1 (MQ=255)
                                      tctGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATg                                 >  1:3079914/1‑70 (MQ=255)
                                       ttacATGTTATCTGGCCGTTTTTGTGGATTTGGGCCACGCTTTACCGCCCGGAGCGGGGATGGGGAATGCa                               <  1:429336/68‑1 (MQ=255)
                                       ctGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTggg                                      >  1:1977202/1‑64 (MQ=255)
                                       ctGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTggg                                      >  1:468772/1‑64 (MQ=255)
                                       ctGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTAt                                  <  1:40474/68‑1 (MQ=255)
                                        tGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTgggg                                     >  1:2762863/1‑64 (MQ=255)
                                        tGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGt                                    >  1:3541952/1‑65 (MQ=255)
                                        tGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGt                                    >  1:3665235/1‑65 (MQ=255)
                                        tGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAg                              >  1:3877369/1‑71 (MQ=255)
                                         gCATGTCATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGt                             >  1:513138/1‑71 (MQ=255)
                                         acATGTTATCTGGCCGTTTTTGTGGATTTGGGCCACGCTTTACCGCCCGGAGCGGGGATGGGGAATGCaa                              <  1:3222498/69‑2 (MQ=255)
                                          cATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAg                              <  1:2177747/69‑1 (MQ=255)
                                          cATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTc                            <  1:490593/71‑1 (MQ=255)
                                          cATGTTATCTGGCCGTTTCTGTGGATCTGGGCCACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAg                              <  1:1728481/69‑1 (MQ=255)
                                            tGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATTCAGTcac                          <  1:532216/71‑1 (MQ=255)
                                             gTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCt                                         <  1:838650/55‑1 (MQ=255)
                                              ttATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTcaca                         <  1:2107844/70‑1 (MQ=255)
                                               tATCTGGCCGTTTTTGTGGATTTGGGCCACGCTTTa                                                          >  1:2957339/1‑36 (MQ=19)
                                               tATCTGGCCGTTTCTGTGGATCTGGGCGACGCACTa                                                          >  1:549532/1‑36 (MQ=255)
                                                 tCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTg                                        <  1:2409752/52‑1 (MQ=255)
                                                  cTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTggg                                      >  1:1266829/1‑53 (MQ=255)
                                                   tGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTggg                                      >  1:1344991/1‑52 (MQ=255)
                                                     gCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCg                  <  1:308330/70‑1 (MQ=255)
                                                       cGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCg               >  1:2185154/1‑71 (MQ=255)
                                                         tttCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGt                      <  1:2909974/62‑1 (MQ=255)
                                                         tttCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGc                   <  1:110587/65‑1 (MQ=255)
                                                         tttCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGc                   <  1:2352216/65‑1 (MQ=255)
                                                           tCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAg            <  1:2412101/70‑1 (MQ=255)
                                                           tCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAg            <  1:2190818/70‑1 (MQ=255)
                                                                 gATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATg                       <  1:3112459/53‑1 (MQ=255)
                                                                 gATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGaaa      <  1:1331031/70‑1 (MQ=255)
                                                                 gATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGAAAg     >  1:3090772/1‑71 (MQ=255)
                                                                 gATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGAAAg     >  1:1316728/1‑71 (MQ=255)
                                                                    cTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGAAAGCGa  >  1:2259260/1‑71 (MQ=255)
                                                                    |                                                                      
GACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAACCGGAGCGTGGCTGGGGTATGCAGTCACATGTTGCGTCGCAGGAGAAAGCGA  >  minE/2569854‑2569992

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: