Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 170,839 | T→G | 16.8% | intergenic (+20/‑550) | 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,839 | 0 | T | G | 16.8% | 71.4 / 8.7 | 41 | intergenic (+20/‑550) | PP_0159/PP_16SA | family III CoA‑transferase/16S ribosomal RNA |
Reads supporting (aligned to +/- strand): ref base T (17/17); new base G (5/2); total (22/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.19e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01 |
GACGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAG > NC_002947/170696‑170985 | gaCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGc > 3:312850/1‑149 (MQ=255) gaCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCGGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGc < 8:350798/149‑1 (MQ=255) tgtaGTAGCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTTATGCGGGGGTGCTCTAGATGAAAAGCGGGGAGGCCTGACGGCCT‑cccc < 8:48584/144‑1 (MQ=255) cGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCGGCGTGATGCGGGGGTGCTCTAGATGAAACGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGttt < 2:48401/148‑1 (MQ=255) cGCCGGGCCCCACCAGTGCTGGGCGTCAATACCGATGGGGTGGGGGGGGATGGGTGGGGGCTGTATGGGGGCGTGTTCCCGGGTCTGCGTTATGCGGGGGGGCTCTATATTAAACGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGttt < 6:226722/148‑1 (MQ=255) caccacTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTa > 5:75485/1‑148 (MQ=255) ctgGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGTGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGt > 2:49398/1‑148 (MQ=255) gCATACCGTAGGGTTGCGGGGGTATCTGTGGGGGCGGGATGCTGGCCAGGTACAGCGTCGGCGTGATGCTGGGGTGCTCTAGATCAAACGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGaa < 8:155332/147‑1 (MQ=37) cggatgtggtgCGGGGGGAAGGGTTGGGGCTGGATGCTGGCGTGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTaaa < 4:312850/143‑1 (MQ=255) cagAAGGGGTGCTGGGGGAAGTGTTGTGGGTGTATGCTGGCTAGGTACAGCGTCTGCGTTATGCTGGGGTGCTATATATTAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTaaa < 5:305932/146‑1 (MQ=38) cTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGGCGCTGTCCCGCCTGCCCCCCTTTGTTTTT‑ATCGCGTTAAAATGTTATTTGAATTTAAAGTTGACGccc > 1:190316/1‑148 (MQ=255) ggggAAGTGTTGGGGCTGGAGGCGGGCGAGGTCCCGCGTCTGCGTGATGCGGGGGGGCTCTAGATCAAACGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCtctc < 2:293620/149‑1 (MQ=38) gatgggtTGGGGGTGGATGCGGGCGAGGTTCAGGGTGTGGTTTATGGTGGGTTTCTCTATATCAAAAGCGGGGTGGCCGTAAGGCGT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGaa < 1:346996/144‑1 (MQ=255) gtTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATcccc > 6:167010/1‑149 (MQ=255) gtTGCGGCTGGATGCTGTCGAGGTACAGCGTCAGCGTGAGGCTGGGGTGCTATAGATCAAAACCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATccc < 7:62939/148‑1 (MQ=255) ggggCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTta > 5:73067/1‑149 (MQ=255) cTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCCGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAAt > 6:344498/1‑149 (MQ=255) aTGCTGGCGAGTTACAGCGTCTGCTGTATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATgcgc < 6:75485/149‑1 (MQ=255) tGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGcccc > 1:86318/1‑148 (MQ=255) tGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCCCGTTGTTTTTAAGCGCGATCAAAAGGAATTTAAATTAAATGTTGACGCCCTCTCGCATACCCTTATAATTCGCCCCa > 4:187978/1‑149 (MQ=255) aCAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGc > 3:99567/1‑148 (MQ=255) cGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGCGGCCGGGAGGCCG‑CCCCGCGTTGTTTTTGAGAGCTATGAAAGGTAATTTGATTTTAAGGTTGCCGCGCGGGCGCAGACCCGTATAATGCGGCCCACTTCCCGCGGCATAGGTAc > 6:21089/1‑149 (MQ=255) cGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAAc > 1:256235/1‑149 (MQ=255) tGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGaaa < 2:140615/148‑1 (MQ=255) tGGGGTGCTCTAGATGAAAAGCGGGGCGGCCGGCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCGCTCGAGTCCCCTTATAATGGGCCCCACTTCCAGCGACATCGGAACGCAAAACt > 5:212266/1‑149 (MQ=25) tGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTg < 6:73067/149‑1 (MQ=255) tctAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCATTGTTTTTCAGCGCTATAAAAGGTAAATGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTACCGCGACATCGGAACGAAAAACACCtt > 1:308769/1‑145 (MQ=25) ctAGATGAAAAGCGGGGAGGCCGGCAGGCCG‑CCCCGCGGTGGTTTTCAGGGCGATCAAAGGGAATTGAAAGTAAAGGGTGCCGCCCGCGCCCATCCCCGTATAAGGCGCCCCACGTCCCGCGGCAGCGGCCCGCAAACCTCCGTgagag > 6:339631/1‑148 (MQ=255) tAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATcccc > 5:373036/1‑97 (MQ=255) tAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATcccc < 6:373036/97‑1 (MQ=255) aTGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCaa < 5:167010/149‑1 (MQ=255) tGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAAt < 5:344498/149‑1 (MQ=255) tGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAAt < 4:99567/149‑1 (MQ=255) gAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGGCCTCTCGAAGCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAAt > 7:176014/1‑148 (MQ=255) cGGGGAGGCCTGACGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTg < 6:116795/133‑1 (MQ=255) cGGGGAGGCCTGACGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTg > 5:116795/1‑133 (MQ=255) ggggAGGCCTGACGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAag > 8:69776/1‑149 (MQ=255) ggggAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCACACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAag > 3:313517/1‑149 (MQ=255) gCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCCGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGtt > 2:298353/1‑148 (MQ=255) cGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAg < 2:86318/149‑1 (MQ=255) cGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAg > 1:377705/1‑149 (MQ=255) | GACGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCT‑CCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAG > NC_002947/170696‑170985 |
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: ‑ |