Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 227,492 C→T 24.8% L96L (CTG→TTG)  crl → DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE227,4920CT24.8% 84.0 / 21.6 68L96L (CTG→TTG) crlDNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (45/6);  new base T (11/6);  total (56/12)
Fisher's exact test for biased strand distribution p-value = 5.96e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

CTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAG  >  minE/227429‑227557
                                                               |                                                                 
cTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAg                                                              <  1:1503086/69‑1 (MQ=255)
cTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAg                                                              <  1:776016/69‑1 (MQ=255)
      gATAAAGCAGGCGACTGGACGAGTGTTCGGATAAAAGATGCCGAAGTGGTTGAACGATTGGAGCACAcc                                                        >  1:436736/1‑69 (MQ=255)
      gATAAAGCAGGCGACTGGACGAGTGTTCGGATAAAAGATGCCGAAGTGGTTGAACGATTGGAGCACAcc                                                        >  1:2675771/1‑69 (MQ=255)
      gATAAAGCAGGCGACTGGACGAGTGTTCGGATAAAAGATGCCGAAGTGGTTGAACGATTGGAGCACAcc                                                        >  1:2408050/1‑69 (MQ=255)
      gATAAAGCAGGCGACTGGACGAGTGTTCGGATAAAAGATGCCGAAGTGGTTGAACGATTGGAGCACAcc                                                        >  1:2339950/1‑69 (MQ=255)
      gATAAAGCAGGCGACTGGACGAGTGTTCGGATAAAAGATGCCGAAGTGGTTGAACGATTGGAGCACAcc                                                        >  1:3212993/1‑69 (MQ=255)
      gATAAAGCAGGCGACTGGACGAGTGTTCGGATAAAAGATGCCGAAGTGGTTGAACGATTGGAGCACAcc                                                        >  1:3296124/1‑69 (MQ=255)
      gATAAAGCAGGCGACTGGACGAGTGTTCGGATAAAAGATGCCGAAGTGGTTGAACGATTGGAGCACAcc                                                        >  1:3300774/1‑69 (MQ=255)
      gATAAAGCAGGCGACTGGACGAGTGTTCGGATAAAAGATGCCGAAGTGGTTGAACGATTGGAGCACAcc                                                        >  1:1369746/1‑69 (MQ=255)
      gATAAAGCAGGCGACTGGACGAGTGTTCGGATAAAAGATGCCGAAGTGGTTGAACGATTGGAGCACAcc                                                        >  1:1365869/1‑69 (MQ=255)
      gATAAAGCAGGCGACTGGACGAGTGTTCGGATAAAAGATGCCGAAGTGGTTGAACGATTGGAGCACAcc                                                        >  1:1279759/1‑69 (MQ=255)
      gATAAAGCAGGCGACTGGACGAGTGTTCGGATAAAAGATGCCGAAGTGGTTGAACGATTGGAGCACAcc                                                        >  1:914328/1‑69 (MQ=255)
                      ggAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTc                                           >  1:788170/1‑66 (MQ=255)
                      ggAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTc                                           >  1:1693349/1‑66 (MQ=255)
                      ggAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTc                                           >  1:2944419/1‑66 (MQ=255)
                      ggAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTc                                           >  1:995538/1‑66 (MQ=255)
                            gtgtTCCGGTAAAAGACACTGAAGTGGTTGAACTACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGc                                  >  1:1695712/1‑69 (MQ=255)
                            gtgtTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGc                                  >  1:567151/1‑69 (MQ=255)
                            gtgtTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGc                                  >  1:475677/1‑69 (MQ=255)
                            gtgtTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGc                                  >  1:3023343/1‑69 (MQ=255)
                            gtgtTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGc                                  >  1:3405390/1‑69 (MQ=255)
                            gtgtTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGc                                  >  1:81580/1‑69 (MQ=255)
                            gtgtTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGc                                  >  1:3170565/1‑69 (MQ=255)
                            gtgtTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGc                                  >  1:1565863/1‑69 (MQ=255)
                            gtgtTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGc                                  >  1:2026157/1‑69 (MQ=255)
                                   ggTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTg                            <  1:2716947/68‑1 (MQ=255)
                                   ggTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTg                            <  1:1887193/68‑1 (MQ=255)
                                   ggTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTg                            <  1:1046640/68‑1 (MQ=255)
                                    gTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTg                            <  1:2693111/67‑1 (MQ=255)
                                         aGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActgct                     >  1:1052905/1‑69 (MQ=255)
                                         aGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActgct                     >  1:864117/1‑69 (MQ=255)
                                         aGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActgct                     >  1:1688300/1‑69 (MQ=255)
                                         aGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActgct                     >  1:3434018/1‑69 (MQ=255)
                                         aGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActgct                     >  1:31800/1‑69 (MQ=255)
                                         aGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActgct                     >  1:3172577/1‑69 (MQ=255)
                                         aGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActgct                     >  1:2736199/1‑69 (MQ=255)
                                         aGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActgc                      >  1:2987621/1‑68 (MQ=255)
                                          gACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActgctg                    >  1:303026/1‑69 (MQ=255)
                                          gACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActgctg                    >  1:3028918/1‑69 (MQ=255)
                                          gACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActgctg                    >  1:920498/1‑69 (MQ=255)
                                          gACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActgctg                    >  1:2099838/1‑69 (MQ=255)
                                          gACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActgctg                    >  1:1665575/1‑69 (MQ=255)
                                          gACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActgctg                    >  1:481535/1‑69 (MQ=255)
                                          gACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActgctg                    >  1:436515/1‑69 (MQ=255)
                                            cacTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTgac                  >  1:3035679/1‑69 (MQ=255)
                                            cacTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTgac                  >  1:2705043/1‑69 (MQ=255)
                                            cacTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTgac                  >  1:116082/1‑69 (MQ=255)
                                            cacTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTgac                  >  1:1334050/1‑69 (MQ=255)
                                            cacTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTgac                  >  1:1619916/1‑69 (MQ=255)
                                            cacTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTgac                  >  1:455622/1‑69 (MQ=255)
                                            cacTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTgac                  >  1:3420277/1‑69 (MQ=255)
                                            cacTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTgac                  >  1:2059207/1‑69 (MQ=255)
                                            cacTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTgac                  >  1:2958733/1‑69 (MQ=255)
                                                gAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGacgacg              >  1:2585842/1‑69 (MQ=255)
                                                gAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGacgacg              >  1:665101/1‑69 (MQ=255)
                                                     ggTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTg                                >  1:706738/1‑46 (MQ=255)
                                                     ggTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGaa         >  1:2375904/1‑69 (MQ=255)
                                                     ggTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGaa         >  1:2990706/1‑69 (MQ=255)
                                                     ggTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGa          >  1:3228813/1‑68 (MQ=255)
                                                     ggTTGAACGACTGGAGCACACCATGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGaa         >  1:584141/1‑69 (MQ=255)
                                                      gTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAAt        >  1:1690027/1‑69 (MQ=255)
                                                            cGATTGGAGCACACCCTGCGTGAATTTCACGAGAAGCTTCGTGAGCTGCTGGCGACGCTGAATCTGAAg  <  1:83764/69‑1 (MQ=255)
                                                            cGATTGGAGCACACCCTGCGTGAATTTCACGAGAAGCTGCGTGAGCTGCTGGCGACGCTGAATCTGAAg  <  1:596310/69‑1 (MQ=255)
                                                            cGATTGGAGCACACCCTGCGTGAATTTCACGAGAAGCTGCGTGAGCTGCTGGCGACGCTGAATCTGAAg  <  1:329984/69‑1 (MQ=255)
                                                            cGATTGGAGCACACCCTGCGTGAATTTCACGAGAAGCTGCGTGAGCTGCTGGCGACGCTGAATCTGAAg  <  1:2715066/69‑1 (MQ=255)
                                                            cGATTGGAGCACACCCTGCGTGAATTTCACGAGAAGCTGCGTGAGCTGCTGGCGACGCTGAATCTGAAg  <  1:2940824/69‑1 (MQ=255)
                                                            cGATTGGAGCACACCCTGCGTGAATTTCACGAGAAGCTGCGTGAGCTGCTGGCGACGCTGAATCTGAAg  <  1:924377/69‑1 (MQ=255)
                                                               |                                                                 
CTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAG  >  minE/227429‑227557

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: