Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,377,077 C→T 32.2% K414K (AAG→AAA mglA ← fused methyl‑galactoside transporter subunitsl and ATP‑binding components of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,377,0770CT32.2% 57.7 / 32.7 68K414K (AAG→AAAmglAfused methyl‑galactoside transporter subunitsl and ATP‑binding components of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base C (11/35);  new base T (18/4);  total (29/39)
Fisher's exact test for biased strand distribution p-value = 1.22e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

CCGCGCGTCGGTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTTTTACCCGCATCG  >  minE/1377012‑1377144
                                                                 |                                                                   
ccGCGCGTCGGTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCttt                                                                  <  1:979615/69‑1 (MQ=255)
ccGCGCGTCGGTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCttt                                                                  <  1:2103738/69‑1 (MQ=255)
ccGCGCGTCGGTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCttt                                                                  <  1:463787/69‑1 (MQ=255)
ccGCGCGTCGGTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCttt                                                                  <  1:191805/69‑1 (MQ=255)
ccGCGCGTCGGTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCttt                                                                  <  1:2297753/69‑1 (MQ=255)
ccGCGCGTCGGTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCttt                                                                  <  1:232550/69‑1 (MQ=255)
ccGCGCGTCGGTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCttt                                                                  <  1:2531446/69‑1 (MQ=255)
    gcgTCGGTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCt                                                              <  1:1627209/69‑1 (MQ=255)
    gcgTCGGTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCt                                                              <  1:3265945/69‑1 (MQ=255)
    gcgTCGGTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCt                                                              <  1:3062594/69‑1 (MQ=255)
    gcgTCGGTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCt                                                              <  1:1909583/69‑1 (MQ=255)
    gcgTCGGTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCt                                                              <  1:3048495/69‑1 (MQ=255)
    gcgTCGGTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCt                                                              <  1:2830221/69‑1 (MQ=255)
    gcgTCGGTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCt                                                              <  1:1986890/69‑1 (MQ=255)
    gcgTCGGTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCt                                                              <  1:457614/69‑1 (MQ=255)
     cgTCGGTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCt                                                              <  1:2193690/68‑1 (MQ=255)
     cgTCGGTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTg                                                             >  1:569497/1‑69 (MQ=255)
      gTCGGTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCttt                                                                  <  1:2266020/63‑1 (MQ=255)
      gTCGGTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCt                                                              <  1:997064/67‑1 (MQ=255)
       tCGGTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGAt                                                           <  1:2740763/69‑1 (MQ=255)
       tCGGTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGAt                                                           <  1:2692987/69‑1 (MQ=255)
       tCGGTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGAt                                                           <  1:994029/69‑1 (MQ=255)
       tCGGTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGAt                                                           <  1:1716496/69‑1 (MQ=255)
       tCGGTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGAt                                                           <  1:2491952/69‑1 (MQ=255)
        cGGTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGAt                                                           <  1:2090635/68‑1 (MQ=255)
         ggTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCtttt                                                                 >  1:659290/1‑61 (MQ=255)
         ggTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGAtt                                                          >  1:456719/1‑68 (MQ=255)
         ggTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGAtt                                                          >  1:2202529/1‑68 (MQ=255)
         ggTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGAtt                                                          >  1:2223942/1‑68 (MQ=255)
         ggTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGAtt                                                          >  1:96956/1‑68 (MQ=255)
         ggTTCATCGAGCATTAATATTTATGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGAtt                                                          >  1:3377670/1‑68 (MQ=255)
            tCATCGAGCATTAATATTTCTGGTTGCGTTAACAGCCAGCGACCAATAATCACTTTTTGCTGATTacca                                                      <  1:1816455/69‑1 (MQ=39)
            tCATCGAGCATTAATATTTCTGGTTGCGTTAACAGCCAGCGACCAATAATCACTTTTTGCTGATTacca                                                      <  1:1073221/69‑1 (MQ=39)
            tCATCGAGCATTAATATTTCTGGTTGCGTTAACAGCCAGCGACCAATAATCACTTTTTGCTGATTacca                                                      <  1:751084/69‑1 (MQ=39)
             cATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCt                                                              <  1:780338/60‑1 (MQ=255)
                     aTTAATATTTCTGGTTGCGTTCACAGCCAGCGACCAATAATCACTTTTTGCTGATTacca                                                      <  1:7612/60‑1 (MQ=37)
                         atatTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCt                                                              <  1:540387/48‑1 (MQ=255)
                            tttCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCaa                                      <  1:3068344/69‑1 (MQ=255)
                            tttCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCaa                                      <  1:2181225/69‑1 (MQ=255)
                                 ggTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCaa                                      <  1:2446003/64‑1 (MQ=255)
                                              gCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTcc                             <  1:558535/60‑1 (MQ=255)
                                              gCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTcc                             <  1:939854/60‑1 (MQ=255)
                                              gCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTcc                             <  1:2746577/60‑1 (MQ=255)
                                              gCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTcc                             <  1:1360392/60‑1 (MQ=255)
                                              gCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTcc                             <  1:137077/60‑1 (MQ=255)
                                              gCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTcc                             <  1:1936562/60‑1 (MQ=255)
                                                cAGCGACCAATAATCACTTTTTGCTGATTACCACCGga                                                 >  1:1879010/1‑38 (MQ=255)
                                                cAGCGACCAATAATCACTTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGCCCCggt                  >  1:2827723/1‑68 (MQ=38)
                                                cAGCGACCAATAATCACTTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGCCCCggt                  >  1:898844/1‑68 (MQ=38)
                                                cAGCGACCAATAATCACTTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGCCCCggt                  >  1:882816/1‑68 (MQ=38)
                                                cAGCGACCAATAATCACTTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGCCCCggt                  >  1:822365/1‑68 (MQ=38)
                                                cAGCGACCAATAATCACTTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGCCCCggt                  >  1:1506733/1‑68 (MQ=38)
                                                cAGCGACCAATAATCACTTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGCCCCggt                  >  1:151383/1‑68 (MQ=38)
                                                cAGCGACCAATAATCACTTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGCCCCggt                  >  1:1653042/1‑68 (MQ=38)
                                                cAGCGACCAATAATCACTTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGCCCCggt                  >  1:1733686/1‑68 (MQ=38)
                                                cAGCGACCAATAATCACTTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGCCCCggt                  >  1:1759104/1‑68 (MQ=38)
                                                cAGCGACCAATAATCACTTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGCCCCggt                  >  1:207090/1‑68 (MQ=38)
                                                cAGCGACCAATAATCACTTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGCCCCggt                  >  1:2082690/1‑68 (MQ=38)
                                                cAGCGACCAATAATCACTTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGCCCCggt                  >  1:2330090/1‑68 (MQ=38)
                                                cAGCGACCAATAATCACTTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGCCCCggt                  >  1:2513362/1‑68 (MQ=38)
                                                cAGCGACCAATAATCACTTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGCCCCggt                  >  1:3253683/1‑68 (MQ=38)
                                                cAGCGACCAATAATCACTTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGCCCCggt                  >  1:3174598/1‑68 (MQ=38)
                                                cAGCGACCAATAATCACTTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGCCCCggt                  >  1:271721/1‑68 (MQ=38)
                                                cAGCGACCAATAATCACTTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGCCCCggt                  >  1:3046076/1‑68 (MQ=38)
                                                            atCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTTTTACCCGc      >  1:2640530/1‑69 (MQ=255)
                                                            atCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTTTTACCCGc      >  1:718668/1‑69 (MQ=255)
                                                            atCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTTTTACCCGc      >  1:2971597/1‑69 (MQ=255)
                                                            atCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTTTTACCCGc      >  1:1359244/1‑69 (MQ=255)
                                                                 cTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTTTTACCCGCATCg  <  1:2325159/68‑1 (MQ=255)
                                                                 |                                                                   
CCGCGCGTCGGTTCATCGAGCATTAATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTTTTACCCGCATCG  >  minE/1377012‑1377144

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: