Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,570,851 C→T 7.6% F252F (TTC→TTT)  yibH → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,570,8510CT 7.6% 237.8 / 7.0 79F252F (TTC→TTT) yibHhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (38/35);  new base T (4/2);  total (42/37)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

CGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGG  >  minE/2570782‑2570919
                                                                     |                                                                    
cGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCa                                                                     >  1:3808903/1‑71 (MQ=255)
 gTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTc                                                                      <  1:632607/69‑1 (MQ=255)
   tACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATc                                                                   >  1:2892101/1‑70 (MQ=255)
     cTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATcccc                                                                <  1:2316874/71‑1 (MQ=255)
     cTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATccc                                                                 <  1:3493113/70‑1 (MQ=255)
     cTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATccc                                                                 <  1:1413914/70‑1 (MQ=255)
     cTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATccc                                                                 <  1:315582/70‑1 (MQ=255)
       cAGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATGGTGTTTATACCCGa                                                              >  1:3094513/1‑71 (MQ=255)
        aGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATGGTGTTTATACCCGat                                                             >  1:1277992/1‑70 (MQ=255)
        aGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATGGTGTTTATACCCGat                                                             >  1:1626206/1‑70 (MQ=255)
        aGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAg                                                             >  1:3257161/1‑71 (MQ=255)
          gTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTc                                                                      >  1:2573784/1‑60 (MQ=255)
          gTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCa                                                           <  1:2080786/71‑1 (MQ=255)
           tACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaa                                                          <  1:3020262/71‑1 (MQ=255)
             cTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaaa                                                        <  1:2419635/71‑1 (MQ=255)
             cTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaa                                                         <  1:70292/70‑1 (MQ=255)
                aTCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACg                                                     <  1:660131/71‑1 (MQ=255)
                 tCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCGAAAACgg                                                    <  1:3952017/71‑1 (MQ=255)
                     cccAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATccc                                                                 <  1:1866897/54‑1 (MQ=255)
                     cccAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCaaa                                                <  1:3843591/71‑1 (MQ=255)
                     cccAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCaa                                                 <  1:2235528/70‑1 (MQ=255)
                      ccAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAAt                                               >  1:2956086/1‑71 (MQ=255)
                        aGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTATTCATCCCCGAGCAAAAACGGCAAATTg                                             <  1:1975554/71‑1 (MQ=255)
                          gTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTc                                           >  1:2704787/1‑71 (MQ=255)
                            aCATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGc                                         >  1:2167269/1‑71 (MQ=255)
                                aCGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCaa                                     >  1:1104172/1‑71 (MQ=255)
                                 cGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCa                                                           >  1:1540175/1‑48 (MQ=255)
                                 cGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAt                                    >  1:1085065/1‑71 (MQ=255)
                                 cGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAt                                    >  1:1774530/1‑71 (MQ=255)
                                  gCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCaa                                                 <  1:2368150/57‑1 (MQ=255)
                                  gCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAtt                                   <  1:3456558/71‑1 (MQ=255)
                                  gCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAtt                                   <  1:3594239/71‑1 (MQ=255)
                                  gCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAtt                                   <  1:722732/71‑1 (MQ=255)
                                   cAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAttt                                  <  1:2847195/71‑1 (MQ=255)
                                    aGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAttt                                  <  1:2036316/70‑1 (MQ=255)
                                    aGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTc                                 >  1:3330861/1‑71 (MQ=255)
                                       tGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaaaa                                                       >  1:1228381/1‑46 (MQ=255)
                                         ccTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaaaa                                                       >  1:1194273/1‑44 (MQ=255)
                                         ccTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaaa                                                        <  1:3066398/43‑1 (MQ=255)
                                         ccTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCa                             >  1:2607976/1‑70 (MQ=255)
                                          cTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaaa                           <  1:1956051/71‑1 (MQ=255)
                                          cTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaa                            >  1:1227510/1‑70 (MQ=255)
                                          cTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaa                            >  1:3238906/1‑70 (MQ=255)
                                           ttGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAAc                                                      <  1:3686817/43‑1 (MQ=255)
                                           ttGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaaaa                          >  1:2131217/1‑71 (MQ=255)
                                           ttGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaaaa                          >  1:369339/1‑71 (MQ=255)
                                            tGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAAc                         >  1:1613731/1‑71 (MQ=255)
                                            tGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAAc                         >  1:3042317/1‑71 (MQ=255)
                                             gCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACt                        >  1:1021290/1‑71 (MQ=255)
                                             gCCGCTACGTCCGGTGATGGTGTTTATACCCGATCAGAAACGACAAATCGTGGCGCAGTTCCGTCAGAACt                        >  1:958001/1‑71 (MQ=255)
                                              ccGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAAtt                                              >  1:1523492/1‑48 (MQ=255)
                                              ccGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTc                       >  1:2138665/1‑71 (MQ=255)
                                              ccGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTc                       <  1:1738223/71‑1 (MQ=255)
                                              ccGCTACGTCCGGTGATGGTGTTTATACCCGATCAGAAACGACAAATCGTGGCGCAGTTCCGTCAGAACTc                       <  1:2482413/71‑1 (MQ=255)
                                              acGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAAtt                                              >  1:983289/2‑48 (MQ=255)
                                               cGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCg                      >  1:441072/1‑71 (MQ=255)
                                               cGCTACGTCCGGTGATGGTGTTTATACCCGATCAGAAACGACAAATCGTGGCGCAGTTCCGTCAGAACTcc                      <  1:2263165/71‑2 (MQ=255)
                                                 cTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAAt                                               >  1:1356220/1‑44 (MQ=255)
                                                 cTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCt                    >  1:1415037/1‑71 (MQ=255)
                                                 cTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCt                    >  1:3393766/1‑71 (MQ=255)
                                                   gCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAt                                    <  1:1683833/53‑1 (MQ=255)
                                                    cGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTg                   <  1:1516537/69‑1 (MQ=255)
                                                    cGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGt                  >  1:1653882/1‑70 (MQ=255)
                                                    cGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGt                  >  1:1658581/1‑70 (MQ=255)
                                                      tCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTAc               <  1:255181/71‑1 (MQ=255)
                                                      tCCGGTGATGGTCTTCATCCCCGAGAAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTAc               >  1:2941404/1‑71 (MQ=255)
                                                       ccGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGc                                         >  1:3330836/1‑44 (MQ=255)
                                                       ccGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTAc               <  1:3767516/70‑1 (MQ=255)
                                                        cGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGt             <  1:2072660/71‑1 (MQ=255)
                                                        cGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGt             <  1:3473696/71‑1 (MQ=255)
                                                        aggTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGt                  <  1:2417400/65‑1 (MQ=255)
                                                          gTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaa                            >  1:2052476/1‑54 (MQ=255)
                                                          gTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCt           >  1:972126/1‑71 (MQ=255)
                                                          gTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCt           >  1:1661046/1‑71 (MQ=255)
                                                            gATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTc            <  1:3912398/68‑1 (MQ=255)
                                                            gATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTc            <  1:1366071/68‑1 (MQ=255)
                                                                gTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAAcc     >  1:2632932/1‑71 (MQ=255)
                                                                 tCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCt    <  1:388802/71‑1 (MQ=255)
                                                                 tCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCt    >  1:1281247/1‑71 (MQ=255)
                                                                 tCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCt    <  1:1480792/71‑1 (MQ=255)
                                                                    tCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTgg  >  1:2152602/1‑70 (MQ=255)
                                                                     |                                                                    
CGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCTGG  >  minE/2570782‑2570919

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: ‑