Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,130,599 C→G 8.6% L102L (CTG→CTC mtr ← tryptophan transporter of high affinity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,130,5990CG 8.6% 229.7 / 9.4 81L102L (CTG→CTCmtrtryptophan transporter of high affinity
Reads supporting (aligned to +/- strand):  ref base C (28/46);  new base G (5/2);  total (33/48)
Fisher's exact test for biased strand distribution p-value = 1.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

TAGTGACATCTCTGCGAAGGTGTGATGCAGAATCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTG  >  minE/2130536‑2130668
                                                               |                                                                     
tAGTGACATCTCTGCGAAGGTGTGATGCAGAATCGAACCACTGGCAGAAATATAGGCATAGGTCAGGata                                                                 >  1:52062/1‑70 (MQ=255)
  gTGACATCTCTGCGAAGGTGTGATGCAGAATCGAACCACTGGCAGAAATATAGGCATAGGTCAGGata                                                                 >  1:3497767/1‑68 (MQ=255)
    gACATCTCTGCGAAGGTGTGATGCAGAATCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATagag                                                            >  1:2327900/1‑71 (MQ=255)
     aCATCTCTGCGAAGGTGTGATGCAGAATCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATagag                                                            >  1:238406/1‑70 (MQ=255)
      cATCTCTGCGAAGGTGTGATGCAGAATCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGca                                                          <  1:829437/71‑1 (MQ=255)
       aTCTCTGCGAAGGTGTGATGCAGAATCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGca                                                          >  1:2079159/1‑70 (MQ=255)
           ctGCGAAGGTGTGATGCAGAATCGAACCACTGGCAGAAATATAGGCATAGGTCa                                                                      >  1:2219800/1‑54 (MQ=255)
            tGCGAAGGTGTGATGCAGAATCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGcaca                                                        <  1:1317160/67‑1 (MQ=255)
            tGCGAAGGTGTGATGCAGAATCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGcaca                                                        <  1:3570034/67‑1 (MQ=255)
            tGCGAAGGTGTGATGCAGAATCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGcaca                                                        <  1:2026017/67‑1 (MQ=255)
               gAAGGTGTGATGCAGAATCGAACCACTGGCAGAAATATAGGCATAGGGCAGGATATaga                                                             <  1:2413942/59‑1 (MQ=255)
                  ggTGTGATGCAGAATCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATg                                               >  1:2699596/1‑70 (MQ=255)
                   gtgtGATGCAGAATCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATgg                                              >  1:1830348/1‑70 (MQ=255)
                     gtgATTCAGAATCGAACCACTGGCAGAAATATAGGCATAGGTCa                                                                      <  1:974314/44‑1 (MQ=255)
                     gtgATGCAGAATCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACaa                                                       <  1:1000448/59‑1 (MQ=255)
                      tgATGCAGAATCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCaa                                                 >  1:1259575/1‑64 (MQ=255)
                      tgATGCAGAATCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAAt                                          >  1:816929/1‑71 (MQ=255)
                      tgAAGCAGAATCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCa                                                  >  1:1606479/1‑63 (MQ=255)
                       gATGCAGAATCGAACCACTGGCAGAAATATAGGCATAGGTCAg                                                                     <  1:901283/43‑1 (MQ=255)
                        aTGCAGAATCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGc                                        >  1:3756300/1‑71 (MQ=255)
                           cAGAATCGAACCACTGGCAGAAATATAGGCATAGGTCAGGa                                                                   >  1:1911112/1‑41 (MQ=255)
                           cAGAATCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACaa                                                       >  1:1656897/1‑53 (MQ=255)
                           cAGAATCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCg                                      <  1:3419028/70‑1 (MQ=255)
                           cAGAATCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCg                                      <  1:2476653/70‑1 (MQ=255)
                           cAGAATCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCg                                      <  1:313305/70‑1 (MQ=255)
                            aGAATCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGtt                                    >  1:3841630/1‑71 (MQ=255)
                             gAATCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGtt                                    >  1:3882045/1‑70 (MQ=255)
                             gAATCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTg                                   <  1:3419420/71‑1 (MQ=255)
                              aaTCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCa                                                  <  1:2064659/55‑1 (MQ=255)
                              aaTCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCa                                                  <  1:1998809/55‑1 (MQ=255)
                               aTCGCACCATTAGCAGAGATATAGGCATAAGTGAGGATATAGAGAACGAAGGCAACGGTAATACCGCTGa                                  >  1:2638131/1‑70 (MQ=255)
                               aTCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGt                                     <  1:220924/67‑1 (MQ=255)
                               aTCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGt                                     <  1:2570390/67‑1 (MQ=255)
                                tCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAAt                                                >  1:615111/1‑55 (MQ=255)
                                tCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGAc                                 >  1:2752857/1‑70 (MQ=255)
                                tCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGAc                                 >  1:1333389/1‑70 (MQ=255)
                                 cGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGaa                                            <  1:2238363/58‑1 (MQ=255)
                                  gAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCAc                              >  1:431476/1‑71 (MQ=255)
                                  gAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCAc                              >  1:3038275/1‑71 (MQ=255)
                                   caCCATTAGCAGAGATATAGGCATAAGTGAGGATATAGAGAACGAAGGCAACGGTAAta                                         >  1:2801671/2‑58 (MQ=255)
                                   caCCATTAGCAGAGATATAGGCATAAGTGAGGATATAGAGAACGAAGGCAACGGTAAta                                         >  1:578192/2‑58 (MQ=255)
                                   aacaacTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATgg                                              <  1:2188696/54‑1 (MQ=255)
                                    aCCATTAGCAGAGATATAGGCATAAGTGAGGATATAGAGAACgaa                                                      >  1:3037448/1‑45 (MQ=255)
                                    aCCATTAGCAGAGATATAGGCATAAGTGAGGATATAGAGAACGAAGGCAACGGTAAtac                                        >  1:2706079/1‑57 (MQ=255)
                                    aCCACTGGCAGAAATATAGGCATAGGTCAGtatata                                                               <  1:3513560/36‑1 (MQ=255)
                                    aCCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGt                            <  1:2963379/71‑1 (MQ=255)
                                    aCCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGt                            <  1:1730270/71‑1 (MQ=255)
                                     ccACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCa                                                  >  1:2968278/1‑48 (MQ=255)
                                     ccACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATg                                               >  1:3844157/1‑51 (MQ=255)
                                     ccACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGc                                        >  1:1487521/1‑58 (MQ=255)
                                      cACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACaaa                                                      <  1:2793215/43‑1 (MQ=255)
                                      cACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGtt                                    <  1:645678/61‑1 (MQ=255)
                                      cACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGtt                                    <  1:2286262/61‑1 (MQ=255)
                                         tGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAg                       >  1:1530889/1‑71 (MQ=255)
                                          ggCAGAAATATAGGCATAGGTCAGGATATAGAGCACaa                                                       <  1:2447768/38‑1 (MQ=255)
                                          ggCAGAAATATAGGCATAGGTCAGGATATAGAGCACaa                                                       <  1:3259477/38‑1 (MQ=255)
                                           gCAGAGATATAGGCATAAGTGAGGATATAGAGAACGAAGGCAACGGTAATACCGCTGATAATGTTCCAggt                     <  1:2258372/71‑3 (MQ=255)
                                            cAGAGATATAGGCATAAGTGAGGATATAGAGAACGAAGGCAACGGTAATACCGCTGataa                               <  1:2207693/60‑4 (MQ=255)
                                            cAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCt                    <  1:631903/71‑1 (MQ=255)
                                            cAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCt                    <  1:3900250/71‑1 (MQ=255)
                                             aGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGtt                                    <  1:315029/54‑1 (MQ=255)
                                             aGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTc                          <  1:1768982/64‑1 (MQ=255)
                                             aGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCtt                   <  1:3786149/71‑1 (MQ=255)
                                             aGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCt                    <  1:476426/70‑1 (MQ=255)
                                               aaaTATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTg                 <  1:2108386/71‑1 (MQ=255)
                                                  tataGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGcc               <  1:1300912/70‑1 (MQ=255)
                                                   ataGGCATAGGTCAGGATATAGAGCACAAAGGCAATgg                                              <  1:1976668/38‑1 (MQ=255)
                                                    taGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCa             <  1:205213/70‑1 (MQ=255)
                                                     aGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTg                                   <  1:159205/47‑1 (MQ=255)
                                                      ggCATACGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGtt                                    <  1:267439/45‑1 (MQ=255)
                                                       gCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCaa         <  1:2968351/71‑1 (MQ=255)
                                                         aTAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGc           <  1:3927449/67‑1 (MQ=255)
                                                         aTAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAAt       <  1:572894/71‑1 (MQ=255)
                                                          tAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATc      >  1:79700/1‑71 (MQ=255)
                                                           aGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCt     >  1:342043/1‑71 (MQ=255)
                                                           aGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCt     >  1:680087/1‑71 (MQ=255)
                                                            ggTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCtt    <  1:2775292/71‑1 (MQ=255)
                                                            ggTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCt     <  1:1548219/70‑1 (MQ=255)
                                                             gTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCttt   <  1:2810825/71‑1 (MQ=255)
                                                             gTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCttt   <  1:2512851/71‑1 (MQ=255)
                                                              tCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTg  <  1:2372374/71‑1 (MQ=255)
                                                               |                                                                     
TAGTGACATCTCTGCGAAGGTGTGATGCAGAATCGAACCACTGGCAGAAATATAGGCATAGGTCAGGATATAGAGCACAAAGGCAATGGAAATGCCGTTGACCACGTTCCAGCCTTTGCCCAGCAAATCTTTG  >  minE/2130536‑2130668

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: