| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,204,124 | 1 | . | G | 21.1% | 41.7 / 15.5 | 19 | coding (805/1326 nt) | zraR | fused DNA‑binding response regulator in two‑component regulatory system with ZraS: response regulator/sigma54 interaction protein |
| Reads supporting (aligned to +/- strand): ref base . (7/8); new base G (1/3); total (8/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.03e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ATAAACGCCGGGAAGGGCGCTTTGTTGAGGCGGACGGCGGCACGCTGTTTCTCGATGAAATTGGCGATATCTCGCCGATGATGCAGGTGCGTCTGCTGCGTGCGATTCAGGAGCGCGAAGTTCAGCGTGTCGGTAGCAACC‑AGATTATCTCGGTTGATGTCCGGCTGATTGCGGCGACCCATCGCGATCTTGCCGCAGAGGTGAATGCCGGGCGTTTTCGCCAGGATCTCTACTATCGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGC > NC_000913/4203984‑4204265 | aTAAACGCCGGGAAGGGCGCGTTGTTGAGGCGGACGGCGGCACGCTGTTTCTCGATGAAATTGGCGATATCTCGCCGATGATGCAGGTGCGTCTGCTGCGTGCGATTCAGGAGCGCGAAGTTCAGCGTGTCGGTAGCAACC‑AGATTAtc < 7:156711/149‑1 (MQ=255) aCGCCGGGAAGGGCGCTTTGTTGAGGCGGACGGCGGCACGCTGTTTCTCGATGAAATTGGCGATATCTCGCCGATGATGCAGGTGCGTCTGCTGCGTGCGATTCAGGAGCGCGAAGTTCAGCGTGTCGGTAGCAACCGAGATAAtctgg > 2:36513/1‑147 (MQ=255) ggAAGGGCGCTTTGTTGAGGCGGACGGCGGCACGCTGTTTCTCGATGAAATTGGCGATATCTCGCCGATGATGCAGGTGCGTCTGCTGCGTGCGATTCAGGAGCGCGAAGTTCAGCGTGTCGGTAGCAACCGAGATAATCT‑GGTTgcta < 1:36513/149‑4 (MQ=255) ggAAGGGCGCTTTGTTGAGGCGGACGGCGGCACGCTGTTTCTCGATGAAATTGGCGATATCTCGCCGATGATGCAGGTGCGTCTGCTGCGTGCGATTCAGGAGCGCGAAGTTCAGCGTGTCGGTAGCAACCGAGATAATCT‑GGTTgcta < 7:129799/149‑4 (MQ=255) gAAGGGCGCTTTGTTGAGGCGGACGGCGGCACGCTGTTTCTCGATGAAATTGGCGATATCTCGCCGATGATGCAGGTGCGTCTGCTGCGTGCGATTCAGGAGCGCGAAGTTCAGCGTGTCGGTAGCAACCGAGATAATCT‑GGTTgctac < 5:78723/149‑5 (MQ=255) gAAGGGCGCTTTGTTGAGGCGGACGGCGGCACGCTGTTTCTCGATGAAATTGGCGATATCTCGCCGATGATGCAGGTGCGTCTGCTGCGTGCGATTCAGGAGCGCGAAGTTCAGCGTGTCGGTAGCAACC‑AGATTATCTCGGTTGATGt > 2:173152/1‑149 (MQ=255) ggCGGACGGCGGCACGCTGTTTCTCGATGAAATTGGCGATATCTCGCCGATGATGCAGGTGCGTCTGCTGCGTGCGATTCAGGAGCGCGAAGTTCAGCGTGTCGGTAGCAACC‑AGATTATCTCGGTTGATGTCCGGCTGATTGCGGCGa < 1:173152/149‑1 (MQ=255) aCGCTGTTTCTCGATGAAATTGGCGATATCTCGCCGATGATGCAGGTGCGTCTGCTGCGTGCGATTCAGGAGCGCGAAGTTCAGCGTGTCGGTAGCAACC‑AGATTATCTCGGTTGATGTCCGGCTGATTGCGGCGACCCATCGCGATCt > 6:79990/1‑149 (MQ=255) gCTGTTTCTCGATGAAATTGGCGATATCTCGCCGATGATGCAGGTGCGTCTGCTGCGTGCGATTCAGGAGCGCGAAGTTCAGCGTGTCGGTAGCAACC‑AGATTATCTCGGTTGATGTCCGGCTGATTGCGGCGACCCATCGCGATCTTg < 5:79990/149‑1 (MQ=255) gATGAAATTGGCGATATCTCGCCGATGATGCAGGTGCGTCTGCTGCGTGCGATTCAGGAGCGCGAAGTTCAGCGTGTCGGTAGCAACC‑AGATTATCTCGGTTGATGTCCGGCTGATTGCGGCGACCCATCGCGATCTTGc < 3:128540/140‑1 (MQ=255) gATGAAATTGGCGATATCTCGCCGATGATGCAGGTGCGTCTGCTGCGTGCGATTCAGGAGCGCGAAGTTCAGCGTGTCGGTAGCAACC‑AGATTATCTCGGTTGATGTCCGGCTGATTGCGGCGACCCATCGCGATCTTGc > 4:128540/1‑140 (MQ=255) tcGCCGATGATGCAGGTGCGTCTGCTGCGTGCGATTCAGGAGCGCGAAGTTCAGCGTGTCGGTAGCAACC‑AGATTATCTCGGTTGATGTCCGGCTGATTGCGGCGACCCATCGCGATCTTGCCGCAGAgg > 5:72775/1‑130 (MQ=255) tcGCCGATGATGCAGGTGCGTCTGCTGCGTGCGATTCAGGAGCGCGAAGTTCAGCGTGTCGGTAGCAACC‑AGATTATCTCGGTTGATGTCCGGCTGATTGCGGCGACCCATCGCGATCTTGCCGCAGAgg < 6:72775/130‑1 (MQ=255) atgatgCAGGTGCGTCTGCTGCGTGCGATTCAGGAGCGCGAAGTTCAGCGTGTCGGTAGCAACC‑AGATTATCTCGGTTGATGTCCGGCTGATTGCGGCGACCCATCGCGATCTTGCCGCAGAGGt > 3:129035/1‑125 (MQ=255) atgatgCAGGTGCGTCTGCTGCGTGCGATTCAGGAGCGCGAAGTTCAGCGTGTCGGTAGCAACC‑AGATTATCTCGGTTGATGTCCGGCTGATTGCGGCGACCCATCGCGATCTTGCCGCAGAGGt < 4:129035/125‑1 (MQ=255) ggTGCGTCTGCTGCGTGCGATTCAGGAGCGCGAAGTTCAGCGTGTCGGTAGCAACC‑AGATTATCTCGGTTGATGTCCGGCTGATTGCGGCGACCCGTCGCGATCTTGCCGCAGAGGTGAATGCCggg > 3:124679/1‑127 (MQ=255) ggTGCGTCTGCTGCGTGCGATTCAGGAGCGCGAAGTTCAGCGTGTCGGTAGCAACC‑AGATTATCTCGGTTGATGTCCGGCTGATTGCGGCGACCCGTCGCGATCTTGCCGCAGAGGTGAATGCCggg < 4:124679/127‑1 (MQ=255) gcgcGAAGTTCAGCGTGTCGGTAGCAACC‑AGATTATCTCGGTTGATGTCCGGCTGATTGCGGCGACCCATCGCGATCTTGCCGCAGAGGTGAATGCCGGGCGTTTTCGCCAGGATCTCTACTATCGCCTGAATGTGGTGGCGATTGAAg > 5:166528/1‑149 (MQ=255) tAGCAACC‑AGATTATCTCGGTTGATGTCCGGCTGATTGCGGCGACCCATCGCGATCTTGCCGCAGAGGTGAATGCCGGGCGTTTTCGCCAGGATCTCTACTATCGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAAcggc < 6:166528/149‑1 (MQ=255) | ATAAACGCCGGGAAGGGCGCTTTGTTGAGGCGGACGGCGGCACGCTGTTTCTCGATGAAATTGGCGATATCTCGCCGATGATGCAGGTGCGTCTGCTGCGTGCGATTCAGGAGCGCGAAGTTCAGCGTGTCGGTAGCAACC‑AGATTATCTCGGTTGATGTCCGGCTGATTGCGGCGACCCATCGCGATCTTGCCGCAGAGGTGAATGCCGGGCGTTTTCGCCAGGATCTCTACTATCGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGC > NC_000913/4203984‑4204265 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |