| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,101,514 | 0 | G | T | 25.7% | 87.8 / 10.4 | 35 | intergenic (+84/‑176) | sodA/kdgT | superoxide dismutase (Mn)/2‑dehydro‑3‑deoxy‑D‑gluconate:H(+) symporter |
| Reads supporting (aligned to +/- strand): ref base G (9/17); new base T (9/0); total (18/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.03e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATT > NC_000913/4101384‑4101650 | tggtgAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTGTGATCAAtt > 2:260985/1‑140 (MQ=255) ggtgAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTGTGATCAAttt > 1:372689/1‑140 (MQ=255) gtgAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTc > 1:417909/1‑140 (MQ=255) gtgAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAAtttt > 1:261793/1‑139 (MQ=255) aaCTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTTGTTTTTTTGGGCTTAAGTTCaaa > 2:155312/1‑140 (MQ=255) aaGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTTTGATCAAATTCAaaaaaaaaaaaa > 2:425116/1‑140 (MQ=255) agcagcGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAaaaaaaaaaaaaa > 2:149820/1‑139 (MQ=255) cACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGat < 1:260985/140‑1 (MQ=255) cACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGat < 1:425116/140‑1 (MQ=255) tGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTTTGATCAATTTTAAAATAAAAAAAAAGATCCGataaaaaaa > 2:542929/1‑138 (MQ=255) gcgAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAAc < 2:417909/140‑1 (MQ=255) gcgAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAAc < 1:495861/140‑1 (MQ=255) aaaaTAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGt < 1:149820/140‑1 (MQ=255) tGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAAATTTAAAAAAAAAAAAATAATCCGATAAAAAAAAAAAAACGTTTCAATTGAtgt > 2:382906/1‑140 (MQ=255) tGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTGTGATCAATTTTAAAAAAAAAAAAATGATCCGATAAAAAAAAAACAACGTTTCAATTGAtgt > 1:469992/1‑140 (MQ=255) ctgcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGAt < 1:382906/140‑1 (MQ=255) ctgcAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATTAATTTTAAAATAAAAAAAAAGATCCGATAAAAAAAAAAAAACGTTTTAATTGATGTGGGTTTTAt > 2:501042/1‑140 (MQ=255) ctgcAATGAGGCGTATAGGCCACATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGAt < 2:424786/140‑1 (MQ=255) aGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTCAAAATAAAAAAAATGATCCGATAAAAATAAAAAAGCGTTTTAATTTATGTGGGTTTGATCACttttt > 2:141868/1‑139 (MQ=255) ggCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTa > 2:469389/1‑139 (MQ=255) ggCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTa < 1:469389/139‑1 (MQ=255) cGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTTTGATCAATTTTAAAAAAAAAAAAAAAATCCGAAAAAAAAAAAAAAACGGTTTAATTGATTTGGGTTTTATCACTttttttt > 2:282900/1‑139 (MQ=255) cGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTTAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTAttt > 1:552399/1‑139 (MQ=255) gCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTGTGATCAAATTTAAAATAAAAACAAAGATCCCATAAAAAAAAAAAAGCGTTTCAATTTATGTGGGTTTTATCACTTTTTTTGATTAaatgaa > 2:244197/1‑140 (MQ=255) tCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCt < 1:282900/140‑1 (MQ=255) aaaaaaTGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAAAAAAAAAAAAGATCCGATAAAAAAAAAAAAGCGTTTTAATTGATGTGGGTTTTATTACTTTTTTTTATTAAATGAAATTCTATCtttt > 2:501192/1‑138 (MQ=255) aaaaTGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGtt < 2:86816/140‑1 (MQ=255) aaTGAACCATCGCCAACGGCGGTGGTTTTTTTTTGATCAATTTCAAAATAAAAAAAATGATCCGATAAAAATAAAAAAGCGTTTTAATTGATGTGGTTTTTATCACTTTTTTTTATTAAATGAATGTCTATCTTCGtttt > 2:196590/1‑139 (MQ=255) tGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCa > 1:374249/1‑140 (MQ=255) aaCCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATc < 1:501042/140‑1 (MQ=255) ccATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCaa < 2:469992/140‑1 (MQ=255) cggcggTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCt < 2:254756/140‑1 (MQ=255) cggcggTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCt < 2:447033/140‑1 (MQ=255) cggcggTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCt < 1:600598/140‑1 (MQ=255) gcggTGGGTTTTTTTTGATCAATTTCAAAAAAAAAAAAATGATCCGATAAAAATAAAAAAGCGTTTCAATTGATGTGGTTTTGATCACTTTTTTTTATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTcc > 2:370406/1‑140 (MQ=255) cggTGGGTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCa > 1:629468/1‑140 (MQ=255) tttttttGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAgg < 1:645722/140‑1 (MQ=255) tttGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAAtt < 1:196590/140‑1 (MQ=255) | TGGTGAACTGGGACGAAGCAGCGGCACGTTTTGCGGCGAAAAAATAATCATTTGCCGCCTGCTGCAATGAGGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATT > NC_000913/4101384‑4101650 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |