Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 454,989 C→G 10.4% L32L (CTG→CTC ybgH ← predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE454,9890CG10.4% 221.2 / 13.2 77L32L (CTG→CTCybgHpredicted transporter
Reads supporting (aligned to +/- strand):  ref base C (33/36);  new base G (4/4);  total (37/40)
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.27e-01

CAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCG  >  minE/454921‑455058
                                                                    |                                                                     
cAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTCGTTGATTGGTGAGATAGAGAATcag                                                                     <  1:2849055/71‑1 (MQ=255)
cAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcag                                                                     >  1:2636942/1‑71 (MQ=255)
cAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcag                                                                     <  1:2689441/71‑1 (MQ=255)
cAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcag                                                                     >  1:1569786/1‑71 (MQ=255)
cAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcag                                                                     <  1:3778093/71‑1 (MQ=255)
     ggCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATca                                                                      >  1:811000/1‑65 (MQ=255)
       cgcTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcagca                                                                   >  1:2224203/1‑66 (MQ=255)
         cTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcag                                                                     <  1:982885/62‑1 (MQ=255)
         cTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcag                                                                     <  1:2444174/62‑1 (MQ=255)
         cTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACg                                                            <  1:2228850/71‑1 (MQ=255)
         cTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACg                                                            <  1:61552/71‑1 (MQ=255)
          tAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcag                                                                     <  1:3327596/61‑1 (MQ=255)
              tAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcag                                                                     <  1:4038403/57‑1 (MQ=255)
              tAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGc                                                       <  1:1913923/71‑1 (MQ=255)
              tAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGc                                                       <  1:4097912/71‑1 (MQ=255)
               aaCTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGcc                                                      >  1:3240464/1‑71 (MQ=255)
                  tCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcag                                                                     <  1:2466773/53‑1 (MQ=255)
                  tCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCg                                                     >  1:3214852/1‑69 (MQ=255)
                  tCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGt                                                    >  1:3580465/1‑70 (MQ=255)
                   cGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTa                                                   >  1:2942984/1‑70 (MQ=255)
                    gTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcagca                                                                   >  1:640421/1‑53 (MQ=255)
                    gTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTaa                                                  <  1:2377919/70‑1 (MQ=255)
                    gTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTa                                                   <  1:2656889/69‑1 (MQ=255)
                    gTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTa                                                   <  1:756171/69‑1 (MQ=255)
                    gTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTa                                                   <  1:136861/69‑1 (MQ=255)
                       ggCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGc                                              >  1:3212056/1‑71 (MQ=255)
                          gTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcagca                                                                   >  1:3908257/1‑47 (MQ=255)
                          gTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcagca                                                                   >  1:1883470/1‑47 (MQ=255)
                           tGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcag                                                                     <  1:3063919/44‑1 (MQ=255)
                           tGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcag                                                                     <  1:1521846/44‑1 (MQ=255)
                             aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcagca                                                                   >  1:3398756/1‑44 (MQ=255)
                             aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATc                                                                       >  1:1175394/1‑40 (MQ=255)
                             aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACg                                                            <  1:1336543/51‑1 (MQ=255)
                             aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACg                                                            <  1:3381634/51‑1 (MQ=255)
                             aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACg                                                            <  1:4046149/51‑1 (MQ=255)
                             aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa                                         >  1:3014661/1‑70 (MQ=255)
                             aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa                                         <  1:700930/70‑1 (MQ=255)
                             aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa                                         <  1:1583992/70‑1 (MQ=255)
                             aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa                                         >  1:272347/1‑70 (MQ=255)
                             aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAat                                        <  1:704828/71‑1 (MQ=255)
                              gTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAata                                       >  1:1330484/1‑71 (MQ=255)
                              gTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAata                                       >  1:4131068/1‑71 (MQ=255)
                               tATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACg                                                            <  1:1607978/49‑1 (MQ=255)
                                 tCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcagca                                                                   >  1:1014040/1‑40 (MQ=255)
                                 tCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATc                                                                       >  1:671265/1‑36 (MQ=255)
                                 tCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACg                                                            <  1:1134709/47‑1 (MQ=255)
                                 tCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTc                                    >  1:1368659/1‑71 (MQ=255)
                                   gTTGTATTTTAGTTGATTGGTGAGATAGAGAATcagca                                                                   >  1:3745011/1‑38 (MQ=255)
                                   gTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTcc                                   >  1:2573458/1‑70 (MQ=255)
                                    ttGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACg                                                            <  1:55879/44‑1 (MQ=255)
                                       tATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGc                                                               <  1:3343843/38‑1 (MQ=255)
                                        aTTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAata                                       >  1:2968897/1‑61 (MQ=255)
                                        aTTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAata                                       >  1:3959543/1‑61 (MQ=255)
                                         ttttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaa                                          >  1:2994322/1‑57 (MQ=255)
                                          tttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTg                           >  1:3699829/1‑71 (MQ=255)
                                           ttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCa                                                          >  1:1472529/1‑39 (MQ=255)
                                           ttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGc                                              >  1:1468493/1‑51 (MQ=255)
                                           ttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTg                           <  1:3326154/70‑1 (MQ=255)
                                           ttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTg                           >  1:1661899/1‑70 (MQ=255)
                                           ttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTg                           <  1:4121510/70‑1 (MQ=255)
                                            tAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTg                           >  1:1307684/1‑69 (MQ=255)
                                              gTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTATAAGCTAAAata                                       <  1:3876936/55‑1 (MQ=255)
                                                 gATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAata                                       >  1:960131/1‑52 (MQ=255)
                                                 gATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAgcgc                     >  1:2511377/1‑70 (MQ=255)
                                                   ttGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTa                                                   <  1:2313457/38‑1 (MQ=255)
                                                   ttGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAAc                  <  1:1794726/71‑1 (MQ=255)
                                                     ggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTaa                                                  >  1:3891473/1‑37 (MQ=37)
                                                      gTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGa                <  1:1915790/70‑1 (MQ=255)
                                                       tGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGa                               <  1:2368093/54‑1 (MQ=255)
                                                       tGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGAc               >  1:1052665/1‑70 (MQ=255)
                                                       tGAGATAGAGAATCAGAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGa                               <  1:206121/53‑1 (MQ=255)
                                                         agaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGata              >  1:3795297/1‑67 (MQ=255)
                                                         agaTAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCATCGACa              <  1:819345/69‑1 (MQ=255)
                                                          gaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGAtata            <  1:3559243/70‑1 (MQ=255)
                                                           aTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACtg                                            <  1:775266/37‑2 (MQ=19)
                                                           aTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACtg                                            <  1:1679394/37‑2 (MQ=19)
                                                           aTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACtg                                            <  1:3218011/37‑2 (MQ=19)
                                                           aTAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAg           >  1:1310295/1‑70 (MQ=255)
                                                           aTAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGt          <  1:1541908/71‑1 (MQ=255)
                                                            tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTcc                                   <  1:2200551/45‑1 (MQ=255)
                                                            tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa         <  1:3613702/71‑1 (MQ=255)
                                                            tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa         <  1:9308/71‑1 (MQ=255)
                                                            tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa         <  1:83469/71‑1 (MQ=255)
                                                             agagAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGAtcgcca                   >  1:352824/1‑58 (MQ=255)
                                                             agagAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAat                                        >  1:2821213/1‑39 (MQ=255)
                                                             agagAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAgc                       >  1:1399834/1‑56 (MQ=255)
                                                             agagAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa         >  1:2988942/1‑70 (MQ=255)
                                                              gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGc                          >  1:3460408/1‑52 (MQ=255)
                                                               agaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa       >  1:581079/1‑70 (MQ=255)
                                                               agaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa       >  1:708108/1‑70 (MQ=255)
                                                                gaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTatatc     >  1:3956993/1‑70 (MQ=255)
                                                                gaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACAt             >  1:822525/1‑63 (MQ=255)
                                                                 aaTCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGAc               >  1:1784888/1‑60 (MQ=255)
                                                                  aTCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAg    <  1:1540939/70‑1 (MQ=255)
                                                                   tCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTc                                    <  1:315200/37‑1 (MQ=255)
                                                                   tCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTc                                    <  1:3549167/37‑1 (MQ=255)
                                                                   tCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTGCCAGATTTGCAGCGCAACGACATAGTAAATAgcg  >  1:1085761/1‑71 (MQ=255)
                                                                   tCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAAt      <  1:3237565/67‑1 (MQ=255)
                                                                   tCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcg  >  1:72838/1‑71 (MQ=255)
                                                                    |                                                                     
CAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCG  >  minE/454921‑455058

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: