Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 106,012 C→G 15.7% A94P (GCC→CCC) ‡ aroP ← aromatic amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE106,0120CG15.7% 137.9 / 15.5 57A94P (GCC→CCC) ‡aroParomatic amino acid transporter
Reads supporting (aligned to +/- strand):  ref base C (28/20);  new base G (2/7);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 7.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

TTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGAC  >  minE/105944‑106081
                                                                    |                                                                     
tttACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGa                                                                     <  1:2831510/71‑1 (MQ=255)
         ggCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaa                                                            <  1:183838/71‑1 (MQ=255)
          gCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCg                                                                >  1:2009484/1‑66 (MQ=255)
           cAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTATCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaaaa                                                          >  1:244246/1‑71 (MQ=255)
           cAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaaaa                                                          >  1:3502966/1‑71 (MQ=255)
           cAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaaaa                                                          >  1:3464949/1‑71 (MQ=255)
           cAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaaaa                                                          >  1:3340757/1‑71 (MQ=255)
           cAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaaaa                                                          >  1:2890682/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)
               cagctcagcCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGc                                                                       >  1:1838204/1‑54 (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)
                   tcagcCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGa                                                                     <  1:1497002/52‑1 (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:4870818/1‑70 (MQ=255)
                                cTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGc                                      >  1:4028542/1‑70 (MQ=255)
                                 tAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAAc                                                                  <  1:1885376/41‑1 (MQ=255)
                                 tAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAAc                                                                  <  1:64571/41‑1 (MQ=255)
                                 tAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCaa                                                            <  1:3970643/47‑1 (MQ=255)
                                 tAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCa                                     >  1:4298674/1‑70 (MQ=255)
                                    caCGAACATTACCCAGTAGTTCCAGCCAGAGAGGAAGCCCGCa                                                             <  1:1986627/42‑1 (MQ=16)
                                    caCGAACATTACCCAGTAGTTCCAGCCAGAGAGGAAGCCCGCa                                                             <  1:2035627/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:3509339/71‑1 (MQ=255)
                                                       ttCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGAc              <  1:361639/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)
                                                                    |                                                                     
TTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGAC  >  minE/105944‑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: