Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 548,043 T→C 37.7% A388A (GCT→GCC 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,0430TC37.7% 27.8 / 32.1 53A388A (GCT→GCCybiTfused predicted transporter subunits and ATP‑binding components of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base T (25/8);  new base C (7/13);  total (32/21)
Fisher's exact test for biased strand distribution p-value = 4.49e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

CTGCAACCGGACAGCGGCACCGTAAAATGGTCTGAGAACGCGCGCATTGGTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATGGATGAGCCAG  >  minE/547984‑548100
                                                           |                                                         
cTGCAACCGGACAGCGGCACCGTAAAATGGTCTGAGAACGCGCGCATTGGTTACTATGCTCAGGACCa                                                   >  1:2477657/1‑68 (MQ=255)
     aCCGGACAGCGGCACCGTAAAATGGTCTGAGAACGCGCGCATTGGTTACTACGCCCAGGACCATGAat                                              >  1:2711408/1‑68 (MQ=255)
     aCCGGACAGCGGCACCGTAAAATGGTCTGAGAACGCGCGCATTGGTTACTACGCCCAGGACCATGAat                                              >  1:996400/1‑68 (MQ=255)
     aCCGGACAGCGGCACCGTAAAATGGTCTGAGAACGCGCGCATTGGTTACTACGCCCAGGACCATGAat                                              >  1:907626/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)
     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:2432672/68‑1 (MQ=255)
                               cTGAGAACGCGCGCATTGGTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCgtg                    <  1:2819744/68‑1 (MQ=255)
                               cTGAGAACGCGCGCATTGGTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCgtg                    <  1:1041959/68‑1 (MQ=255)
                               cTGAGAACGCGCGCATTGGTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCgtg                    <  1:2507351/68‑1 (MQ=255)
                                tGAGAACGCGCGCATTGGTTACTACGCCCAGGACCATGAATATGAGTTCGaa                                   <  1:2224282/52‑1 (MQ=255)
                                tGAGAACGCGCGCATTGGTTACTACGCCCAGGACCATGAATATGAGTTCGaa                                   <  1:1515548/52‑1 (MQ=255)
                                tGAGAACGCGCGCATTGGTTACTACGCCCAGGACCATGAATATGAGTTCGaa                                   <  1:685610/52‑1 (MQ=255)
                                tGAGAACGCGCGCATTGGTTACTACGCCCAGGACCATGAATATGAGTTCGaa                                   <  1:3363391/52‑1 (MQ=255)
                                tGAGAACGCGCGCATTGGTTACTACGCCCAGGACCATGAATATGAGTTCGaa                                   <  1:1840407/52‑1 (MQ=255)
                                tGAGAACGCGCGCATTGGTTACTACGCCCAGGACCATGAATATGAGTTCGaa                                   <  1:2673709/52‑1 (MQ=255)
                                tGAGAACGCGCGCATTGGTTACTACGCCCAGGACCATGAATATGAGTTCGaa                                   <  1:526454/52‑1 (MQ=255)
                                tGAGAACGCGCGCATTGGTTACTACGCCCAGGACCATGAATATGAGTTCGaa                                   <  1:496673/52‑1 (MQ=255)
                                tGAGAACGCGCGCATTGGTTACTACGCCCAGGACCATGAATATGAGTTCGaa                                   <  1:946176/52‑1 (MQ=255)
                                  agaACGCGCGCATTGGTTACTACGCCCAGTACCATGAATATGAGTTCGaa                                   <  1:2755516/50‑1 (MQ=25)
                                  agaACGCGCGCATTGGTTACTACGCCCAGGACCATGAATATGAGTTCGaa                                   <  1:153281/50‑1 (MQ=255)
                                  agaACGCGCGCATTGGTTACTACGCCCAGGACCATGAATATGAGTTCGaa                                   <  1:2842147/50‑1 (MQ=255)
                                       gcgcgcATTGGTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATgg           <  1:804114/69‑1 (MQ=255)
                                        cgcgcATTGGTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATgg           >  1:3216270/1‑68 (MQ=255)
                                        cgcgcATTGGTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATgg           >  1:3177855/1‑68 (MQ=255)
                                        cgcgcATTGGTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATgg           >  1:344198/1‑68 (MQ=255)
                                        cgcgcATTGGTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATgg           >  1:1397716/1‑68 (MQ=255)
                                        cgcgcATTGGTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATgg           >  1:2336788/1‑68 (MQ=255)
                                        cgcgcATTGGTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATgg           >  1:583532/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:1223567/1‑69 (MQ=255)
                                                ggTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATGGATGAGCCAg  >  1:449855/1‑69 (MQ=255)
                                                ggTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATGGATGAGCCAg  >  1:2701799/1‑69 (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:1817240/1‑69 (MQ=255)
                                                ggTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATGGATGAGCCAg  >  1:1295238/1‑69 (MQ=255)
                                                ggTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATGGATGAGCCAg  >  1:91978/1‑69 (MQ=255)
                                                          cTCAGGACCATGAATATGAGTTTGAAAATGATCTGACCGTGTTCg                >  1:234811/1‑45 (MQ=255)
                                                          cTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCg                >  1:3097152/1‑45 (MQ=255)
                                                          cTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCg                >  1:977920/1‑45 (MQ=255)
                                                          cTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCg                >  1:939211/1‑45 (MQ=255)
                                                          cTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCg                >  1:2709711/1‑45 (MQ=255)
                                                          cTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCg                >  1:303606/1‑45 (MQ=255)
                                                          cTCAGGACAACGAATATGAGTTTGAAAATGATCTGACCGTGTTCg                >  1:3335119/1‑45 (MQ=255)
                                                          atCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCg                >  1:1967808/2‑45 (MQ=255)
                                                           |                                                         
CTGCAACCGGACAGCGGCACCGTAAAATGGTCTGAGAACGCGCGCATTGGTTACTATGCTCAGGACCACGAATATGAGTTTGAAAATGATCTGACCGTGTTCGAATGGATGAGCCAG  >  minE/547984‑548100

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: