Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 2,620,218 | T→A | 16.5% | intergenic (‑58/+28) | uraA ← / ← upp | uracil permease/uracil phosphoribosyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,620,218 | 0 | T | A | 16.5% | 98.6 / 11.2 | 36 | intergenic (‑58/+28) | uraA/upp | uracil permease/uracil phosphoribosyltransferase |
Reads supporting (aligned to +/- strand): ref base T (14/16); new base A (4/2); total (18/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.58e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01 |
ACCAAACATGGCGAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAA‑AGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCGGGTGCGCT > NC_000913/2620071‑2620360 | aCCAAACATGGCGAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTa > 1:146782/1‑149 (MQ=255) tGGCGAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTTTAAG‑AGTCGGc < 4:323834/149‑1 (MQ=255) tGGCGAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAA‑AGTCGGc < 2:146782/149‑1 (MQ=255) gAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAA‑AGTCGGCTTTa > 1:178445/1‑149 (MQ=255) ttGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAA‑AGTCGGCTTTAATTAttttt > 1:262677/1‑149 (MQ=255) cGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAA‑AGTCGGCTTTAATTATTTTTATTCTTTATTTc < 1:256789/149‑1 (MQ=255) gggATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAA‑AGTCGGCTTTAATTATTTTTATTCTTTATTTCg > 4:124670/1‑149 (MQ=255) gggATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAA‑AGTCGGCTTTAATTATTTTTATTCTTTATTTCg < 2:262677/149‑1 (MQ=255) tCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAA‑AGTCGGCTTTAATTATTTTTATTCTTTATTTCGTa < 4:201937/144‑1 (MQ=255) tCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAA‑AGTCGGCTTTAATTATTTTTATTCTTTATTTCGTa > 3:201937/1‑144 (MQ=255) tCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAA‑AGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCaaa < 4:615373/149‑1 (MQ=255) tAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAA‑AGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGAtt > 4:49144/1‑149 (MQ=255) ggCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAA‑AGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCAc < 3:49144/149‑1 (MQ=255) cGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAA‑AGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCg > 2:131647/1‑149 (MQ=255) ctcacCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGGCTTTAAG‑AGTCGGCtttttt > 1:402596/1‑96 (MQ=255) ctcacCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGGCTTTAAG‑AGTCGGCtttttt < 2:402596/99‑4 (MQ=255) acCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAA‑AGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGc > 1:514475/1‑149 (MQ=255) ccccGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAA‑AGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGcc < 1:131647/149‑1 (MQ=255) aCGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAA‑AGTCGGCTTTAATTATTTTTATTCTTTATTTCGt < 2:232255/106‑1 (MQ=255) aCGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAA‑AGTCGGCTTTAATTATTTTTATTCTTTATTTCGt > 1:232255/1‑106 (MQ=255) aCGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAA‑AGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATg < 3:124670/149‑1 (MQ=255) tctGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAA‑AGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCt < 1:341735/149‑1 (MQ=255) atGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAA‑AGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATc < 2:178445/149‑1 (MQ=255) gtTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAA‑AGTCGGCGTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGTCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATg > 1:289770/1‑149 (MQ=255) ttATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAA‑AGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGtt < 3:447175/128‑1 (MQ=255) ttATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAA‑AGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGtt > 4:447175/1‑128 (MQ=255) cgcTTTACTCAAAAAAAAGCCGACTCTTAA‑AGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCgg < 3:461969/144‑1 (MQ=255) cgcTTTACTCAAAAAAAAGCCGACTCTTAA‑AGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCgg > 4:461969/1‑144 (MQ=255) tCAAAAAAAAGCCGACTCTTAA‑AGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGa > 3:380218/1‑123 (MQ=255) tCAAAAAAAAGCCGACTCTTAA‑AGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGa < 4:380218/123‑1 (MQ=255) aaaGCCGACTTTAAG‑AGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCggg > 2:30294/1‑149 (MQ=255) aaGCCGACTTTAAG‑AGTCGGCTTTAATTCTTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCGGGt > 1:339686/1‑149 (MQ=255) aaGCCGACT‑TTAAGAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCGGGt < 1:30294/149‑1 (MQ=255) gCCGACTTTAAG‑AGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCGGGTgc > 3:465424/1‑149 (MQ=255) cGACTCTTAA‑AGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCGGGTgcgc < 2:514475/149‑1 (MQ=255) gACTCTTAA‑AGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCGGGTGCGCt > 4:286108/1‑149 (MQ=255) | ACCAAACATGGCGAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAA‑AGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCGGGTGCGCT > NC_000913/2620071‑2620360 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |