Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,248,436 T→C 66.1% G279G (GGA→GGG)  ruvB ← ATP‑dependent DNA helicase, component of RuvABC resolvasome

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,248,4360TC66.1% 39.6 / 26.1 53G279G (GGA→GGG) ruvBATP‑dependent DNA helicase, component of RuvABC resolvasome
Reads supporting (aligned to +/- strand):  ref base T (6/12);  new base C (30/5);  total (36/17)
Fisher's exact test for biased strand distribution p-value = 3.41e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

GTTTCACGTTCTTCGCCAATGGCTGCCGCCAGGTTATCCAGACCTACAGGTCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGAAACCTTCAGCATC  >  minE/1248386‑1248504
                                                  |                                                                    
gTTTCACGTTCTTCGCCAATGGCTGCCGCCAGGTTATCCAGACCTACAGGTCCACCAAAGAACTTATCg                                                    >  1:103916/1‑69 (MQ=255)
          cttcGCCAATGGCTGCCGCCAGGTTATCCAGACCTACAGGTCCACCAAAGAACTTATCGATTACCGcc                                           >  1:552473/1‑68 (MQ=255)
          cttcGCCAATGGCTGCCGCCAGGTTATCCAGACCTACAGGTCCACCAAAGAACTTATCGATTACCGcc                                           >  1:118898/1‑68 (MQ=255)
                  aTGGCTGCCGCCAGGTTATCCAGACCTACAGGTCCACCAAAGAACTTATCGAt                                                  >  1:2784318/1‑53 (MQ=255)
                  aTGGCTGCCGCCAGGTTATCCAGACCTACAGGTCCACCAAAGAACTTATCGAt                                                  >  1:1626930/1‑53 (MQ=255)
                  aTGGCTGCCGCCAGCTTATCCAGACCTACAGGTCCACCAAAGAACTTATCGAt                                                  >  1:81738/1‑53 (MQ=255)
                     gCTGCCGCCAGGTTATCCAGACCTACAGGTCCACCAAAGAACTTATc                                                     <  1:2872487/47‑1 (MQ=255)
                      cTGCCGCCAGGTTATCCAGACCTACAGGTCCACCAAAGAACTTATc                                                     <  1:1717581/46‑1 (MQ=255)
                         ccgccAGGTTATCCAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCgg                            >  1:425629/1‑68 (MQ=255)
                         ccgccAGGTTATCCAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCgg                            >  1:516421/1‑68 (MQ=255)
                         ccgccAGGTTATCCAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCgg                            >  1:662687/1‑68 (MQ=255)
                         ccgccAGGTTATCCAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCgg                            >  1:362832/1‑68 (MQ=255)
                         ccgccAGGTTATCCAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCgg                            >  1:784987/1‑68 (MQ=255)
                         ccgccAGGTTATCCAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCgg                            >  1:2804603/1‑68 (MQ=255)
                         ccgccAGGTTATCCAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCgg                            >  1:2059076/1‑68 (MQ=255)
                         ccgccAGGTTATCCAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCgg                            >  1:1930380/1‑68 (MQ=255)
                         ccgccAGGTTATCCAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCgg                            >  1:485973/1‑68 (MQ=255)
                         ccgccAGGTTATCCAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCgg                            >  1:3257700/1‑68 (MQ=255)
                         ccgccAGGTTATCCAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCgg                            >  1:1633426/1‑68 (MQ=255)
                         ccgccAGGTTATCCAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCgg                            >  1:3308060/1‑68 (MQ=255)
                         ccgccAGGTTATCCAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCg                             >  1:2265954/1‑67 (MQ=255)
                         ccgccAGGTTATCCAGACCTACCGGCCCACCAAAGAACTAATCGATTACCGCCAGCAACAATTTGCgg                            >  1:2556752/1‑68 (MQ=255)
                                     ccAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAAttt                                >  1:36410/1‑52 (MQ=255)
                                     ccAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa               >  1:2983956/1‑69 (MQ=255)
                                     ccAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa               >  1:388313/1‑69 (MQ=255)
                                     ccAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa               >  1:3391706/1‑69 (MQ=255)
                                     ccAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa               >  1:3390850/1‑69 (MQ=255)
                                     ccAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa               >  1:2767501/1‑69 (MQ=255)
                                     ccAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa               >  1:3360846/1‑69 (MQ=255)
                                     ccAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa               >  1:1043090/1‑69 (MQ=255)
                                     ccAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa               >  1:2064971/1‑69 (MQ=255)
                                     ccAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa               >  1:1083724/1‑69 (MQ=255)
                                     ccAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa               >  1:1723396/1‑69 (MQ=255)
                                     ccAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa               >  1:2701282/1‑69 (MQ=255)
                                     ccAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa               >  1:1302741/1‑69 (MQ=255)
                                     ccAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa               >  1:1535007/1‑69 (MQ=255)
                                     ccAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa               >  1:1417199/1‑69 (MQ=255)
                                     ccAGACCTACCGGCCCACCAAAGAACTTATAGATAACCGCCAGcaac                                     >  1:1412103/1‑47 (MQ=37)
                                      cAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGaa              <  1:2614013/69‑1 (MQ=255)
                                      cAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGaa              <  1:2242019/69‑1 (MQ=255)
                                      cAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGaa              <  1:2158973/69‑1 (MQ=255)
                                      cAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGaa              <  1:2965184/69‑1 (MQ=255)
                                      cAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGaa              <  1:1957110/69‑1 (MQ=255)
                                      cAGACCTACAGGTCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa               <  1:3037999/68‑1 (MQ=255)
                                      cAGACCTACAGGTCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa               <  1:1415028/68‑1 (MQ=255)
                                      cAGACCTACAGGTCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa               <  1:1633789/68‑1 (MQ=255)
                                      cAGACCTACAGGTCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa               <  1:2823590/68‑1 (MQ=255)
                                        gACCTACAGGTCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGaaa             <  1:503262/68‑1 (MQ=255)
                                        gACCTACAGGTCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGaaa             <  1:973018/68‑1 (MQ=255)
                                        gACCTACAGGTCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGaaa             <  1:2912393/68‑1 (MQ=255)
                                        gACCTACAGGTCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGaaa             <  1:2729972/68‑1 (MQ=255)
                                        gACCTACAGGTCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGaaa             <  1:1186453/68‑1 (MQ=255)
                                                ggTCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGAAACCTTCAGc     <  1:1111674/68‑1 (MQ=255)
                                                  tCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGAAACCTTCAGCATc  >  1:3315122/1‑69 (MQ=255)
                                                  tCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGAAACCTTCAGCATc  >  1:3015802/1‑69 (MQ=255)
                                                  tCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGAAACCTTCAGCATc  >  1:2763677/1‑69 (MQ=255)
                                                  tCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGAAACCTTCAGCATc  >  1:617666/1‑69 (MQ=255)
                                                  tCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGAAACCTTCAGCATc  >  1:2390383/1‑69 (MQ=255)
                                                  tCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGAAACCTTCAGCATc  >  1:221484/1‑69 (MQ=255)
                                                  |                                                                    
GTTTCACGTTCTTCGCCAATGGCTGCCGCCAGGTTATCCAGACCTACAGGTCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGAAACCTTCAGCATC  >  minE/1248386‑1248504

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: ‑