| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 4,619,393 | C→T | L389L (CTG→TTG) | deoA → | thymidine phosphorylase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,619,393 | 0 | C | T | 95.5% | 58.5 / ‑6.5 | 22 | L389L (CTG→TTG) | deoA | thymidine phosphorylase |
| Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (9/12); minor base A (1/0); total (10/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.55e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01 | |||||||||||
AAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTC > NC_000913/4619160‑4619638 | aaGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTGCGGTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGAGACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGAc < 2:237023/251‑1 (MQ=255) atGGTCGTATGGTAGCGGCACAAAAAGGCCCAGCCGTCGTCGTTGAGAACTACGCGAAGGATCGGCCGACAGCGATGCTGACTAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCAGCTGACACCAGTGATGACAGCGTCGGCTTGTCTGATAAGGCGCGTTTGGGCGACCAGGTAGACGGTc < 2:249418/248‑1 (MQ=255) tGGTCGTAAGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGCACTACGCGAAGTATCTGCCGCCAGCGATGCTGACGAACGCAGTCTATGCTGCAACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGAAATAGGCGGCGGTCGTCGTCCGACATCTGACACCATTGATTACAGCGTCGGCTTTACTGCTATGGCCCGTTTGGGCGACCAGATAGACggac > 1:172634/1‑246 (MQ=255) ggCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGaa < 1:133100/251‑1 (MQ=255) cGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGGAaacagct > 1:69913/1‑248 (MQ=255) gATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCCTTGATTACAGCGTCGGCTTTACTGATATGGCGCGGTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAATGCGAAAACAGCCGGCAGGAAGGGGCGaaaaggggggaagcggggc > 2:135155/1‑235 (MQ=255) cGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGGCGAAAACAGCTGGCAGGAAGCGGCGAAAACGGTGAAAGCGGGAATTaaaatt > 2:70286/1‑248 (MQ=255) tGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGAAGAAAAAAGCTGGCAGGAAGCGGGGAAAGCGGTGGAAGCGGCAATTAAACTTGCCGATAAAGAACCgg > 1:177147/1‑251 (MQ=255) gAAATGGATTCCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCg < 2:83814/251‑1 (MQ=255) ccccGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAg < 1:1020/250‑1 (MQ=255) cgggcTGGGGATGGCAGTGGTTGCAATGGGCAGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGTCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGATGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCTGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGGCTAGCGCCGTATCAGCGAAt < 2:51015/248‑1 (MQ=255) gcgcTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGGCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa < 1:132595/251‑1 (MQ=255) tGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATCTGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGAGAAAGACGAAAACAGCTGGCAGGGAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACgg > 1:157100/1‑251 (MQ=255) ggggATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGGTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACgg < 1:182724/250‑1 (MQ=255) cAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCCCGCGAAAGACg > 1:81677/1‑135 (MQ=255) cAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACg < 2:81677/135‑1 (MQ=255) gCAATGGGCGGCTGTCGTCGTCAGGCATCTGACACCATTGAGTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGGCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCa < 2:146923/250‑1 (MQ=255) cGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGAAGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGGGAAAACGGCAATTAAAATTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCCTTTAAATTTGAAGCACTgg > 1:143201/1‑250 (MQ=255) tCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGGAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAACACCCGAAAGCACACCAACTGTCTATCGCCGTATCAACCAATAAAGGTATACTGATCTGATCATTTAAATTTTAAGCACTGAgg > 1:24114/1‑250 (MQ=255) acacCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGTTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACat < 1:16774/251‑1 (MQ=255) gttACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGTGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACAAACTGTCTATCGCCGTATCAGCGAATAACGGTATACCGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAACCGTGCATTTATTATGGt < 1:60269/250‑1 (MQ=255) cgcgTATGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACCCGAAAGACGAAAACCACTGGCAGGAAGCGGCGAAAGCGGGGAAAGCGGCAATTAAAATTGCCGATAAAGCACCGGAAAGCACACCAACCGTCTCTCGCCGTATCAGCGAAAAACGGTAAACTGATCTGATCCTTTAAATTTTAAGCACTGAGGACGGGGGACATATGAAACGGGCATTTCTTATGGGGCTGGGATCATTc > 2:115239/1‑251 (MQ=255) | AAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTC > NC_000913/4619160‑4619638 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |