Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 703,834 A→C 5.3% F271V (TTC→GTC)  ymcA ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE703,8340AC 5.3% 198.7 / 4.1 76F271V (TTC→GTC) ymcAconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (45/27);  new base C (2/2);  total (47/29)
Fisher's exact test for biased strand distribution p-value = 6.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

ATATGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGCTGGCAGCGCGATAAAC  >  minE/703769‑703903
                                                                 |                                                                     
atatGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCa                                                                   <  1:877951/70‑1 (MQ=255)
atatGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCa                                                                   <  1:4355681/70‑1 (MQ=255)
  atGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCaaa                                                                 >  1:5579130/1‑70 (MQ=255)
  atGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCaaa                                                                 >  1:2177605/1‑70 (MQ=255)
      cGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATc                                                             >  1:999812/1‑70 (MQ=255)
      cGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATc                                                             >  1:1110835/1‑70 (MQ=255)
      cGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCa                                                            <  1:1662944/71‑1 (MQ=255)
       ggTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCaa                                                           >  1:109960/1‑71 (MQ=255)
       ggTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCaa                                                           >  1:3381006/1‑71 (MQ=255)
              gCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACg                                                         <  1:5654407/66‑1 (MQ=255)
               cGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTAc                                                   <  1:2031798/71‑1 (MQ=255)
                 tATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGc                                                                     <  1:2895887/51‑1 (MQ=255)
                 tATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACgt                                                        <  1:2771650/64‑1 (MQ=255)
                 tATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACgt                                                        <  1:2692703/64‑1 (MQ=255)
                   tCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACcgcg                                               >  1:5338065/1‑71 (MQ=255)
                     gcgcAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGtt                                             <  1:2274816/71‑1 (MQ=255)
                     gcgcAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGtt                                             <  1:4107375/71‑1 (MQ=255)
                        cagggcagggCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCAt                                          <  1:4726849/71‑1 (MQ=255)
                         agggcagggcGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCAt                                          <  1:1171639/70‑1 (MQ=255)
                          gggcagggcGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATaa                                        >  1:1789026/1‑71 (MQ=255)
                          gggcagggcGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATaa                                        >  1:1433550/1‑71 (MQ=255)
                            gacgggcGCAGATCGTTAAAGTTGGTGCGCAACGTAACGCCAAACATAAAGGTGTTGCCACGTTCATAGCt                                      <  1:1149429/68‑1 (MQ=255)
                             cagggcGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACt                                      >  1:5634674/1‑70 (MQ=255)
                             acgggcGCAGATCGTTAAAGTTGGTGCGCAACGTAACGCCAAACATAAAGgtg                                                       >  1:3028354/3‑53 (MQ=255)
                              agggcGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACt                                      >  1:2925019/1‑69 (MQ=255)
                               gggcGCAGATCGTTAAAGTTGGTGCGCAACGTAACGCCAAACATaaagg                                                         <  1:3648392/49‑3 (MQ=255)
                                 gcgcAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTAcc                                                  >  1:3580591/1‑54 (MQ=255)
                                 gcgcAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTAcc                                                  >  1:2101668/1‑54 (MQ=255)
                                 gcgcAGATCGTTAAAGTTGGTGCGCAACGTAACGCCAAACATAAAGGTGTTGCCACGTTCATAGCTAAGGt                                 >  1:5142751/1‑71 (MQ=255)
                                   gcAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTg                               <  1:4706880/71‑1 (MQ=255)
                                    cAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACgtg                                                       <  1:3158008/46‑1 (MQ=255)
                                     aGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTgag                             <  1:3768694/71‑1 (MQ=255)
                                      gATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGtt                                                     <  1:2110083/46‑1 (MQ=255)
                                      gATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGtt                                >  1:1829395/1‑67 (MQ=255)
                                      gATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGtt                                >  1:2571267/1‑67 (MQ=255)
                                       aTCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGtt                                >  1:185409/1‑66 (MQ=255)
                                       aTCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACAGCGTTCATAACTCAGGtt                                >  1:5658889/1‑66 (MQ=255)
                                         cGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCaa                                                           >  1:1864379/1‑37 (MQ=255)
                                         cGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCaa                                                           >  1:628646/1‑37 (MQ=255)
                                         cGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAAc                                                          >  1:3437851/1‑38 (MQ=255)
                                         cGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAAc                                                          >  1:921970/1‑38 (MQ=255)
                                         cGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAgg                                  >  1:80547/1‑62 (MQ=255)
                                         cGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAgg                                  >  1:1087217/1‑62 (MQ=255)
                                         cGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAgg                                  >  1:2300236/1‑62 (MQ=255)
                                         cGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAgc                                  >  1:76829/1‑61 (MQ=255)
                                         cGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCt                         >  1:4514070/1‑71 (MQ=255)
                                          gTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTga                              >  1:4265829/1‑65 (MQ=255)
                                          gTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTga                              >  1:3238848/1‑65 (MQ=255)
                                            tGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGcc                      >  1:3713938/1‑71 (MQ=255)
                                            tGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGcc                      >  1:3591559/1‑71 (MQ=255)
                                             gAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTga                              >  1:2545410/1‑62 (MQ=255)
                                             gAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTga                              >  1:644927/1‑62 (MQ=255)
                                             gAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTga                              >  1:5049547/1‑62 (MQ=255)
                                                  tGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCAt                                          >  1:2886105/1‑45 (MQ=255)
                                                  tGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCAt                                          >  1:3567367/1‑45 (MQ=255)
                                                  tGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGCTg                >  1:3031926/1‑71 (MQ=255)
                                                       cGTAACGTGAAGCCAAACATCAACTTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGCTGGCAgc           >  1:2481980/1‑71 (MQ=255)
                                                       cGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGCTGGCAgc           >  1:1838828/1‑71 (MQ=255)
                                                         tAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACt                                      <  1:5390072/42‑1 (MQ=255)
                                                         tAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGCTGGCAgcgc         <  1:3595754/71‑1 (MQ=255)
                                                         tAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGCTGGCAgcg          >  1:4762669/1‑70 (MQ=255)
                                                          aaCGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACt                                      <  1:5655753/41‑1 (MQ=255)
                                                           aCGTGAAGCCAAACATCAACGTGTTACCGCGTTCAt                                          >  1:3557265/1‑36 (MQ=255)
                                                           aCGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAgg                                  >  1:2509495/1‑44 (MQ=255)
                                                           aCGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGc                  >  1:4430949/1‑60 (MQ=255)
                                                           aCGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGCTGGCAGCGCGa       >  1:5025772/1‑71 (MQ=255)
                                                           aCGAGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAg                   >  1:2075463/1‑59 (MQ=255)
                                                             gTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGCTGGCAGCGCtata     <  1:4420229/71‑1 (MQ=255)
                                                             gTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGCTGGCAGCGCGATa     <  1:4395318/71‑1 (MQ=255)
                                                             gTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGCTGGCAGCGCGATa     <  1:599358/71‑1 (MQ=255)
                                                             gTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGCTGGCAGCGCGATa     <  1:1676637/71‑1 (MQ=255)
                                                             gTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGCTGGCAGCGCGATa     <  1:877608/71‑1 (MQ=255)
                                                             gTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGCTGGCAGCGCGATa     <  1:4676954/71‑1 (MQ=255)
                                                              tGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGCTGGCAGCGCGATaa    <  1:3476720/71‑1 (MQ=255)
                                                              tGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGCTGGCAGCGCGATaa    <  1:853949/71‑1 (MQ=255)
                                                               gAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGCTGGCAGCGCGATaaa   >  1:254318/1‑71 (MQ=255)
                                                                aaGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGCTGGCAGCGCGATaaa   <  1:4267328/70‑1 (MQ=255)
                                                                aaGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGCTGGCAGCGCGATaaa   <  1:239660/70‑1 (MQ=255)
                                                                aaGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGCTGGCAGCGCGATAAAc  >  1:3654422/1‑71 (MQ=255)
                                                                 |                                                                     
ATATGCCGGTTTTGGCGTATCGCGCAGGGCAGGGCGCAGATCGTTGAAATTGGTCCGTAACGTGAAGCCAAACATCAACGTGTTACCGCGTTCATAACTCAGGTTGAGATCTGCCCAGCTGGCAGCGCGATAAAC  >  minE/703769‑703903

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: