Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,248,433 A→C 67.3% P280P (CCT→CCG 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,4330AC67.3% 56.4 / 31.8 52P280P (CCT→CCGruvBATP‑dependent DNA helicase, component of RuvABC resolvasome
Reads supporting (aligned to +/- strand):  ref base A (6/11);  new base C (30/5);  total (36/16)
Fisher's exact test for biased strand distribution p-value = 4.20e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

GTTTCACGTTCTTCGCCAATGGCTGCCGCCAGGTTATCCAGACCTACAGGTCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGAAA  >  minE/1248386‑1248493
                                               |                                                            
gTTTCACGTTCTTCGCCAATGGCTGCCGCCAGGTTATCCAGACCTACAGGTCCACCAAAGAACTTATCg                                         >  1:103916/1‑69 (MQ=255)
          cttcGCCAATGGCTGCCGCCAGGTTATCCAGACCTACAGGTCCACCAAAGAACTTATCGATTACCGcc                                >  1:118898/1‑68 (MQ=255)
          cttcGCCAATGGCTGCCGCCAGGTTATCCAGACCTACAGGTCCACCAAAGAACTTATCGATTACCGcc                                >  1:552473/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:2804603/1‑68 (MQ=255)
                         ccgccAGGTTATCCAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCgg                 >  1:3257700/1‑68 (MQ=255)
                         ccgccAGGTTATCCAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCgg                 >  1:3308060/1‑68 (MQ=255)
                         ccgccAGGTTATCCAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCgg                 >  1:2059076/1‑68 (MQ=255)
                         ccgccAGGTTATCCAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCgg                 >  1:362832/1‑68 (MQ=255)
                         ccgccAGGTTATCCAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCgg                 >  1:1930380/1‑68 (MQ=255)
                         ccgccAGGTTATCCAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCgg                 >  1:1633426/1‑68 (MQ=255)
                         ccgccAGGTTATCCAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCgg                 >  1:425629/1‑68 (MQ=255)
                         ccgccAGGTTATCCAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCgg                 >  1:485973/1‑68 (MQ=255)
                         ccgccAGGTTATCCAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCgg                 >  1:516421/1‑68 (MQ=255)
                         ccgccAGGTTATCCAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCgg                 >  1:662687/1‑68 (MQ=255)
                         ccgccAGGTTATCCAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCgg                 >  1:784987/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:388313/1‑69 (MQ=255)
                                     ccAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa    >  1:3391706/1‑69 (MQ=255)
                                     ccAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa    >  1:3390850/1‑69 (MQ=255)
                                     ccAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa    >  1:3360846/1‑69 (MQ=255)
                                     ccAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa    >  1:2983956/1‑69 (MQ=255)
                                     ccAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa    >  1:1043090/1‑69 (MQ=255)
                                     ccAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa    >  1:2767501/1‑69 (MQ=255)
                                     ccAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa    >  1:1083724/1‑69 (MQ=255)
                                     ccAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa    >  1:1302741/1‑69 (MQ=255)
                                     ccAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa    >  1:1417199/1‑69 (MQ=255)
                                     ccAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa    >  1:2701282/1‑69 (MQ=255)
                                     ccAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa    >  1:1535007/1‑69 (MQ=255)
                                     ccAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa    >  1:1723396/1‑69 (MQ=255)
                                     ccAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa    >  1:2064971/1‑69 (MQ=255)
                                     ccAGACCTACCGGCCCACCAAAGAACTTATAGATAACCGCCAGcaac                          >  1:1412103/1‑47 (MQ=37)
                                      cAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGaa   <  1:1957110/69‑1 (MQ=255)
                                      cAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGaa   <  1:2158973/69‑1 (MQ=255)
                                      cAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGaa   <  1:2242019/69‑1 (MQ=255)
                                      cAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGaa   <  1:2614013/69‑1 (MQ=255)
                                      cAGACCTACCGGCCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGaa   <  1:2965184/69‑1 (MQ=255)
                                      cAGACCTACAGGTCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa    <  1:2823590/68‑1 (MQ=255)
                                      cAGACCTACAGGTCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa    <  1:1633789/68‑1 (MQ=255)
                                      cAGACCTACAGGTCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa    <  1:3037999/68‑1 (MQ=255)
                                      cAGACCTACAGGTCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGa    <  1:1415028/68‑1 (MQ=255)
                                        gACCTACAGGTCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGaaa  <  1:503262/68‑1 (MQ=255)
                                        gACCTACAGGTCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGaaa  <  1:2729972/68‑1 (MQ=255)
                                        gACCTACAGGTCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGaaa  <  1:1186453/68‑1 (MQ=255)
                                        gACCTACAGGTCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGaaa  <  1:2912393/68‑1 (MQ=255)
                                        gACCTACAGGTCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGaaa  <  1:973018/68‑1 (MQ=255)
                                               |                                                            
GTTTCACGTTCTTCGCCAATGGCTGCCGCCAGGTTATCCAGACCTACAGGTCCACCAAAGAACTTATCGATTACCGCCAGCAACAATTTGCGGTCCATATAATCGAAA  >  minE/1248386‑1248493

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: