Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,613,318 A→G 7.0% F353S (TTC→TCC)  eutH ← predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,613,3180AG 7.0% 209.6 / 6.5 71F353S (TTC→TCC) eutHpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (37/29);  new base G (2/3);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

CCGACAATCATCGGGAAGATCATGGCGTTCATGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGACTTTGCCGCGGGTATCCATCTGCTTCATCATGCCG  >  minE/1613249‑1613386
                                                                     |                                                                    
ccGACAATCATCGGGAATATCATGGCGTTCATGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGa                                                                      >  1:3978735/1‑70 (MQ=255)
ccGACAATCATCGGGAAGATCATGGCGTTCATGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGa                                                                      >  1:3230653/1‑70 (MQ=255)
ccGACAATCATCGGGAAGATCATGGCGTTCATGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGa                                                                      >  1:764681/1‑70 (MQ=255)
ccGACAATCATCGGGAAGATCATGGCGTTCATGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGa                                                                      <  1:3578205/70‑1 (MQ=255)
ccGACAATCATCGGGAAGATCATGGCGTTCATGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGa                                                                      >  1:364070/1‑70 (MQ=255)
    cAATCATCGGGAAGATCATGGCGTTCATGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAagca                                                                 <  1:2954110/71‑1 (MQ=255)
     aaTCATCGGGAAGATCATGGCGTTCATGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAag                                                                   >  1:4403211/1‑68 (MQ=255)
         atcGGGAAGATCATGGCGTTCATGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGa                                                            <  1:34999/71‑1 (MQ=255)
            gggAAGATCATGGCGTTCATGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGaaa                                                          <  1:601292/70‑1 (MQ=255)
                aGATCATGGCGTTCATGTTGGCAGCGGCGAGGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCg                                                     <  1:4358247/71‑1 (MQ=255)
                aGATCATGGCGTTCATGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCg                                                     <  1:3170280/71‑1 (MQ=255)
                aGATCATGGCGTTCATGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCg                                                     <  1:29239/71‑1 (MQ=255)
                  aTCATGGCGTTCATGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGa                                                            >  1:4619239/1‑62 (MQ=255)
                  aTCATGGCGTTCATGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGa                                                            >  1:4235575/1‑62 (MQ=255)
                  aTCATGGCGTTCATGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGa                                                            >  1:4663427/1‑62 (MQ=255)
                  aTCATGGCGTTCATGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGa                                                            >  1:4212955/1‑62 (MQ=255)
                    cATGGCGTTCATGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAgg                                                 >  1:1870585/1‑71 (MQ=255)
                      tGGCGTTCATGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGgcg                                               >  1:3184626/1‑71 (MQ=255)
                        gCGTTCATGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGgcg                                               >  1:1777509/1‑69 (MQ=255)
                        gCGTTCATGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGgcg                                               >  1:3383079/1‑69 (MQ=255)
                           ttCATGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGtt                                          <  1:1813397/71‑1 (MQ=255)
                           ttCATGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGtt                                          <  1:3236263/71‑1 (MQ=255)
                           ttCATGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGtt                                          <  1:4633175/71‑1 (MQ=255)
                             cATGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGaaa                                                                    >  1:3544519/1‑43 (MQ=255)
                             cATGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGaaa                                                                    >  1:1673388/1‑43 (MQ=255)
                             cATGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGaaa                                                                    >  1:2309559/1‑43 (MQ=255)
                             cATGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGaaa                                                                    >  1:85991/1‑43 (MQ=255)
                             cATGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGGAAGCAGCGGAAACGGCGAAGGCGCAGTtga                                        >  1:2028695/1‑71 (MQ=255)
                              aTGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGgcg                                               <  1:1367794/63‑1 (MQ=255)
                              aTGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTtgat                                       >  1:2004945/1‑71 (MQ=255)
                              aTGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTtgat                                       >  1:1990977/1‑71 (MQ=255)
                              aTGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTtga                                        >  1:4105387/1‑70 (MQ=255)
                              aTGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTtga                                        <  1:2641312/70‑1 (MQ=255)
                                 ttGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGGAAGCAGCGGAAACGGCGAAGGCGCAGTtgatga                                     <  1:4184952/70‑1 (MQ=255)
                                 ttGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACg                                                        >  1:890587/1‑51 (MQ=255)
                                 ttGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGa                                                    >  1:2026665/1‑55 (MQ=255)
                                 ttGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGa                                                    >  1:2338373/1‑55 (MQ=255)
                                 ttGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGa                                                    >  1:2659576/1‑55 (MQ=255)
                                  tGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGACt                                   >  1:4195015/1‑71 (MQ=255)
                                  tGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGACt                                   >  1:880223/1‑71 (MQ=255)
                                  tGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGACt                                   >  1:4482328/1‑71 (MQ=255)
                                      aGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTtgat                                       <  1:3586298/63‑1 (MQ=255)
                                       gcggcgAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGAc                                    <  1:1003063/65‑1 (MQ=255)
                                       gcggcgAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGAc                                    <  1:1273993/65‑1 (MQ=255)
                                       gcggcgAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGAc                                    <  1:3506369/65‑1 (MQ=255)
                                        cggcgAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGtt                                          <  1:2549887/58‑1 (MQ=255)
                                        cggcgAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGtt                                          <  1:4178327/58‑1 (MQ=255)
                                        cggcgAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGtt                                          <  1:682950/58‑1 (MQ=255)
                                        cggcgAAGCCTAAGTGGGCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTtgat                                       <  1:1140925/61‑1 (MQ=255)
                                         ggcgAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGACTTTGCcgc                            >  1:1503552/1‑71 (MQ=255)
                                          gcgAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGACTTTGCcgcg                           >  1:650608/1‑71 (MQ=255)
                                              aGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGACTTTGCCGCGGGTa                       >  1:1269144/1‑71 (MQ=255)
                                              aGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGACTTTGCCGCGGGTa                       >  1:1325086/1‑71 (MQ=255)
                                                  tAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGaa                                                   <  1:2489488/39‑1 (MQ=255)
                                                  tAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGACTTTGCCGCGGGTATCCa                   >  1:1350836/1‑71 (MQ=255)
                                                    aGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGACTTTGCCGCGGGTATCCa                   <  1:3490095/69‑1 (MQ=255)
                                                    aGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGACTTTGCCGCGGGTATCCAt                  <  1:1323048/70‑1 (MQ=255)
                                                     gTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGa                                                    <  1:320304/35‑1 (MQ=255)
                                                       cgTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGACTTTGCCGCGGGTa                       >  1:2976175/2‑62 (MQ=255)
                                                          cGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGtt                                          <  1:4139895/40‑1 (MQ=255)
                                                           gCCCAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGACTTTGCCGCGGGTATCCATCTGCTtca          >  1:4348282/1‑71 (MQ=255)
                                                            cccAGGGCGAAAGCAGCGGAAACGGCGAAGGTGCAGTtga                                        >  1:2532848/1‑40 (MQ=255)
                                                            cccAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTtga                                        >  1:201562/1‑40 (MQ=255)
                                                             ccAGGGCGGAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGACTTTGCCGCGGGTATCCATCTGCTtcat         <  1:4075282/70‑1 (MQ=255)
                                                             ccAGGGCGGAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGACTTTGCCGCGGGTATCCATCTGCTtcat         <  1:1143531/70‑1 (MQ=255)
                                                             ccAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGACTTTGCcg                             >  1:4697004/1‑50 (MQ=255)
                                                             ccAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGACTTTGCGGCGGGTATCCATCTGCTtca          >  1:3112066/1‑69 (MQ=255)
                                                               aGGGCGGAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGACTTTGCCGCGGGTATCCATCTGCTtcatca       >  1:3318732/1‑70 (MQ=255)
                                                               aGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGACTTTGCCGCGGGTATCCATCTGCTtcatca       >  1:3857181/1‑70 (MQ=255)
                                                               aGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGACTTTGCCGCGGGTATCCATCTGCTTcatcat      >  1:3925957/1‑71 (MQ=255)
                                                               aGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGACTTTGCCGCGGGTATCCATCTGCTTcatcat      >  1:3673506/1‑71 (MQ=255)
                                                               aGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGACTTTGCCGCGGGTATCCATCTGCTTcatcat      >  1:2676068/1‑71 (MQ=255)
                                                               aGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGACTTTGCCGCGGGTATCCATCTGCTTcatcat      >  1:1113294/1‑71 (MQ=255)
                                                                gggCGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGAc                                    <  1:3466278/40‑1 (MQ=255)
                                                                 ggCGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGACTTTGCCGCGGGTATCCATCTGCTTCATCATGc    >  1:4259493/1‑71 (MQ=255)
                                                                 ggCGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGACTTTGCCGCGGGTATCCATCTGCTTCATCATGc    >  1:2529882/1‑71 (MQ=255)
                                                                   cGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGACTTTGCCGCGGGTATCCATCTGCTTCATCATGCCg  <  1:4283595/71‑1 (MQ=255)
                                                                   cGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGACTTTGCCGCGGGTATCCATCTGCTTCATCATGCCg  <  1:3370965/71‑1 (MQ=255)
                                                                   cGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGACTTTGCCGCGGGTATCCATCTGCTTCATCATGCCg  <  1:1933686/71‑1 (MQ=255)
                                                                   cGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGACTTTGCCGCGGGTATCCATCTGCTTCATCATGCCg  <  1:1339516/71‑1 (MQ=255)
                                                                     |                                                                    
CCGACAATCATCGGGAAGATCATGGCGTTCATGTTGGCAGCGGCGAAGCCTAAGTGGTCGCCCAGGGCGAAAGCAGCGGAAACGGCGAAGGCGCAGTTGATGACTTTGCCGCGGGTATCCATCTGCTTCATCATGCCG  >  minE/1613249‑1613386

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: