Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,334,224 | 0 | T | A | 52.6% | 13.5 / 27.4 | 19 | D1142V (GAT→GTT) | yfaL | AIDA‑I family autotransporter adhesin YfaL/EhaC |
Reads supporting (aligned to +/- strand): ref base T (0/9); new base A (10/0); total (10/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.08e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCCATTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATCACCACACCACGCCCCGGTAACCACTGATACCCCGCCTC > NZ_CP009273/2334141‑2334309 | gCCATTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCct < 1:40669/90‑1 (MQ=255) gCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTGAAATCATCCTGCTGCAc < 2:744658/90‑1 (MQ=255) gCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCAc < 2:477256/90‑1 (MQ=255) aGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGGAAAAACATCCCGCTGCACGCCCCGATaaaac > 2:387391/1‑88 (MQ=255) aGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGGAAAAAAATACTGCTGCCCGCCCCGATaaaac > 2:477537/1‑88 (MQ=255) atatCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCAcctgcgcctgc < 2:360446/90‑1 (MQ=255) tcaCCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCa < 1:647875/90‑1 (MQ=255) aCCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAAACATCCTGCTGCACGCCCCGATAAATCACCTGCGCCCGCGGGTCAAt > 2:199672/1‑90 (MQ=255) cTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAAACATCCTGCTGCACGCCCCGATAAAACACCTGCGCCTGCGGTTCAAacac > 2:647875/1‑90 (MQ=255) cTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAAACATCCTGCTGCACGCCCCGATAAAACACCCGCGCCCGCGGGTCAATcac > 2:346596/1‑90 (MQ=255) cTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAAACATCCCGCTGCCCGCCCCGATAAAAAACCCGCGCCCGCGGGTTAATcac > 2:150896/1‑90 (MQ=255) gCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAat > 2:296933/1‑39 (MQ=255) gCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAat < 1:296933/39‑1 (MQ=255) gCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAat < 1:113634/39‑1 (MQ=255) gCGATTGTGACACGCGCGCACGGTTAGCGGCGGGAAAAAAATCCTCCTGCCCGCCCCGATAAATTACCTGCGCCCTCGGGGCAATcacca > 1:461765/1‑90 (MQ=255) cGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATcaccac < 2:585806/90‑1 (MQ=255) cacGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATCACCACACCACGccc < 1:199672/90‑1 (MQ=255) gcgcgcACGGTTAGCGGCGGTAAAAACATCCTGCTGCACGCCCCGATAAAACACCCGCGCCCGCGGGTCAATCACCACACCAcccccccg > 2:334989/1‑88 (MQ=255) gcggTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATCACCACACCACGCCCCGGTAACCACTGATAcccc < 2:589100/90‑1 (MQ=255) ggaaaaACATCCTGCTGCACCCCCCGATAAAACACCCGCGCCCGCGGGTCAAACACCCCACCCCCCCCCCGCAACCACTGATACcccccc > 2:458265/3‑90 (MQ=255) gTAAAAACATCCTTCTGCACGCCCCGATAAAACACCCGCGCCCGCGGGTCAATCACCACACCCCCCCCCCGGAAACACCGATACcccccc > 2:322235/1‑90 (MQ=255) aaaaTCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATCACCACACCACGCCCCGGTAACCACTGATACCCCGCCTc < 1:458265/90‑1 (MQ=255) | GCCATTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCACCTGCGCCTGCGGTTCAATCACCACACCACGCCCCGGTAACCACTGATACCCCGCCTC > NZ_CP009273/2334141‑2334309 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |