Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 524,180 | T→G | 20.3% | intergenic (+47/‑769) | tyrS → / → PP_16SC | tyrosine‑‑tRNA ligase/16S ribosomal RNA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 524,180 | 0 | T | G | 20.3% | 73.6 / 2.7 | 47 | intergenic (+47/‑769) | tyrS/PP_16SC | tyrosine‑‑tRNA ligase/16S ribosomal RNA |
Reads supporting (aligned to +/- strand): ref base T (8/25); new base G (6/2); total (20/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.20e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGAGGCCGGTAGGCC‑TCCCCGC‑‑TTTGT‑TTTT‑GGGG‑TT‑TGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACC > NC_002947/524034‑524327 | tGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGAGGCCGGTAGGCC‑Tcc < 3:8523/149‑1 (MQ=255) tGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGAGGCCGGTAGGCC‑Acc > 2:5496/1‑149 (MQ=255) aGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGAGGCCGGTAGGCC‑GCCCCGC‑‑TTTGT‑TTTT‑g < 3:228548/149‑1 (MQ=255) ttCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGGGGGCCGGT‑GGGC‑GCCCCCC‑‑CTTGT‑TTTT‑TGGG‑tt > 1:129986/1‑149 (MQ=255) ttCATGTTTGTCCTAGGTGCTACCCATGTGTGCCAACCGGGCAAGAAGTCGTTAGGCCGGGTTACCCTCAAGGCTGAGTGATTCCCTGGATGGTGTAGCTATATAGAAGC‑GGGGATGGCCGTAGGCCAACCCCCT‑‑TTTGT‑TTTT‑ggctt > 7:79262/1‑146 (MQ=255) cATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGCGGCCGGTAGGCC‑TCCCCGC‑‑GTTGTGTTTT‑GGGG‑ttg > 5:236872/1‑148 (MQ=255) cATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGAGGCCGGTAGGCC‑TCCCCGC‑‑TTTGT‑TTTT‑GGGG‑tt < 1:28504/147‑1 (MQ=255) cATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGAGGCCGGTAAGCC‑AAACCCC‑‑ATTGT‑TTTT‑TGTG‑tt‑tt > 6:103914/1‑148 (MQ=255) cTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGGGGGCGGGTGGGC‑GGCCCGC‑‑GGTGT‑TTttggggg‑tt‑tgggggtggtg > 2:79300/1‑145 (MQ=255) cTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGGGGGCGGGGTGGC‑GGCCCGC‑‑GTTTT‑TTTT‑TTGG‑TT‑ttgtgggttgtg > 1:291965/1‑143 (MQ=255) cTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGAGGCCGGTAGGCC‑TACCCGC‑‑GTTGT‑TTTT‑TGGG‑TT‑TGGTGGTTGGtg > 2:27155/1‑149 (MQ=255) cTACCCTTGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTTCCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATTTAGAAGC‑GGGGAGGCCGGTAGGCC‑TCCCCGC‑‑TTTGT‑TTTT‑GGGG‑TT‑TGGTGGTTGCTGTTGTagg < 1:107129/149‑1 (MQ=255) cTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGAGGCCGGTTGGCC‑ACCCCGC‑‑GTTGT‑TTTT‑TGGG‑TT‑TTGTGGTTTCagttgtttg > 2:67914/1‑146 (MQ=255) cTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGAGGCCGGTAGGCC‑TCCCCGC‑‑TTTGT‑TTTT‑GGGG‑TT‑TGGTGGTTGCTGTTGTagg < 3:91446/149‑1 (MQ=255) cccATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGAGGCCGGTAGGCC‑TCCCCGC‑‑TTTGT‑TTTT‑GGGG‑TT‑TGGTGGTTGCTGTTGTaggg < 3:88891/147‑1 (MQ=255) gtgtgCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGGGGGCGGGTGGCC‑GCCCCGC‑‑GTTGT‑TTTT‑TGGG‑TT‑TGGTGGTTGGGGTTGTtgggggggggg > 8:184796/1‑149 (MQ=255) tgtgCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGAGGCCGGTAGGCC‑TCCCCGC‑‑TTTGT‑TTTT‑GGGG‑TT‑TGGTGGTTGCTGTTGTagggggaggggg < 1:5496/149‑1 (MQ=255) gggCATGAAGTGGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATTGTGTAGCTTTATAGAAGC‑GGGGAGGCCGGTAGGCC‑TCCCCGC‑‑TTTGT‑TTTT‑g < 7:11892/106‑1 (MQ=38) gggCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGAGGCCGGTAGGCC‑TCCCCGC‑‑TTTGT‑TTTT‑gg < 1:17439/107‑1 (MQ=255) ggCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGCGGCCGGTTGGCC‑GGCCCGG‑‑GTTTT‑TTTT‑GGGGTTT‑TGGGTGTTG‑TGTTGTTGGGGGAGGGGGCTGtttttt > 5:150455/1‑145 (MQ=255) ggCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGAGGCCGGTAGGCC‑TCCCCGC‑‑TTTGT‑TTTT‑GGGG‑TT‑TGGTGGTTGCTGTTGTAGGGGGAGGGGGCGGTTTTgaga < 2:69671/149‑1 (MQ=255) ggCAAGAAGTCGTTTGGCCGGGCTACCCTCAAGGCTGAGGGATTGCCTGGATGGGGTAGCTATATAGAAGC‑GGGGAGGCCGGTAGGCC‑TCCCCGC‑‑TTTGT‑TTTT‑GGGG‑TT‑TGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTgaga < 8:306025/149‑1 (MQ=255) gCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGGGGCGGGGTGTCC‑CCCCCCCCATTTTT‑TTTT‑GGGT‑TTGTGGTGGTGGTTTGTTTGGGGGGGGGGGGGTttttatg > 6:25859/1‑146 (MQ=255) tgtagTCGTTTGGCGGGGTTACCGTCACGTCTTTGTGATTGCCTTGTTGGTGTAGCTATATAGAAGC‑GGGGAGGCCGGTAGGCC‑TCCCCGC‑‑TTTGT‑TTTT‑GGGG‑TT‑TGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGgt < 5:18006/146‑1 (MQ=37) tctggtGGCCGGGTTAGCCTCACGGCTGTCTGATTGCCTGGATTGTGTCGCTATATAGAAGC‑GGGGAGGCCGCTTGGCC‑TCCCCGC‑‑TTTGT‑TTTT‑GGGG‑TT‑TGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTgg < 8:235086/144‑1 (MQ=37) tGGCCGGGTTACCCTCAAGGCTGTGTGATTGCCTGGATGGTGTTGCTATTTAGTAGC‑GGGGCGGCCGGTAGGCC‑GCCCCGC‑‑TTTGT‑TTTT‑GGGG‑TT‑TGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGt < 5:25859/148‑1 (MQ=38) ggCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGAGGCCGGTAGGCC‑TCCCCGC‑‑TTTGT‑TTTT‑GGGG‑TT‑TGGTGGTTGCTGTTGTAGGGGGTGGGGGCTGTTTTGAGAGGGTGTTGGCGGTg < 1:27155/148‑1 (MQ=255) gcggggTTACCCTCAAGGCGGTGTTTTTGGCTGGATGGTTTTGCTATTTTGATGC‑GGGGCGGCGGGTAGGCC‑TCCCCGC‑‑TTTGT‑TTTT‑GGGG‑TT‑TGGGGTTTGGAGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTgag < 1:67914/146‑1 (MQ=25) gCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGAGGCCGGTAGGCC‑TCCCCGC‑‑TTTGT‑TTTT‑GGGG‑TT‑TGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTga < 5:103914/148‑1 (MQ=255) cGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGAGGCCGGTAGGCC‑TCCCCGC‑‑TTTGT‑TTTT‑GGGG‑TT‑TGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGAt < 1:197052/149‑1 (MQ=255) aaGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGAGGCCGGTTGGCC‑TCCCCGC‑‑TTTGT‑TTTT‑TGGG‑TT‑TTGTGGTTGCTGTTGTATGGGGAGGGGGCTGTTTTTAGAGGGTGTTTGGGGTGAGATGTGTACGCcac > 5:328381/1‑148 (MQ=255) tggCTGAGTGATTGCCTGGAGGTGGTCGCTATATAGAAGC‑GGGGAGGCCGGTAGGCC‑TCCCCGC‑‑TTTGT‑TTTT‑GGGG‑TT‑TGGTGGGTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACaa < 5:64123/148‑1 (MQ=25) gCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGGGGGCCGGAGGGCG‑CCCCCCG‑‑TTTTT‑TTTT‑GGGG‑TT‑TTG‑GTTTGGGGTTTTTGGGGGGGGGGGCGGGTTTTTGTGGGTGTTTGGGGGGGGATGTGTGCGGCACaaa > 6:64123/1‑147 (MQ=255) gAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGAGGCCGGTAGGCC‑TCCCCGC‑‑TTTGT‑TTTT‑GGGG‑TT‑TGGTGGTTGCt < 5:64759/85‑1 (MQ=255) tGCCTGGATGGGGTAGCTATATAGACGC‑GGGGCGGCCGGTAGGCC‑TCCCCGC‑‑TTTGT‑TTTT‑GGGG‑TT‑TGGGGGGTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAgcgc < 7:77396/149‑1 (MQ=25) cTGGTTGGTGTCGCTATATAGAAGC‑GGGGAGGCGGGTAGGCC‑TCCCCGC‑‑TTTGT‑TTTT‑GGGG‑TT‑TGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGAt < 1:131000/121‑1 (MQ=255) gATGGTGTAGCTATATAGAAGC‑GGGGAGGCCGGTAGGCC‑TCCCCGC‑‑TTTGT‑TTTT‑GGGG‑TT‑TGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGCGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCtgt < 1:229552/147‑1 (MQ=255) aaGC‑GGGGAGGCCGGTAGGCC‑GCCCCCC‑‑GTTGT‑TTTT‑TGGG‑TT‑TGGTGGTTGGTGTTGTAGGGGGAGGGGGCTGTTTTTAGAGGGTGTTGGCGGTGAGATGTTTCCGCCACAACCCATGCAGCGCCAGTGAGATAGAGCCCCGCCcccg > 2:167748/1‑149 (MQ=255) c‑GGGGAGGCCGGTAGGCC‑TCCCCGC‑‑TTTGT‑TTTT‑GGGG‑TT‑TGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTcc > 5:126637/1‑149 (MQ=255) gCCGGTAGGCC‑TCCCCGC‑‑TTTGT‑TTTT‑GGGG‑TT‑TGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCtt > 8:25339/1‑147 (MQ=255) cGGTAGGCC‑TCCCCGC‑‑TTTGT‑TTTT‑GGGG‑TT‑TGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTgaga < 7:184796/149‑1 (MQ=255) cGGTAGGCC‑TCCCCGC‑‑TTTGT‑TTTT‑GGGG‑TT‑TGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTgaga < 6:328381/149‑1 (MQ=255) cGGTAGGCC‑TCCCCGC‑‑TTTGT‑TTTT‑GGGG‑TT‑TGGTGGGTGCTGGTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTgaga < 1:79300/149‑1 (MQ=255) ggTAGGCC‑TCCCCGC‑‑TTTGT‑TTTT‑GGGG‑TT‑TGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTgaga > 2:217272/1‑148 (MQ=255) gCC‑TCCCCGC‑‑TTTGT‑TTTT‑GGGG‑TT‑TGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCa > 6:74515/1‑149 (MQ=255) cc‑TCCCCGC‑‑TTTGT‑TTTT‑GGGG‑TG‑TGGTGGTTGCTGTTGTAGGGGGTGGGGGCTGTTTTGTGAGGGTGTTGGCGGTTAGATGTGTCCGGCAGAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACTTTCag < 8:79262/149‑2 (MQ=25) c‑TCCCCGC‑‑TTTGT‑TTTT‑GGGG‑TT‑TGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCAcc > 6:632/1‑149 (MQ=255) | TGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGAGGCCGGTAGGCC‑TCCCCGC‑‑TTTGT‑TTTT‑GGGG‑TT‑TGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACC > NC_002947/524034‑524327 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |