Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,027,446 A→G 30.8% N134S (AAC→AGC)  malX → fused maltose and glucose‑specific PTS enzyme IIBC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,027,4460AG30.8% 30.2 / 18.6 26N134S (AAC→AGC) malXfused maltose and glucose‑specific PTS enzyme IIBC components
Reads supporting (aligned to +/- strand):  ref base A (9/9);  new base G (4/4);  total (13/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGTAAACTTCTGGTTGACCAATAAAGGCATTCTGCCAACCACGGATGCCGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGATC  >  minE/1027376‑1027489
                                                                      |                                           
gcgGTAAACTTCTGGTTGACCAATAAAGGCATTCTGCCAACCACGGATGCCGCGGTTCTGAAAGCCAATaa                                             >  1:1738975/1‑71 (MQ=255)
       aCTTCTGGTTGACCAATAAAGGCATTCTGCCAACCACGGATGCCGCGGTTCTGAAAGCCAATAACATCCag                                      <  1:209211/71‑1 (MQ=255)
          tCTGGTTGACCAATAAAGGCATTCTGCCAACCACGGATGCCGCGGTTCTGAAAGCCAATAACATCCagag                                    >  1:1052178/1‑70 (MQ=255)
           cTGGTTGACCAATAAAGGCATTCTGCCAACCACGGATGCCGCGGTTCTGAAAGCCAATAGCATCCAGAGCa                                  >  1:287880/1‑71 (MQ=255)
           cTGGTTGACCAATAAAGGCATTCTGCCAACCACGGATGCCGCGGTTCTGAAAGCCAATAACATCCagag                                    <  1:2061924/69‑1 (MQ=255)
                 gACCAATAAAGGCATTCTGCCAACCACGGATGCCGCGGTTCTGAAAGCCAATAGCATCCAGAGCATTCtt                             <  1:1731545/70‑1 (MQ=255)
                      aTAAAGGCATTCTGCCAACCACGGATGCCGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATc                       >  1:602019/1‑71 (MQ=255)
                      aTAAAGGCATTCTGCCAACCACGGATGCCGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATc                       <  1:1760587/71‑1 (MQ=255)
                      aTAAAGGCATTCTGCCAACCACGGATGCCGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATc                       <  1:737154/71‑1 (MQ=255)
                       tAAAGGCATTCTGCCAACCACGGATGCCGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATc                       <  1:718484/70‑1 (MQ=255)
                         aaGGCATTCTGCCAACCACGGATGCCGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATCCa                     <  1:1506200/70‑1 (MQ=255)
                         aaGGCATTCTGCCAACCACGGATGCCGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATCCa                     <  1:746430/70‑1 (MQ=255)
                            gCATTCTGCCAACCACGGATGCCGCGGTTCTGAAAGCCAATAGCATCCAGAGCATTCTTGGGATCCAGtcg                 <  1:1962751/71‑1 (MQ=255)
                               ttCTGCCAACCACGGATTCCGCGGTTCTGAAAGCCAATAGCATCCAGAGCATTCTTGGGATCCa                     <  1:1882941/64‑1 (MQ=255)
                                 cTGCCAACCACGGATGCCGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATCCAGtcgatcga            >  1:1841893/1‑71 (MQ=255)
                                 cTGCCAACCACGGATGCCGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATCCAGtcgatcga            >  1:1485260/1‑71 (MQ=255)
                                 cTGCCAACCACGGATGCCGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATCCAGtcgatcg             >  1:473532/1‑70 (MQ=255)
                                 cTGCCAACCACGGATGCCGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATCCAGtcg                 <  1:1272354/66‑1 (MQ=255)
                                         cACGGATGCCGCGGTTCTGAAAGCCAATAGCATCCAGAGCATTCTTGGGATCCAGTCGATCGATa          >  1:1052549/1‑65 (MQ=255)
                                         cACGGATGCCGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGa    >  1:131600/1‑71 (MQ=255)
                                          aCGGATGCCGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATCCAGtcgat               >  1:1739958/1‑59 (MQ=255)
                                           cGGATGCCGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGATc  <  1:805618/71‑1 (MQ=255)
                                                 ccGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCt                              >  1:1890070/1‑37 (MQ=255)
                                                       ttCTGAAAGCCAATAGCATCCAGAGCATTCTTGGGATc                       <  1:1418467/38‑1 (MQ=255)
                                                           gAAAGCCAATAGCATCCAGAGCATTCTTGGGATCCAGtcgatcta            >  1:194858/1‑43 (MQ=255)
                                                           gAAAGCCAATAGCATCCAGAGCATTCTTGGGATCCAGtcgatcga            >  1:637755/1‑45 (MQ=255)
                                                                  aaTAACATCCAGAGCATTCTTTGGATCCAGTCGATCGATAc         >  1:344359/1‑41 (MQ=255)
                                                                     aaCATCCAGAGCATTCTTGGGATCCAGtcgatcga            >  1:797074/1‑35 (MQ=255)
                                                                      |                                           
GCGGTAAACTTCTGGTTGACCAATAAAGGCATTCTGCCAACCACGGATGCCGCGGTTCTGAAAGCCAATAACATCCAGAGCATTCTTGGGATCCAGTCGATCGATACCGGGATC  >  minE/1027376‑1027489

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: