Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,708,913 G→A 78.1% E54E (GAG→GAA yhfZ → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,708,9130GA78.1% 113.9 / 29.9 64E54E (GAG→GAAyhfZconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (7/7);  new base A (30/20);  total (37/27)
Fisher's exact test for biased strand distribution p-value = 5.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

AGATGAGCTGGCAAATGAATGCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAGGAATA  >  minE/2708847‑2708975
                                                                  |                                                              
aGATGAGCTGGCAAATGAATGCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAATc                                                              >  1:4853117/1‑69 (MQ=255)
aGATGAGCTGGCAAATGAATGCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAATc                                                              >  1:4782177/1‑69 (MQ=255)
aGATGAGCTGGCAAATGAATGCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAATc                                                              >  1:1432258/1‑69 (MQ=255)
 gATGAGCTGGCAAATGAATGCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAATCAAg                                                           >  1:3417358/1‑71 (MQ=255)
        tGGCAAATGAATGCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAATCAAGCgg                                                        <  1:847574/67‑1 (MQ=255)
         ggCAAATGAATGCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGa                                                       <  1:1448394/67‑1 (MQ=255)
             aatgaatgCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACg                                               <  1:2943138/71‑1 (MQ=255)
               tgaatgCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGa                                                       >  1:3459342/1‑61 (MQ=255)
                 aatgCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGAtt                                           >  1:3209483/1‑71 (MQ=255)
                 aatgCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGAtt                                           >  1:3036705/1‑71 (MQ=255)
                  atgCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTg                                          <  1:4978589/71‑1 (MQ=255)
                   tgCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTg                                          <  1:1501590/70‑1 (MQ=255)
                   tgCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTg                                          <  1:5660691/70‑1 (MQ=255)
                      cGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAAc                                       <  1:5201614/70‑1 (MQ=255)
                       gTTCATCCGTTGGCCTGACGCAGGCCGCTTTGAAAACGCTGGAATCAAGCgg                                                        <  1:295796/52‑1 (MQ=255)
                          cATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGt                                  >  1:1564634/1‑71 (MQ=255)
                          cATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGa                                                   >  1:1278612/1‑54 (MQ=255)
                          cATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGt                                  >  1:4033568/1‑71 (MQ=255)
                          cATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGt                                  >  1:1989860/1‑71 (MQ=255)
                           aTCCTTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGt                                  <  1:3333137/70‑1 (MQ=255)
                                ttGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGgcgc                            <  1:1010979/71‑1 (MQ=255)
                                ttGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGgcgc                            <  1:5140493/71‑1 (MQ=255)
                                ttGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGgcgc                            <  1:544019/71‑1 (MQ=255)
                                ttGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGgcg                             >  1:1410223/1‑70 (MQ=255)
                                 tGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGgcgc                            >  1:2489090/1‑70 (MQ=255)
                                 tGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGgcgc                            >  1:5332886/1‑70 (MQ=255)
                                  ggCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCaa                          >  1:2724693/1‑71 (MQ=255)
                                  ggCCTGACGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGt                                  >  1:1410442/1‑63 (MQ=255)
                                  ggCCTGACGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGt                                  >  1:1404947/1‑63 (MQ=255)
                                   gCCTGACGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAAt                         <  1:3925667/71‑1 (MQ=255)
                                    ccTGACGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATg                        <  1:5284244/71‑1 (MQ=255)
                                        aCGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGaa                                        >  1:3828186/1‑51 (MQ=255)
                                        aCGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAg                    >  1:4499871/1‑71 (MQ=255)
                                         cGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGt                   >  1:4935379/1‑71 (MQ=255)
                                         cGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGt                   >  1:203992/1‑71 (MQ=255)
                                          gCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGtt                  <  1:1168167/71‑1 (MQ=255)
                                          gCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGt                   >  1:1385710/1‑70 (MQ=255)
                                          gCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGt                   >  1:3402161/1‑70 (MQ=255)
                                          gCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGt                   >  1:2416082/1‑70 (MQ=255)
                                               ccGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTg             >  1:1183657/1‑71 (MQ=255)
                                               ccGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTg             >  1:5609510/1‑71 (MQ=255)
                                               ccGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTg             <  1:3074273/71‑1 (MQ=255)
                                                cgcgTTGAAAACGCTGGAATCAAGCGGAGCGATACGGAt                                            >  1:2138272/1‑39 (MQ=255)
                                                 gcgTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGAtt                                           >  1:2002797/1‑39 (MQ=255)
                                                 gcgTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAg                    >  1:3847029/1‑62 (MQ=255)
                                                 gcgTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTgg            <  1:4010145/70‑1 (MQ=255)
                                                  cgTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTc          <  1:2477216/71‑1 (MQ=255)
                                                     tGAAAACGCTGGAATCAAGCTGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTgg            >  1:5043449/1‑66 (MQ=255)
                                                     tGAAAACGCTGGAATCAAGCGGATCGATACGGATTGAACGCCGTGGgcgc                            <  1:609743/50‑1 (MQ=255)
                                                     tGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGaa                                        >  1:637766/1‑38 (MQ=255)
                                                     tGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGgcgc                            <  1:560475/50‑1 (MQ=255)
                                                      gAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGt                   >  1:3028858/1‑58 (MQ=255)
                                                      gAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTAt                >  1:4562422/1‑61 (MQ=255)
                                                      gAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTc          >  1:381545/1‑67 (MQ=255)
                                                      gAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAg       <  1:4815048/70‑1 (MQ=255)
                                                      gAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAg       <  1:3661627/70‑1 (MQ=255)
                                                      gAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAg       <  1:4808684/70‑1 (MQ=255)
                                                       aaaaCGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAgg      <  1:3347439/70‑1 (MQ=255)
                                                          aCGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAGGAata  <  1:1983704/71‑1 (MQ=255)
                                                          aCGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAGGAata  <  1:1965070/71‑1 (MQ=255)
                                                          aCGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAGGAata  <  1:5573913/71‑1 (MQ=255)
                                                           cGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAga      >  1:1499004/1‑65 (MQ=255)
                                                           cGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAga      >  1:2567475/1‑65 (MQ=255)
                                                            gCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTAt                >  1:4489624/1‑55 (MQ=255)
                                                                  |                                                              
AGATGAGCTGGCAAATGAATGCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAGGAATA  >  minE/2708847‑2708975

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: