Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,570,848 C→G 7.4% V251V (GTC→GTG)  yibH → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,570,8480CG 7.4% 236.3 / 7.0 81V251V (GTC→GTG) yibHhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (38/36);  new base G (4/2);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

GCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCT  >  minE/2570778‑2570917
                                                                      |                                                                     
gCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTc                                                                       >  1:3838290/1‑71 (MQ=255)
 cTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCt                                                                      <  1:2534870/71‑1 (MQ=255)
    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:315582/70‑1 (MQ=255)
         cTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATccc                                                               <  1:3493113/70‑1 (MQ=255)
         cTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATccc                                                               <  1:1413914/70‑1 (MQ=255)
           cAGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATGGTGTTTATACCCGa                                                            >  1:3094513/1‑71 (MQ=255)
            aGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATGGTGTTTATACCCGat                                                           >  1:1626206/1‑70 (MQ=255)
            aGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATGGTGTTTATACCCGat                                                           >  1:1277992/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:722732/71‑1 (MQ=255)
                                      gCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAtt                                 <  1:3456558/71‑1 (MQ=255)
                                      gCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAtt                                 <  1:3594239/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:369339/1‑71 (MQ=255)
                                               ttGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaaaa                        >  1:2131217/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:3393766/1‑71 (MQ=255)
                                                     cTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCt                  >  1:1415037/1‑71 (MQ=255)
                                                       gCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAAt                                  <  1:1683833/53‑1 (MQ=255)
                                                        cGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTg                 <  1:1516537/69‑1 (MQ=255)
                                                        cGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGt                >  1:1658581/1‑70 (MQ=255)
                                                        cGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGt                >  1:1653882/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:3473696/71‑1 (MQ=255)
                                                            cGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGt           <  1:2072660/71‑1 (MQ=255)
                                                            aggTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGt                <  1:2417400/65‑1 (MQ=255)
                                                              gTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCaa                          >  1:2052476/1‑54 (MQ=255)
                                                              gTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCt         >  1:1661046/1‑71 (MQ=255)
                                                              gTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCt         >  1:972126/1‑71 (MQ=255)
                                                                gATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTc          <  1:1366071/68‑1 (MQ=255)
                                                                gATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTc          <  1:3912398/68‑1 (MQ=255)
                                                                    gTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAAcc   >  1:2632932/1‑71 (MQ=255)
                                                                     tCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCt  >  1:1281247/1‑71 (MQ=255)
                                                                     tCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCt  <  1:388802/71‑1 (MQ=255)
                                                                     tCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCt  <  1:1480792/71‑1 (MQ=255)
                                                                      |                                                                     
GCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGGCAAATTGTCGCCCAATTTCGGCAAAACTCGCTGTTACGTCTGAAACCT  >  minE/2570778‑2570917

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