Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,524,808 C→A 52.0% V189L (GTA→TTA)  yfcJ ← predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,524,8080CA52.0% ‑1.4 / 90.5 73V189L (GTA→TTA) yfcJpredicted transporter
Reads supporting (aligned to +/- strand):  ref base C (9/26);  new base A (15/23);  total (24/49)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.08e-01

GCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCA  >  minE/1524759‑1524872
                                                 |                                                                
gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTgtggtg                                               <  1:456956/69‑1 (MQ=255)
gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTgtggtg                                               <  1:2327818/69‑1 (MQ=255)
gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTAAgg                                                                >  1:1836110/1‑52 (MQ=255)
gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTAAGGGTa                                                             >  1:341945/1‑55 (MQ=255)
gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTAAGGGTAATGCCATTgtggtg                                               >  1:643693/1‑69 (MQ=255)
gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTAAGGGTAATGCCATTgtggtg                                               >  1:2740386/1‑69 (MQ=255)
gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTAAGGGTAATGCCATTgtggtg                                               >  1:1955249/1‑69 (MQ=255)
gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTAAGGGTAATGCCATTgtggtg                                               >  1:74986/1‑69 (MQ=255)
gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTAAGGGTAATGCCATTgtggtg                                               >  1:1467661/1‑69 (MQ=255)
gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTAAGGGTAATGCCATTgtggtg                                               >  1:818854/1‑69 (MQ=255)
gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTAAGGGTAATGCCATTgtggtg                                               >  1:1376443/1‑69 (MQ=255)
gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTAAGGGTAATGCCATTgtggtg                                               >  1:329928/1‑69 (MQ=255)
gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTAAGGGTAATGCCATTgtggtg                                               >  1:3407431/1‑69 (MQ=255)
gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTAAGGGTAATGCCATTgtggtg                                               >  1:2225111/1‑69 (MQ=255)
gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTAAGGGTAATGCCATTgtggtg                                               >  1:876530/1‑69 (MQ=255)
gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTAAGGGTAATGCCATTgtggtg                                               >  1:523422/1‑69 (MQ=255)
gCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTAAGGGTAATGCCATTgtggtg                                               >  1:2294517/1‑69 (MQ=255)
            aCTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGc                                    >  1:1719329/1‑68 (MQ=255)
            aCTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGc                                    >  1:573995/1‑68 (MQ=255)
            aCTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGc                                    >  1:543424/1‑68 (MQ=255)
            aCTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGc                                    >  1:2736595/1‑68 (MQ=255)
                tGCGCACTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTgtggt                                                >  1:868366/1‑52 (MQ=255)
                tGCGCACTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCaa                               >  1:886128/1‑69 (MQ=255)
                tGCGCACTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCaa                               >  1:1159716/1‑69 (MQ=255)
                tGCGCACTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCaa                               >  1:315549/1‑69 (MQ=255)
                tGCGCACTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCaa                               >  1:1217630/1‑69 (MQ=255)
                      cTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGt                         <  1:3237713/69‑1 (MQ=255)
                      cTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGt                         <  1:1115368/69‑1 (MQ=255)
                      cTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGt                         <  1:2271000/69‑1 (MQ=255)
                      cTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGt                         <  1:2162214/69‑1 (MQ=255)
                      cTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGt                         <  1:2313088/69‑1 (MQ=255)
                      cTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGt                         <  1:2054488/69‑1 (MQ=255)
                      cTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGt                         <  1:464076/69‑1 (MQ=255)
                      cTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGt                         <  1:3072062/69‑1 (MQ=255)
                      cTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGt                         <  1:2404191/69‑1 (MQ=255)
                      cTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGt                         <  1:690444/69‑1 (MQ=255)
                      cTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGt                         <  1:301375/69‑1 (MQ=255)
                      cTGTGCCGTAACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGt                         <  1:2038910/69‑1 (MQ=255)
                        gtgCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCaaa                              <  1:193489/62‑1 (MQ=255)
                        gtgCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCaaa                              <  1:122541/62‑1 (MQ=255)
                             gTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGTAATGGc                   <  1:1703169/68‑1 (MQ=255)
                             gTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGTAATGGc                   <  1:36994/68‑1 (MQ=255)
                             gTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGTAATGGc                   <  1:1168355/68‑1 (MQ=255)
                             gTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGTAATGGc                   <  1:858895/68‑1 (MQ=255)
                             gTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGTAATGGc                   <  1:1268736/68‑1 (MQ=255)
                             gTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGTAATGGc                   <  1:1496220/68‑1 (MQ=255)
                             gTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGTAATGGc                   <  1:1673372/68‑1 (MQ=255)
                             gTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGTAATGGc                   <  1:2931143/68‑1 (MQ=255)
                             gTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGTAATGGc                   <  1:1783030/68‑1 (MQ=255)
                              ttACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGTAATGGc                   <  1:2539593/67‑1 (MQ=255)
                                             aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:745721/69‑1 (MQ=255)
                                             aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:557377/69‑1 (MQ=255)
                                             aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:706657/69‑1 (MQ=255)
                                             aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:658395/69‑1 (MQ=255)
                                             aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:2274082/69‑1 (MQ=255)
                                             aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:2964725/69‑1 (MQ=255)
                                             aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:1175901/69‑1 (MQ=255)
                                             aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:1210673/69‑1 (MQ=255)
                                             aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:1224998/69‑1 (MQ=255)
                                             aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:145531/69‑1 (MQ=255)
                                             aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:1552615/69‑1 (MQ=255)
                                             aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:1694601/69‑1 (MQ=255)
                                             aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:1755998/69‑1 (MQ=255)
                                             aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:185659/69‑1 (MQ=255)
                                             aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:1977183/69‑1 (MQ=255)
                                             aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:2062857/69‑1 (MQ=255)
                                             aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:2091464/69‑1 (MQ=255)
                                             aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:2105658/69‑1 (MQ=255)
                                             aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:226710/69‑1 (MQ=255)
                                             aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:2379799/69‑1 (MQ=255)
                                             aGTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:45183/69‑1 (MQ=255)
                                              gTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:1012181/68‑1 (MQ=255)
                                              gTAAGGGTAATGCCATTGTGGTGAGCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCa  <  1:1104699/68‑1 (MQ=255)
                                                 |                                                                
GCCAGGGCCGGTACTTTGCGCACTGTGCCGTTACAGGCCCACGCCAGTACGGGTAATACCATTGTGGTGATCGCCAGTGCGGCAAAACCGTAATGGCTATGAATCAACAGGCCA  >  minE/1524759‑1524872

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