Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 455,055 A→C 19.5% A10A (GCT→GCG ybgH ← predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE455,0550AC19.5% 115.6 / 38.4 72A10A (GCT→GCGybgHpredicted transporter
Reads supporting (aligned to +/- strand):  ref base A (30/28);  new base C (7/7);  total (37/35)
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.41e-01

AGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATAGT  >  minE/454990‑455121
                                                                 |                                                                  
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:742297/70‑1 (MQ=255)
    aggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGc                                                           <  1:2253888/70‑1 (MQ=255)
     ggCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCt                                                          <  1:1058675/71‑1 (MQ=255)
     ggCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCt                                                          <  1:1765940/71‑1 (MQ=255)
      gCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtg                                                         <  1:523346/71‑1 (MQ=255)
      gCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCt                                                          >  1:925881/1‑70 (MQ=255)
       cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                                                        <  1:3820351/71‑1 (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:2396610/71‑1 (MQ=255)
         cGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGa                                                      <  1:1320841/71‑1 (MQ=255)
          gCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGAt                                                     >  1:1978969/1‑71 (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:3627767/70‑1 (MQ=255)
               ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgt                                                 <  1:844890/70‑1 (MQ=255)
               ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgt                                                 <  1:970811/70‑1 (MQ=255)
               ccGTAAAAGCTAAAATATTCCCAGATTTGAAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgt                                                 <  1:2950957/70‑1 (MQ=255)
                  tAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGttt                                             <  1:1840787/71‑1 (MQ=255)
                  tAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGagg                                                    <  1:2582739/64‑3 (MQ=255)
                    aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGttttca                                          <  1:3945306/71‑4 (MQ=255)
                    aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGttttca                                          >  1:975180/1‑68 (MQ=255)
                     aaGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTAtt                                          <  1:1835985/71‑1 (MQ=255)
                      aGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                                                        <  1:3140237/56‑1 (MQ=255)
                         tAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                                                        <  1:3642785/53‑1 (MQ=255)
                         tAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTa                                            <  1:2425680/65‑1 (MQ=255)
                          aaaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCgg                                                            <  1:2976604/48‑1 (MQ=255)
                          aaaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCgg                                                            <  1:3761966/48‑1 (MQ=255)
                           aaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGa                                                      >  1:768312/1‑53 (MQ=255)
                           aaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGa                                                      >  1:3120921/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)
                                                gcAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAAc               >  1:1308397/1‑71 (MQ=255)
                                                  aaCGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGt             >  1:197753/1‑71 (MQ=255)
                                                  aaCGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGt             >  1:316874/1‑71 (MQ=255)
                                                  aaCGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGt             >  1:854809/1‑71 (MQ=255)
                                                    cGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAAtt                                >  1:2443820/1‑50 (MQ=255)
                                                    cGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTAATCATAATAATTCCACAATTGACGGCa                 >  1:2142190/1‑65 (MQ=255)
                                                     gACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCAc          >  1:149303/1‑71 (MQ=255)
                                                        aTAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGt        >  1:1375216/1‑70 (MQ=255)
                                                         tAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTaa      <  1:2576594/71‑1 (MQ=255)
                                                         tAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTaa      <  1:1664833/71‑1 (MQ=255)
                                                            tAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATa    <  1:3792915/70‑1 (MQ=255)
                                                            tAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATa    <  1:591364/70‑1 (MQ=255)
                                                            tAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATAg   >  1:2831583/1‑71 (MQ=255)
                                                            tAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATAg   <  1:855418/71‑1 (MQ=255)
                                                            tAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATAg   >  1:1353590/1‑71 (MQ=255)
                                                             aaaTAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATAGt  >  1:1821709/1‑71 (MQ=255)
                                                              aaTAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGa                      >  1:3001603/1‑50 (MQ=255)
                                                              aaTAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATAGt  >  1:2918522/1‑70 (MQ=255)
                                                                tAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTc                               <  1:1720520/39‑1 (MQ=255)
                                                                 aGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGt        >  1:35047/1‑61 (MQ=255)
                                                                 |                                                                  
AGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATAGT  >  minE/454990‑455121

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: