Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 2,066,253 | A→G | intergenic (‑489/+52) | cobU ← / ← insH‑6 | cobinamide‑P guanylyltransferase/cobinamide kinase/CP4‑44 prophage; IS5 transposase and trans‑activator |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,066,253 | 0 | A | G | 94.7% | 36.5 / ‑4.1 | 19 | intergenic (‑489/+52) | cobU/insH‑6 | cobinamide‑P guanylyltransferase/cobinamide kinase/CP4‑44 prophage; IS5 transposase and trans‑activator |
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (10/8); total (10/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.74e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01 |
GAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAG > NC_000913/2066126‑2066301 | gAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGGATCAg < 1:527703/133‑1 (MQ=255) aaTTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGGATCAGc < 2:567141/133‑1 (MQ=255) ttATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGGATCAGCCt > 2:569746/1‑133 (MQ=255) cGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGGATCAGCCTATTTa > 2:15842/1‑133 (MQ=255) gtggCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGGATCAGCCTATTTAGACCGtt > 1:337627/1‑133 (MQ=255) aCTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGGATCAGCCTATTTAGACCGTTTCTTCg < 2:337627/133‑1 (MQ=255) aattAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGGATCAGCCTATTTAGACCGTTTCTTAGCCAttt > 2:407127/4‑133 (MQ=0) gggAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGGATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGc > 1:429504/1‑133 (MQ=17) gggAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGGATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGc > 1:558747/1‑133 (MQ=18) ggAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTAATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGGATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCg > 2:308076/1‑133 (MQ=18) aaGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGGATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGtt > 1:494019/1‑133 (MQ=17) gCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGGATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATccc < 1:83432/133‑1 (MQ=17) cGAATAAGTGGGGAACTTATTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATcccc < 1:295300/133‑1 (MQ=17) cGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGGATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATcccc < 1:387908/133‑1 (MQ=17) cGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGGATCAGCCTATTTAGACCGTTTCTTAGCCATTTAAGGCGTTATcccc < 1:407127/133‑1 (MQ=17) gAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGGATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCa < 2:558747/133‑1 (MQ=17) gAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGGATCAGCCTATTTAAACCGTTTCTTCGCCATTTAAGGCGTTATcc > 1:43605/1‑130 (MQ=17) gAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGGATCAGCCTATTTAAACCGTTTCTTCGCCATTTAAGGCGTTATcc < 2:43605/130‑1 (MQ=12) aaTAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGGATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAg > 2:24874/1‑133 (MQ=17) | GAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAG > NC_000913/2066126‑2066301 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |