Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 106,043 A→G 16.6% A83A (GCT→GCC aroP ← aromatic amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE106,0430AG16.6% 88.4 / 11.4 42A83A (GCT→GCCaroParomatic amino acid transporter
Reads supporting (aligned to +/- strand):  ref base A (21/14);  new base G (2/5);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

AACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATG  >  minE/105974‑106106
                                                                     |                                                               
aaCTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAg                                                                >  1:725031/1‑71 (MQ=255)
  cTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGc                                                               >  1:4870818/1‑70 (MQ=255)
  cTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGc                                                               >  1:4028542/1‑70 (MQ=255)
   tAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCa                                                              >  1:4298674/1‑70 (MQ=255)
      aaCGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGt                                                          <  1:1876714/71‑1 (MQ=255)
       aCGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTg                                                         >  1:2065810/1‑71 (MQ=255)
       aCGAACATTACCCAGTAGTTCCAGCCAGAGAGGAAGCCCGCAAACGGTCCCCAGTATTTATAGGCaa                                                             <  1:4883786/67‑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)
              gTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAg                                                  <  1:2718070/71‑1 (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:1183856/1‑57 (MQ=255)
                           ccagccagAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTaaa                                                   >  1:5387778/1‑57 (MQ=255)
                           ccagccagAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAg                                     >  1:1694569/1‑71 (MQ=255)
                           ccagccagAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAg                                     >  1:994333/1‑71 (MQ=255)
                           ccagccagAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAg                                     >  1:4309736/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)
                                                aaaaCTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGaccac                        >  1:5197546/1‑63 (MQ=255)
                                                    cTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAg            <  1:5270303/71‑1 (MQ=255)
                                                      tccccAGTATTTATAGGCAAAGTGGGCAAATGAACCGGATACCGGCTCCTCAACCACCATTTCGCCAAGCt          <  1:5349701/70‑1 (MQ=255)
                                                      tccccAGTATTTATAGGCAAAGTGGGCAAATGAACCGGATACCGGCTCCTCAACCACCATTTCGCCAAGCt          <  1:1085254/70‑1 (MQ=255)
                                                      gCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCt          >  1:2700314/1‑71 (MQ=255)
                                                           aGTATTTATAGGCAAAGTGGGCAAATGAACCGGATACCGGCTCCTCAACCACCATTTCGCCAAGCTGGCg      >  1:1605763/1‑70 (MQ=255)
                                                           aGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGc     <  1:606703/71‑1 (MQ=255)
                                                            gTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCa                 >  1:805067/1‑58 (MQ=255)
                                                             tATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGttcttc                              <  1:2509002/44‑1 (MQ=255)
                                                              aTTTATAGGCAAAGTGGGCAAATGAACCGGATACCGGCTCCTCAACCACCATTTCGCCAAGCTGGCGCAt   >  1:1419409/1‑70 (MQ=255)
                                                               tttATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTg         <  1:856727/63‑1 (MQ=255)
                                                               tttATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATg  <  1:5494395/70‑1 (MQ=255)
                                                                     |                                                               
AACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATG  >  minE/105974‑106106

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: