Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 3,214,968 | (T)6→7 | coding (506/507 nt) | mug ← | G/U mismatch‑specific DNA glycosylase, xanthine DNA glycosylase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,214,962 | 1 | . | T | 95.8% | 77.3 / ‑3.2 | 24 | intergenic (+74/+5) | rpoD/mug | RNA polymerase, sigma 70 (sigma D) factor/G/U mismatch‑specific DNA glycosylase, xanthine DNA glycosylase |
Reads supporting (aligned to +/- strand): ref base . (0/1); new base T (11/12); total (11/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01 |
ACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCG‑TT‑TTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGA > NC_000913/3214837‑3215091 | aCTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTACTCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTTGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTT‑TTTTATCGcc < 2:706085/139‑1 (MQ=255) aCTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTT‑TTTTATCGcc < 1:177351/139‑1 (MQ=255) aCTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTT‑TTTTATCGcc < 1:949917/139‑1 (MQ=255) tGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTT‑TTTTATCGCCCa < 1:1028900/139‑1 (MQ=255) gAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGTCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTT‑TTTTATCGCCCACGCACTACCa < 1:412592/139‑1 (MQ=255) gAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTT‑TTTTATCGCCCACGCACTACCa > 2:677384/1‑139 (MQ=255) cGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTT‑TTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTcgcg > 1:452913/1‑139 (MQ=255) gCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTT‑TTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATAc < 1:103644/139‑1 (MQ=255) gACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTT‑TTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAAc > 2:701491/1‑139 (MQ=255) gATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTT‑TTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAg > 1:1026802/1‑139 (MQ=255) gCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTT‑TTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTa > 2:394124/1‑139 (MQ=255) tACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTT‑TTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGt < 2:132794/113‑1 (MQ=255) tACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTT‑TTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGt > 1:132794/1‑113 (MQ=255) cgccgcACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTT‑TTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGc < 2:127244/139‑1 (MQ=255) cgccgcACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTT‑TTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGc < 2:452913/139‑1 (MQ=255) cgccgcACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTT‑TTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGc < 1:251143/139‑1 (MQ=255) cACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTT‑TTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGa > 2:17767/1‑139 (MQ=255) aCCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTT‑TTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGAt < 2:1055068/139‑1 (MQ=255) ggCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTT‑TTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTgg > 2:643006/1‑139 (MQ=255) caATCGCCACCTTTTCGGTGGCG‑TTATTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGt < 2:717234/139‑1 (MQ=255) cGCCACCTTTTCGGTGGCGTTT‑TTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGaa < 2:334108/139‑1 (MQ=255) cGCCACCTTTTCGGTGGCGTTT‑TTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGaa > 2:1023265/1‑139 (MQ=255) ccACCTTTTCGGTGGCGTTT‑TTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAAcc > 2:199104/1‑139 (MQ=255) tCGGTGGCGTTT‑TTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGa > 2:566191/1‑139 (MQ=255) | ACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCG‑TT‑TTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGA > NC_000913/3214837‑3215091 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |