Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 2,025,196 +C coding (167/786 nt) fliR → flagellar export pore protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,025,1921.C97.6% 147.7 / ‑3.0 41coding (163/786 nt)fliRflagellar export pore protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (0/40);  total (0/41)
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.95e-01

GATGATCACGTTCGCCATTGCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCTTT  >  W3110S.gb/2025173‑2025236
                    |                                             
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:413408/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:209567/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:2129233/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:2162701/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:2444473/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:2543295/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:294915/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:298338/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:339290/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:409720/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:1109982/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:437623/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:475323/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:546210/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:552821/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:583019/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:627524/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:738524/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:865575/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:908931/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:1487924/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:1021156/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:1173397/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:1231510/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:1288437/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:1320430/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:1354754/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:1453813/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:1503405/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:1529075/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:1531192/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:1723107/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:1848967/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:1959660/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:1978005/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:2037714/65‑1 (MQ=255)
gatgatCACGTTCGCCATTGCGCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:1339949/65‑1 (MQ=255)
 atgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:8184/64‑1 (MQ=255)
 atgatCACGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:2056747/64‑1 (MQ=255)
        cGTTCGCCATTGCCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:2064042/57‑1 (MQ=255)
           tCGCCATTGCCCCCAGCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCttt  <  1:1423530/54‑1 (MQ=255)
                    |                                             
GATGATCACGTTCGCCATTGCCCCATCATTACCTGCCAACGATGTTCCTGTTTTTTCGTTCTTT  >  W3110S.gb/2025173‑2025236

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