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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE455,0360GT30.0% 148.9 / 49.1 83L17M (CTG→ATG) ‡ybgHpredicted transporter
Reads supporting (aligned to +/- strand):  ref base G (22/36);  new base T (9/16);  total (31/52)
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.28e-01

ATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGG  >  minE/454971‑455104
                                                                 |                                                                    
atgggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGAtcgc                                                                   <  1:2110966/66‑1 (MQ=255)
atgggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGAtcgc                                                                   <  1:594407/66‑1 (MQ=255)
 ttGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCaa                                                                 >  1:3571166/1‑70 (MQ=255)
  tGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCaa                                                                 >  1:2071687/1‑69 (MQ=255)
   ggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACg                                                               <  1:3413146/70‑1 (MQ=255)
   ggTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACg                                                               <  1:1588269/70‑1 (MQ=255)
     tGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGata                                                            <  1:479848/71‑3 (MQ=255)
      gagaTAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACa                                                            >  1:2183757/1‑70 (MQ=255)
       agaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGAtata                                                          <  1:1031011/71‑1 (MQ=255)
         aTAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGt                                                        <  1:2289002/71‑1 (MQ=255)
          tAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTa                                                       >  1:1461850/1‑71 (MQ=255)
          tAGAGAATCAGCAGGACACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa                                                       <  1:220749/71‑1 (MQ=255)
           agagAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAgcgc                                                                   <  1:2383786/58‑1 (MQ=255)
            gagaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTata                                                     <  1:2013045/71‑3 (MQ=255)
            gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                     <  1:1765787/71‑1 (MQ=255)
            gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                     <  1:3595261/71‑1 (MQ=255)
            gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                     <  1:1366534/71‑1 (MQ=255)
             agaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGAtata                                                          <  1:3048127/65‑1 (MQ=255)
              gaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTatatc                                                   >  1:665482/1‑70 (MQ=255)
              gaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTatatc                                                   <  1:3408650/71‑2 (MQ=255)
              gaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAAt                                                    <  1:1390922/70‑1 (MQ=255)
                   agcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAg                                                                      >  1:2611027/1‑47 (MQ=255)
                   agcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgcg                                              <  1:954568/71‑1 (MQ=255)
                      agGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCgg                                           >  1:574366/1‑71 (MQ=255)
                       gggCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGcgcgcg                                            >  1:3215920/1‑69 (MQ=255)
                       gggCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGcgcgcg                                            >  1:1412661/1‑69 (MQ=255)
                       aggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGc                                          <  1:2253888/70‑1 (MQ=255)
                       aggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGc                                          <  1:742297/70‑1 (MQ=255)
                        ggCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCt                                         <  1:1058675/71‑1 (MQ=255)
                        ggCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCt                                         <  1:1765940/71‑1 (MQ=255)
                         gCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACg                                                               <  1:3558903/48‑1 (MQ=255)
                         gCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtg                                        <  1:523346/71‑1 (MQ=255)
                         gCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCt                                         >  1:925881/1‑70 (MQ=255)
                          cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAgc                                                                     >  1:2350882/1‑41 (MQ=255)
                          cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAg                                                                      >  1:198889/1‑40 (MQ=255)
                          cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                                       <  1:3820351/71‑1 (MQ=255)
                           aCGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGAc                                                             >  1:2373244/1‑48 (MQ=255)
                           aCGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCtgt                                       >  1:2978068/1‑70 (MQ=255)
                            cGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgc                                                 >  1:1591018/1‑59 (MQ=255)
                            cGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCtgt                                       >  1:2146594/1‑69 (MQ=255)
                            cGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGa                                     <  1:1320841/71‑1 (MQ=255)
                            cGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGa                                     <  1:2396610/71‑1 (MQ=255)
                             gCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGAt                                    >  1:1978969/1‑71 (MQ=255)
                             gCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTa                                                       <  1:2916924/52‑1 (MQ=255)
                             gCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTa                                                       <  1:1194229/52‑1 (MQ=255)
                             gCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGag                                    >  1:1627916/1‑70 (MQ=255)
                             gCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGag                                    >  1:3258357/1‑70 (MQ=255)
                             gCATGCAGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGag                                    >  1:513700/1‑70 (MQ=255)
                               aTGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGc                                  >  1:954552/1‑71 (MQ=255)
                                tGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGag                                    <  1:3130043/68‑2 (MQ=255)
                                  ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgtg                               >  1:3201955/1‑71 (MQ=255)
                                  ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgt                                <  1:970811/70‑1 (MQ=255)
                                  ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgt                                <  1:844890/70‑1 (MQ=255)
                                  ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgt                                <  1:3627767/70‑1 (MQ=255)
                                  ccGTAAAAGCTAAAATATTCCCAGATTTGAAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgt                                <  1:2950957/70‑1 (MQ=255)
                                   cGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAAt                                                    <  1:428715/49‑1 (MQ=255)
                                    gTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGatat                                                           <  1:1236708/41‑1 (MQ=255)
                                     tAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACa                                                            <  1:3026772/39‑1 (MQ=255)
                                     tAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACa                                                            <  1:3207931/39‑1 (MQ=255)
                                     tAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGttt                            <  1:1840787/71‑1 (MQ=255)
                                     tAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGagg                                   <  1:2582739/64‑3 (MQ=255)
                                       aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGA‑GGGTGttttca                         >  1:975180/1‑68 (MQ=255)
                                       aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGA‑GGGTGttttca                         <  1:3945306/71‑4 (MQ=255)
                                        aaGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTAtt                         <  1:1835985/71‑1 (MQ=255)
                                         aGCTAAAATATTCCCAGATTTGCAGCGCAACGACATa                                                          >  1:1545148/1‑37 (MQ=255)
                                         aGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                                       <  1:3140237/56‑1 (MQ=255)
                                            tAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                                       <  1:3642785/53‑1 (MQ=255)
                                            tAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTa                           <  1:2425680/65‑1 (MQ=255)
                                             aaaaTATTCCCAGATTTGCAGCGCAACGACATAGTa                                                       <  1:2826580/36‑1 (MQ=255)
                                             aaaaTATTCCCAGATTTGCAGCGCAACGACATAGTa                                                       <  1:834973/36‑1 (MQ=255)
                                             aaaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCgg                                           <  1:3761966/48‑1 (MQ=255)
                                             aaaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCgg                                           <  1:2976604/48‑1 (MQ=255)
                                              aaaTATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                     <  1:3947282/37‑1 (MQ=255)
                                              aaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGa                                     >  1:3120921/1‑53 (MQ=255)
                                              aaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGa                                     >  1:768312/1‑53 (MQ=255)
                                                atatTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTa                           <  1:2360379/61‑1 (MQ=255)
                                                atatTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCAtaataa                 <  1:1649166/71‑1 (MQ=255)
                                                 tatTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                                       >  1:1749513/1‑48 (MQ=255)
                                                  atTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAAtt               >  1:1736829/1‑71 (MQ=255)
                                                   ttCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTc              >  1:3255694/1‑71 (MQ=255)
                                                    tCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTc              >  1:2045534/1‑70 (MQ=255)
                                                      ccAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCcac           <  1:3367121/71‑1 (MQ=255)
                                                      ccAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCca            <  1:2367512/70‑1 (MQ=255)
                                                        aGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCcaca          >  1:3313784/1‑70 (MQ=255)
                                                        aGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACaa         <  1:775624/71‑1 (MQ=255)
                                                         gATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGATTATTCATAATAATTCCACAAt        >  1:2068334/1‑71 (MQ=255)
                                                              gCAGCGCATCGACATAGTAAATAGCGCTCGGCTGTGATGCGTGTTTa                           <  1:3225894/47‑1 (MQ=39)
                                                               cAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCcaca          >  1:2811476/1‑63 (MQ=255)
                                                                aGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACgg  >  1:2086334/1‑70 (MQ=255)
                                                                 |                                                                    
ATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGG  >  minE/454971‑455104

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: ‑