Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 106,019 G→C 17.3% A91A (GCC→GCG aroP ← aromatic amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE106,0190GC17.3% 119.7 / 15.2 52A91A (GCC→GCGaroParomatic amino acid transporter
Reads supporting (aligned to +/- strand):  ref base G (27/16);  new base C (2/7);  total (29/23)
Fisher's exact test for biased strand distribution p-value = 6.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGAC  >  minE/105953‑106081
                                                                  |                                                              
ggCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaa                                                            <  1:183838/71‑1 (MQ=255)
 gCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCg                                                                >  1:2009484/1‑66 (MQ=255)
  cAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTATCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaaaa                                                          >  1:244246/1‑71 (MQ=255)
  cAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaaaa                                                          >  1:2890682/1‑71 (MQ=255)
  cAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaaaa                                                          >  1:3340757/1‑71 (MQ=255)
  cAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaaaa                                                          >  1:3464949/1‑71 (MQ=255)
  cAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaaaa                                                          >  1:3502966/1‑71 (MQ=255)
   ggTCAGCTCTGCCATTCCCACCAGCACGAACATTACCCAGTAGTTCCAGCCAGAGAGGAAGCCCGCaaacg                                                         <  1:1526286/70‑3 (MQ=255)
    gTCAGCTCTGCCATTCCCACCAGCACGAACATTACCCAGTAGTTCCAGCCAGAGAGGAAGCCCGCaaacgg                                                        >  1:3006078/1‑68 (MQ=255)
    gTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGc                                                              <  1:3549408/65‑1 (MQ=255)
      cagctcagcCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTg                                                       <  1:1434863/70‑1 (MQ=255)
       agctcagcCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGcc                                                     >  1:2651705/1‑71 (MQ=255)
        gctcagcCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaaaa                                                          >  1:489770/1‑65 (MQ=255)
        gctcagcCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaaaa                                                          >  1:1616064/1‑65 (MQ=255)
            agcCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAg                                                 <  1:2289323/70‑1 (MQ=255)
                aTGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTAttt                                            <  1:1976010/71‑1 (MQ=255)
                   gCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTtat                                          >  1:1507965/1‑70 (MQ=255)
                     aaCTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAg                                       >  1:725031/1‑71 (MQ=255)
                       cTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGc                                      >  1:4028542/1‑70 (MQ=255)
                       cTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGc                                      >  1:4870818/1‑70 (MQ=255)
                        tAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaa                                                            <  1:3970643/47‑1 (MQ=255)
                        tAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCa                                     >  1:4298674/1‑70 (MQ=255)
                           caCGAACATTACCCAGTAGTTCCAGCCAGAGAGGAAGCCCGCa                                                             <  1:2035627/42‑1 (MQ=16)
                           caCGAACATTACCCAGTAGTTCCAGCCAGAGAGGAAGCCCGCa                                                             <  1:1986627/42‑1 (MQ=16)
                           aaCGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGt                                 <  1:1876714/71‑1 (MQ=255)
                            aCGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTg                                >  1:2065810/1‑71 (MQ=255)
                            aCGAACATTACCCAGTAGTTCCAGCCAGAGAGGAAGCCCGCAAACGGTCCCCAGTATTTATAGGCaa                                    <  1:4883786/67‑1 (MQ=255)
                             cGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGccc                                                    <  1:1543649/50‑1 (MQ=255)
                                aCATTACCCAGTAGTTCCAGCCAGAGAGGAAGCCCGCAAACGGTCCCCAGTATTTATAGGCAAAGTgggca                            <  1:1351023/71‑4 (MQ=255)
                                aCAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCa                                     <  1:1130490/62‑1 (MQ=255)
                                aCAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCa                                     <  1:5635643/62‑1 (MQ=255)
                                 cAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTaa                           <  1:4489955/71‑1 (MQ=255)
                                  attACCCAGTAGTTCCAGCCAGAGAGGAAGCCCGCAAACGGTCCCCAGTATTtat                                          <  1:336187/53‑1 (MQ=255)
                                   gTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAg                         <  1:2718070/71‑1 (MQ=255)
                                     aCCCAGTAGTTCCAGCCAGAGAGGAAGCCCGCAAACGGTCCCCAGTATTta                                           >  1:3805356/1‑51 (MQ=255)
                                      cccAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAg                                  >  1:4906245/1‑59 (MQ=255)
                                      cccAGTAGTTCCAGCCAGAGAGGAAGCCCGCAAACGGTCCCCAGTATTTATAGGCa                                     <  1:4990146/56‑1 (MQ=255)
                                        cAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGcc                    >  1:2739150/1‑71 (MQ=255)
                                         agtagtTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAg                                  >  1:3061920/1‑56 (MQ=255)
                                              ttCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCa                                     <  1:1213210/48‑1 (MQ=255)
                                              ttCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGAc              <  1:361639/71‑1 (MQ=255)
                                              ttCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGAc              <  1:3509339/71‑1 (MQ=255)
                                               tCCAGCCAGAGGCGAAACCGGCTAAACTGCCCCAGTATTTATAAGCAAAGTGGCTaaa                          >  1:2531373/1‑58 (MQ=255)
                                                ccagccagAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTaaa                          >  1:5387778/1‑57 (MQ=255)
                                                ccagccagAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTaaa                          >  1:1183856/1‑57 (MQ=255)
                                                ccagccagAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAg            >  1:4309736/1‑71 (MQ=255)
                                                ccagccagAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAg            >  1:1694569/1‑71 (MQ=255)
                                                ccagccagAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAg            >  1:994333/1‑71 (MQ=255)
                                                     cagAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGttct       >  1:1899357/1‑71 (MQ=255)
                                                      agagGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAg                                  >  1:963747/1‑43 (MQ=255)
                                                       gagGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTa                            <  1:3860395/48‑1 (MQ=255)
                                                          gCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGac  >  1:4109909/1‑71 (MQ=255)
                                                          gCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGac  >  1:2519690/1‑71 (MQ=255)
                                                                  |                                                              
GGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGAC  >  minE/105953‑106081

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: