Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 170,841 | A→C | 9.1% | intergenic (+22/‑548) | PP_0159 → / → PP_16SA | family III CoA‑transferase/16S ribosomal RNA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 170,841 | 0 | A | C | 9.1% | 138.7 / 7.8 | 55 | intergenic (+22/‑548) | PP_0159/PP_16SA | family III CoA‑transferase/16S ribosomal RNA |
Reads supporting (aligned to +/- strand): ref base A (28/22); new base C (3/2); total (31/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.78e-02 |
GAGACGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAGGCA > NC_002947/170694‑170988 | gagaCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGAGCTCTA‑GATGAAAAGCGGGGAGGCCGTCAg < 1:31590/149‑1 (MQ=255) gaCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGCCGTCAGGc < 8:230049/149‑1 (MQ=255) gCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGCGGCCGTCAGGCCTc > 1:172278/1‑149 (MQ=255) gCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGGGTTGGGGCTGGATGCTGGCGAGGTACAGCGGCTGCGTGAGGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGCCGGCAGGCCTc > 5:289559/1‑149 (MQ=255) cAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGCCGTCAGGCCTccc < 4:215131/149‑1 (MQ=255) aGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGGGGTGCGGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCTTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGCCTGACGGCCTcccc < 8:288625/149‑1 (MQ=255) aGTGCAGTATCCCCGGCCGCCACCACTGCGGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGGGTTGGGGCGGGATGCGGGCGCGGTACAGCGTCTGCGTGATGCTGGGGGGCTCTA‑GATGAAACGCGGGGAGGCCTGACGGCCTcccc < 4:65354/149‑1 (MQ=255) gtaGTTTCGCGGGGCGCACCCAGTGGTGGGGGAGCATACCGATGGGGTGCTGGGGTATGTGTGGGGGCTGGATGCTGTCGAGGTCCAGCGTCGGGGTGTGGCTGGGGTGCGCTA‑GATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCt < 8:197504/147‑1 (MQ=255) ggAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCt > 2:126095/1‑149 (MQ=255) caccacTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTa > 3:103255/1‑148 (MQ=255) accccggctgGGCGAGCATACCGCCGTGGGGCTGGGGGAAGTGGTGGGGCTGGAGGCGGGCGAGGTACAGCGTCTGCGTGATGGTGGGGGGCGCTA‑GCTGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATc < 6:289559/143‑1 (MQ=37) aCCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCGATCAAAGGTAATTGAAAttt > 3:240489/1‑146 (MQ=255) ctaaGTGGTGGTGGGGTATGTGTGGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTTATGCTGGGTTGCTCTA‑TATCAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAgg < 5:203582/147‑1 (MQ=38) cGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGttttt > 3:229882/1‑118 (MQ=255) tgCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGcc > 1:51937/1‑149 (MQ=255) cTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGTCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGccc > 5:113374/1‑148 (MQ=255) ggggAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCGTTGTTTTTAAGCGGTATGAAAGGTAATTGAAATTTAAGGTTGACGCCctcgc > 4:152885/1‑147 (MQ=255) ggggAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGCCGTAAGGCCAACCCGCTTTTTTTTTAAGCGCTATAAAAGGTAATTGAAATTAAAGGTTTACGCCCtgtc > 8:56011/1‑149 (MQ=255) ggAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGGGGGCGGCAGGCCTCTCCGCGTTGTTTTTGAGCGCGATAAAAGGTTATTGGAATTAAAGGTTTAAGCCCtctct > 3:10383/1‑147 (MQ=255) ggAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCtctc < 8:107704/147‑1 (MQ=255) ggAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCg > 7:107704/1‑148 (MQ=255) gAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCaaa < 2:259286/115‑1 (MQ=255) gAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCaaa > 1:259286/1‑115 (MQ=255) agggtTGGGGCTGGAGGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGTGGCTCTA‑TATTAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCctct < 4:327261/140‑1 (MQ=39) aGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCctct > 3:327261/1‑143 (MQ=255) tgtTGGGGCTGTAGGCTGCCGAGGTACAGGGTCTGGGTGAGGCGGGGGTGCTCTA‑GATCAAACGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATcc < 4:240489/148‑1 (MQ=38) ttGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCt > 3:344152/1‑149 (MQ=255) ggggCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCATTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTTAAGGTTGAAGCGCTCTCGAATCCCCTta > 4:239520/1‑149 (MQ=255) gCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTtata > 2:65402/1‑148 (MQ=255) gCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTtata < 1:65402/148‑1 (MQ=255) aTGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATgcgc < 4:103255/149‑1 (MQ=255) tGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGcc > 7:113316/1‑149 (MQ=255) cTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGCGGCCGTCAGGCCTCCCCGCGTTGTTTTTGAGCGCTTTCAAAGGTAATTGAAATTAAAGGTTGACGCCCTATCGAATCCCCTTATAATGCGccc > 3:61032/1‑148 (MQ=38) tGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGCCGTCAGGCCACCCCGCTTTGTTTTTCAGCGCTATAAAAGGTAATTGAAATTTAACGTTGAAGCCCGCTCGAATCCCCTTATAATGCGccccc > 8:325115/1‑148 (MQ=39) gCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACt < 4:344152/149‑1 (MQ=255) aCAGCGTCTGCGTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATAAAAGGTAATTTTAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGAGCCCCACTTCCAGc > 3:62192/1‑148 (MQ=255) cGTCTTCGTGATGGGGGGGTGC‑CTATTAAGAAAAGCGGTGAGGGCGTGAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACAt < 7:79794/149‑1 (MQ=21) gttaTGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGc < 1:341015/134‑1 (MQ=255) gTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCGTTGTTTTTAAGGGCTATAAAAGGTAATTTAAATTAAAGGTTGACGCCCGCTCGAATGCCCTTATAATGAGCCCCACATCCAga > 2:341015/1‑136 (MQ=38) tGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTATCGAATCCCCGTATAATGCGCCCCACTTCCAGCGACATCGGAACg > 5:7547/1‑148 (MQ=255) aTGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGGCGTCAGGCCTACCCGCTTTGTTTTTCAGCGCTTTCAAAGGTAATTGTAATTAAAGGTTGACGGCCTGGGGCATGCCCGTATAATGCGGCCCACTTCGCGCGCCATCGGAACgta > 6:124868/1‑146 (MQ=255) tGCTCTA‑GATTAAAAGCGGGCAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTg < 3:152885/149‑1 (MQ=255) ctctA‑GATGAAAAGCGGGGAGGCCGTCATGCCGACCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGCATCCCCGTAAAATGAGCCCCCCTTACAGCGACATCGGAACGAAAAACTCCtt > 3:324003/1‑146 (MQ=25) tA‑GATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGAt > 1:250525/1‑148 (MQ=255) gATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCa < 2:30319/149‑1 (MQ=255) gATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTGTGTTTTTCACCGCAATCAAAGGTAATTGAAAATAAAGGTTTACGCCCACTCGAATCCCCATATAATGCGCCCCACTTCCAGCGACATAGGAACGAAAAACAACTTGAGAttg > 8:183597/1‑147 (MQ=21) aaacaaaGCGGGGAGGCCTGACGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGCCGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCaa > 6:265872/5‑149 (MQ=18) aTGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCaa < 2:250525/149‑1 (MQ=255) ttaaaaGGGGGCAGGCCGTCAGGCCGCCCCGCTTTTTTCATCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCaa < 1:126095/146‑1 (MQ=18) tGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAAt < 7:11399/149‑1 (MQ=255) aaaGCGGGGAGGCCTGACGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCG‑CCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGt > 3:66772/1‑149 (MQ=18) aaaGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGa < 1:114361/148‑1 (MQ=255) aaGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGt < 3:239520/149‑1 (MQ=255) aaGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGAAGCCCTCTCGAATCCCCTTATAATGCACCCAACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGt > 2:246827/1‑149 (MQ=255) ggggAGGCCTGACGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAag > 4:76201/1‑149 (MQ=255) aGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAGGCa > 8:271753/1‑148 (MQ=255) | GAGACGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTA‑GATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAGGCA > NC_002947/170694‑170988 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |