Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 548,052 C→T 36.7% H391H (CAC→CAT ybiT → fused predicted transporter subunits and ATP‑binding components of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE548,0520CT36.7% 34.8 / 31.9 57H391H (CAC→CATybiTfused predicted transporter subunits and ATP‑binding components of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base C (28/8);  new base T (8/13);  total (36/21)
Fisher's exact test for biased strand distribution p-value = 4.32e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGGACAGCGGCACCGTAAAATGGTCTGAGAACGCGCGCATTGGTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATGGATGAGCCAGTGGAAGCAGGAAGGCG  >  minE/547989‑548116
                                                               |                                                                
aCCGGACAGCGGCACCGTAAAATGGTCTGAGAACGCGCGCATTGGTTACTACGCCCAGGACCATGAat                                                              >  1:996400/1‑68 (MQ=255)
aCCGGACAGCGGCACCGTAAAATGGTCTGAGAACGCGCGCATTGGTTACTACGCCCAGGACCATGAat                                                              >  1:2711408/1‑68 (MQ=255)
aCCGGACAGCGGCACCGTAAAATGGTCTGAGAACGCGCGCATTGGTTACTACGCCCAGGACCATGAat                                                              >  1:2386573/1‑68 (MQ=255)
aCCGGACAGCGGCACCGTAAAATGGTCTGAGAACGCGCGCATTGGTTACTACGCCCAGGACCATGAat                                                              >  1:2332840/1‑68 (MQ=255)
aCCGGACAGCGGCACCGTAAAATGGTCTGAGAACGCGCGCATTGGTTACTACGCCCAGGACCATGAat                                                              >  1:2199537/1‑68 (MQ=255)
aCCGGACAGCGGCACCGTAAAATGGTCTGAGAACGCGCGCATTGGTTACTACGCCCAGGACCATGAat                                                              >  1:907626/1‑68 (MQ=255)
aCCGGACAGCGGCACCGTAAAATGGGCTGAGAACGCGCGCATTGGTTACTACGCCCAGGACCATGaa                                                               >  1:1777992/1‑67 (MQ=255)
          ggCACCGTAAAATGGTCTGAGAACGCGCGCATTGGTTACTATGCTCAGGACCACGAATATGAGTTTGaa                                                   <  1:1578945/69‑1 (MQ=255)
           gCACCGTAAAATGGTCTGAGAACGCGCGCATTGGTTACTATGCTCAGGACCACGAATATGAGTTTGaa                                                   <  1:147478/68‑1 (MQ=255)
               cGTAAAATGGTCTGAGAACGCGCGCATTGGTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGa                                              <  1:565340/69‑1 (MQ=255)
                          cTGAGAACGCGCGCATTGGTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCgtg                                    <  1:2507351/68‑1 (MQ=255)
                          cTGAGAACGCGCGCATTGGTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCgtg                                    <  1:2432672/68‑1 (MQ=255)
                          cTGAGAACGCGCGCATTGGTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCgtg                                    <  1:1041959/68‑1 (MQ=255)
                          cTGAGAACGCGCGCATTGGTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCgtg                                    <  1:2819744/68‑1 (MQ=255)
                           tGAGAACGCGCGCATTGGTTACTACGCCCAGGACCATGAATATGAGTTCGaa                                                   <  1:496673/52‑1 (MQ=255)
                           tGAGAACGCGCGCATTGGTTACTACGCCCAGGACCATGAATATGAGTTCGaa                                                   <  1:2673709/52‑1 (MQ=255)
                           tGAGAACGCGCGCATTGGTTACTACGCCCAGGACCATGAATATGAGTTCGaa                                                   <  1:685610/52‑1 (MQ=255)
                           tGAGAACGCGCGCATTGGTTACTACGCCCAGGACCATGAATATGAGTTCGaa                                                   <  1:946176/52‑1 (MQ=255)
                           tGAGAACGCGCGCATTGGTTACTACGCCCAGGACCATGAATATGAGTTCGaa                                                   <  1:2224282/52‑1 (MQ=255)
                           tGAGAACGCGCGCATTGGTTACTACGCCCAGGACCATGAATATGAGTTCGaa                                                   <  1:3363391/52‑1 (MQ=255)
                           tGAGAACGCGCGCATTGGTTACTACGCCCAGGACCATGAATATGAGTTCGaa                                                   <  1:1840407/52‑1 (MQ=255)
                           tGAGAACGCGCGCATTGGTTACTACGCCCAGGACCATGAATATGAGTTCGaa                                                   <  1:526454/52‑1 (MQ=255)
                           tGAGAACGCGCGCATTGGTTACTACGCCCAGGACCATGAATATGAGTTCGaa                                                   <  1:1515548/52‑1 (MQ=255)
                             agaACGCGCGCATTGGTTACTACGCCCAGTACCATGAATATGAGTTCGaa                                                   <  1:2755516/50‑1 (MQ=25)
                             agaACGCGCGCATTGGTTACTACGCCCAGGACCATGAATATGAGTTCGaa                                                   <  1:2842147/50‑1 (MQ=255)
                             agaACGCGCGCATTGGTTACTACGCCCAGGACCATGAATATGAGTTCGaa                                                   <  1:153281/50‑1 (MQ=255)
                                  gcgcgcATTGGTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATgg                           <  1:804114/69‑1 (MQ=255)
                                   cgcgcATTGGTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATgg                           >  1:3177855/1‑68 (MQ=255)
                                   cgcgcATTGGTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATgg                           >  1:3216270/1‑68 (MQ=255)
                                   cgcgcATTGGTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATgg                           >  1:344198/1‑68 (MQ=255)
                                   cgcgcATTGGTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATgg                           >  1:583532/1‑68 (MQ=255)
                                   cgcgcATTGGTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATgg                           >  1:2336788/1‑68 (MQ=255)
                                   cgcgcATTGGTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATgg                           >  1:1397716/1‑68 (MQ=255)
                                       cATTGGTTACTACGCCCAGGACCATGAATATGAGTTCGaa                                                   <  1:1722129/40‑1 (MQ=255)
                                           ggTTACTATGCTCAGGACCACGAATATTAGTTTGAAAATGATCTGAcc                                       >  1:1803526/1‑48 (MQ=255)
                                           ggTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCg                                      >  1:2319037/1‑49 (MQ=255)
                                           ggTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATGGATGAGCCAg                  >  1:574017/1‑69 (MQ=255)
                                           ggTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATGGATGAGCCAg                  >  1:713353/1‑69 (MQ=255)
                                           ggTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATGGATGAGCCAg                  >  1:78299/1‑69 (MQ=255)
                                           ggTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATGGATGAGCCAg                  >  1:91978/1‑69 (MQ=255)
                                           ggTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATGGATGAGCCAg                  >  1:2701799/1‑69 (MQ=255)
                                           ggTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATGGATGAGCCAg                  >  1:449855/1‑69 (MQ=255)
                                           ggTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATGGATGAGCCAg                  >  1:1817240/1‑69 (MQ=255)
                                           ggTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATGGATGAGCCAg                  >  1:1295238/1‑69 (MQ=255)
                                           ggTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATGGATGAGCCAg                  >  1:1223567/1‑69 (MQ=255)
                                                     cTCAGGACCATGAATATGAGTTTGAAAATGATCTGACCGTGTTCg                                >  1:234811/1‑45 (MQ=255)
                                                     cTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCg                                >  1:3097152/1‑45 (MQ=255)
                                                     cTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCg                                >  1:303606/1‑45 (MQ=255)
                                                     cTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCg                                >  1:2709711/1‑45 (MQ=255)
                                                     cTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCg                                >  1:939211/1‑45 (MQ=255)
                                                     cTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCg                                >  1:977920/1‑45 (MQ=255)
                                                     cTCAGGACAACGAATATGAGTTTGAAAATGATCTGACCGTGTTCg                                >  1:3335119/1‑45 (MQ=255)
                                                     atCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCg                                >  1:1967808/2‑45 (MQ=255)
                                                           aCCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATGGATGAGCCAGTGGAAGCAGGAAGGCg  >  1:1013128/1‑69 (MQ=255)
                                                           aCCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATGGATGAGCCAGTGGAAGCAGGAAGGCg  >  1:1716708/1‑69 (MQ=255)
                                                           aCCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATGGATGAGCCAGTGGAAGCAGGAAGGCg  >  1:1061148/1‑69 (MQ=255)
                                                           aCCACGAATATGAGTTTGAAAAAGATCTGACCGTGTTCGAATGGATGAGCCAGt                 >  1:1015325/1‑54 (MQ=255)
                                                               |                                                                
ACCGGACAGCGGCACCGTAAAATGGTCTGAGAACGCGCGCATTGGTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATGGATGAGCCAGTGGAAGCAGGAAGGCG  >  minE/547989‑548116

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: