Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 3,510,121 | +G | intergenic (‑358/+151) | yhfZ ← / ← trpS | hypothetical protein/tryptophan tRNA synthetase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,510,120 | 1 | . | G | 97.9% | 120.4 / ‑3.2 | 48 | intergenic (‑357/+152) | yhfZ/trpS | hypothetical protein/tryptophan tRNA synthetase |
Reads supporting (aligned to +/- strand): ref base . (1/0); new base G (30/17); total (31/17) | |||||||||||
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 |
CGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCC‑GTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAAT > NC_000913/3510045‑3510190 | cGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGCTTcccc > 1:2172209/1‑76 (MQ=255) gaTGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCC‑ggt > 2:1268740/1‑74 (MQ=255) gtgGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCC‑Ggtttc > 2:4540380/1‑74 (MQ=255) gAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCACATGAGGCTACAAAAACCGGGTTTCCCCGGTTTCaga > 1:1067084/1‑76 (MQ=255) gAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCaga > 1:1029524/1‑76 (MQ=255) aaaCGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCagag < 1:1879432/76‑1 (MQ=255) aCGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGCTTCCCCGGGTTCAGAGtg > 2:546900/1‑76 (MQ=255) cGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGtga > 2:1048062/1‑76 (MQ=255) ggCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGtgat > 2:766509/1‑76 (MQ=255) gCGGGGCGCTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGtgatg < 2:583485/76‑1 (MQ=255) gCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGtgatg > 2:4093752/1‑76 (MQ=255) gAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTgatgat > 2:4640633/1‑76 (MQ=255) gAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTgatgat > 2:4068618/1‑76 (MQ=255) aGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATa > 2:3732385/1‑76 (MQ=255) ggCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATaa > 1:4434097/1‑76 (MQ=255) cACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATaaaa < 1:4652511/76‑1 (MQ=255) cTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGc > 1:1351460/1‑76 (MQ=255) cTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGc > 1:4613175/1‑76 (MQ=255) caTTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAAt < 1:701873/76‑1 (MQ=255) aTTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAAtt > 2:333318/1‑76 (MQ=255) ttttCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGcc < 1:4683714/76‑1 (MQ=255) tCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGa < 1:657590/76‑1 (MQ=255) cGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGAt > 1:580697/1‑76 (MQ=255) attattGCAAATGAGGATGAAAAAACCGGGTTGCCCCGGTTGCAGAGTGATGATAAAAGCAAAATTGCCTGATgcg < 2:3802728/76‑1 (MQ=37) attattGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATgcg > 2:2054224/1‑76 (MQ=255) ttattGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATgcgc > 1:4118409/1‑76 (MQ=255) tattGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCt > 1:4265171/1‑76 (MQ=255) tGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACg > 1:2837180/1‑76 (MQ=255) tcAAATGAGGATAAAACAGCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGc < 2:3630702/75‑1 (MQ=37) cAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCt > 2:919152/1‑76 (MQ=255) aaaTGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCtt > 2:2577049/1‑76 (MQ=255) aaaTGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCtt > 1:4738690/1‑76 (MQ=255) aaTGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTa > 1:3313981/1‑76 (MQ=255) aTGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTAt < 1:983510/76‑1 (MQ=255) gAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCa < 1:4357753/76‑1 (MQ=255) aGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAg < 1:39656/76‑1 (MQ=255) tgATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAgg > 2:3103627/2‑76 (MQ=255) ggATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAgg > 1:1747757/1‑76 (MQ=255) aTAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGcc > 1:4000050/1‑76 (MQ=255) tAAAAAAACAGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCGTATCAGGCCt > 2:3040015/1‑76 (MQ=25) aaaaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTAc < 1:4487176/76‑1 (MQ=255) aaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACAtt < 2:3572964/76‑1 (MQ=255) aaaCCCGGTTTCCCCGGGTTCAGGGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACAtt < 2:314274/76‑1 (MQ=18) aaCCGTGTTTCCCCGGTTTCAGAGTGATTATAAAAGCACCATTGCCTGATGCGCTACTCTTCTCAGGCCTACAttt < 1:546900/76‑1 (MQ=18) aCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTc > 2:1137967/1‑76 (MQ=255) aCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTc < 1:4338653/76‑1 (MQ=255) ccGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTcc > 1:1109351/1‑76 (MQ=255) ggTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTg > 2:2400616/1‑76 (MQ=255) tCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAat < 2:513309/76‑1 (MQ=255) tCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAat > 1:1135835/1‑76 (MQ=255) | CGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCC‑GTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAAT > NC_000913/3510045‑3510190 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |