Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,294,637 | 0 | T | A | 53.8% | 8.7 / 20.0 | 26 | T652T (ACT→ACA) | PP_2893 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (2/10); new base A (14/0); total (16/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.24e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.23e-02 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACGTTGCAGGCCAACGGCAAGGGCGACATTCGCCAACAGGCCTTGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCACTGGCCGAATCCAGGCCGGTGGCCAGGAGTGGCAGTTGCAGGCGCCTGCACGCCTGCAGCGGCTGGCCAGCGGCCAGCTCGACTTCGGGGCCCACTGCTGGCGCTCCGGCCAGGCAAGCCTGTGT > NC_002947/3294503‑3294760 | aCGTTGCAGGCCAACGGCAAGGGCGACATTCGCCAACAGGCCTTGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGactggcccctc > 1:428270/1‑135 (MQ=255) gCAGGCCAACGGCAAGGGCGACATTCGCCAACAGGCCTTGACCCTGGCTCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACAGGCCACAGGCCgc > 1:979256/1‑135 (MQ=255) gCAGGCCAACGGCAAGGGCGACATTCGCCAACAGGCCTTGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCCCAGGCCgc > 2:687601/1‑135 (MQ=255) gCAGGCCAACGGCAAGGGCGACATTCGCCAACAGGCCTTGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCCCAGGCCGa > 1:557644/1‑136 (MQ=255) aGGCCAACGGCAAGGGCGACATTCGCCAACAGGCCTTGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCCCAGGCCgcaa > 2:705684/1‑133 (MQ=255) aGGCCAACGGCAAGGGCGACATTCGCCAACAGGCCTTGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCACAGGCCGaaa > 2:679304/1‑135 (MQ=255) ccAACGGCAAGGGCGACATTCGCCAACAGGCCTTGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCACTGGCCGAATCCa < 2:606457/136‑1 (MQ=255) aaCGGCAAGGGCGACATTCGCCAACAGGCCTTGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCACTGGCCGAATCCAgg < 2:236439/136‑1 (MQ=255) aaCGGCAAGGGCGACATTCGCCAACAGGCCTTGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCACTGGCCGAATCCAgg < 2:162693/136‑1 (MQ=255) aaCGGCAAGGGCGACATTCGCCAACAGGCCTTGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCACTGGCCGAATCCAgg < 2:686377/136‑1 (MQ=255) cAAGGGCGACATTCGCCAACAGGCCTTGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCCCTGGGCGCATACCGGcccgg > 1:49406/1‑133 (MQ=255) cAAGGGCGACATTCGCCAACAGGCCTTGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTCGCCTGGATGGCCAGTTGATCAAAGGTGACTGGCGCGGCCGAGAGGCCACAGGCCGAATCCAGGCCggg > 2:231665/1‑135 (MQ=255) gCGACATTCGCCAACAGGCCTTGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACATGCCACACGCCGCAACCCCGCCGGTgggcc > 2:400002/1‑133 (MQ=255) gCCTTGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACGGGCCCCAGGCCGAATACACGCCCGGGGCCAGGAGTGGCAGTTTCAgg > 1:320683/1‑136 (MQ=255) gtGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGAGGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCACTGGCCGAATCCAGGCCGGTGGCCAGGAGTGGCAGTTGCAGgcgc < 1:715816/135‑1 (MQ=255) gACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGAGAGGCCACAGGCCGAATACACGCCGGTGGGCCGGCGTGGCAGTTTCAGGcgccc > 1:199021/1‑135 (MQ=255) gcgcGGGACGGCCCGCAGCTCACGCGCGACCTTGGCCTGGAGGGCCAGTTGAACAATGGTGCCTGGCGCGGCCGACTGGCCACTGGCCGAATCCTGGCCGGTGGCCTGGAGTGGCTGTTGCAGGcgcctgcacgcc < 2:591378/136‑1 (MQ=255) cgcTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGAGGGGCCCCAGGGCGCATACCGGGCGGGGGGCCGGCGTGGCAGTTGGCGGcgccggcccgccc > 1:728770/1‑135 (MQ=255) gcTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCCCAGGCCGCATCCCGGCCGGGGGGCAGGCGTGGCAGGTGCAGGcgccggcccgccgg > 1:282786/1‑134 (MQ=255) ccGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCACTGGCCGAATCCAGGCCGGTGGCCAGGAGTGGCAGTTGCAGGCGCCTGCACGCCTGCAGCGGCTgg < 1:761031/136‑1 (MQ=255) cGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCTACTGGCCGAATCCAGGCCGGTGGCCAGGAGTGGCAGTTGCAGGCGCCTGCACGCCTGCAGCGGCTggc > 1:310568/1‑136 (MQ=255) cGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCACAGGCCGAATCCAGGCCGGTGGCCAGGAGTGGCAGTTGCAGGCGCCTGCACGCCTGCAGCGGCTggc > 2:1005955/1‑136 (MQ=255) cGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCACAGGGCGAATACACGCCGGGGGGCAGGAGTGGGAGTTTGCGGCGCCTGCCCGCCGGGCGCGGCGGGCCCGCGGGCCGCTc > 1:890581/1‑136 (MQ=255) ttGAACAACGGTGACGGGCGCGGCCGACTGCCCACTGGCCGAATCCCGCCCGGTGGCCAGTAGTGGCAGTTGCAGGCGCCTGCACGCCTGCAGCGGCTGGCCAGCGGCCAGCTCGACTTCGGGGCCCACTGCTGgc < 1:705684/136‑1 (MQ=255) ttGAACAAAGGTGACTGGCGCGGCCGACTGGCCACTGGCCGAAACCAGGCCGGTGGCCAGGAGTGGCAGTTGCAGGCGCCTGCACGCCTGCAGCGGCTGGCCAGCGGCCAGCTCGACTTCGGGGCCCACTGCTGgc < 1:679304/136‑1 (MQ=255) tGACTGGCGCGGCCGACGGGCCACAGGCCGAATACACGCCCGTGGCCAGGAGTGGCAGTTTCAGGCGCCTGCACGCCGGCAGCCGCAGGCCCGCCGCCACCGCGACTTCGGGGCCCCCATCAGGCGCTACCGcccc > 2:547571/1‑134 (MQ=255) ccGACGGGCCACTGGCCGAAGCCAGGCCGGGGGCCAGGAGTGGCAGTTGCAGGCGCCTGCACGCCTGCAGCGGCTGGCCAGCGGCCAGCTCGACTTCGGGGCCCACTGCTGGCGCTCCGGCCAGGCAAGCCtgtgt < 1:440585/136‑1 (MQ=255) | ACGTTGCAGGCCAACGGCAAGGGCGACATTCGCCAACAGGCCTTGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCACTGGCCGAATCCAGGCCGGTGGCCAGGAGTGGCAGTTGCAGGCGCCTGCACGCCTGCAGCGGCTGGCCAGCGGCCAGCTCGACTTCGGGGCCCACTGCTGGCGCTCCGGCCAGGCAAGCCTGTGT > NC_002947/3294503‑3294760 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |