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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE106,0940CT72.7% 75.5 / 29.0 73L66L (CTG→CTAaroParomatic amino acid transporter
Reads supporting (aligned to +/- strand):  ref base C (11/9);  new base T (19/34);  total (30/43)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

CCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGC  >  minE/106031‑106157
                                                               |                                                               
ccAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTg                                                            <  1:1232346/69‑1 (MQ=255)
            aGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTg                                                            <  1:2528664/57‑1 (MQ=255)
              cAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCTAGCTGACGCATGATCAGaa                                              <  1:2350698/69‑1 (MQ=255)
              cAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCTAGCTGACGCATGATCAGaa                                              <  1:3212050/69‑1 (MQ=255)
              cAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCTAGCTGACGCATGATCAGaa                                              <  1:2071360/69‑1 (MQ=255)
              cAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCTAGCTGACGCATGATCAGaa                                              <  1:1533180/69‑1 (MQ=255)
              cAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCTAGCTGACGCATGATCAGaa                                              <  1:1463506/69‑1 (MQ=255)
              cAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCTAGCTGACGCATGATCAGaa                                              <  1:212345/69‑1 (MQ=255)
              cAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCTAGCTGACGCATGATCAGaa                                              <  1:290953/69‑1 (MQ=255)
              cAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCTAGCTGACGCATGATCAGaa                                              <  1:554447/69‑1 (MQ=255)
              cAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCTAGCTGACGCATGATCAGaa                                              <  1:1240875/69‑1 (MQ=255)
              cAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCTAGCTGACGCATGATCAGaa                                              <  1:2491077/69‑1 (MQ=255)
              cAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCTAGCTGACGCATGATCAGaa                                              <  1:781421/69‑1 (MQ=255)
              cAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCTAGCTGACGCATGATCAGaa                                              <  1:1160444/69‑1 (MQ=255)
              cAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCTAGCTGACGCATGATCAGaa                                              <  1:2317411/69‑1 (MQ=255)
                   tGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGc                                          >  1:127260/1‑68 (MQ=255)
                    ggCTAAAGGAGCCTGCGACAGGTTCTTCGCCCACCATTTCACCCAGCTGACGCATGATCAGAAAGGc                                          >  1:1680934/1‑67 (MQ=255)
                    ggCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCg                                         >  1:2721557/1‑68 (MQ=255)
                    ggCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCg                                         >  1:3323967/1‑68 (MQ=255)
                    ggCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCg                                         >  1:1488725/1‑68 (MQ=255)
                    ggCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCg                                         >  1:934510/1‑68 (MQ=255)
                    ggCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCg                                         >  1:2096631/1‑68 (MQ=255)
                    ggCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCg                                         >  1:2951654/1‑68 (MQ=255)
                    ggCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCg                                         >  1:1922506/1‑68 (MQ=255)
                          aGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTg                                                            <  1:1838021/43‑1 (MQ=255)
                           ggAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAAcc                                 >  1:977258/1‑69 (MQ=255)
                           ggAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAAcc                                 >  1:975581/1‑69 (MQ=255)
                              gCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGc                              <  1:321839/69‑1 (MQ=255)
                              gCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGc                              <  1:3238415/69‑1 (MQ=255)
                              gCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGc                              <  1:2865807/69‑1 (MQ=255)
                              gCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGc                              <  1:555357/69‑1 (MQ=255)
                              gCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGc                              <  1:1189702/69‑1 (MQ=255)
                              gCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGc                              <  1:2402817/69‑1 (MQ=255)
                                                 accaccATTTCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:3056900/68‑1 (MQ=255)
                                                 accaccATTTCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:2946373/68‑1 (MQ=255)
                                                 accaccATTTCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:3243939/68‑1 (MQ=255)
                                                 accaccATTTCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:3278875/68‑1 (MQ=255)
                                                 accaccATTTCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:457659/68‑1 (MQ=255)
                                                 accaccATTTCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:272721/68‑1 (MQ=255)
                                                 accaccATTTCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:59390/68‑1 (MQ=255)
                                                 accaccATTTCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:1053930/68‑1 (MQ=255)
                                                 accaccATTTCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:1218553/68‑1 (MQ=255)
                                                 accaccATTTCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:1829678/68‑1 (MQ=255)
                                                 accaccATTTCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:2371888/68‑1 (MQ=255)
                                                 accaccATTTCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:1963366/68‑1 (MQ=255)
                                                 accaccATTTCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:1443058/68‑1 (MQ=255)
                                                 accaccATTTCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:1417561/68‑1 (MQ=255)
                                                 accaccATTTCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:1369079/68‑1 (MQ=255)
                                                 accaccATTTCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:2210392/68‑1 (MQ=255)
                                                 accaccATTTCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:1196817/68‑1 (MQ=255)
                                                 accaccATTTCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:225792/68‑1 (MQ=255)
                                                 accaccATTTCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:1023595/68‑1 (MQ=255)
                                                 accaccATTTCACCTAGCTGACGCATGATCAGAAAGGAGATAAAACCAGCAATGGCGTAAACCAGGAt            <  1:571299/68‑1 (MQ=255)
                                                   caccaTTTCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:960183/66‑1 (MQ=255)
                                                          tCACCTAGCTGACGCATGATCAGAAGGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTgg   >  1:3048116/1‑68 (MQ=255)
                                                          tCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            >  1:3352979/1‑59 (MQ=255)
                                                          tCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTgg   >  1:1185890/1‑68 (MQ=255)
                                                          tCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGc  >  1:67060/1‑69 (MQ=255)
                                                          tCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGc  >  1:3234856/1‑69 (MQ=255)
                                                          tCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGc  >  1:1136193/1‑69 (MQ=255)
                                                          tCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGc  >  1:1713474/1‑69 (MQ=255)
                                                          tCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGc  >  1:1736841/1‑69 (MQ=255)
                                                          tCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGc  >  1:1099553/1‑69 (MQ=255)
                                                          tCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGc  >  1:3140701/1‑69 (MQ=255)
                                                          tCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGc  >  1:2317031/1‑69 (MQ=255)
                                                          tCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGc  >  1:2465475/1‑69 (MQ=255)
                                                          tCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGc  >  1:246858/1‑69 (MQ=255)
                                                          tCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGc  >  1:2530115/1‑69 (MQ=255)
                                                          tCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGc  >  1:2543282/1‑69 (MQ=255)
                                                          tCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGc  >  1:2207139/1‑69 (MQ=255)
                                                          tCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGc  >  1:290155/1‑69 (MQ=255)
                                                          tCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGc  >  1:2195819/1‑69 (MQ=255)
                                                          tCACCTAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGc  >  1:2928722/1‑69 (MQ=255)
                                                               cAGCTGACGCATGTTCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:2903647/54‑1 (MQ=255)
                                                               cAGCTGACGCATGATCAGAAAGGCTATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:773306/54‑1 (MQ=255)
                                                               cAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:3370666/54‑1 (MQ=255)
                                                               cAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:534067/54‑1 (MQ=255)
                                                               cAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:1269864/54‑1 (MQ=255)
                                                               cAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:2233073/54‑1 (MQ=255)
                                                               cAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:2705719/54‑1 (MQ=255)
                                                               cAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:802740/54‑1 (MQ=255)
                                                               cAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:831257/54‑1 (MQ=255)
                                                               cAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:867151/54‑1 (MQ=255)
                                                               cAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:2247268/54‑1 (MQ=255)
                                                               cAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:2086176/54‑1 (MQ=255)
                                                               cAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:2067030/54‑1 (MQ=255)
                                                               cAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGAt            <  1:1806251/54‑1 (MQ=255)
                                                               |                                                               
CCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGC  >  minE/106031‑106157

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: