| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 2,653,367 | G→A | 36.8% | intergenic (+28/‑488) | sseA → / → ryfA | 3‑mercaptopyruvate sulfurtransferase/novel sRNA, function unknown |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,653,367 | 0 | G | A | 36.8% | 10.1 / 17.5 | 19 | intergenic (+28/‑488) | sseA/ryfA | 3‑mercaptopyruvate sulfurtransferase/novel sRNA, function unknown |
| Reads supporting (aligned to +/- strand): ref base G (6/6); new base A (4/3); total (10/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-03 | |||||||||||
TTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCA > NC_000913/2653233‑2653489 | ttttGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGtt > 4:832107/1‑147 (MQ=255)ttttGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGtt < 3:832107/147‑1 (MQ=255) ccAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGaa > 2:806157/1‑108 (MQ=255) ccAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGaa < 1:806157/108‑1 (MQ=255) aaCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTg > 3:599724/1‑109 (MQ=255) aaCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGAAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTg < 4:599724/109‑1 (MQ=255) gAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGtt < 2:775596/85‑1 (MQ=255) gAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGtt > 1:775596/1‑85 (MQ=255) tttACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAAttt > 4:206616/1‑58 (MQ=255) tttACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAAttt < 3:206616/58‑1 (MQ=255) aTAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAg < 4:575511/90‑1 (MQ=255) aTAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCCGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAg > 3:575511/1‑90 (MQ=255) aagcgccaAATTCGGCGGTTTTTTATTGCCAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTca < 3:831782/89‑1 (MQ=37) aaCCGCCAAATTCGGCGGTTTTTTATTGCCAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTca > 4:831782/1‑92 (MQ=37) aCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCa < 1:947307/82‑1 (MQ=255) aCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCa > 2:947307/1‑82 (MQ=255) aCCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCtttt > 1:188210/1‑124 (MQ=255) aCCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCtttt < 2:188210/124‑1 (MQ=255) aCCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCa < 1:238520/129‑1 (MQ=255) aCCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCa > 2:238520/1‑129 (MQ=255) ccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTAc > 2:202752/1‑66 (MQ=255) ccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTAc < 1:202752/66‑1 (MQ=255) ccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCaaa > 1:751229/1‑94 (MQ=255) ccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCaaa < 2:751229/94‑1 (MQ=255) ccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCtttt > 1:175802/1‑120 (MQ=255) ccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCtttt < 2:175802/120‑1 (MQ=255) | TTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCA > NC_000913/2653233‑2653489 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |