Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 838,010:1 | +C | 13.8% | intergenic (‑43/+37) | pth ← / ← rplY | peptidyl‑tRNA hydrolase/50S ribosomal protein L25/general stress protein Ctc |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 838,010 | 1 | . | C | 13.8% | 100.8 / 12.4 | 29 | intergenic (‑43/+37) | pth/rplY | peptidyl‑tRNA hydrolase/50S ribosomal protein L25/general stress protein Ctc |
Reads supporting (aligned to +/- strand): ref base . (11/14); new base C (2/2); total (13/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.94e-01 |
GTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGC‑GGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGCCT > NC_002947/837877‑838146 | gTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGC‑GGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCG‑GCACTGTAACCCGATc > 1:587123/1‑150 (MQ=255) gTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGC‑GGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCG‑GCACTGTAACCCGATc < 2:587123/150‑1 (MQ=255) aaCGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGC‑GGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTATAGTGCGCGCG‑GCACTTTAACCCGATCAACg > 2:291099/1‑150 (MQ=255) tgatATGCCGGTTCTTTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGC‑GGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGAt < 1:200542/147‑1 (MQ=255) gCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGC‑GGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCg > 2:690427/1‑150 (MQ=255) gggTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCCACGATAAACAGGCCGGCAGGCCCTGGCA‑GAGATCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTACGCGGGGCTTACTGGGCa > 2:442250/1‑150 (MQ=255) gggTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGC‑GGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCa > 1:151732/1‑150 (MQ=255) gTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGC‑GGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAgc < 2:146639/150‑1 (MQ=255) ggCCGGGGTTACCCAGGCCAACGATCAACTGGATGGC‑GGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGa > 1:219060/1‑150 (MQ=255) ggggTTACCCAGGCCAACGATCAACTGGATGGC‑GGTCACGTCAGGGGCTCTCCCTTGGAGTTGCTGGGTTACAGTGCGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCaa > 1:382782/1‑150 (MQ=255) cgccAGGCCAACGATCAACTGGATGGC‑GGGCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTg < 2:232515/148‑1 (MQ=255) aCCCAGGCCAACTATCAACGCGATGGC‑GGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTg < 2:606768/150‑1 (MQ=255) aCCCAGGCCAACGATCAACTGGATGGC‑GGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTg < 1:690427/150‑1 (MQ=255) ccAGGCCAACGATCAACTGGATGGC‑GGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGa > 1:240269/1‑150 (MQ=255) cAACGATCAACTGGATGGC‑GGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGa < 2:240269/150‑1 (MQ=255) tCAACTGGATGGC‑GGTCACGTCAGGGGCTCTTCCTTGGAGTTTGTGGGTTACAGTGCGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTGCTGGAGCAACGCGTGGAGCGTGAACGTTg < 2:382782/150‑1 (MQ=255) tCAACTGGATGGC‑GGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTg < 2:219060/150‑1 (MQ=255) tGGATGGC‑GGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGC‑CGCGCGCACTGTAAcc < 2:535743/68‑1 (MQ=25) tGGATGGC‑GGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGC‑CGCGCGCACTGTAAcc > 1:535743/1‑68 (MQ=25) gC‑GGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCttt < 2:151732/150‑1 (MQ=255) gTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGc > 2:506176/1‑72 (MQ=255) gTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGc < 1:506176/72‑1 (MQ=255) gggCTCGTCCTTGGAGTTGGTGGGTTACAGTGCGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGGCATCACCgggg < 1:442250/150‑1 (MQ=255) ttCCTTGGAGTTGGTGGGTTACAGTGCGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCaga < 1:90841/150‑1 (MQ=255) ttCCTTGGAGTTGGTGGGTTACAGTGCGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCaga > 2:90841/1‑150 (MQ=255) gTTACAGTGCGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGcc > 2:505361/1‑150 (MQ=255) gTTACAGTGC‑CGCGCGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCg > 2:478745/1‑58 (MQ=21) gTTACAGTGC‑CGCGCGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCg < 1:478745/58‑1 (MQ=21) ttACAGTGCGCGCG‑GCACTGTACCCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTTGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGCCt < 1:291099/150‑1 (MQ=255) | GTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGC‑GGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGCCT > NC_002947/837877‑838146 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |