Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,307,281 G→A 6.4% D105D (GAC→GAT yeeF ← predicted amino‑acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,307,2810GA 6.4% 233.2 / 3.9 78D105D (GAC→GATyeeFpredicted amino‑acid transporter
Reads supporting (aligned to +/- strand):  ref base G (31/42);  new base A (3/2);  total (34/44)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

GGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGT  >  minE/1307219‑1307348
                                                              |                                                                   
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)
      agagCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGagaaga                                                       >  1:3644884/1‑71 (MQ=255)
           ttCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTc                                                                   <  1:4417825/54‑1 (MQ=255)
           ttCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTc                                                                   <  1:4890419/54‑1 (MQ=255)
               aaaTAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACccac                                              <  1:4898365/71‑1 (MQ=255)
                gaTAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGAc                                                      >  1:910253/2‑62 (MQ=255)
                aaTAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACccac                                              >  1:3863930/1‑70 (MQ=255)
                 aTAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAgcgat                                                         <  1:2430133/58‑2 (MQ=255)
                  tAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCagag                                                           <  1:3037931/55‑1 (MQ=255)
                  tAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCagag                                                           <  1:2871970/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)
                    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:1059404/43‑1 (MQ=255)
                                tAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGagaa                                                         <  1:3408543/43‑1 (MQ=255)
                                tAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGagaa                                                         <  1:4748720/43‑1 (MQ=255)
                                tAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGagaa                                                         <  1:3742541/43‑1 (MQ=255)
                                tAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGagaa                                                         <  1:4259053/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:392577/70‑1 (MQ=255)
                                         gATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATg                     <  1:1677791/70‑1 (MQ=255)
                                         gATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATg                     >  1:1177148/1‑70 (MQ=255)
                                            catcGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGa                                   >  1:4463388/1‑53 (MQ=255)
                                            catcGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGa                                   >  1:4126299/1‑53 (MQ=255)
                                            catcGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGa                                   >  1:4278524/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:3162498/70‑1 (MQ=255)
                                                ggCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCt              <  1:2195351/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:3806044/1‑69 (MQ=255)
                                                       aCAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGt        >  1:4094206/1‑69 (MQ=255)
                                                       aCAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGt        >  1:5150024/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:748109/1‑71 (MQ=255)
                                                          gATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAgtg   <  1:913097/71‑1 (MQ=255)
                                                          gATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGAGTAAgt    <  1:5576616/70‑1 (MQ=255)
                                                           aTAGTCGAGCAGAGAAGACCAACGCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAgtgt  <  1:4062448/71‑1 (MQ=255)
                                                           aTAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTc                              <  1:2875575/43‑1 (MQ=255)
                                                           aTAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTc                              <  1:5528301/43‑1 (MQ=255)
                                                           aTAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTc                              <  1:4057712/43‑1 (MQ=255)
                                                           aTAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAgtg   >  1:382116/1‑70 (MQ=255)
                                                              gTCGAGCAGAGATGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAgtg   <  1:4155727/67‑1 (MQ=255)
                                                              |                                                                   
GGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGT  >  minE/1307219‑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: