Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 228,196 | 0 | T | G | 54.5% | 8.0 / 10.3 | 13 | R196R (CGT→CGG) | yafD | endonuclease/exonuclease/phosphatase family protein |
Reads supporting (aligned to +/- strand): ref base T (2/3); new base G (5/1); total (9/4) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.42e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.05e-02 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTG > NZ_CP009273/228114‑228284 | aGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGGTTTAATGCCTGGGGGCGGAGAAGGa > 1:312058/1‑90 (MQ=255) cAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTAATGCCTGGGGGCGGAGAAGGGTg > 2:321795/1‑90 (MQ=255) aCTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAAcgc < 2:676355/90‑1 (MQ=255) aTTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGTTTTAATTACTGGGGCCGGAGGAGGGAGAAAGCGGTttat > 1:130652/1‑90 (MQ=255) tCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTAATGCCTGGGGCCGTAGAAGGATGAAAGCGGTATATTGCTTTgc > 2:64818/1‑90 (MQ=255) caccacAGCGGCCCGGTCATTATGGCGGGGGGTTTCAATGCCTTGGGCCGCAGGAGGAGGAACGCGGTATTTTGCTTTTCTCGGGGAATg > 1:597347/1‑90 (MQ=255) caccacAGCGGCCCGGTCATTATGGCGGGGGATTTTAATTCCCGGGGCCGGAGAAGGGTGAACGCGGTATATTGCTTTTCGCGGGgaaat > 2:453878/1‑88 (MQ=255) ccacAGCGGCCCGGTCATTATGGCGGGGGATTTTAATGCCTGGGGGCCGAGAAGGGGGAAAGCGGTTTTTTTCTTTTCGCGGGGAAAGTc > 2:334912/1‑90 (MQ=255) caGCGGCCCGGTCATTATGGCGGGGGATTTTAATGCCTGGGGGCGCAGAAGGAAGAAAGCGTTATATTGCTTTGTGCGGGGAAAGTCGcg > 2:387787/1‑89 (MQ=255) caGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCt < 1:231367/90‑1 (MQ=255) cGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAgg < 2:641357/90‑1 (MQ=255) cGGTCAATAAGGCCGGAGATTTCAATGCCTGGAGCCGGAAAAAGATGAACCCGTTTTATCCCCTTGCGCGGGAAATGTCGCTGCGCCAgg < 2:457783/90‑1 (MQ=255) gTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTg > 1:214799/1‑90 (MQ=255) | AGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTG > NZ_CP009273/228114‑228284 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |