Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,356,182 Δ1 bp 100% coding (444/1035 nt) yehA ← predicted fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,356,1820A.89.5% 130.9 / 11.9 38coding (444/1035 nt)yehApredicted fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base A (4/0);  new base . (0/34);  total (4/34)
Fisher's exact test for biased strand distribution p-value = 1.35e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCA  >  minE/1356139‑1356221
                                           |                                       
aTGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAg                   >  1:848996/1‑66 (MQ=255)
aTGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGc                  >  1:149631/1‑67 (MQ=255)
   aCTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCAc                >  1:1599960/1‑66 (MQ=255)
   aCTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCAc                >  1:484881/1‑66 (MQ=255)
                 aaaCTCAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:168155/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTTTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:1077931/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:82697/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:75417/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:751429/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:662981/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:860383/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:8965/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:432680/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:1025413/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:274002/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:228573/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:973272/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:1941633/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:1928271/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:1891591/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:1874318/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:1491977/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:1099255/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:1109684/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:1121858/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:1241015/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:1287280/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:1403830/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:1444673/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:1864583/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:1544016/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:1551045/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:1611019/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:1649535/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:1722546/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:1753272/65‑1 (MQ=255)
                 aaaCACAGTAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:1780047/65‑1 (MQ=255)
                         tAGAAATTGTACTGCGGCCTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa  <  1:792160/57‑1 (MQ=255)
                                           |                                       
ATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCA  >  minE/1356139‑1356221

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: