Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,645,878 C→A 29.2% L49L (CTG→CTT yhgN ← predicted antibiotic transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,645,8780CA29.2% 88.8 / 35.4 65L49L (CTG→CTTyhgNpredicted antibiotic transporter
Reads supporting (aligned to +/- strand):  ref base C (3/43);  new base A (3/16);  total (6/59)
Fisher's exact test for biased strand distribution p-value = 3.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.28e-01

CTAAGAAATGCCAGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATGATTGCCCGCCGTCTTTTCGGTTCAGTATGTTTCAGT  >  minE/2645828‑2645944
                                                  |                                                                  
cTAAGAAATGCCAGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAAc                                                  <  1:2870715/69‑1 (MQ=255)
  aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                                                                  <  1:1803429/51‑1 (MQ=38)
  aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                                                                  <  1:2255677/51‑1 (MQ=38)
  aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                                                                  <  1:525392/51‑1 (MQ=38)
  aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                                                                  <  1:2151246/51‑1 (MQ=38)
  aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                                                                  <  1:2417073/51‑1 (MQ=38)
  aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                                                                  <  1:521450/51‑1 (MQ=38)
  aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                                                                  <  1:1003002/51‑1 (MQ=38)
  aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                                                                  <  1:867879/51‑1 (MQ=38)
  aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                                                                  <  1:375713/51‑1 (MQ=38)
  aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                                                                  <  1:1763779/51‑1 (MQ=38)
  aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                                                                  <  1:1560081/51‑1 (MQ=38)
  aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                                                                  <  1:369667/51‑1 (MQ=38)
  aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                                                                  <  1:347122/51‑1 (MQ=38)
  aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                                                                  <  1:3023577/51‑1 (MQ=38)
  aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                                                                  <  1:1117783/51‑1 (MQ=38)
  aaGAAATGCCAGAATTTTCTCGCCCGCAAACAGGAACACCAGCATCACAAg                                                                  <  1:3059658/51‑1 (MQ=38)
        tGCCAGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCTCAc                                          <  1:1588184/69‑1 (MQ=255)
        tGCCAGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCAc                                          >  1:1053818/1‑69 (MQ=255)
        tGCCAGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCAc                                          <  1:2722648/69‑1 (MQ=255)
        tGCCAGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCAc                                          <  1:2536556/69‑1 (MQ=255)
        tGCCAGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCAc                                          <  1:2329313/69‑1 (MQ=255)
        tGCCAGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCAc                                          >  1:231812/1‑69 (MQ=255)
        tGCCAGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCAc                                          <  1:3398187/69‑1 (MQ=255)
        tGCCAGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCAc                                          <  1:2121346/69‑1 (MQ=255)
        tGCCAGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCAc                                          <  1:2125378/69‑1 (MQ=255)
        tGCCAGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCAc                                          <  1:973679/69‑1 (MQ=255)
            ggAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg                                      <  1:2058965/68‑1 (MQ=255)
            aGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg                                      <  1:3323271/69‑1 (MQ=255)
            aGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg                                      <  1:3101180/69‑1 (MQ=255)
            aGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg                                      <  1:3406497/69‑1 (MQ=255)
            aGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg                                      <  1:1214504/69‑1 (MQ=255)
            aGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg                                      <  1:2964101/69‑1 (MQ=255)
            aGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg                                      <  1:872645/69‑1 (MQ=255)
            aGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg                                      <  1:1632783/69‑1 (MQ=255)
            aGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg                                      <  1:1764084/69‑1 (MQ=255)
            aGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg                                      <  1:868751/69‑1 (MQ=255)
            aGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg                                      <  1:1897031/69‑1 (MQ=255)
            aGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg                                      <  1:2648695/69‑1 (MQ=255)
            aGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg                                      <  1:2627273/69‑1 (MQ=255)
            aGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg                                      <  1:623404/69‑1 (MQ=255)
            aGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg                                      <  1:606147/69‑1 (MQ=255)
            aGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg                                      <  1:244963/69‑1 (MQ=255)
             gAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATg                                      <  1:3195410/68‑1 (MQ=255)
                                     acagCAGCATCACAAGGAGAGCAATAAGCAACTCTCGCACCATGATTGCCCGCCGTCTTTTCGGTTCAg             >  1:976148/1‑69 (MQ=255)
                                     acacCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATGATTGCCCGCCGTCTTTTCGGTTCAg             <  1:2714924/69‑1 (MQ=255)
                                     acacCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATGATTGCCCGCCGTCTTTTCGGTTCAg             <  1:106236/69‑1 (MQ=255)
                                     acacCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATGATTGCCCGCCGTCTTTTCGGTTCAg             <  1:1300576/69‑1 (MQ=255)
                                     acacCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATGATTGCCCGCCGTCTTTTCGGTTCAg             <  1:2984691/69‑1 (MQ=255)
                                     acacCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATGATTGCCCGCCGTCTTTTCGGTTCAg             <  1:2618057/69‑1 (MQ=255)
                                     acacCAGCATCACAAGGAGAGCAATAAGCAACTCTCGCACCATGATTGCCCGCCGTCTTTTCGGTTCAg             >  1:2207607/1‑69 (MQ=255)
                                     acacCAGCATCACAAGGAGAGCAATAAGCAACTCTCGCACCATGATTGCCCGCCGTCTTTTCGGTTCAg             >  1:2015742/1‑69 (MQ=255)
                                       acCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATGATTGCCCGCCGTCTTTTCGGTTCAGTa           <  1:3350615/69‑1 (MQ=255)
                                       acCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATGATTGCCCGCCGTCTTTTCGGTTCAGTa           <  1:1242653/69‑1 (MQ=255)
                                       acCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATGATTGCCCGCCGTCTTTTCGGTTCAGTa           <  1:1809248/69‑1 (MQ=255)
                                       acCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATGATTGCCCGCCGTCTTTTCGGTTCAGTa           <  1:1826556/69‑1 (MQ=255)
                                       acCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATGATTGCCCGCCGTCTTTTCGGTTCAGTa           <  1:676515/69‑1 (MQ=255)
                                       acCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATGATTGCCCGCCGTCTTTTCGGTTCAGTa           <  1:671975/69‑1 (MQ=255)
                                       acCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATGATTGCCCGCCGTCTTTTCGGTTCAGTa           <  1:2832863/69‑1 (MQ=255)
                                       acCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATGATTGCCCGCCGTCTTTTCGGTTCAGTa           <  1:289999/69‑1 (MQ=255)
                                       acCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATGATTGCCCGCCGTCTTTTCGGTTCAGTa           <  1:3289122/69‑1 (MQ=255)
                                          aGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATGATTGCCCGCCGTCTTTTCGGTTCAGTATGt        <  1:1975825/69‑1 (MQ=255)
                                          aGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATGATTGCCCGCCGTCTTTTCGGTTCAGTATGt        <  1:3050767/69‑1 (MQ=255)
                                          aGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATGATTGCCCGCCGTCTTTTCGGTTCAGTATGt        <  1:3128129/69‑1 (MQ=255)
                                                aCCAGGAGAGCAATAAGCAACTCTCGCACCATGATTGCCCGCCGTCTTTTCGGTTCAGTATGTTTCAGt  >  1:978037/1‑69 (MQ=255)
                                                  |                                                                  
CTAAGAAATGCCAGAATTTTCTCACCCGCAAACAGGAACACCAGCATCACCAGGAGAGCAATAAGCAACTCTCGCACCATGATTGCCCGCCGTCTTTTCGGTTCAGTATGTTTCAGT  >  minE/2645828‑2645944

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