Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,450,298 A→G 22.2% W78R (TGG→CGG)  yifE ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,450,2980AG22.2% 82.2 / 24.0 45W78R (TGG→CGG) yifEconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (19/16);  new base G (6/4);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.28e-01

ACCTGCGGTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGTTTCTCTTCCTC  >  minE/2450228‑2450367
                                                                      |                                                                     
aCCTGCGGTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACca                                                                       <  1:288825/71‑1 (MQ=255)
   tGCGGTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCgcac                                                                    >  1:1983996/1‑71 (MQ=255)
      ggTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCACAcgcg                                                                 <  1:1224786/71‑1 (MQ=255)
      ggTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCACAcgcg                                                                 <  1:1227252/71‑1 (MQ=255)
       gTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGt                                                                <  1:1742132/71‑1 (MQ=255)
         ttACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTc                                                              >  1:865472/1‑71 (MQ=255)
         ttACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTc                                                              >  1:1643955/1‑71 (MQ=255)
          tACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCACACGCGTTCt                                                             >  1:1889395/1‑71 (MQ=255)
            ccgccgGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGCCCACACGCGTTCTGc                                                           >  1:635406/1‑71 (MQ=255)
                 ggAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTg                                                      >  1:295147/1‑71 (MQ=255)
                    aaGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCa                                                   >  1:1423847/1‑71 (MQ=255)
                        gtgtGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACcac                                                                      >  1:907518/1‑48 (MQ=255)
                           tgAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTc                                            <  1:2117600/71‑1 (MQ=255)
                           tgAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTc                                            <  1:1451253/71‑1 (MQ=255)
                           tgAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCACACGCGt                                                                >  1:842358/1‑51 (MQ=255)
                           tgAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCACACGCGTTCTGCTTc                                                        >  1:61730/1‑59 (MQ=255)
                            gAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCa                                           >  1:143374/1‑71 (MQ=255)
                             aaaaCGTTTTGGACGCTTAATACGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACTGGCTCAc                                          >  1:2126602/1‑71 (MQ=255)
                             aaaaCGTTTTGGACGCTTAATACGCGTCATATACTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCAc                                          >  1:442841/1‑71 (MQ=255)
                             aaaaCGTTTTGGACGCTTAATACGCGTCATATACTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCAc                                          >  1:796128/1‑71 (MQ=255)
                             aaaaCGTTTTGGACGCTTAATACGCGTCATATACTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCAc                                          >  1:184510/1‑71 (MQ=255)
                               aaCGTTTTGGACGCTTAATACGCGTCATATACTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCACGt                                        <  1:327112/71‑1 (MQ=255)
                                       ggACGCTTAATACGCGTCATATACTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCAc                                 <  1:1703472/70‑1 (MQ=255)
                                        gACGCTTAATACGCGTCATATACTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACgg                               <  1:740765/71‑1 (MQ=255)
                                         aCGCTTAATACGCGTCATATACTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGc                              <  1:2163454/71‑1 (MQ=255)
                                         aCGCTTAATACGCGTCATATACTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGc                              <  1:301457/71‑1 (MQ=255)
                                          cGCTTAATACGCGTCATATACTTGGACCACAcgcg                                                                 <  1:1234381/35‑1 (MQ=255)
                                             ttAATACGCGTCATATACTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCACGt                                        >  1:187892/1‑57 (MQ=255)
                                             ttAATACGCGTCATATACTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCACGt                                        >  1:266182/1‑57 (MQ=255)
                                              tAATACGCGTCATATACTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCAc                                 <  1:954258/63‑1 (MQ=255)
                                               aaTACGCGTCATATACTTGGACCACACGCGTTCTGCTTCTGTc                                                    <  1:226029/43‑1 (MQ=255)
                                                aTACGCGTCATATACTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGc                       >  1:970498/1‑71 (MQ=255)
                                                aTACGCGTCATATACTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGc                       <  1:1021187/71‑1 (MQ=255)
                                                aTACGCGTCATATACTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGc                       <  1:354822/71‑1 (MQ=255)
                                                aTACGCGTCATATACTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGc                       >  1:1933489/1‑71 (MQ=255)
                                                 tACGCGTCATATACTTGGACCACACGCGTTCTGCTTCTGt                                                     <  1:883447/40‑1 (MQ=255)
                                                  aCGCGTCATATACTTGGACCGCACGCGTTCTGCTTCTGTCACt                                                 >  1:1948201/1‑43 (MQ=255)
                                                         atatACTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCAc                                 >  1:638525/1‑52 (MQ=255)
                                                          tataCTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGt             >  1:741195/1‑71 (MQ=255)
                                                           ataCTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCa                             <  1:612592/54‑1 (MQ=255)
                                                            taCTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCg                                     >  1:19166/1‑45 (MQ=255)
                                                            taCTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCg                                     >  1:812470/1‑45 (MQ=255)
                                                            taCTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCg                                     >  1:230607/1‑45 (MQ=255)
                                                            taCTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGttt           >  1:1941665/1‑71 (MQ=255)
                                                            taCTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGtt            >  1:263402/1‑70 (MQ=255)
                                                               ttGGACCACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGt             <  1:2128979/66‑1 (MQ=255)
                                                                 ggACCGCACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGTTTctct       <  1:48376/70‑1 (MQ=255)
                                                                 ggACCACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGTTtctc        <  1:150615/69‑1 (MQ=255)
                                                                     cacaCGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCa                             <  1:399534/44‑1 (MQ=255)
                                                                     cacaCGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCa                             <  1:1576301/44‑1 (MQ=255)
                                                                     cacaCGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCa                             <  1:363503/44‑1 (MQ=255)
                                                                     cacaCGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCa                             <  1:1532006/44‑1 (MQ=255)
                                                                     cacaCGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGTTTCTCTTCCTc  <  1:1756620/71‑1 (MQ=255)
                                                                      |                                                                     
ACCTGCGGTTTACCGCCGGAAAGGGTGTGAAAACGTTTTGGACGCTTAATACGCGTCATATACTTGGACCACACGCGTTCTGCTTCTGTCACTGGCTCACGTTCGCCACGGCATACTGCTACGAAGAGTTTCTCTTCCTC  >  minE/2450228‑2450367

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: