| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,121,432 | 0 | T | G | 52.9% | 9.8 / 18.1 | 17 | D19A (GAC→GCC) | dinI | DNA damage‑inducible protein I |
| Reads supporting (aligned to +/- strand): ref base T (1/7); new base G (9/0); total (10/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.11e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-03 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGG‑AAAGTTCGC‑CAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCG‑TACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCG > NC_000913/1121194‑1121669 | agtaaAATAAAAAAGCCGGGGCGACCCGTCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTGATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGG‑AAAGTTCGC‑CAGACAGGGCGTCAATAGCCCCAg < 2:235918/251‑1 (MQ=255) ggCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGG‑GAAGTTCGC‑CAGCCAGGGCGGCAATAACCCCCGCCGGCAATGGAGGAGATTTTGCTATgg > 1:77717/1‑251 (MQ=255) cACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGG‑AAAGGTCGC‑CAGCCCGGGCGGCAAAAAACCCAGCCGGGCATGGAGAAGGTTTTGCTATGGGGACTTCCAttt > 1:40329/1‑250 (MQ=255) cTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGGAAAGT‑CGCGCGACCAGGG‑GGCAATAGACCCCGCTGGCAATGGGGAAGGTTTTGCTTTgggggagtcacctccc > 1:237564/1‑237 (MQ=255) gCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGG‑AAAGTTTGC‑CCGCCCGAGGGTCAATATCCCCACCTGGCGATTGAGAAGTTTTTTCCATTGTGACTTCAACTCGCataataatccccc > 1:145302/1‑244 (MQ=255) tGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGG‑GAAGTTTGC‑CAGCCAAGGGGGCAATAACCCCCGCTGGCAATGGAGAAGATTTCTCTATTGTGACTTCAATTTGCATAATAGCCCCCtct > 1:170243/1‑249 (MQ=255) gACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGG‑AAAGTTCGC‑CAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGt < 2:40329/250‑1 (MQ=255) cATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAAGTACTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGG‑AAAGTTCGC‑CAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTgg < 1:325938/251‑1 (MQ=255) aCGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAAAACGGCGGG‑AAAGTTCGC‑CAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTAAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTa < 2:257333/251‑1 (MQ=255) tGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGG‑AAAGATTGC‑CAGCCCGGGCGGCCATAACCCCCGCTGGCCATGGGGAAGGTTTTTCTATTGGGGCTTCCCTTCGCATAATAACCCCCCCTTGGATATACTGGTTTTTTATACAGGTAAAAAAAc > 1:19527/1‑250 (MQ=255) ttcttTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGG‑GAAGATCGC‑CAGCCAGGGGGGCAATAACCACCCCCGGGAAGGGAGAAGGTTTTTCTCTGGTGGGTTCAATTTGCATAATAGCCCCCTGTTTAATATACTGGTTATTTATAAAGGTAAAATAAACTAATTAAAAAAGGAAGATACGATTTTTAt > 2:41244/1‑251 (MQ=255) tGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGG‑AAAGTTCGC‑CAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCG‑TACCg < 1:247973/215‑1 (MQ=255) aTCAGGAAACGCATACTGAATACGGCGGG‑AAAGGTTGC‑CAGCCCGGGGGGCAAAAGCCCCCGCCGGCAAAGGGGAAGGTTTTGCTATGGTGGCTTCAATTCCCATAATAAACCCCCGTTTAATTTAATGGTTTTTTTTAAAAGAAAAATAACCTAATGAAAACAGGAAGCTAAGATTTTTTTTGTTTTAAGGACCAGCGGAACCGTTTCCCCGCCG‑GGAATGGGACCCGCGTGtcctccatccgttgggtg > 1:261568/1‑235 (MQ=255) tCAGGAAACGCATACTGAATACGGCGGG‑GAAAATCGC‑CACCCAGGGCGGCAATAGCCCCAGCTGGGAATTGAGAAGTTTTCGCTATGGTGACTTCAATTCGCCTAATAGCCCCCCTTTTTATATACTGGTTTTTTATACAAGGAAAATAACCTAATTACAAAAGGAAGCTACGATTTTTTTTTTTTTACGGACCAGCG‑GACCCTTTTCCCCGCGAGGAATGGCACCACCGTGTCATCCGTCCGCGCGATGc > 1:3670/1‑251 (MQ=255) cAGGAAACGCATACTGAATACGGCGGG‑AAAGTTTGC‑CCACCAGGGCGGCAATAAACCCACCTGGCAATGGAGAAGTTTTCGCTATGGGGACTTCAATTCGCATAATAACCCCCTTTTGAATATACTGGTTTTTTTTAAAGGTCAAATAACCTCATGACAAAAGGCAGCTACGATTTTTTTTTTTTTACGGACCACCG‑GACCGTTTCCCCCGCGGGGAATGGGCCCACCGGGTCATCAGTCAGTGCGATGCt > 1:325931/1‑251 (MQ=255) cGC‑CAGCCAGGGCGTCAATAGCCCCAGCTTGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCG‑TACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGtt < 2:203038/250‑1 (MQ=255) gc‑cagccaGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCG‑TACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCg < 2:19527/251‑1 (MQ=255) | AGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGG‑AAAGTTCGC‑CAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCG‑TACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCG > NC_000913/1121194‑1121669 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |