Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,334,224 | 0 | T | A | 56.0% | 4.3 / 26.4 | 25 | D1142V (GAT→GTT) | yfaL | AIDA‑I family autotransporter adhesin YfaL/EhaC |
Reads supporting (aligned to +/- strand): ref base T (4/7); new base A (14/0); total (18/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.86e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.48e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGCCATTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATCACCACACCACGCCCCGGTAACCACTGATACCCCG > NZ_CP009273/2334140‑2334305 | cGCCATTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATcc < 2:805562/90‑1 (MQ=255) gCCATTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCct < 1:594137/90‑1 (MQ=255) aTTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAatcat > 1:195593/1‑84 (MQ=255) aTTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAatcat < 2:195593/84‑1 (MQ=255) aTTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAAACATCccgct > 2:23628/1‑90 (MQ=255) aTTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGGAAAAACATCctgct > 2:23224/1‑90 (MQ=255) aTTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGGAAAAACATCctgct > 2:197290/1‑90 (MQ=255) aTTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGGAAAAAAAACctgct > 2:194100/1‑90 (MQ=255) tAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTg < 1:23224/90‑1 (MQ=255) tAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTg < 1:845600/90‑1 (MQ=255) aaCCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAAACATCCCGCTGCCCCCCCCGAt > 2:88411/1‑90 (MQ=255) aGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAAACATCCTGCTGCACGCCCTGATAAATc > 1:942982/1‑90 (MQ=255) aGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAAACATCCTGCTGCACGCCCTGATAAATc > 1:204198/1‑90 (MQ=255) aGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAAACATCCCGCTGCACGCCCCGATAAAta > 1:905158/1‑89 (MQ=255) aGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAAAAATCCCCCTGCCCCCCCCGATAAATc > 2:59728/1‑90 (MQ=255) aGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGGAAAAACATCCTGCTGCCCGCCCCGATaaaac > 1:907449/1‑88 (MQ=255) tcatcaCCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTCCTGCACGCCCCGGTAAATCACCCGCGCCCGCggg > 2:686232/1‑89 (MQ=255) catcaCCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGtt < 1:747575/90‑1 (MQ=255) aCCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGAAAAAAACCTCTGCCCCACGCCCCGTTAAATAACCCCCGCCCGCCGGTCaaa > 1:880124/1‑89 (MQ=255) ccTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAAACATCCTGCTGCCCGCCCCGATAAAACACCCGCGCCCGCGGGTCAATca > 2:522755/1‑90 (MQ=255) gtgACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATcacca > 1:315111/1‑84 (MQ=255) gtgACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATcacca < 2:315111/84‑1 (MQ=255) gtgACACGCGCGCACGGTTAGCGGCGGTAAAAACATCCTGCTGCCCCCCCCGGTAAATCACCCGCGCCTGCGGGTCAATCaccaccccac > 1:592800/1‑90 (MQ=255) gcgcgcACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGGTCAATCACCACACCACGcccccg > 1:840123/1‑88 (MQ=255) gcgcACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATCACCACACCACGCCCCGGTa < 2:592800/90‑1 (MQ=255) cgcACGGTTAGCGGCGGTAAAAACATCCTGCTGCACGCCCCGATAAATCACCTGCGCCTGCGGGTCAATCACCACACCACGCCCCCGTaa > 1:134186/1‑90 (MQ=255) cggTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATCACCACACCACCCCCCGGTAACCACTGATAccccc > 1:254519/1‑89 (MQ=255) | CGCCATTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATCACCACACCACGCCCCGGTAACCACTGATACCCCG > NZ_CP009273/2334140‑2334305 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |