Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,307,251 T→C 7.3% L115L (CTA→CTG) ‡ yeeF ← predicted amino‑acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,307,2510TC 7.3% 233.7 / 5.9 82L115L (CTA→CTG) ‡yeeFpredicted amino‑acid transporter
Reads supporting (aligned to +/- strand):  ref base T (36/40);  new base C (3/3);  total (39/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCG  >  minE/1307182‑1307321
                                                                     |                                                                      
caccaGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTa                                                                       <  1:1247791/71‑1 (MQ=255)
  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:2437463/71‑1 (MQ=255)
           ccacAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGatc                                                            <  1:4672811/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)
              caAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTa                                                                       >  1:2517701/1‑57 (MQ=255)
              caAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTa                                                                       >  1:2510987/1‑57 (MQ=255)
              caAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTa                                                                       >  1:5696434/1‑57 (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)
                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:5514547/1‑49 (MQ=255)
                          atgggatggAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAAt                                                                   >  1:2036669/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:3017427/71‑1 (MQ=255)
                                      gCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCaga                                 <  1:4528372/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:2235700/55‑1 (MQ=255)
                                                       tAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCagag                                <  1:3037931/55‑1 (MQ=255)
                                                       tAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCagag                                <  1:2871970/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:3742541/43‑1 (MQ=255)
                                                                     tAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGagaa                              <  1:3408543/43‑1 (MQ=255)
                                                                     tAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGagaa                              <  1:1059404/43‑1 (MQ=255)
                                                                     tAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCAACCATAAAGCCGACAGTCg  <  1:5383184/71‑1 (MQ=255)
                                                                     |                                                                      
CACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCG  >  minE/1307182‑1307321

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: