Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 564,054 C→A D283Y (GAT→TAT)  ybjL ← predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE564,0540CA82.5% 72.6 / 19.9 40D283Y (GAT→TAT) ybjLpredicted transporter
Reads supporting (aligned to +/- strand):  ref base C (0/7);  new base A (0/33);  total (0/40)
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.89e-01

ACTTCTTTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACG  >  minE/564026‑564089
                            |                                   
aCTTCTTTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCaa    <  1:857582/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:1586760/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:1066337/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:1202207/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:307588/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:305394/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:292898/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:1476546/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:1766024/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:161898/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:1768585/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:1777840/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:1783770/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:1819126/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:1850728/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:198415/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:326476/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:585242/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:661910/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:717218/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:783881/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:998592/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:1700560/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:1016268/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:1040269/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:1049113/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:1116566/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:1182818/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:1191366/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:1205462/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:1260444/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:1347498/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:1360197/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:1429792/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:1530177/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:1543107/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:155391/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:1007311/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACg  <  1:1710150/62‑1 (MQ=255)
  ttcttTACCGTTACGGAAGCTGGGATAGAGTCGGGCATGGGCGGCGGGATAGCCTACCAACg  <  1:1009660/62‑1 (MQ=255)
                            |                                   
ACTTCTTTACCGTTACGGAAGCTGGGATCGAGTCGGGCATGGGCGTCGGGATAGCCTACCAACG  >  minE/564026‑564089

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: