Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,307,293 A→T 9.5% S101S (TCT→TCA yeeF ← predicted amino‑acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,307,2930AT 9.5% 169.7 / 9.2 63S101S (TCT→TCAyeeFpredicted amino‑acid transporter
Reads supporting (aligned to +/- strand):  ref base A (28/29);  new base T (3/3);  total (31/32)
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.48e-01

AGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGT  >  minE/1307225‑1307348
                                                                    |                                                       
agagCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGagaaga                                                       >  1:3644884/1‑71 (MQ=255)
         aaaTAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACccac                                              <  1:4898365/71‑1 (MQ=255)
          gaTAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGAc                                                      >  1:910253/2‑62 (MQ=255)
          aaTAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACccac                                              >  1:3863930/1‑70 (MQ=255)
             aGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCAt                                           >  1:1470460/1‑70 (MQ=255)
             aGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCa                                                    >  1:2110271/1‑61 (MQ=255)
             aGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCAt                                           >  1:1544400/1‑70 (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)
                              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:3691516/71‑1 (MQ=255)
                                 ttGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAAt                      <  1:1330337/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:1677791/70‑1 (MQ=255)
                                   gATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATg                     >  1:1177148/1‑70 (MQ=255)
                                   gATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATg                     <  1:392577/70‑1 (MQ=255)
                                      catcGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGa                                   >  1:4463388/1‑53 (MQ=255)
                                      catcGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGa                                   >  1:4278524/1‑53 (MQ=255)
                                      catcGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGa                                   >  1:4126299/1‑53 (MQ=255)
                                      catcGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACa                                 <  1:4916529/55‑1 (MQ=255)
                                        tcGGCGCGAACAGATAGTCGAGCAGAGAAGCCCAACCCCCCATAAAGCCGACAGTCGGGCTaa                       >  1:4688343/1‑63 (MQ=255)
                                        tcGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCTACCATAAAGCCGACAGt                               >  1:1236954/1‑55 (MQ=255)
                                        tcGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCgg                            >  1:949094/1‑58 (MQ=255)
                                        tcGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCgg                            >  1:2563448/1‑58 (MQ=255)
                                        tcGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTc               <  1:1068143/71‑1 (MQ=255)
                                        tcGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTc               >  1:3012123/1‑71 (MQ=255)
                                          ggCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCCCCATaaa                                        >  1:4536053/1‑44 (MQ=255)
                                          ggCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCt              <  1:2195351/70‑1 (MQ=255)
                                          ggCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCt              <  1:3162498/70‑1 (MQ=255)
                                          ggCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTg             >  1:4553048/1‑71 (MQ=255)
                                           gcgcGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCg                                    <  1:1780573/47‑1 (MQ=255)
                                                 aCAGATAGTCGAGCAGATAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGt        >  1:3487932/1‑69 (MQ=255)
                                                 aCAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGt        >  1:4094206/1‑69 (MQ=255)
                                                 aCAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGt        >  1:5150024/1‑69 (MQ=255)
                                                 aCAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGt        >  1:3806044/1‑69 (MQ=255)
                                                 aCAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTaa      <  1:1115919/71‑1 (MQ=255)
                                                  cAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTc                              >  1:1711269/1‑46 (MQ=255)
                                                    gATAGTCGAGCAGAGATGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAgtg   <  1:3934290/71‑1 (MQ=255)
                                                    gATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTa                        <  1:1989448/50‑1 (MQ=255)
                                                    gATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGt        <  1:560737/66‑1 (MQ=255)
                                                    gATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAgtg   <  1:913097/71‑1 (MQ=255)
                                                    gATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAgtg   >  1:748109/1‑71 (MQ=255)
                                                    gATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGAGTAAgt    <  1:5576616/70‑1 (MQ=255)
                                                     aTAGTCGAGCAGAGAAGACCAACGCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAgtgt  <  1:4062448/71‑1 (MQ=255)
                                                     aTAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTc                              <  1:5528301/43‑1 (MQ=255)
                                                     aTAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTc                              <  1:4057712/43‑1 (MQ=255)
                                                     aTAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTc                              <  1:2875575/43‑1 (MQ=255)
                                                     aTAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAgtg   >  1:382116/1‑70 (MQ=255)
                                                        gTCGAGCAGAGATGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAgtg   <  1:4155727/67‑1 (MQ=255)
                                                                gagaAGACCAACCCACCATAAAGCCGACAGTCGGGCt                         <  1:1795077/37‑1 (MQ=255)
                                                                gagaAGACCAACCCACCATAAAGCCGACAGTCGGGCt                         <  1:3771317/37‑1 (MQ=255)
                                                                 agaagaCCAACCCACCATAAAGCCGACAGTCGGGCt                         <  1:4256360/36‑1 (MQ=255)
                                                                 agaagaCCAACCCACCATAAAGCCGACAGTCGGGCt                         <  1:2773797/36‑1 (MQ=255)
                                                                 agaagaCCAACCCACCATAAAGCCGACAGTCGGGCTAATGGAttt                <  1:3824735/45‑1 (MQ=255)
                                                                    |                                                       
AGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGT  >  minE/1307225‑1307348

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: