Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,569,871 G→T 9.5% A47A (GCG→GCT yibI → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,569,8710GT 9.5% 231.1 / 5.8 72A47A (GCG→GCTyibIpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (35/30);  new base T (3/4);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAAC  >  minE/2569801‑2569941
                                                                      |                                                                      
cgATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATAcggcg                                                                        <  1:1529839/71‑1 (MQ=255)
     cATGAGATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTg                                                                   >  1:2952413/1‑71 (MQ=255)
     cATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCt                                                                    <  1:3340647/70‑1 (MQ=255)
     cATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTg                                                                   >  1:3311068/1‑71 (MQ=255)
     cATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTg                                                                   >  1:1367069/1‑71 (MQ=255)
       tGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTggg                                                                 >  1:356224/1‑71 (MQ=255)
       tGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTggg                                                                 <  1:2255860/71‑1 (MQ=255)
        gACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGt                                                                >  1:2772236/1‑71 (MQ=255)
         aCATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGt                                                                >  1:4112618/1‑70 (MQ=255)
            tCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGt                                                                <  1:193167/67‑1 (MQ=255)
               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:1917940/1‑70 (MQ=255)
                     tGATTGCCAAAAAACGTAATCATCCTCATGCCGACGCCATTCATGTTGCTGGTTGGGTGAGTCTTTTTAc                                                    >  1:1228502/1‑70 (MQ=255)
                       aTTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTa                                                     <  1:3277672/67‑1 (MQ=255)
                        ttGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTg                                                >  1:1638252/1‑71 (MQ=255)
                        ttGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTg                                                <  1:1446602/71‑1 (MQ=255)
                        ttGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTg                                                <  1:514100/71‑1 (MQ=255)
                          gCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATAcggcg                                                                        >  1:2166476/1‑45 (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:1575045/1‑71 (MQ=255)
                                               cATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCtgtg                         >  1:1727044/1‑71 (MQ=255)
                                               cATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCtgtg                         >  1:3910992/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:2201061/1‑68 (MQ=255)
                                                      aCGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATc                     >  1:3326168/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:1827599/71‑1 (MQ=255)
                                                               aTACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGctc         >  1:650214/1‑71 (MQ=255)
                                                                tACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTgg                        <  1:1251511/55‑1 (MQ=255)
                                                                  cggcggGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCtgt                          >  1:2853821/1‑51 (MQ=255)
                                                                  cggcggGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCtgt                          >  1:3424403/1‑51 (MQ=255)
                                                                  cggcggGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCtgt                          >  1:3263133/1‑51 (MQ=255)
                                                                     cggGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGct          <  1:1857656/64‑1 (MQ=255)
                                                                     cggGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGct          <  1:1708824/64‑1 (MQ=255)
                                                                      gggCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAAc  <  1:222475/71‑1 (MQ=255)
                                                                      |                                                                      
CGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTCTATCAAC  >  minE/2569801‑2569941

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: