Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,673,569 | 0 | A | C | 63.2% | 5.4 / 15.0 | 19 | T296P (ACC→CCC) | ydgH | DUF1471 family protein YdgH |
Reads supporting (aligned to +/- strand): ref base A (6/1); new base C (0/12); total (6/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.58e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.51e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCT > NZ_CP009273/1673493‑1673655 | cTACGGAAAAATGCCCGAAGTCTCCTATCAGGTTGCGAAACGGGCCGCGAAGAAAGGTGCAAAGTATTACCACACCCCCCGCCAGTGGCa < 2:7602/90‑1 (MQ=255) accATGACGGACTTCTCCTCTCGGGTGGGGAAGCGGCGCGCGAAAAAAGTTGCTAATTTTTCCCACCCCCCCCGCCAGTGGCAGGAACGt < 1:494470/88‑1 (MQ=255) aaCATGCCCAAAGTCTCCTATCGGGTTGGAAAACGTCCCGAGAAAAAAGGTGCAAATTTTTCCCACACCCCCCGCCAGTGGCAGGAACGt < 1:150352/90‑1 (MQ=255) tgcccGAAGTCTCCTATCAGTTTGCGAAACGTGCCGCGAAAAAAGGTGCAAATTATTCCCACATCCCCCGCCAGTGGCAGGAACGTGGta < 2:353922/87‑1 (MQ=255) gACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGtaa > 1:560626/1‑90 (MQ=255) tATCGGGTTCCGAAACGCGCGGAGAAAAAAGGTGCTAAGTTTCCCCCCACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTc < 1:83604/90‑1 (MQ=255) tATCAGGTGGCGAAACGTGCCGAAAAGAAGGGTGCTAAGTTTTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTc < 1:552196/90‑1 (MQ=255) ggTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGGAATAACCTGACCGTCAgcgc > 1:234349/1‑90 (MQ=255) ggTTGCAAAGCGTGCCGGGAAGAAGGGTGCAAAGTTTTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc < 1:414045/90‑1 (MQ=255) ggTTGCAAACCGGGCCGCGAAAAAAGGTGCTAATTATTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc < 2:212819/90‑1 (MQ=255) aaaCGTGCCGAGAAGAAGGGTGCTAAGTTTTCCCCCCTCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTg < 1:12374/90‑1 (MQ=255) cgagaagaaGGGTGCTAAGTTTTCCCCCCCCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTataaata < 1:494725/87‑1 (MQ=255) cgagaaaaaAGTGGCGAAGTATTCCCCCATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTataaata < 1:172726/87‑1 (MQ=255) aaaaGGTGATAAGTATTCCCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTAtt < 2:620705/89‑1 (MQ=255) cTAAGTATTCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACt < 2:341896/90‑1 (MQ=255) ttACCACATCACCCGCCAGTGGCAGGAACGTGGGAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATaa > 1:341896/1‑90 (MQ=255) ttACCACATCACCCGCCAGTGGCAGGAACGTGGGAATAACCTGACCGGCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATaa > 2:474953/1‑90 (MQ=255) acaTCACCCGCCAGTGGCAGGAACGGGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAggcgg > 1:55072/1‑90 (MQ=255) aTCACCCGCCAGTGGCAGGAACGTGGGAATAAACTGACCGTCAGCGCAGATCTGTATAAATAAGTATTCAGTTACTCCTATAATGCGGCt > 1:209912/1‑90 (MQ=255) | CTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCT > NZ_CP009273/1673493‑1673655 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |