Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 106,094 C→A 12.1% L66L (CTG→CTT aroP ← aromatic amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE106,0940CA12.1% 118.8 / 4.1 41L66L (CTG→CTTaroParomatic amino acid transporter
Reads supporting (aligned to +/- strand):  ref base C (17/19);  new base A (2/3);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.78e-01

CTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCG  >  minE/106026‑106163
                                                                    |                                                                     
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)
          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:452202/71‑1 (MQ=255)
                                     cTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATg                                <  1:815917/71‑1 (MQ=255)
                                     cTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATg                                <  1:2768478/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:89526/1‑38 (MQ=255)
                                                              ttCACCCAGCTGACGCATGATCAGAAAGGCGATAAAAc                                        >  1:4069261/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)
                                                                   ccAGCTGACGCATGATCAGAAAGGCGATAAAACCAg                                     >  1:2647180/1‑36 (MQ=255)
                                                                   ccAGCTGACGCATGATCAGAAAGGCGATAAAACCAg                                     >  1:1595452/1‑36 (MQ=255)
                                                                   ccAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAAt                                 >  1:3754345/1‑40 (MQ=255)
                                                                   ccAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCg  <  1:1093594/71‑1 (MQ=255)
                                                                    |                                                                     
CTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCG  >  minE/106026‑106163

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: