Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,307,253 G→A 5.9% L115L (CTA→TTA) ‡ yeeF ← predicted amino‑acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,307,2530GA 5.9% 254.2 / 3.7 84L115L (CTA→TTA) ‡yeeFpredicted amino‑acid transporter
Reads supporting (aligned to +/- strand):  ref base G (33/46);  new base A (3/2);  total (36/48)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

CCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGG  >  minE/1307184‑1307322
                                                                     |                                                                     
ccaGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAg                                                                       >  1:590542/1‑70 (MQ=255)
ccaGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAg                                                                       >  1:3653919/1‑70 (MQ=255)
 caGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGaa                                                                     <  1:2889049/71‑1 (MQ=255)
   gcgcCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATg                                                                   <  1:209712/71‑1 (MQ=255)
     gcCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGtt                                                                 <  1:3336601/71‑1 (MQ=255)
       caccacAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGa                                                               <  1:4515339/71‑1 (MQ=255)
         ccacAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGatc                                                             <  1:4933914/71‑1 (MQ=255)
         ccacAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGatc                                                             <  1:4672811/71‑1 (MQ=255)
         ccacAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGatc                                                             <  1:2437463/71‑1 (MQ=255)
         ccacAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGat                                                              >  1:111523/1‑70 (MQ=255)
          cacaAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGtt                                                                 <  1:707316/66‑1 (MQ=255)
          cacaAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGatca                                                            <  1:1716307/71‑1 (MQ=255)
            caAATATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGatcatc                                                          >  1:3027411/1‑71 (MQ=255)
            caAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGatcatc                                                          <  1:2850238/71‑1 (MQ=255)
            caAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGatcatc                                                          <  1:1724083/71‑1 (MQ=255)
            caAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGatcatc                                                          <  1:2894500/71‑1 (MQ=255)
            caAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGatcatc                                                          >  1:2303092/1‑71 (MQ=255)
            caAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGa                                                               >  1:886771/1‑66 (MQ=255)
            caAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGa                                                               >  1:5324891/1‑66 (MQ=255)
              aaCATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATg                                                                   <  1:1467974/60‑1 (MQ=255)
              aaCATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATg                                                                   <  1:4721457/60‑1 (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:4097693/1‑68 (MQ=255)
                    cACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGgcgc                                                     >  1:1986030/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)
                        atgggatggAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGaa                                                                     >  1:767571/1‑48 (MQ=255)
                        atgggatggAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAAt                                                                    >  1:2036669/1‑49 (MQ=255)
                        atgggatggAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAAt                                                                    >  1:5514547/1‑49 (MQ=255)
                        atgggatggAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGa                                              <  1:5291204/71‑1 (MQ=255)
                        atgggatggAAGGCAACAGAGCTTCAAAATAGATTTTCGCCAGTAGAAt                                                                    >  1:3868982/1‑49 (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)
                                  aGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAAt                                                                    <  1:666152/39‑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:2235700/55‑1 (MQ=255)
                                                     tAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCagag                                 <  1:3037931/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:3742541/43‑1 (MQ=255)
                                                                   tAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGagaa                               <  1:4259053/43‑1 (MQ=255)
                                                                   tAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGagaa                               <  1:3408543/43‑1 (MQ=255)
                                                                   tAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGagaa                               <  1:4748720/43‑1 (MQ=255)
                                                                   tAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGagaa                               <  1:1059404/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)
                                                                     |                                                                     
CCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGG  >  minE/1307184‑1307322

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: