| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 3,834,811 | 0 | T | G | 59.3% | 28.4 / 26.9 | 27 | D115E (GAT→GAG) | rfaP | kinase that phosphorylates core heptose of lipopolysaccharide |
| Reads supporting (aligned to +/- strand): ref base T (8/3); new base G (16/0); total (24/3) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.64e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGGGTTTGGCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGATCTCACTCCCACAATTAGCCTTGAAGATTATTG > W3110S.gb/3834751‑3834843 | tGGGTTTGGCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGATCTCa > 1:952219/1‑65 (MQ=255)tGGGTTTGGCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGATCTCa > 1:1521133/1‑65 (MQ=255)tGGGTTTGGCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGATCTCa > 1:614886/1‑65 (MQ=255)tGGGTTTGGCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGATCTCa > 1:2464331/1‑65 (MQ=255)tGGGTTTGGCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGATCTCa > 1:2414669/1‑65 (MQ=255)tGGGTTTGGCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGATCTCa > 1:2322374/1‑65 (MQ=255)tGGGTTTGGCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGATCTCa > 1:2050631/1‑65 (MQ=255)tGGGTTTGGCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGATCTCa > 1:1443964/1‑65 (MQ=255) ttGGCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGAGCacacgtc > 1:2015699/1‑61 (MQ=255) ttGGCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGAGCacacgtc > 1:2434124/1‑61 (MQ=255) ttGGCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGAGCacacgtc > 1:2418538/1‑61 (MQ=255) ttGGCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGAGCacacgtc > 1:2394054/1‑61 (MQ=255) ttGGCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGAGCacacgtc > 1:2152420/1‑61 (MQ=255) ttGGCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGAGCacacgtc > 1:2064175/1‑61 (MQ=255) ttGGCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGAGCacacgtc > 1:1913639/1‑61 (MQ=255) ttGGCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGAGCacacgtc > 1:1803910/1‑61 (MQ=255) ttGGCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGAGCacacgtc > 1:1758014/1‑61 (MQ=255) ttGGCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGAGCacacgtc > 1:1675135/1‑61 (MQ=255) ttGGCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGAGCacacgtc > 1:1570557/1‑61 (MQ=255) ttGGCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGAGCacacgtc > 1:1560843/1‑61 (MQ=255) ttGGCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGAGCacacgtc > 1:1535823/1‑61 (MQ=255) ttGGCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGAGCacacgtc > 1:1519992/1‑61 (MQ=255) ttGGCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGAGCacacgt > 1:1807192/1‑61 (MQ=255) ttGGCGAAAAAGGGTTAAATCCATTAAATCGCGCATCATTTATTATTACCGAAGAGCacacgtc > 1:489299/1‑61 (MQ=255) ttAACTCGCGCATCATTTATTATTACCGAAGATCTCACTCCCACAATTAGCCTTGAAGATTATtg < 1:190461/65‑1 (MQ=255) ttAACTCGCGCATCATTTATTATTACCGAAGATCTCACTCCCACAATTAGCCTTGAAGATTATtg < 1:2179425/65‑1 (MQ=255) ttAACTCGCGCATCATTTATTATTACCGAAGATCTCACTCCCACAATTAGCCTTGAAGATTATtg < 1:1571968/65‑1 (MQ=255) | TGGGTTTGGCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGATCTCACTCCCACAATTAGCCTTGAAGATTATTG > W3110S.gb/3834751‑3834843 |
| 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: ‑ |