Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,229,092 C→A 29.4% A23A (GCC→GCA holE → DNA polymerase III, theta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,229,0920CA29.4% 108.1 / 42.0 78A23A (GCC→GCAholEDNA polymerase III, theta subunit
Reads supporting (aligned to +/- strand):  ref base C (28/27);  new base A (17/6);  total (45/33)
Fisher's exact test for biased strand distribution p-value = 8.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACA  >  minE/1229025‑1229154
                                                                   |                                                              
tGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCg                                                               >  1:2069001/1‑69 (MQ=255)
  cTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCgg                                                              <  1:3218497/68‑1 (MQ=255)
   tGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCAggtg                                                            >  1:407488/1‑67 (MQ=255)
   tGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCAgggg                                                            >  1:2317843/1‑69 (MQ=255)
   tGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCAgggg                                                            >  1:2137566/1‑69 (MQ=255)
   tGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCAgggg                                                            >  1:218123/1‑69 (MQ=255)
   tGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCAgggg                                                            >  1:2106473/1‑69 (MQ=255)
   tGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCAgggg                                                            >  1:3227181/1‑69 (MQ=255)
   tGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCAgggg                                                            >  1:191628/1‑69 (MQ=255)
   tGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCAgggg                                                            >  1:2207742/1‑69 (MQ=255)
   tGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCAgggg                                                            >  1:488795/1‑69 (MQ=255)
   tGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCAgggg                                                            >  1:1702278/1‑69 (MQ=255)
   tGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCAgggg                                                            >  1:3117324/1‑69 (MQ=255)
   tGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCAgggg                                                            >  1:2677547/1‑69 (MQ=255)
   tGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCAgggg                                                            >  1:1310251/1‑69 (MQ=255)
   tGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCAgggg                                                            >  1:2895097/1‑69 (MQ=255)
   tGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCAgggg                                                            >  1:2605125/1‑69 (MQ=255)
   tGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCAgggg                                                            >  1:1486743/1‑69 (MQ=255)
   tGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCAggg                                                             >  1:557311/1‑68 (MQ=255)
        aaTCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGc                                                        <  1:2192628/68‑1 (MQ=255)
         aTCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCgg                                                              >  1:2761736/1‑61 (MQ=255)
         aTCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCgg                                                              >  1:744858/1‑61 (MQ=255)
         aTCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCgg                                                              >  1:2505750/1‑61 (MQ=255)
         aTCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCgg                                                              >  1:1836795/1‑61 (MQ=255)
         aTCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCgg                                                              >  1:1725357/1‑61 (MQ=255)
          tCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCa                                                       >  1:156044/1‑67 (MQ=255)
          tCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCa                                                       >  1:2966883/1‑67 (MQ=255)
          tCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCAt                                                      >  1:2508612/1‑68 (MQ=255)
          tCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCAt                                                      >  1:3142940/1‑68 (MQ=255)
          tCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCAt                                                      >  1:2828545/1‑68 (MQ=255)
          tCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCAt                                                      >  1:978343/1‑68 (MQ=255)
          tCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCAt                                                      >  1:637297/1‑68 (MQ=255)
          tCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCAt                                                      >  1:54946/1‑68 (MQ=255)
           cTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCAtt                                                     >  1:953203/1‑68 (MQ=255)
           cTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCAtt                                                     >  1:1375960/1‑68 (MQ=255)
           cTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCAtt                                                     >  1:3275063/1‑68 (MQ=255)
           cTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCAtt                                                     >  1:2263290/1‑68 (MQ=255)
            tGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCAtt                                                     >  1:1265878/1‑67 (MQ=255)
                      ggATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTa                                         <  1:2147979/69‑1 (MQ=255)
                      ggATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTa                                         <  1:1495508/69‑1 (MQ=255)
                      ggATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTa                                         <  1:1574416/69‑1 (MQ=255)
                      ggATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTa                                         <  1:1714358/69‑1 (MQ=255)
                      ggATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTa                                         <  1:2102030/69‑1 (MQ=255)
                      ggATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTa                                         <  1:2108889/69‑1 (MQ=255)
                      ggATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTa                                         <  1:2166230/69‑1 (MQ=255)
                      ggATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTa                                         <  1:3219747/69‑1 (MQ=255)
                         tCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACaa                                      <  1:3044639/69‑1 (MQ=255)
                         tCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACaa                                      <  1:324405/69‑1 (MQ=255)
                         tCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACaa                                      <  1:2357917/69‑1 (MQ=255)
                         tCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACaa                                      <  1:2609444/69‑1 (MQ=255)
                         tCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACaa                                      <  1:1145160/69‑1 (MQ=255)
                                  aaTGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGt                             <  1:1522111/69‑1 (MQ=255)
                                  aaTGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGt                             <  1:429893/69‑1 (MQ=255)
                                  aaTGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGt                             <  1:705541/69‑1 (MQ=255)
                                  aaTGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGt                             <  1:3231023/69‑1 (MQ=255)
                                  aaTGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGt                             <  1:3116103/69‑1 (MQ=255)
                                  aaTGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGt                             <  1:2454587/69‑1 (MQ=255)
                                  aaTGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGCTGGCATTTAAAGAACGCTACAATATGCCGGt                             <  1:2182630/69‑1 (MQ=255)
                                        tAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCg                        >  1:2086463/1‑68 (MQ=255)
                                        tAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCg                        >  1:1476654/1‑68 (MQ=255)
                                                   tCGATTTGGCGGTGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGAt                          >  1:684988/1‑55 (MQ=255)
                                                   tCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGAt                          >  1:2134639/1‑55 (MQ=255)
                                                   tCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGAt                          >  1:2129850/1‑55 (MQ=255)
                                                   tCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGAt                          >  1:688093/1‑55 (MQ=255)
                                                   tCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGAt                          >  1:805406/1‑55 (MQ=255)
                                                   tCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGAt                          >  1:456290/1‑55 (MQ=255)
                                                   tCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTg            >  1:2062939/1‑69 (MQ=255)
                                                       ttttGCGGCGGCAGGGGTGGCATTTAAAGAACGCTACAATATGCCAGTGATCGCTGAAGCGGTTGAACg        <  1:1878639/69‑1 (MQ=255)
                                                       tttGGCGGCGGCAGGGGTGGCATTTAAAGAACGCTACAATATGCCAGTGATCGCTGAAGCGGTTGAACg        <  1:1513308/69‑1 (MQ=255)
                                                       tttGGCGGCGGCAGGGGTGGCATTTAAAGAACGCTACAATATGCCAGTGATCGCTGAAGCGGTTGAACg        <  1:1226604/69‑1 (MQ=255)
                                                       tttGGCGGCGGCAGGGGTGGCATTTAAAGAACGCTACAATATGCCAGTGATCGCTGAAGCGGTTGAACg        <  1:1695332/69‑1 (MQ=255)
                                                       tttGGCGGCGGCAGGGGTGGCATTTAAAGAACGCTACAATATGCCAGTGATCGCTGAAGCGGTTGAACg        <  1:1887150/69‑1 (MQ=255)
                                                        ttGGCGGCGGCAGGGGTGGCATTTAAAGAACGCTACAATATGCCAGTGATCGCTGAAGCGGTTGAACg        <  1:131114/68‑1 (MQ=255)
                                                         tGGCGGCGGCCTGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTg      <  1:32672/69‑1 (MQ=255)
                                                             gtcggcCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACa  <  1:2987055/67‑1 (MQ=255)
                                                             ggcggcCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGt                             <  1:2820244/42‑1 (MQ=255)
                                                             ggcggcCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACa  <  1:1811515/69‑1 (MQ=255)
                                                             ggcggcCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACa  <  1:2483815/69‑1 (MQ=255)
                                                                   |                                                              
TGCTGAAGAATCTGGCTAAACTGGATCAAACAGAAATGGATAAAGTGAATGTCGATTTGGCGGCGGCCGGGGTGGCATTTAAAGAACGCTACAATATGCCGGTGATCGCTGAAGCGGTTGAACGTGAACA  >  minE/1229025‑1229154

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: