Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,336,785 | 0 | A | C | 55.6% | 14.0 / 19.3 | 18 | G288G (GGT→GGG) | yfaL | AIDA‑I family autotransporter adhesin YfaL/EhaC |
Reads supporting (aligned to +/- strand): ref base A (6/2); new base C (0/10); total (6/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.51e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCTCCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATC > NZ_CP009273/2336696‑2336874 | gCTCCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCa > 1:159957/1‑90 (MQ=255) gCTCCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCa > 2:218/1‑90 (MQ=255) gTATTGCAAAAACCCAGCGTTTTTCCGTCGGAAATATAAAAGGTAACTTCCCTTAAGCCGAGATAAATAAAGCCCCCCGCCGCAGaggaa < 1:256290/90‑1 (MQ=255) ttcccAACACCCAGCTTTTTCCCGGCGGAAATTACAAAGGAAACTTCACTTAAGCCGAGATAAAAAACCCCCCCCGCCGCAGaggaagga < 1:728974/87‑1 (MQ=255) gCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGAcc > 2:265366/1‑90 (MQ=255) ccAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCa < 2:541476/90‑1 (MQ=255) cAATAACCAGCTTTTTTCCGTCGGAAATATAAAAGGTAACTTCACTTAAGCCGAGATAAAAAAAGCCCCCCGCCGCAGAGGAAGGACCAt < 1:577211/90‑1 (MQ=255) cAATAACCAGCGTTTTTCCGTCGGAAATATAAAAGGTAACTTCCCTTAAGCCGAGATACAAAAAGCCCCCCGCCGCAGAGGAAGGACCAt < 2:335853/90‑1 (MQ=255) aCCAGCTTTTTCCCGCCGGCAATATAAAAGGTAACTTCCCTTAAGCCGAAATAAAAAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCa < 1:507748/90‑1 (MQ=255) aGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAg > 2:912697/1‑90 (MQ=255) aGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAg > 2:675275/1‑90 (MQ=255) cccgccgGCAATAAAAAAGGTAATTTCACTAACGCCGAAATACATAACGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCtgctgc < 2:750748/89‑1 (MQ=255) taaaaGGTACCTTCACTTAAGCCGAAATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCa < 1:74763/88‑1 (MQ=255) aaGGTACCTTCACTTAAGCCGAGATACATAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCg < 2:235075/90‑1 (MQ=255) aGGTACCTCCACTAAAGCCGAGATACAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGa < 2:268754/90‑1 (MQ=255) gaaaCTCCACTAACGCCGAGATACAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACt < 1:565007/88‑1 (MQ=255) accTTCACTTAAGCCGAGATACAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTaa < 2:211405/88‑1 (MQ=255) aaCTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTaa > 1:1133147/1‑90 (MQ=255) ttAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGccct > 1:268754/1‑88 (MQ=255) ccACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTg > 1:691192/1‑87 (MQ=255) aCCCGCCGCAGAGGAAGGACCATCTCCATAGCctgc > 2:1121117/1‑36 (MQ=255) aCCCGCCGCAGAGGAAGGACCATCTCCATAGCctgc < 1:1121117/36‑1 (MQ=255) aCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGtcatc > 1:1442427/1‑90 (MQ=255) aCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGtcatc > 1:1329730/1‑90 (MQ=255) | GCTCCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATC > NZ_CP009273/2336696‑2336874 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |