Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,336,785 | 0 | A | C | 54.5% | 18.9 / 25.6 | 22 | G288G (GGT→GGG) | yfaL | AIDA‑I family autotransporter adhesin YfaL/EhaC |
Reads supporting (aligned to +/- strand): ref base A (8/2); new base C (0/12); total (8/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.41e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCTCCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATC > NZ_CP009273/2336696‑2336874 | gCTCCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCa > 1:566015/1‑90 (MQ=255) ccGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCAccc > 1:464043/1‑90 (MQ=255) ttGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGaggaagga > 2:82500/1‑90 (MQ=255) gCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGAcc > 2:213518/1‑90 (MQ=255) cAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTACCTTCATTTAAGCCGAGATACAAAACGCCCCCCGCCGCAGAGGAAGGACCAt < 2:34435/90‑1 (MQ=255) aCCAGCTTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCa < 2:61802/90‑1 (MQ=255) ggctttttcCCGCCGGAAATAAAAAAGGTACCTTCACTAAAGCCGAAATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAg < 1:370378/86‑1 (MQ=255) aGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAg > 1:591046/1‑90 (MQ=255) ttttCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACAAAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCctgc < 2:101295/90‑1 (MQ=255) ttttCCGCCGCAAATAACAAAGTTACCTCCCCTAAAGCCGAGATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCctgc < 2:419997/90‑1 (MQ=255) tttcCCGTCGGCAATTACAAGGGAACCTTCTCTAAACCCGAGATAAAAAAACCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCctgc < 2:591046/90‑1 (MQ=255) tttcCCGCCGGAAATAAAAAAGGTAACTCCACTAAAGCCGAGATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCctgc < 1:82500/90‑1 (MQ=255) cGTCGCAAATATCAAAGGAACCCTCGCTTAAGCGGGGTCACTTAAGGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTgcgc < 1:318330/90‑1 (MQ=255) gTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTgcgcc > 2:370378/1‑89 (MQ=255) cAAAGGTAACTTCACTTAAGCCGAAATAAAAAACCCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCAt < 2:414216/90‑1 (MQ=255) aaaGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGaggaagga > 1:460237/1‑55 (MQ=255) aaaGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGaggaagga < 2:460237/55‑1 (MQ=255) gTAACTCCACTAAACCCGAGATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACt < 2:162819/90‑1 (MQ=255) aaCTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTaa > 1:582822/1‑90 (MQ=255) aaCTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTaa > 2:458297/1‑90 (MQ=255) ttCACTAAAGCCGAAAAAAAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAacac < 1:253066/90‑1 (MQ=255) ccGAAATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTc < 1:183578/90‑1 (MQ=255) aaaaCAAAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGc < 1:503708/87‑1 (MQ=255) aCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGtcatc > 2:129438/1‑90 (MQ=255) | GCTCCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATC > NZ_CP009273/2336696‑2336874 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |