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.9% | 14.4 / 24.5 | 19 | G288G (GGT→GGG) | yfaL | AIDA‑I family autotransporter adhesin YfaL/EhaC |
Reads supporting (aligned to +/- strand): ref base A (8/0); new base C (0/11); total (8/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.32e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCTCCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAA‑CTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATC > NZ_CP009273/2336696‑2336874 | gCTCCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAA‑CTTCACTTAAGCCGAGATACATAAAGCCa > 1:667813/1‑90 (MQ=255) ttGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAA‑CTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGaggaagga > 1:531621/1‑90 (MQ=255) gcaaaTACCCAGCTTTTTTCCGCCGGAAATACAAAAGGAAA‑CTCCACTAAAGCCGAAAAAAAAAAACCCCCCCGCCGCAGAGGAAGGAcc < 2:416928/87‑1 (MQ=255) gcaaaTAACCAGCTTTTTTCCGTCGGAAATATAAAAGGTAA‑CTTCACTTAAGCCGAGATACAAAAACCCCCCCGCCGCAGAGGAAGGAcc < 2:28926/87‑1 (MQ=255) gCCAATAACCAGCTTTTTCCCGTCGGCAATATAAAAGGTAA‑CTTCACTTAAGCCGAGATAAAAAAAGCCCCCCGCCGCAGAGGAAGGAcc < 1:728995/90‑1 (MQ=255) gCCAATAACCAGCGTTTTTCCGTCGGCAATATAAAAGGTAA‑CTTCTCTTAAGCCGAGATAAAAAAACCCCCCCGCCGCAGAGGAAGGAcc < 1:851423/90‑1 (MQ=255) ccAATACCCAGCTTTTTTCCGTCGGCAATATAAAAGGTAA‑CTTCCCTTAAGCCGAGATAAAAAAAGCCCCCCGCCGCAGAGGAAGGACCa < 1:469376/90‑1 (MQ=255) ccAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAA‑CTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCa > 2:410611/1‑90 (MQ=255) caacaaCCAGCTTTTTTCCGTCGGCAATATCAAAGGTAA‑CTTCACTTAAGCCGAGATACATAAAGCCCCCCGCCGCAGAGGAAGGACCAt < 2:787612/90‑1 (MQ=255) aaTAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAA‑CTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCAtc > 2:330211/1‑90 (MQ=255) tAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAA‑CTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCAtctc > 2:307736/1‑90 (MQ=255) aGCGTTTTTCCGTCGGCAATATCAAAGGTAA‑CTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAg > 1:309993/1‑90 (MQ=255) aGCGTTTTTCCGTCGGCAATATCAAAGGTAA‑CTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAg > 2:73323/1‑90 (MQ=255) tttcCCGCCGACAATA‑CAAAAGTAACCTCCACTAAACCCGAAATAAAAAACCCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCctgc < 2:667813/90‑1 (MQ=255) cGGCAATAAAAAAGGAAA‑CTTCACTAAAGCCGAAATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTAt < 2:41188/90‑1 (MQ=255) cAAAGGTAA‑CTTCACTAAAGCCGAAATACAAAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCAt < 1:330211/90‑1 (MQ=255) aGGTAA‑CTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGa > 2:842410/1‑90 (MQ=255) gCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGtt > 1:226619/1‑90 (MQ=255) ccGAGATACAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTc < 2:535657/90‑1 (MQ=255) aaGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTg < 1:739389/90‑1 (MQ=255) aCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGtcatc > 1:747801/1‑90 (MQ=255) | GCTCCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAA‑CTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATC > 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 < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |