Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,336,785 | 0 | A | C | 57.1% | 27.6 / 36.6 | 29 | G288G (GGT→GGG) | yfaL | AIDA‑I family autotransporter adhesin YfaL/EhaC |
Reads supporting (aligned to +/- strand): ref base A (11/1); new base C (0/16); total (11/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.59e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-06 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCTCCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTT‑CACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATC > NZ_CP009273/2336696‑2336874 | gCTCCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTT‑CACTTAAGCCGAGATACATAAAGCCa > 1:746679/1‑90 (MQ=255) gTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTT‑CACTTAAGCCGAGATACATAAAGCCACccgc > 1:414785/1‑90 (MQ=255) tGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTT‑CACTTAAGCCGAGATACATAAAGCCCCCCGCCGCAGagga < 2:338802/90‑1 (MQ=255) gTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTT‑CACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGaggaa > 2:581184/1‑90 (MQ=255) gTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTT‑CACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGaggaa > 2:502617/1‑90 (MQ=255) gTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTT‑CACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGaggaa > 2:48035/1‑90 (MQ=255) tATTGCAAATACCCAGCTTTTTCCCGGCGGAAATAAAAAAGGTACTTT‑CACTTAAGCCGAGATAAAAAAACCCCCCCGCCGCAGaggaag < 1:833716/90‑1 (MQ=255) ccAAAAACCCGCGTTTTCCCCTCGGCAATATAAAAGGTAACTT‑CACTTAAGCCGAGATACAAAAACCCCCCCGCCGCAGAGGAAGGACCa < 1:581184/90‑1 (MQ=255) aaTAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTT‑CACTTAAGCCGAGATACATAAAGCCACCCGCCGCa < 2:50315/76‑1 (MQ=255) aaTAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTT‑CACTTAAGCCGAGATACATAAAGCCACCCGCCGCa > 1:50315/1‑76 (MQ=255) aGCGTTTTTCCGTCGGCAATATCAAAGGTAACTT‑CACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAg > 1:659288/1‑90 (MQ=255) aGCGTTTTTCCGTCGGCAATATCAAAGGTAACTT‑CACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAg > 1:783561/1‑90 (MQ=255) tttcCCGTCGGAAAAAAAAAAGGAACCTC‑CACTAAAGCCGAAATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCctgc < 1:344744/90‑1 (MQ=255) tttcCCGTCGCAAATATCAAAGG‑AACTTCCACTAACGCAAAAATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCctgc < 2:746679/90‑1 (MQ=255) cccGTCGGCAATAAAAAAGGTAACTT‑CACTAAAGCCGAGATACAAAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCtgctgc < 2:868597/89‑1 (MQ=255) cgccGGCAATCACAAAGGTACCTT‑CACTTAAGCCGAAATACAAAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTgcgc < 1:640709/87‑1 (MQ=255) ccGGAAATACAAAAGGTAACTT‑CACTAAAGCCGAAATACAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTa < 1:444200/89‑1 (MQ=255) aaaTAACAAAGTAACCTT‑CCCTAAAGCCGAAATAAAAAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGt < 1:502617/89‑1 (MQ=255) aaaaGGTACCTC‑CCCTAAACCCGAAATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCAt < 1:344097/89‑1 (MQ=255) aaaaGGTAACTC‑CACTAAAGCCGAAATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCAt < 1:298236/89‑1 (MQ=255) aaCTT‑CACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTaa > 2:311967/1‑90 (MQ=255) tccccTTAAGCCGAGATAAAAAAGCCCGCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAacac < 1:497949/86‑1 (MQ=255) aCTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGcc > 1:698981/1‑90 (MQ=255) ctaaaGCCGAGAAACAAAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCt < 2:698981/87‑1 (MQ=255) cTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCct > 1:494361/1‑57 (MQ=255) cGAAATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCt < 2:659288/90‑1 (MQ=255) gagaTACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTg > 2:757083/1‑90 (MQ=255) gaaaTACAAAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTg < 2:202324/87‑1 (MQ=255) aaaaCCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGc < 2:305953/89‑1 (MQ=255) aaaaCCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGc < 2:6905/89‑1 (MQ=255) aCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGtcatc > 1:326858/1‑90 (MQ=255) | GCTCCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTT‑CACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATC > NZ_CP009273/2336696‑2336874 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |