Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,078,064 G→A 44.9% A151A (GCC→GCT)  uxaA ← altronate hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,078,0640GA44.9% 7.4 / 37.3 49A151A (GCC→GCT) uxaAaltronate hydrolase
Reads supporting (aligned to +/- strand):  ref base G (11/16);  new base A (5/17);  total (16/33)
Fisher's exact test for biased strand distribution p-value = 2.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.38e-01

GCTGTGAGCAGCCGTAGGTGTGGCTGAAGAGGAACACGCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGACACAGCCC  >  minE/2078012‑2078127
                                                    |                                                               
gCTGTGAGCAGCCGTAGGTGTGGCTGAAGAGGAACACGCCGTCAGTACCTTCGGCGTTGTTGGTCTCtt                                                 <  1:3106581/69‑1 (MQ=255)
                ggTGTGGCTGAAGAGGAACACGCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTg                                 >  1:2567235/1‑69 (MQ=255)
                  tgtgGCTGAAGAGGAACACGCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTgg                                <  1:2947421/68‑1 (MQ=255)
                   gtgGCTGAAGAGGAACACGCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTgg                                <  1:709345/67‑1 (MQ=255)
                              ggAACACGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGa                   >  1:3159527/1‑69 (MQ=255)
                              ggAACACGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGa                   >  1:1685572/1‑69 (MQ=255)
                              ggAACACGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGa                   >  1:1852938/1‑69 (MQ=255)
                              ggAACACGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGa                   >  1:2084722/1‑69 (MQ=255)
                              ggAACACGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGa                   >  1:3329225/1‑69 (MQ=255)
                               gAACACGCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTgg                                <  1:1841839/55‑1 (MQ=255)
                                 acacGCGGTCTGTACCTTCAGCGTTGTTGGTCTCTTTCTGGAAACGATTTTGGATCTGCCGCGCGATGc                <  1:3044091/69‑1 (MQ=255)
                                 acacGCCGTCGGTGCCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc                <  1:1798262/69‑1 (MQ=255)
                                 acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc                <  1:2067543/69‑1 (MQ=255)
                                 acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc                <  1:2901926/69‑1 (MQ=255)
                                 acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc                <  1:1039965/69‑1 (MQ=255)
                                 acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc                <  1:2997281/69‑1 (MQ=255)
                                 acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc                <  1:2737380/69‑1 (MQ=255)
                                 acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc                <  1:867023/69‑1 (MQ=255)
                                 acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc                <  1:248017/69‑1 (MQ=255)
                                 acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc                <  1:3058572/69‑1 (MQ=255)
                                 acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc                <  1:728283/69‑1 (MQ=255)
                                 acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc                <  1:1422785/69‑1 (MQ=255)
                                 acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc                <  1:1862837/69‑1 (MQ=255)
                                 acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc                <  1:1727676/69‑1 (MQ=255)
                                 acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc                <  1:818916/69‑1 (MQ=255)
                                 acacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTCGATCTGCCGCGCGATGc                <  1:180504/69‑1 (MQ=255)
                                  cacGCCGTCGGTACCTTCAGCGTTGTTGGTCTCTTTCAGGAAACGATTTTGGATCTGCCGCGCGATGc                <  1:604527/68‑1 (MQ=255)
                                  cacGCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGcc               >  1:3428281/1‑69 (MQ=255)
                                     gCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTg                           >  1:955038/1‑54 (MQ=255)
                                     gCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCg              >  1:998354/1‑67 (MQ=255)
                                     gCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGtt            >  1:1663247/1‑69 (MQ=255)
                                     gCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGtt            >  1:1506129/1‑69 (MQ=255)
                                     gCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGt             >  1:846910/1‑68 (MQ=255)
                                     gCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGt             >  1:1911962/1‑68 (MQ=255)
                                       cGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGa          <  1:835379/69‑1 (MQ=255)
                                       cGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGa          <  1:2281934/69‑1 (MQ=255)
                                       cGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGa          <  1:1505307/69‑1 (MQ=255)
                                       cGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGa          <  1:3360240/69‑1 (MQ=255)
                                           aGTACGTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCcgcgc                     <  1:278673/54‑1 (MQ=255)
                                           aGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCcgcgc                     <  1:407657/54‑1 (MQ=255)
                                           aGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCcgcgc                     <  1:3164099/54‑1 (MQ=255)
                                           aGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCcgcgc                     <  1:274742/54‑1 (MQ=255)
                                           aGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCcgcgc                     <  1:2696541/54‑1 (MQ=255)
                                           aGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCcgcgc                     <  1:2612671/54‑1 (MQ=255)
                                           aGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCcgcgc                     <  1:1114633/54‑1 (MQ=255)
                                           aGTACCTTCGGCGTT‑TTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCcgcgc                     <  1:759826/53‑1 (MQ=255)
                                               ccTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGACACAGccc  >  1:3149105/1‑69 (MQ=255)
                                               ccTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGACACAGccc  >  1:939995/1‑69 (MQ=255)
                                               ccTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGACACAGccc  >  1:1162840/1‑69 (MQ=255)
                                                    |                                                               
GCTGTGAGCAGCCGTAGGTGTGGCTGAAGAGGAACACGCCGTCAGTACCTTCGGCGTTGTTGGTCTCTTTCAGAAAACGGTTCTGGATCTGCCGCGCGATGCCGTTGACACAGCCC  >  minE/2078012‑2078127

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