Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,203,776 T→C 73.2% A211A (GCA→GCG aaeA ← p‑hydroxybenzoic acid efflux system component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,203,7760TC73.2% 70.3 / 22.4 63A211A (GCA→GCGaaeAp‑hydroxybenzoic acid efflux system component
Reads supporting (aligned to +/- strand):  ref base T (7/10);  new base C (34/12);  total (41/22)
Fisher's exact test for biased strand distribution p-value = 2.04e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

CCTGCGGCAACACTATCAACAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCTGCACGATACCCCGGACGCACCCCTTCCAGCTTAGTTTCTTCCATATAGGCCAGTACATAGAAGG  >  minE/2203714‑2203840
                                                              |                                                                
ccTGCGGCAACACTATCAACAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCTGCACGa                                                            <  1:3426934/69‑1 (MQ=255)
ccTGCGGCAACACTATCAACAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCTGCACGa                                                            <  1:605284/69‑1 (MQ=255)
  tGCGGCAACACTATCAACAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCTGCACGATa                                                          >  1:496241/1‑69 (MQ=255)
  tGCGGCAACACTATCAACAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCTGCACGATa                                                          >  1:1957403/1‑69 (MQ=255)
         acacTATCAATAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCCGCACGATACCCAggg                                                   >  1:2368976/1‑68 (MQ=255)
         acacTATCAACAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCCGCAc                                                              >  1:1823661/1‑58 (MQ=255)
         acacTATCAACAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCCGCACGATACCCAggg                                                   >  1:990030/1‑68 (MQ=255)
         acacTATCAACAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCCGCACGATACCCAggg                                                   >  1:1746253/1‑68 (MQ=255)
         acacTATCAACAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCCGCACGATACCCAggg                                                   >  1:1647840/1‑68 (MQ=255)
         acacTATCAACAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCCGCACGATACCCAggg                                                   >  1:14608/1‑68 (MQ=255)
         acacTATCAACAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCCGCACGATACCCAggg                                                   >  1:58528/1‑68 (MQ=255)
         acacTATCAACAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCCGCACGATACCCAggg                                                   >  1:1092337/1‑68 (MQ=255)
         acacTATCAACAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCCGCACGATACCCAggg                                                   >  1:689237/1‑68 (MQ=255)
         acacTATCAACAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCCGCACGATACCCAggg                                                   >  1:3051332/1‑68 (MQ=255)
              aTCAACAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCTGCACGATACCCCGGACGc                                                >  1:1272594/1‑67 (MQ=255)
              aTCAACAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCTGCACGATACCCCGGACGc                                                >  1:1730831/1‑67 (MQ=255)
              aTCAACAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCTGCACGATACCCCGGACGc                                                >  1:1794730/1‑67 (MQ=255)
              aTCAACAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCTGCACGATACCCCGGACGc                                                >  1:368113/1‑67 (MQ=255)
                   cAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCt                                          <  1:2586016/68‑1 (MQ=255)
                   cAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCt                                          <  1:1906910/68‑1 (MQ=255)
                   cAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCt                                          <  1:1846626/68‑1 (MQ=255)
                   cAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCt                                          <  1:508569/68‑1 (MQ=255)
                   cAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCt                                          <  1:1800230/68‑1 (MQ=255)
                   cAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCt                                          <  1:1799692/68‑1 (MQ=255)
                   cAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCt                                          <  1:1602812/68‑1 (MQ=255)
                   cAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCt                                          <  1:1583028/68‑1 (MQ=255)
                   cAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCt                                          <  1:153517/68‑1 (MQ=255)
                   cAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCt                                          <  1:3399794/68‑1 (MQ=255)
                   cAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCt                                          <  1:2716515/68‑1 (MQ=255)
                   cAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCt                                          <  1:1099978/68‑1 (MQ=255)
                                   tttGTTACTGCCAAGCGGCGTGATCTCTGCACGTTACCCCGGACGCACCCCTTCCAGCTTAGTttcttc                         <  1:2062339/69‑1 (MQ=255)
                                   tttGTTACTGCCAAGCGGCGTGATCTCTGCACGATACCCCGGACGCACCCCTTCCAGCTTAGTttcttc                         <  1:2812274/69‑1 (MQ=255)
                                       ttACTGCCAAGCGGCGTGATCTCTGCACGATACCCCGGACGCACCCCTTCCAGCTTAGTTTCTTCCa                       <  1:2959027/67‑1 (MQ=255)
                                       ttACTGCCAAGCGGCGTGATCTCTGCACGATACCCCGGACGCACCCCTTCCAGCTTAGTTTCTTCCa                       <  1:2849249/67‑1 (MQ=255)
                                       ttACTGCCAAGCGGCGTGATCTCTGCACGATACCCCGGACGCACCCCTTCCAGCTTAGTTTCTTCCa                       <  1:1600577/67‑1 (MQ=255)
                                       ttACTGCCAAGCGGCGTGATCTCTGCACGATACCCCGGACGCACCCCTTCCAGCTTAGTTTCTTCCa                       <  1:585671/67‑1 (MQ=255)
                                       ttACTGCCAAGCGGCGTGATCTCTGCACGATACCCCGGACGCACCCCTTCCAGCTTAGTTTCTTCCa                       <  1:1241754/67‑1 (MQ=255)
                                        tACTGCCAAGCGGCGTGATCTCTGCACGATACCCCGGACGCACCCCTTCCAGCTTAGTTTCTTCCa                       <  1:1254944/66‑1 (MQ=255)
                                          cTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCCCCCCTTCCAGCTTAGTTTCTTCCata                     >  1:3066180/1‑66 (MQ=255)
                                          cTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCt                                          >  1:105882/1‑45 (MQ=38)
                                          cTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCata                     >  1:47851/1‑66 (MQ=255)
                                          cTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCata                     >  1:464206/1‑66 (MQ=255)
                                          cTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCata                     >  1:514251/1‑66 (MQ=255)
                                          cTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCata                     >  1:375299/1‑66 (MQ=255)
                                          cTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCata                     >  1:568058/1‑66 (MQ=255)
                                          cTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCata                     >  1:365302/1‑66 (MQ=255)
                                          cTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCata                     >  1:3068435/1‑66 (MQ=255)
                                          cTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCata                     >  1:971643/1‑66 (MQ=255)
                                          cTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCata                     >  1:1007569/1‑66 (MQ=255)
                                          cTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCata                     >  1:998234/1‑66 (MQ=255)
                                          cTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCata                     >  1:229633/1‑66 (MQ=255)
                                          cTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCata                     >  1:1108314/1‑66 (MQ=255)
                                          cTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCata                     >  1:118479/1‑66 (MQ=255)
                                          cTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCata                     >  1:1261897/1‑66 (MQ=255)
                                          cTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCata                     >  1:1379276/1‑66 (MQ=255)
                                          cTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCata                     >  1:1458196/1‑66 (MQ=255)
                                          cTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCata                     >  1:2181733/1‑66 (MQ=255)
                                          cTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCata                     >  1:2469986/1‑66 (MQ=255)
                                          cTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCata                     >  1:2625635/1‑66 (MQ=255)
                                          cTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCata                     >  1:2804025/1‑66 (MQ=255)
                                          cTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCata                     >  1:2640797/1‑66 (MQ=255)
                                          cTGCCAAGCGGCGTGATCTCCGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCat                      >  1:1114209/1‑65 (MQ=255)
                                                         aTCTCTGCACGATACCCCGGACGCACCCCTTCCAGCTTAGTTTCTTCCATATAGGCCAGTACATAGa     >  1:2718967/1‑67 (MQ=255)
                                                           ctccGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCATATAGGCCAGTACATAGAAgg  <  1:2085181/68‑1 (MQ=255)
                                                           ctccGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCATATAGGCCAGTACATAGAAgg  <  1:2881522/68‑1 (MQ=255)
                                                           ctccGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCATATAGGCCAGTACATAGAAgg  <  1:1118419/68‑1 (MQ=255)
                                                           ctccGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCATATAGGCCAGTACATAGAAgg  <  1:333694/68‑1 (MQ=255)
                                                           ctccGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCATATAGGCCAGTACATAGAAgg  <  1:179602/68‑1 (MQ=255)
                                                           ctccGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCATATAGGCCAGTACATAGAAgg  <  1:261790/68‑1 (MQ=255)
                                                           ctccGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCATATAGGCCAGTACATAGAAgg  <  1:2151997/68‑1 (MQ=255)
                                                           ctccGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCATATAGGCCAGTACATAGAAgg  <  1:2177981/68‑1 (MQ=255)
                                                           ctccGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCATATAGGCCAGTACATAGAAgg  <  1:3354799/68‑1 (MQ=255)
                                                           ctccGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCATATAGGCCAGTACATAGAAgg  <  1:2890130/68‑1 (MQ=255)
                                                           ctccGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCATATAGGCCAGTACATAGAAgg  <  1:2582768/68‑1 (MQ=255)
                                                           ctccGCACGATACCCAGGGCGCACCCCTTCCAGCTTAGTTTCTTCCATATAGGCCAGTACATAGAAgg  <  1:2615564/68‑1 (MQ=255)
                                                              |                                                                
CCTGCGGCAACACTATCAACAGTCCCTTTCAGCACTTTGTTACTGCCAAGCGGCGTGATCTCTGCACGATACCCCGGACGCACCCCTTCCAGCTTAGTTTCTTCCATATAGGCCAGTACATAGAAGG  >  minE/2203714‑2203840

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: