Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 455,036 G→T 30.3% L17M (CTG→ATG) ‡ ybgH ← predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE455,0360GT30.3% 195.8 / 69.9 112L17M (CTG→ATG) ‡ybgHpredicted transporter
Reads supporting (aligned to +/- strand):  ref base G (39/39);  new base T (11/23);  total (50/62)
Fisher's exact test for biased strand distribution p-value = 1.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 123 total aligned reads displayed.

GTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGC  >  minE/454967‑455105
                                                                     |                                                                     
 ttGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAgcg                                                                     >  1:2372739/1‑71 (MQ=255)
 ttGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAg                                                                       <  1:2972474/69‑1 (MQ=255)
 ttGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAg                                                                       <  1:1774744/69‑1 (MQ=255)
  tGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAgcg                                                                     >  1:1587338/1‑70 (MQ=255)
   gATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAgcgc                                                                    >  1:596092/1‑70 (MQ=255)
   gATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAgcgc                                                                    >  1:1619090/1‑70 (MQ=255)
     ttGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAAc                                                                 >  1:2081188/1‑71 (MQ=255)
     tgggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCAc                                                                 <  1:3862299/69‑1 (MQ=255)
     tgggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCAc                                                                 <  1:3353582/69‑1 (MQ=255)
         tGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACa                                                             <  1:1790594/71‑1 (MQ=255)
           agaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGatat                                                            <  1:3291639/70‑1 (MQ=255)
           agaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGatat                                                            <  1:3771237/70‑1 (MQ=255)
            gaTAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAg                                                          >  1:3082007/1‑71 (MQ=255)
            gaTAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAg                                                          <  1:265496/71‑1 (MQ=255)
             aTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGAtata                                                           <  1:2225637/69‑1 (MQ=255)
             aTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGAtata                                                           >  1:3039509/1‑69 (MQ=255)
              tAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGAtata                                                           <  1:5071689/68‑1 (MQ=255)
              tAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGt                                                         <  1:4606702/70‑1 (MQ=255)
              tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCTCAACGACATAGTa                                                        <  1:1553546/71‑1 (MQ=255)
              tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa                                                        <  1:5403146/71‑1 (MQ=255)
              tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTACCAGATTTGCAGCGCAACGACATAGTa                                                        <  1:4030605/71‑1 (MQ=255)
               agagAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTa                                                        >  1:136720/1‑70 (MQ=255)
               agagAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa                                                        <  1:4914839/70‑1 (MQ=255)
                gagaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGata                                                             >  1:2970162/1‑62 (MQ=255)
                gagaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGata                                                             >  1:3129678/1‑62 (MQ=255)
                gagaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGt                                                         >  1:1350989/1‑68 (MQ=255)
                gagaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGt                                                         >  1:4002479/1‑68 (MQ=255)
                gagaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTata                                                      <  1:700401/71‑3 (MQ=255)
                gagaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTat                                                       <  1:5468723/70‑2 (MQ=255)
                gagaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTat                                                       <  1:700044/70‑2 (MQ=255)
                gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                      <  1:4845549/71‑1 (MQ=255)
                gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                      <  1:2283969/71‑1 (MQ=255)
                gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                      <  1:2087903/71‑1 (MQ=255)
                gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                      <  1:3352644/71‑1 (MQ=255)
                gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                      >  1:328063/1‑71 (MQ=255)
                gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                      <  1:3796277/71‑1 (MQ=255)
                gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                      <  1:4334378/71‑1 (MQ=255)
                 agaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                      >  1:4696071/1‑70 (MQ=255)
                   aaTCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATa                                                           >  1:4607506/1‑63 (MQ=255)
                   aaTCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATa                                                           >  1:1265823/1‑63 (MQ=255)
                      cagcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgc                                                >  1:4318729/1‑71 (MQ=255)
                      cagcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgc                                                >  1:2903474/1‑71 (MQ=255)
                       agcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgcg                                               >  1:1768082/1‑71 (MQ=255)
                        gcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                      >  1:494033/1‑63 (MQ=255)
                         taaggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATcgcgcgcg                                             >  1:4027102/4‑71 (MQ=255)
                         cagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGcgcgcg                                             >  1:925877/1‑71 (MQ=255)
                          agGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACa                                                             <  1:283032/54‑1 (MQ=255)
                          agGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCgg                                            >  1:1929373/1‑71 (MQ=255)
                          aaggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCgg                                            >  1:1049272/3‑71 (MQ=255)
                          aaggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCgg                                            >  1:7948/3‑71 (MQ=255)
                             gCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtg                                         <  1:3150876/71‑1 (MQ=255)
                              cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                                        <  1:4503315/71‑1 (MQ=255)
                              cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtg                                         <  1:1933962/70‑1 (MQ=255)
                               aCGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAg                                                                       >  1:5562898/1‑39 (MQ=255)
                               aCGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAg                                                                       >  1:1352737/1‑39 (MQ=255)
                                cGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGa                                      <  1:70855/71‑1 (MQ=255)
                                cGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGa                                      <  1:294826/71‑1 (MQ=255)
                                 gCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGag                                     >  1:2119683/1‑70 (MQ=255)
                                  cATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAAc                                                                 >  1:3208739/1‑42 (MQ=255)
                                  cATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGTTGTGATg                                    >  1:3128724/1‑71 (MQ=255)
                                  cATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAg                                                          <  1:3572438/49‑1 (MQ=255)
                                  cATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGagg                                    <  1:121246/71‑3 (MQ=255)
                                    tGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa                                                        <  1:4810159/49‑1 (MQ=255)
                                    tGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTatat                                                     <  1:4975880/52‑1 (MQ=255)
                                    tGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGagg                                    <  1:1998365/69‑3 (MQ=255)
                                      ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCa                                                                   >  1:4173295/1‑36 (MQ=255)
                                      ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCa                                                                   >  1:4221607/1‑36 (MQ=255)
                                      ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCa                                                                   >  1:4484290/1‑36 (MQ=255)
                                      ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgtg                                >  1:442478/1‑71 (MQ=255)
                                         tAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATa                                                           <  1:4168569/41‑1 (MQ=255)
                                         tAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGt                                                         <  1:941405/43‑1 (MQ=255)
                                         tAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGt                                                         <  1:1091625/43‑1 (MQ=255)
                                         tAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                      <  1:4982778/46‑1 (MQ=255)
                                           aaaaCTGAAGTACTCCCAGATTTGTATCGCCACGATATAGTATATCGCGCGCGGCTGTGag                                     >  1:3834926/1‑60 (MQ=255)
                                           aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGttttc                           <  1:3652356/70‑3 (MQ=255)
                                           aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGttttc                           <  1:1798406/70‑3 (MQ=255)
                                           aaaGCTAAAATATTCCCAGATTTGCAGCGCAACGACAt                                                            >  1:776803/1‑38 (MQ=255)
                                           aaaGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTAt                           >  1:3473302/1‑71 (MQ=255)
                                              gCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGc                                           <  1:2532398/52‑1 (MQ=255)
                                                tAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                                        <  1:515501/53‑1 (MQ=255)
                                                 aaaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTa                            >  1:4483838/1‑64 (MQ=255)
                                                  aaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTa                            >  1:5443618/1‑63 (MQ=255)
                                                     tatTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGtt                              <  1:1532654/58‑1 (MQ=255)
                                                       ttCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTc               <  1:4033631/71‑1 (MQ=255)
                                                          ccAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCcac            <  1:5599699/71‑1 (MQ=255)
                                                          ccAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCca             >  1:3020669/1‑70 (MQ=255)
                                                          ccAGATTTGCAGAGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCca             >  1:1175136/1‑70 (MQ=255)
                                                           cAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCcaca           <  1:1510223/71‑1 (MQ=255)
                                                            aGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGa                                      >  1:3796198/1‑43 (MQ=255)
                                                            aGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGa                                      >  1:2925206/1‑43 (MQ=255)
                                                             gATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAAt         >  1:5553338/1‑71 (MQ=255)
                                                              aTTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAAtt        <  1:5153035/71‑1 (MQ=255)
                                                              aTTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAAtt        >  1:1924153/1‑71 (MQ=255)
                                                                ttGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGa      >  1:2460373/1‑71 (MQ=255)
                                                                ttGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGa      >  1:5031569/1‑71 (MQ=255)
                                                                ttGCAGAGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGa      >  1:737611/1‑71 (MQ=255)
                                                                  gCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGa      <  1:5020296/69‑1 (MQ=255)
                                                                   cAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACgg   >  1:3548632/1‑71 (MQ=255)
                                                                    aGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGc  <  1:1001850/71‑1 (MQ=255)
                                                                     gcgcAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGc  <  1:4602369/70‑1 (MQ=255)
                                                                     |                                                                     
GTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGC  >  minE/454967‑455105

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: