Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,156,434 C→T 39.6% P138P (CCG→CCA 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,4340CT39.6% 29.3 / 40.4 63P138P (CCG→CCAobgEGTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base C (14/24);  new base T (10/15);  total (24/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTT  >  minE/2156379‑2156493
                                                       |                                                           
tCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTaa                                                >  1:2706274/1‑69 (MQ=255)
tCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTaa                                                >  1:727381/1‑69 (MQ=255)
tCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTaa                                                >  1:1871465/1‑69 (MQ=255)
tCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTaa                                                >  1:1870657/1‑69 (MQ=255)
tCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGGTGGTTTTCTGCCGCGGTGTACGGTTaa                                                >  1:2492928/1‑69 (MQ=255)
  aGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGtt                                                  <  1:3076579/65‑1 (MQ=255)
  aGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGtt                                                  <  1:2314784/65‑1 (MQ=255)
  aGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGtt                                                  <  1:2441776/65‑1 (MQ=255)
  aGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGtt                                                  <  1:1914522/65‑1 (MQ=255)
  aGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGtt                                                  <  1:283867/65‑1 (MQ=255)
  aGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGtt                                                  <  1:2868152/65‑1 (MQ=255)
  aGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGtt                                                  <  1:180459/65‑1 (MQ=255)
  aGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGtt                                                  <  1:3373689/65‑1 (MQ=255)
  aGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGtt                                                  <  1:3396198/65‑1 (MQ=255)
  aGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGtt                                                  <  1:1518374/65‑1 (MQ=255)
  aGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGtt                                                  <  1:987447/65‑1 (MQ=255)
  aGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGtt                                                  <  1:101426/65‑1 (MQ=255)
    cTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACgg                                             <  1:2139456/68‑1 (MQ=255)
    cTCCAGCAGCAGCTCACGCTTATCGCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGtt                                                  >  1:1496113/1‑63 (MQ=255)
    cTCCAGCAGCAGCTCACGCTTATCGCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAActg                                             >  1:1049433/1‑66 (MQ=255)
    cTCCAGCAGCAGCTCACGCTTATCGCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAActg                                             >  1:3408871/1‑66 (MQ=255)
    cTCCAGCAGCAGCTCACGCTTATCGCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAActg                                             >  1:525244/1‑66 (MQ=255)
    cTCCAGCAGCAGCTCACGCTTATCGCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAActg                                             >  1:735507/1‑66 (MQ=255)
    cTCCAGCAGCAGCTCACGCTTATCGCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAActg                                             >  1:760635/1‑66 (MQ=255)
    cTCCAGCAGCAGCTCACGCTTATCGCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAActg                                             >  1:989330/1‑66 (MQ=255)
    cTCCAGCAGCAGCTCACGCTTATCGCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAActg                                             >  1:2424796/1‑66 (MQ=255)
    cTCCAGCAGCAGCTCACGCTTATCGCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAActg                                             >  1:1867218/1‑66 (MQ=255)
    cTCCAGCAGCAGCTCACGCTTATCGCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAAct                                              >  1:2654424/1‑66 (MQ=255)
             cagcTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGaa                                   >  1:553358/1‑69 (MQ=255)
             cagcTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGaa                                   >  1:3296405/1‑69 (MQ=255)
             cagcTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGCTGTACGGTTAACGGACGATTTGaa                                   >  1:3188072/1‑69 (MQ=255)
              agcTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGaaa                                  >  1:2267726/1‑69 (MQ=255)
                  cgcgCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACggg                              <  1:898731/69‑1 (MQ=255)
                  cgcgCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACggg                              <  1:1316631/69‑1 (MQ=255)
                           cGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGcc                       <  1:1250715/67‑1 (MQ=255)
                           cGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGcc                       <  1:2073098/67‑1 (MQ=255)
                            gCCCTGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAACTGACGACTTGAAACGGGTATTGCCCa                     <  1:2327599/68‑1 (MQ=255)
                            gCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAACTGACGACTTGAAACGGGTATTGCCCa                     <  1:1075850/68‑1 (MQ=255)
                            gCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAACTGACGACTTGAAACGGGTATTGCCCa                     <  1:764190/68‑1 (MQ=255)
                            gCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAACTGACGACTTGAAACGGGTATTGCCCa                     <  1:650685/68‑1 (MQ=255)
                            gCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAACTGACGACTTGAAACGGGTATTGCCCa                     <  1:369676/68‑1 (MQ=255)
                            gCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAACTGACGACTTGAAACGGGTATTGCCCa                     <  1:3339155/68‑1 (MQ=255)
                            gCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAACTGACGACTTGAAACGGGTATTGCCCa                     <  1:1753333/68‑1 (MQ=255)
                            gCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAACTGACGACTTGAAACGGGTATTGCCCa                     <  1:3007254/68‑1 (MQ=255)
                            gCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAACTGACGACTTGAAACGGGTATTGCCCa                     <  1:2327774/68‑1 (MQ=255)
                            gCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAACTGACGACTTGAAACGGGTATTGCCCa                     <  1:2424056/68‑1 (MQ=255)
                            gCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAACTGACGACTTGAAACGGGTATTGCCCa                     <  1:2194750/68‑1 (MQ=255)
                            gCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAACTGACGACTTGAAACGGGTATTGCCCa                     <  1:2470731/68‑1 (MQ=255)
                            gCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAACTGACGACTTGAAACGGGTATTGCCCa                     <  1:2029854/68‑1 (MQ=255)
                            gCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAACTGACGACTTGAAACGGGTATTGCCCa                     <  1:193438/68‑1 (MQ=255)
                             cccGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAACTGACGACTTGAAACGGGTATTGCCCa                     <  1:3427035/67‑1 (MQ=255)
                                 gCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGt               <  1:3064786/69‑1 (MQ=255)
                                 gCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGt               <  1:3355/69‑1 (MQ=255)
                                           ggTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCg                >  1:1513432/1‑58 (MQ=255)
                                           ggTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAgccgcc     >  1:2429258/1‑69 (MQ=255)
                                           ggTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAgccgcc     >  1:1537801/1‑69 (MQ=255)
                                           ggTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAgccgcc     >  1:1547316/1‑69 (MQ=255)
                                           ggTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAgccgc      >  1:1355741/1‑68 (MQ=255)
                                              tttCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCtt  <  1:404278/69‑1 (MQ=255)
                                              tttCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCtt  <  1:1910712/69‑1 (MQ=255)
                                              tttCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCtt  <  1:1352837/69‑1 (MQ=255)
                                              tttCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCtt  <  1:838898/69‑1 (MQ=255)
                                              tttCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCtt  <  1:1093401/69‑1 (MQ=255)
                                                       |                                                           
TCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTT  >  minE/2156379‑2156493

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: