Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,156,632 A→G 71.2% N72N (AAT→AAC obgE ← GTPase involved in cell partioning and DNA repair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,156,6320AG71.2% 78.3 / 32.0 76N72N (AAT→AACobgEGTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base A (4/18);  new base G (10/44);  total (14/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCACTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTG  >  minE/2156576‑2156686
                                                        |                                                      
ggCACTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACcgcg                                             <  1:462197/68‑1 (MQ=255)
               aCGTCTTTACCGCGCTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGAt                             >  1:2566396/1‑69 (MQ=255)
               aCGTCTTTACCGCGCTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGAt                             >  1:3376024/1‑69 (MQ=255)
               aCGTCTTTACCGCGCTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGAt                             >  1:2378122/1‑69 (MQ=255)
               aCGTCTTTACCGCGCTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGAt                             >  1:1876314/1‑69 (MQ=255)
               aCGTCTTTACCGCGCTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGAt                             >  1:3195900/1‑69 (MQ=255)
               aCGTCTTTACCGCGCTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGAt                             >  1:681309/1‑69 (MQ=255)
               aCGTCTTTACCGCGCTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGAt                             >  1:2933592/1‑69 (MQ=255)
               aCGTCTTTACCGCGCTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGAt                             >  1:834479/1‑69 (MQ=255)
               aCGTCTTTACCGCGCTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGAt                             >  1:1325927/1‑69 (MQ=255)
               aCGTCTTTACCGCGCTTACCGGTACAGTCTAGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGAt                             >  1:2150443/1‑69 (MQ=255)
                           cgcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCTCGCTCTGCACGGAAAGATTTTTCAAAACGa                 <  1:468834/69‑1 (MQ=255)
                           cgcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGa                 <  1:287199/69‑1 (MQ=255)
                           cgcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGa                 <  1:2841064/69‑1 (MQ=255)
                           cgcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGa                 <  1:3066834/69‑1 (MQ=255)
                           cgcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGa                 <  1:3285618/69‑1 (MQ=255)
                           cgcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGa                 <  1:2956248/69‑1 (MQ=255)
                           cgcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGa                 <  1:3304085/69‑1 (MQ=255)
                           cgcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGa                 <  1:2408350/69‑1 (MQ=255)
                           cgcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGa                 <  1:2386232/69‑1 (MQ=255)
                           cgcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGa                 <  1:299117/69‑1 (MQ=255)
                           cgcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGa                 <  1:3363124/69‑1 (MQ=255)
                           cgcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGa                 <  1:3123098/69‑1 (MQ=255)
                           cgcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGa                 <  1:3413873/69‑1 (MQ=255)
                           cgcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGa                 <  1:2171045/69‑1 (MQ=255)
                           cgcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGa                 <  1:872315/69‑1 (MQ=255)
                           cgcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGa                 <  1:2062348/69‑1 (MQ=255)
                           cgcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGa                 <  1:453127/69‑1 (MQ=255)
                           cgcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGa                 <  1:857380/69‑1 (MQ=255)
                           cgcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGa                 <  1:531082/69‑1 (MQ=255)
                           cgcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGa                 <  1:85667/69‑1 (MQ=255)
                           cgcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGa                 <  1:1716074/69‑1 (MQ=255)
                            gcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGAGCGCGCTCTGCACGGa                                  <  1:529021/51‑1 (MQ=255)
                            gcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGa                                  <  1:907398/51‑1 (MQ=255)
                            gcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGa                                  <  1:807703/51‑1 (MQ=255)
                            gcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGa                                  <  1:766222/51‑1 (MQ=255)
                            gcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGa                                  <  1:3136117/51‑1 (MQ=255)
                            gcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGa                                  <  1:635810/51‑1 (MQ=255)
                            gcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGa                                  <  1:3242669/51‑1 (MQ=255)
                            gcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGa                                  <  1:576003/51‑1 (MQ=255)
                            gcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGa                                  <  1:3346998/51‑1 (MQ=255)
                            gcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGa                                  <  1:418049/51‑1 (MQ=255)
                            gcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGa                                  <  1:3430808/51‑1 (MQ=255)
                            gcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGa                                  <  1:2995929/51‑1 (MQ=255)
                            gcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGa                                  <  1:1175521/51‑1 (MQ=255)
                            gcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGa                                  <  1:2004016/51‑1 (MQ=255)
                            gcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGa                                  <  1:2051852/51‑1 (MQ=255)
                            gcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGa                                  <  1:2118457/51‑1 (MQ=255)
                            gcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGa                                  <  1:2362669/51‑1 (MQ=255)
                            gcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGa                                  <  1:2468055/51‑1 (MQ=255)
                            gcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGa                                  <  1:2468964/51‑1 (MQ=255)
                            gcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGa                                  <  1:2679223/51‑1 (MQ=255)
                            gcTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTACACGGa                                  <  1:1640716/51‑1 (MQ=255)
                            gcTTACCGGTACAGTCTCGGCTTGCACAGTTCTGACCGCGCTCTGCACGGa                                  <  1:777451/51‑1 (MQ=255)
                             cTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGa                                  <  1:1572463/50‑1 (MQ=255)
                                   ggTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGa          <  1:3337216/68‑1 (MQ=255)
                                   ggTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGa          <  1:168127/68‑1 (MQ=255)
                                   ggTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGa          <  1:1788567/68‑1 (MQ=255)
                                   ggTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGa          <  1:1797808/68‑1 (MQ=255)
                                   ggTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGa          <  1:1801780/68‑1 (MQ=255)
                                   ggTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGa          <  1:1881706/68‑1 (MQ=255)
                                   ggTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGa          <  1:2766193/68‑1 (MQ=255)
                                   ggTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGa          <  1:2050122/68‑1 (MQ=255)
                                   ggTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGa          <  1:353625/68‑1 (MQ=255)
                                     tACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGa          <  1:2370680/66‑1 (MQ=255)
                                       cAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAAc                   <  1:1836803/55‑1 (MQ=255)
                                       cAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAAc                   <  1:669564/55‑1 (MQ=255)
                                       cAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAAc                   <  1:2140309/55‑1 (MQ=255)
                                       cAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAAc                   <  1:2160614/55‑1 (MQ=255)
                                       cAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAAc                   <  1:2433504/55‑1 (MQ=255)
                                       cAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAAc                   <  1:3120594/55‑1 (MQ=255)
                                        aGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAAc                   <  1:1368742/54‑1 (MQ=255)
                                            gcgGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCgtg  >  1:2713830/1‑67 (MQ=255)
                                            gcgGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCgtg  >  1:411519/1‑67 (MQ=255)
                                            gcgGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCgtg  >  1:1210869/1‑67 (MQ=255)
                                            gcgGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCgtg  >  1:117145/1‑67 (MQ=255)
                                                        |                                                      
GGCACTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTG  >  minE/2156576‑2156686

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: