Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 4,265,670 | T→G | 5.9% | F453V (TTT→GTT) | dnaB → | replicative DNA helicase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,265,670 | 0 | T | G | 5.9% | 245.0 / 7.0 | 68 | F453V (TTT→GTT) | dnaB | replicative DNA helicase |
Reads supporting (aligned to +/- strand): ref base T (31/33); new base G (2/2); total (33/35) | |||||||||||
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.97e-01 |
ATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGT > NC_000913/4265535‑4265798 | aTCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAAt > 1:512678/1‑148 (MQ=255) cGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATgg < 3:27731/148‑1 (MQ=255) aGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTcgcg < 4:880957/148‑1 (MQ=255) gATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCt < 2:512678/148‑1 (MQ=255) gCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAAc > 2:933710/1‑129 (MQ=255) gCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAAc < 1:933710/129‑1 (MQ=255) tCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTcg < 1:97090/131‑1 (MQ=255) tCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTcg > 2:97090/1‑131 (MQ=255) tCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCGTTAAAGGTc > 4:143056/1‑123 (MQ=255) tCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCGTTAAAGGTc < 3:143056/123‑1 (MQ=255) ttATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGcc > 3:523200/1‑148 (MQ=255) atctatcGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGc < 4:523200/148‑1 (MQ=255) tctatcGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCa > 1:735267/1‑148 (MQ=255) cGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGt > 4:77536/1‑115 (MQ=255) cGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTacg > 2:122153/1‑148 (MQ=255) cGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTCAACGGTCAATGGt < 3:77536/115‑1 (MQ=255) atgaGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTacgacga > 3:543286/1‑148 (MQ=255) tgaGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTacgacgac > 4:724180/1‑148 (MQ=255) gaGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTacgacgacg < 2:735267/148‑1 (MQ=255) ggTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGa < 1:1020820/117‑1 (MQ=255) ggTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGa > 2:1020820/1‑117 (MQ=255) tgtATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAGACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCg > 2:328767/1‑125 (MQ=255) tgtATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTTACCTTTAACGGTCAATGGTCGTGCTTCGACAACTATGCg < 1:328767/125‑1 (MQ=255) tgtATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGaata < 3:724180/148‑1 (MQ=255) cACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGaataataat < 1:122153/148‑1 (MQ=255) cGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATGGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGaataa > 4:663857/1‑142 (MQ=255) cGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGaataa < 3:663857/142‑1 (MQ=255) aaaaCAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAattatt < 4:543286/148‑1 (MQ=255) gTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTAcgacgacga < 4:901875/130‑1 (MQ=255) gTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTAcgacgacga > 3:901875/1‑130 (MQ=255) tGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCgg > 3:87481/1‑110 (MQ=255) tGATTTAAAAGGCATCGCAGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCgg < 4:87481/110‑1 (MQ=255) gATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGTGACGGTACGCCTGACCTTTAACgg > 2:1110540/1‑80 (MQ=255) gATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACgg < 1:1110540/80‑1 (MQ=255) ttAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATg < 1:241780/103‑1 (MQ=255) ttAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATg > 2:241780/1‑103 (MQ=255) aaGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTcgcg < 4:828693/85‑1 (MQ=255) aaGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTcgcg > 3:828693/1‑85 (MQ=255) cATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCa > 1:516201/1‑72 (MQ=255) cATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCa < 2:516201/72‑1 (MQ=255) cgcgGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCGTTaaag < 4:741837/65‑3 (MQ=255) cgcgGAAATTATTATCGGTAAACAACGTAACGGCCCAATCAGGACGGTACGCCTGACCGTTaaag > 3:741837/1‑63 (MQ=255) ggAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAatt > 3:1114503/1‑120 (MQ=255) ggAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAatt < 4:1114503/120‑1 (MQ=255) tatCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGc < 1:157016/89‑1 (MQ=255) tatCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGc > 2:157016/1‑89 (MQ=255) tAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGc > 1:4831/1‑75 (MQ=255) tAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGc < 2:4831/75‑1 (MQ=255) tAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTa < 1:288310/111‑1 (MQ=255) tAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTa > 2:288310/1‑111 (MQ=255) acaacGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCggg < 3:639589/75‑1 (MQ=255) acaacGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCggg > 4:639589/1‑75 (MQ=255) acaacGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGc < 3:860806/77‑1 (MQ=255) acaacGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGc > 4:860806/1‑77 (MQ=255) tAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCa > 1:189644/1‑148 (MQ=255) aaCGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCgg < 1:1020403/67‑1 (MQ=255) aaCGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCgg > 2:1020403/1‑67 (MQ=255) gCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACtt > 4:941414/1‑124 (MQ=255) gCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACtt < 3:941414/124‑1 (MQ=255) ccAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGc > 4:1095790/1‑67 (MQ=255) ccAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGc < 3:1095790/67‑1 (MQ=255) cAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCaa > 3:681558/1‑148 (MQ=255) cGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAac < 4:627615/127‑1 (MQ=255) cGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAac > 3:627615/1‑127 (MQ=255) cGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACtgt < 4:681558/148‑1 (MQ=255) cGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACtgt < 2:189644/148‑1 (MQ=255) gACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCaa < 1:551722/96‑1 (MQ=255) gACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCaa > 2:551722/1‑96 (MQ=255) | ATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGT > NC_000913/4265535‑4265798 |
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: ‑ |