Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,683,516 G→A 16.2% R486R (CGC→CGT yfhM ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,683,5160GA16.2% 102.4 / 8.5 61R486R (CGC→CGTyfhMconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (30/21);  new base A (5/5);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01

ACGGAGAATTTCACTTCATCTTTCGGCGTTAGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGATACTGATTATCGCCCGTGTTAGCG  >  minE/1683451‑1683582
                                                                 |                                                                  
acGGAGAATTTCACTTCATCTTTCGGCGTTAGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATg                                                                    >  1:1131661/1‑66 (MQ=255)
acGGAGAATTTCACTTCATCTTTCGGCGTTAGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATg                                                                    >  1:2673204/1‑66 (MQ=255)
acGGAGAATTTCACTTCATCTTTCGGCGTTAGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATg                                                                    >  1:1172914/1‑66 (MQ=255)
acGGAGAATTTCACTTCATCTTTCGGCGTTAGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATg                                                                    >  1:2712788/1‑66 (MQ=255)
acGGAGAATTTCACTTCATCTTTCGGCGTTAGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATg                                                                    >  1:2376981/1‑66 (MQ=255)
                              tgCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTc                                                <  1:302640/55‑1 (MQ=255)
                              aGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATgcg                                                                  >  1:335309/1‑38 (MQ=255)
                              aGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTc                                                <  1:1024534/56‑1 (MQ=255)
                              aGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTc                                                <  1:353002/56‑1 (MQ=255)
                              aGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTc                                                <  1:2019522/56‑1 (MQ=255)
                              aGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTc                                                <  1:554467/56‑1 (MQ=255)
                              aGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTc                                                <  1:564667/56‑1 (MQ=255)
                              aGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTc                                                <  1:68278/56‑1 (MQ=255)
                              aGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTc                                                <  1:239105/56‑1 (MQ=255)
                              aGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAAt                                     >  1:1854131/1‑67 (MQ=255)
                              aGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATcc                                   >  1:2030423/1‑69 (MQ=255)
                              aGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATcc                                   >  1:1692343/1‑69 (MQ=255)
                              aGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATcc                                   >  1:2014224/1‑69 (MQ=255)
                              aGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATcc                                   >  1:1807382/1‑69 (MQ=255)
                              aGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATcc                                   >  1:3267368/1‑69 (MQ=255)
                              aGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATcc                                   >  1:1747016/1‑69 (MQ=255)
                              aGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATcc                                   >  1:1444542/1‑69 (MQ=255)
                              aGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATcc                                   >  1:472036/1‑69 (MQ=255)
                              aGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATcc                                   >  1:1265058/1‑69 (MQ=255)
                              aGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATcc                                   >  1:682067/1‑69 (MQ=255)
                              aGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATcc                                   >  1:809347/1‑69 (MQ=255)
                              aGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATcc                                   >  1:892674/1‑69 (MQ=255)
                              aGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATcc                                   >  1:974456/1‑69 (MQ=255)
                              aGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGACATAAAATCTTCGACGTGGAAATcc                                   >  1:211726/1‑69 (MQ=255)
                              aGCGGGGTTTTCTCACCGGTCAGATTCACCGCCATGCGCTCTGGCATAAAATCTTc                                                <  1:2769336/56‑1 (MQ=255)
                                 ggggTTTTCTCACCGGTCAGATTCAGCGCCATACGCTCTGGCATAAAATCTTCGACGTGGAAATCCcac                                >  1:312921/1‑69 (MQ=255)
                                 ggggTTTTCTCACCGGTCAGATTCAGCGCCATACGCTCTGGCATAAAATCTTCGACGTGGAAATCCcac                                >  1:2579507/1‑69 (MQ=255)
                                 ggggTTTTCTCACCGGTCAGATTCAGCGCCATACGCTCTGGCATAAAATCTTCGACGTGGAAATCCcac                                >  1:682245/1‑69 (MQ=255)
                                 ggggTTTTCTCACCGGTCAGATTCAGCGCCATACGCTCTGGCATAAAATCTTCGACGTGGAAATCCcac                                >  1:1136886/1‑69 (MQ=255)
                                 ggggTTTTCTCACCGGTCAGATTCAGCGCCATACGCTCTGGCATAAAATCTTCGACGTGGAAATCCcac                                >  1:1149327/1‑69 (MQ=255)
                                       ttCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATCCcaca                               >  1:3410787/1‑64 (MQ=255)
                                       ttCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGa                          >  1:2554210/1‑69 (MQ=255)
                                       ttCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGa                          >  1:2378919/1‑69 (MQ=255)
                                       ttCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGa                          >  1:1334253/1‑69 (MQ=255)
                                       ttCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGa                          >  1:1654152/1‑69 (MQ=255)
                                                  cAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATcc                                   <  1:1180345/49‑1 (MQ=255)
                                                  cAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATcc                                   <  1:2967322/49‑1 (MQ=255)
                                                  cAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATcc                                   <  1:512755/49‑1 (MQ=255)
                                                  cAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGATACTGATTATc               >  1:902758/1‑69 (MQ=255)
                                                  cAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGATACTGATTATc               >  1:1258687/1‑69 (MQ=255)
                                                  cAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGATACTGATTATc               >  1:672441/1‑69 (MQ=255)
                                                  cAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGATACTGATTATc               >  1:391191/1‑69 (MQ=255)
                                                     aTTCAGCGCCATACGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGATACTGGTTATCGcc            <  1:2614284/69‑1 (MQ=255)
                                                     aTTCAGCGCCATACGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGATACTGGTTATCGcc            <  1:539612/69‑1 (MQ=255)
                                                     aTTCAGCGCCATACGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGATACTGGTTATCGcc            <  1:1755498/69‑1 (MQ=255)
                                                     aTTCAGCGCCATACGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGATACTGGTTATCGcc            <  1:3250546/69‑1 (MQ=255)
                                                     aTTCAGCGCCATACGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGATACTGGTTATCGcc            <  1:2346910/69‑1 (MQ=255)
                                                       tCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGATACTGATTATCGCCCg          <  1:1591397/69‑1 (MQ=255)
                                                        cAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGATACTGATTATCGCCCg          <  1:1408709/68‑1 (MQ=255)
                                                           cgcCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGATACTGATTATCGCCCGTGtt      <  1:1029146/69‑1 (MQ=255)
                                                           cgcCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGATACTGATTATCGCCCGTGtt      <  1:965285/69‑1 (MQ=255)
                                                           cgcCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGATACTGATTATCGCCCGTGtt      <  1:2653625/69‑1 (MQ=255)
                                                           cgcCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGATACTGATTATCGCCCGTGtt      <  1:2529474/69‑1 (MQ=255)
                                                           cgcCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGATACTGATTATCGCCCGTGtt      <  1:2250846/69‑1 (MQ=255)
                                                           cgcCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGATACTGATTATCGCCCGTGtt      <  1:1354410/69‑1 (MQ=255)
                                                           cgcCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGATACTGATTATCGCCCGTGtt      <  1:1794654/69‑1 (MQ=255)
                                                                tGCGCTCTGGCATAAATTCTTCGACGTGGAAATCCCACATCCGATACTGATTATCGCCCGTGTTAgcg  >  1:2873815/1‑68 (MQ=255)
                                                                tGCGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGATACTGATTATCGCCCGTGTTAgcg  >  1:27564/1‑68 (MQ=255)
                                                                tGCGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGATACTGATTATCGCCCGTGTTAgcg  >  1:2676600/1‑68 (MQ=255)
                                                                tGCGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGATACTGATTATCGCCCGTGTTAgcg  >  1:767240/1‑68 (MQ=255)
                                                                tGCGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGATACTGATTATCGCCCGTGTTAgcg  >  1:827304/1‑68 (MQ=255)
                                                                tGCGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGATACTGATTATCGCCCGTGTTAgcg  >  1:1133893/1‑68 (MQ=255)
                                                                 |                                                                  
ACGGAGAATTTCACTTCATCTTTCGGCGTTAGCGGGGTTTTCTCACCGGTCAGATTCAGCGCCATGCGCTCTGGCATAAAATCTTCGACGTGGAAATCCCACATCCGATACTGATTATCGCCCGTGTTAGCG  >  minE/1683451‑1683582

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: