Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,307,260 G→A 6.1% I112I (ATC→ATT yeeF ← predicted amino‑acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,307,2600GA 6.1% 242.3 / 3.8 81I112I (ATC→ATTyeeFpredicted amino‑acid transporter
Reads supporting (aligned to +/- strand):  ref base G (32/44);  new base A (3/2);  total (35/46)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.56e-01

CACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATG  >  minE/1307191‑1307329
                                                                     |                                                                     
caccacAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGa                                                                      <  1:4515339/71‑1 (MQ=255)
  ccacAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGatc                                                                    <  1:4933914/71‑1 (MQ=255)
  ccacAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGatc                                                                    <  1:2437463/71‑1 (MQ=255)
  ccacAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGatc                                                                    <  1:4672811/71‑1 (MQ=255)
  ccacAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGat                                                                     >  1:111523/1‑70 (MQ=255)
   cacaAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGatca                                                                   <  1:1716307/71‑1 (MQ=255)
     caAATATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGatcatc                                                                 >  1:3027411/1‑71 (MQ=255)
     caAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGatcatc                                                                 <  1:1724083/71‑1 (MQ=255)
     caAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGatcatc                                                                 <  1:2894500/71‑1 (MQ=255)
     caAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGatcatc                                                                 <  1:2850238/71‑1 (MQ=255)
     caAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGatcatc                                                                 >  1:2303092/1‑71 (MQ=255)
     caAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGa                                                                      >  1:5324891/1‑66 (MQ=255)
     caAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGa                                                                      >  1:886771/1‑66 (MQ=255)
        aCATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGgc                                                              >  1:2406952/1‑71 (MQ=255)
         cATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGa                                                                      <  1:973296/62‑1 (MQ=255)
         cATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGgcg                                                             >  1:5364665/1‑71 (MQ=255)
             cACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGgcgc                                                            >  1:1986030/1‑68 (MQ=255)
             cACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGgcgc                                                            >  1:4097693/1‑68 (MQ=255)
             cACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGaa                                                         >  1:2977472/1‑71 (MQ=255)
              aCGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGaa                                                         >  1:2693746/1‑70 (MQ=255)
               cGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGaa                                                         >  1:4482169/1‑69 (MQ=255)
                gatgggatggAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGaa                                                         <  1:2586918/68‑1 (MQ=255)
                 atgggatggAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGa                                                     <  1:5291204/71‑1 (MQ=255)
                  tgggatggAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGa                                                     >  1:1643169/1‑70 (MQ=255)
                  tggaatggaaGGCACCAGAGCTTCAAAATAGATTTTCGCCAGCAGAATGTTGATCATCGGCGCGAACa                                                       <  1:1818989/68‑1 (MQ=255)
                   gggatggAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGatcatc                                                                 <  1:2901160/57‑1 (MQ=255)
                    ggatggAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGat                                                                     >  1:2596483/1‑52 (MQ=255)
                       tggAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGatc                                                                    <  1:2964775/50‑1 (MQ=255)
                       tggAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGat                                                                     <  1:913706/49‑1 (MQ=255)
                            ggCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGa                                                     <  1:1346579/60‑1 (MQ=255)
                            ggCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGa                                                     <  1:1315512/60‑1 (MQ=255)
                            ggCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAg                                             >  1:4504357/1‑68 (MQ=255)
                             gCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCaga                                         <  1:4528372/71‑1 (MQ=255)
                             gCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCaga                                         <  1:3017427/71‑1 (MQ=255)
                              cACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCagag                                        <  1:1412215/71‑1 (MQ=255)
                                ccAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCagag                                        <  1:3000696/69‑1 (MQ=255)
                                  agagCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGgc                                                              >  1:1375809/1‑45 (MQ=255)
                                  agagCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGagaaga                                    >  1:3644884/1‑71 (MQ=255)
                                    agCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAg                                                      <  1:948535/51‑1 (MQ=255)
                                    agCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATa                                                   >  1:1827460/1‑54 (MQ=255)
                                      cTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGgcg                                                             <  1:62154/42‑1 (MQ=255)
                                       ttCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTc                                                <  1:4417825/54‑1 (MQ=255)
                                       ttCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTc                                                <  1:4890419/54‑1 (MQ=255)
                                        tCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGcgcg                                                           >  1:1330570/1‑42 (MQ=255)
                                           aaaTAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACccac                           <  1:4898365/71‑1 (MQ=255)
                                            gaTAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGAc                                   >  1:910253/2‑62 (MQ=255)
                                            aaTAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACccac                           >  1:3863930/1‑70 (MQ=255)
                                             aTAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATa                                                   >  1:1880476/1‑45 (MQ=255)
                                             aTAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAgcgat                                      <  1:2430133/58‑2 (MQ=255)
                                              tAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCagag                                        <  1:2871970/55‑1 (MQ=255)
                                              tAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCagag                                        <  1:3037931/55‑1 (MQ=255)
                                              tAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCagag                                        <  1:2235700/55‑1 (MQ=255)
                                               aGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCAt                        >  1:1470460/1‑70 (MQ=255)
                                               aGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCa                                 >  1:2110271/1‑61 (MQ=255)
                                               aGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCAt                        >  1:1544400/1‑70 (MQ=255)
                                                gATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATa                                                   >  1:996585/1‑42 (MQ=255)
                                                gATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATa                                                   >  1:1789257/1‑42 (MQ=255)
                                                gATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATa                       <  1:4106017/70‑1 (MQ=255)
                                                gATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATa                       <  1:5051139/70‑1 (MQ=255)
                                                 aTTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATaaa                     <  1:3428018/71‑1 (MQ=255)
                                                  ttttCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAg                    <  1:763865/71‑1 (MQ=255)
                                                   tttCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACccacc                          >  1:5242652/1‑64 (MQ=255)
                                                      cGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCg                 <  1:2248072/70‑1 (MQ=255)
                                                           gTAGAATGTTGATCATCGGCGCGCACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTc           >  1:105633/1‑71 (MQ=255)
                                                           gTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTc           >  1:667762/1‑71 (MQ=255)
                                                            tAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGagaa                                      <  1:4748720/43‑1 (MQ=255)
                                                            tAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGagaa                                      <  1:4259053/43‑1 (MQ=255)
                                                            tAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGagaa                                      <  1:1059404/43‑1 (MQ=255)
                                                            tAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGagaa                                      <  1:3742541/43‑1 (MQ=255)
                                                            tAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGagaa                                      <  1:3408543/43‑1 (MQ=255)
                                                            tAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCAACCATAAAGCCGACAGTCg          <  1:5383184/71‑1 (MQ=255)
                                                             aGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAACGcc                  >  1:4503115/1‑62 (MQ=255)
                                                             aGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCgg         >  1:2587812/1‑71 (MQ=255)
                                                                aTGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCt      <  1:2625261/71‑1 (MQ=255)
                                                                aTGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCt      <  1:2251210/71‑1 (MQ=255)
                                                                  gTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTaa    >  1:2875201/1‑71 (MQ=255)
                                                                   ttGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTaa    <  1:4386829/70‑1 (MQ=255)
                                                                   ttGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAAt   <  1:1330337/71‑1 (MQ=255)
                                                                   ttGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAAt   <  1:3691516/71‑1 (MQ=255)
                                                                    tGATCATCGGCGCGAACAGATAGTCGAGCAGagaa                                      >  1:1076955/1‑35 (MQ=255)
                                                                    tGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGt            >  1:4004448/1‑61 (MQ=255)
                                                                     gATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAAccc                             >  1:4920409/1‑43 (MQ=255)
                                                                     gATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATg  <  1:392577/70‑1 (MQ=255)
                                                                     gATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATg  <  1:1677791/70‑1 (MQ=255)
                                                                     gATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATg  >  1:1177148/1‑70 (MQ=255)
                                                                     |                                                                     
CACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATG  >  minE/1307191‑1307329

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: