Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,336,785 | 0 | A | C | 60.0% | 9.3 / 19.6 | 20 | G288G (GGT→GGG) | yfaL | AIDA‑I family autotransporter adhesin YfaL/EhaC |
Reads supporting (aligned to +/- strand): ref base A (7/1); new base C (0/12); total (7/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.03e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.43e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCAT > NZ_CP009273/2336699‑2336873 | ccGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCAccc > 2:228327/1‑90 (MQ=255) gTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCCCCCGCCGCAGaggaa < 2:158830/90‑1 (MQ=255) tGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGAc > 1:263841/1‑90 (MQ=255) gcaaaTAACCAGCGTTTTTCCGTCGGAAATATCAAAGGTAACTTCACTTAAGCCGAGATACAAAAAGCCCCCCGCCGCAGAGGAAGGAcc < 2:296881/87‑1 (MQ=255) gCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGAcc > 2:289244/1‑90 (MQ=255) cAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCCCCCGCCGCAGAGGAAGGACCAt < 2:52810/90‑1 (MQ=255) aaaaaCCAGCGTTTTCCCGCCGCCAATAACAAGGGTACCTCCCCTTAAGCCGAAATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCAtc < 1:264168/87‑1 (MQ=255) aaTACCCAGCTTTTTTCCGCCGGCAATATAAAAGGTAACTTCACTAAAGCCGAGATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCAtc < 2:483810/90‑1 (MQ=255) aaTAACCACCTTTTTCCCGTCGGAAATAACAAAGGTAACTTCCCTAAAGCCGAAAAACAAAAACCCCCCCGCCGCAGAGGAAGGACCAtc < 1:80011/90‑1 (MQ=255) cccAGCTTTTTCCCGCCGGAAATATCAAAGGTAACTTCACTTAAGCCGAAATACAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCa < 2:361548/89‑1 (MQ=255) aGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAg > 1:347092/1‑90 (MQ=255) gCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGc > 1:343238/1‑90 (MQ=255) ttttCCGCCGGCAATACAAAAGGTAACTTCACTAAAGCCGAGATACAAAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCctgc < 2:552727/90‑1 (MQ=255) tttcCCGTCGGCAATATAAAAGGTAACTTCCCTAAAGCCGAGATAAAAAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCctgc < 1:101513/90‑1 (MQ=255) cttcccGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCctgc < 2:245899/86‑1 (MQ=255) cccGTCGCCAATAAAAAAGGTAACTTCACTTAAGCGGAGATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCtgctgc < 1:289244/89‑1 (MQ=255) gacaCTTCATTAAACCCGAGAAACATAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACt < 1:122749/87‑1 (MQ=255) aTACAAAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCt < 2:359399/90‑1 (MQ=255) cACCCGCCGCGGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGtcat > 2:429310/1‑90 (MQ=255) cACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGtcat > 2:99625/1‑90 (MQ=255) | CCGTCATTCTCTGTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCAT > NZ_CP009273/2336699‑2336873 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |