Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 106,091 A→G 12.8% G67G (GGT→GGC aroP ← aromatic amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE106,0910AG12.8% 94.5 / 5.5 39G67G (GGT→GGCaroParomatic amino acid transporter
Reads supporting (aligned to +/- strand):  ref base A (17/17);  new base G (2/3);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTG  >  minE/106026‑106155
                                                                 |                                                                
cTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAg                                                             <  1:5270303/71‑1 (MQ=255)
  tccccAGTATTTATAGGCAAAGTGGGCAAATGAACCGGATACCGGCTCCTCAACCACCATTTCGCCAAGCt                                                           <  1:1085254/70‑1 (MQ=255)
  tccccAGTATTTATAGGCAAAGTGGGCAAATGAACCGGATACCGGCTCCTCAACCACCATTTCGCCAAGCt                                                           <  1:5349701/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)
          aTTTATAGGCAAAGTGGGCAAATGAACCGGATACCGGCTCCTCAACCACCATTTCGCCAAGCTGGCGCAt                                                    >  1:1419409/1‑70 (MQ=255)
           tttATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTg                                                          <  1:856727/63‑1 (MQ=255)
           tttATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATg                                                   <  1:5494395/70‑1 (MQ=255)
                    aaaGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCt                                                           >  1:3789774/1‑53 (MQ=255)
                    aaaGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAgg                                         >  1:4151564/1‑71 (MQ=255)
                       gTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTTACCCa                                                              <  1:2182475/47‑1 (MQ=255)
                        tGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGAt                                     <  1:178533/71‑1 (MQ=255)
                        tGGCTAAAGGAGCATGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGAt                                     <  1:482131/71‑1 (MQ=255)
                              aaGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAAc                                >  1:3476152/1‑70 (MQ=255)
                               aGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCa                              <  1:2706303/71‑1 (MQ=255)
                                 gAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGc                            <  1:3640472/71‑1 (MQ=255)
                                   gCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCaa                          >  1:3768620/1‑71 (MQ=255)
                                     cTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATg                        <  1:2768478/71‑1 (MQ=255)
                                     cTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATg                        <  1:815917/71‑1 (MQ=255)
                                     cTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATg                        <  1:452202/71‑1 (MQ=255)
                                     cTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATg                        <  1:168769/71‑1 (MQ=255)
                                      tGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGc                            <  1:2659560/66‑1 (MQ=255)
                                      tGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATgg                       >  1:924761/1‑71 (MQ=255)
                                             ggTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAAcc                >  1:4568779/1‑71 (MQ=255)
                                                tcttcGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAAcc                               >  1:3464436/1‑53 (MQ=255)
                                                tcttcGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAAc                                >  1:2033309/1‑52 (MQ=255)
                                                tcttcGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAg             >  1:705734/1‑71 (MQ=255)
                                                tcttcGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAg             >  1:846298/1‑71 (MQ=255)
                                                 cttcGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAg             <  1:4017139/70‑1 (MQ=255)
                                                  ttcGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGa           >  1:3414702/1‑71 (MQ=255)
                                                     aaccaccATTTCGCCAAGCTGGCGCATAATCAGGAAAGc                                        <  1:3470610/38‑1 (MQ=12)
                                                        caccaTTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGt                    <  1:1100088/56‑1 (MQ=255)
                                                             tttCACCCAGCTGACGCATGATCAGAAAGGCGATaa                                   <  1:1120233/36‑1 (MQ=255)
                                                              ttCACCCAGCTGACGCATGATCAGAAAGGCGATaaaa                                 >  1:4942523/1‑37 (MQ=255)
                                                              ttCACCCAGCTGACGCATGATCAGAAAGGCGATAAAAc                                >  1:4069261/1‑38 (MQ=255)
                                                              ttCACCCAGCTGACGCATGATCAGAAAGGCGATAAAAc                                >  1:89526/1‑38 (MQ=255)
                                                               tCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGTATAATCCCTg  >  1:3089723/1‑67 (MQ=255)
                                                               tCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATcc     >  1:2811649/1‑64 (MQ=255)
                                                                 aCCCATCTGACGCATGATCAGAAAGGCGATAAAAcc                               <  1:953404/36‑1 (MQ=255)
                                                                 aCCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAAt       <  1:2547607/60‑1 (MQ=255)
                                                                 |                                                                
CTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTG  >  minE/106026‑106155

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: