Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,569,867 C→T 9.7% T46M (ACG→ATG) ‡ yibI → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,569,8670CT 9.7% 217.4 / 10.9 72T46M (ACG→ATG) ‡yibIpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (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 = 7.85e-01

TCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTC  >  minE/2569798‑2569934
                                                                     |                                                                   
tcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATAcg                                                                    <  1:3889651/71‑1 (MQ=255)
 cgcgATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATAcgg                                                                   >  1:3130024/1‑71 (MQ=255)
  gcgATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATAcggc                                                                  >  1:2063193/1‑71 (MQ=255)
  gcgATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATAcggc                                                                  >  1:2343356/1‑71 (MQ=255)
   cgATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATAcggcg                                                                 <  1:1529839/71‑1 (MQ=255)
        cATGAGATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTg                                                            >  1:2952413/1‑71 (MQ=255)
        cATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCt                                                             <  1:3340647/70‑1 (MQ=255)
        cATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTg                                                            >  1:1367069/1‑71 (MQ=255)
        cATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTg                                                            >  1:3311068/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:1446602/71‑1 (MQ=255)
                           ttGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTg                                         >  1:1638252/1‑71 (MQ=255)
                           ttGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTg                                         <  1:514100/71‑1 (MQ=255)
                            tGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATAc                                                                     >  1:2306928/1‑42 (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:3910992/1‑71 (MQ=255)
                                                  cATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCtgtg                  >  1:1575045/1‑71 (MQ=255)
                                                  cATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCtgtg                  >  1:1727044/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:2339690/1‑71 (MQ=255)
                                                      ccGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATc              >  1:210277/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:2201061/1‑68 (MQ=255)
                                                         aCGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATc              >  1:3326168/1‑68 (MQ=255)
                                                         aCGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATc              >  1:1413131/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:3525466/1‑71 (MQ=255)
                                                              aTTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGAc      >  1:4000188/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:3424403/1‑51 (MQ=255)
                                                                     cggcggGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCtgt                   >  1:2853821/1‑51 (MQ=255)
                                                                     cggcggGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCtgt                   >  1:3263133/1‑51 (MQ=255)
                                                                     |                                                                   
TCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCTGTTTACTCTGCATGTTATCTGGCCGTTTCTGTGGATCTGGGCGACGCTC  >  minE/2569798‑2569934

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: