Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP009273 | 1,032,559 | C→T | 9.4% | R119W (CGG→TGG) | hyaF → | protein involved in nickel incorporation into hydrogenase‑1 proteins |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP009273 | 1,032,559 | 0 | C | T | 9.4% | 106.1 / 2.5 | 43 | R119W (CGG→TGG) | hyaF | protein involved in nickel incorporation into hydrogenase‑1 proteins |
Reads supporting (aligned to +/- strand): ref base C (21/18); new base T (2/2); total (23/20) | |||||||||||
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.61e-01 |
TGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCGGGTGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGGCGTTGTGGCAGGCGGCAACGCAAAATCTCTTGCCGACAGATTCGCTGTTACCGCCGCCCATTGATGGCCTGATGAATGGCC > CP009273/1032413‑1032706 | tGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCgg < 2:462840/149‑1 (MQ=255) gCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTACCAATACGCTTCTCGGGGAAGGGGAAGTGGGAGTGCGCATTCATCGTGCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCGGGTgcg < 2:205946/149‑1 (MQ=255) aaTCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCGGGTGCGCAGACGTCGCGGCGa < 1:495353/149‑1 (MQ=255) tgtcAACGCCGACGATCGTCACTTTATCAACGCGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGGAATGCAGGAGGCGATCTTCTGCGGATTATGGCGGGTGCGCAGACGTCGCGGCGaaa < 1:665034/147‑1 (MQ=255) cAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCGGGTGCGCAGGCGTCGCGGCGACAAGt > 1:415843/1‑149 (MQ=255) cAACGCCGACGATCGTCACATTATCAATACGCTTCACGGGGAAAGGGAACTGTCAGAGCGAAAACAGCAGGCTAACGACAGTGAAAGTGAAAAACAGTAAGCAATCAACTGCGGATAATGGCGGGTGCGCAGACGTCGCGACGAATAga > 1:571001/1‑148 (MQ=255) aaCGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCGGGTGGGCAGACGTCGCGGGGGAAAGtt > 1:50425/1‑149 (MQ=255) ctggggtccGTGACTTTTTTGATTCGCTTGTCGGGGAGGGGGAAGTGTCAGTGCGCATTCAGGAGGCTGACGACAGTGAAAGTGAAATACAGGGGGCGATCTTCTGCGGATTATGGCGGGTGCGCAGACGTCGCGGTGAAAAGTTGCTg < 1:343559/141‑1 (MQ=255) gacgaTCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGGGGATTATGGCGGGTGCGGAGACGTCGGCGCGGAAAGTtggtggg > 1:343775/1‑148 (MQ=255) gacgaTCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTCTGGCGGGTGCGCAGACGGCGCGGCGAAAAGTTGCTgga > 2:702979/1‑149 (MQ=255) gacgaTCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATGATGGCGGGTGCGGCGAAGTCGCGGCGAAAAGTTGCTgta > 2:173434/1‑147 (MQ=255) cTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCGGGTGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTgg < 2:742867/149‑1 (MQ=255) cTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCGGGGGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTgg > 2:684896/1‑149 (MQ=255) aTCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCGGGTGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAggc < 2:293776/149‑1 (MQ=255) aatccgcgcctcGGGGAAGGGGAAGTGTCAGTGCGGATTCCGCAGGCTGAAGACAGTGAAACTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCGGGTGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAggctgg < 1:421895/140‑1 (MQ=255) aCGCTTCTCGGGTAAGGGGAAGTTTCAGTGCTCCTTCAGCAGGCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCGGGTGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAGgctggctg < 2:688367/149‑1 (MQ=255) aCGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCGGGTGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAGgctggctg < 2:777727/149‑1 (MQ=255) ggaagTGCCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCGGGTGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGGCGTtgt < 2:311892/149‑1 (MQ=255) aagTGTCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGTGGGTGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGGCGTTGTgg > 1:48644/1‑149 (MQ=255) aagTGTCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGTGGGTGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGGCGTTGTgg > 1:48560/1‑149 (MQ=255) gtCAGTGCGCATTCAGCAGGCTGAAGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCGGGTGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGGCGTTGTGGCAgg > 2:440721/1‑149 (MQ=255) agcagGCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCGGGTGCGCAGACGTAGCGGCGAAAAGCCGCTAGAGGACCCACTGCAGGCTGGCTGCACGCGGCCGGCGTGGTGGCAtgcggggaggcgagat > 1:387142/1‑138 (MQ=255) agGCTGACGACAGTGAAAGTGAACTACAGGAGGCGATCTTCTGCGGATTATGGCGGGTGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGGCGTTGTGGCAGGCGGCAACGCAAAAtctc < 2:47217/149‑1 (MQ=255) agGCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCGGGTGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGGCGTTGTGGCAGGCGGCAACGCAAAAtctc < 2:10812/149‑1 (MQ=255) agGCTGACGACAGGGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCGGGTGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGGCGTTGTGGCAGGCGGCAACGCAAAAtctc < 2:171793/149‑1 (MQ=255) gCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCGGGTGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGGCGTTGTGGCAGGCGGCAACGCAAAATCTCtt > 2:576336/1‑149 (MQ=255) gCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCGGGTGCGCAGACGTCGCAGCGTAAAGTTGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGGCGGTGGGGCCGGCGGCAACGGAAAAtcttct > 2:735501/1‑146 (MQ=255) tGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCGGGTGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGGCGTTGTGGCAGGCGGCAACGCAAAATCTCTTGc > 2:153409/1‑149 (MQ=255) tGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCGGGTGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAGGCTGGCTGCGCACCGCTGGCGTTGTGGCAGGCGGCAACGCAAAATCTCTTGc > 2:255345/1‑149 (MQ=255) tgaaagtgaaaTACAGGAGGCGATGTTCTGCGGATTATGGCGGGTGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGGCGTTGTGGCAGGCGGCAACGCAAAATCTCTTGCCGACAGAtt > 2:644849/1‑149 (MQ=255) gacagtgaaaTACAGGAGGCGGTCTTCTGCGGATTATGGCGGGTGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAGGCTGGGTGCGCGGCGCTGGCGTTCTGGCTGGTGGCAACGCAAACTCTCTTGCGGACAGATTc < 1:273399/146‑1 (MQ=255) agtgaaaTACAGGAGGGGATCTTCTGCGGATTATGGCGGGTGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGGCGTTGTGGCAGGCGGCAACGCAAAATCTCTTGCCGACAGATTCGCt > 2:549450/1‑149 (MQ=255) aaatcgcggaggCGGTCTTCTGCGGATTCTGGCGGGTGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGGCGTTGTGGCAGGCGGCAACGCAAAATCTCTTGCCGACAGATTCGCTGTTa < 1:303170/142‑1 (MQ=255) aggaggCGATCTTCTGCGGATTATGGCGGGTGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGGCGTTGTGGCAGGCGGCAACGCAAAATCTCTTGCCGACAGATTCGCTGGTAccgccg > 1:793975/1‑149 (MQ=255) aggaggCGATCTTCTGCGGATTATGGCGGGTGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGGCGGTGTGGCAGGCGGCAACGCAAAATCTCTTGCCGACAGATTCGCTGTTGccgccg > 1:708261/1‑149 (MQ=255) gggcaTCTTCTGCTGATTATGGCTGGTGCGCAGACGTCGTGGCGAAAAGGTGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGTCGTTGTGGCAGGCGGCAACGCAAAATCTCTTGCCGACAGATTCGCTGTTACCGCCGCCCa < 2:609918/145‑1 (MQ=255) cGATCTTCTGCGGATTATGGTGGGTGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGGCGTTGTGGCAGGCGGCAACGCAAAATCTCTTGCCGACAGATTCGCTGTTACCGCCGCCCAtt < 2:48644/149‑1 (MQ=255) cGATCTTCTGCGGATTATGGTGGGTGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGGCGTTGTGGCAGGCGGCAACGCAAAATCTCTTGCCGACAGATTCGCTGTTACCGCCGCCCAtt < 2:48560/149‑1 (MQ=255) tctGCGGATTATGGCGGGTGCGCAGACGTCGCGGCGGAAAGTTGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGGCGGTGGGGCAGGCGGCAACGCAAAATCTCTTGCCGACAGATTCGCTGGTACCGCCTCCCATTGATGGc > 1:628129/1‑149 (MQ=255) tctGCGGATTATGGCGGGTGCGCAGACGTCGCGGCGAAAAGTTGCCGGAGGACAAACCGGAGGCTGGCGGCGCGCCGCTGGCGTTGTGGCAGGCGGCAACGCAAAATCTCTTGCCGACAGATTCCCCGGTACCGCCGCCCCTTGATgcc > 1:71932/1‑147 (MQ=255) gCGGATTATGGCGGGTGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGGCGTTGTGGCAGGCGGCAACGCCAAATCTCTTGCCGACCGATTCGCTGTTACCGCCGCCCATTGATGGCCtg > 2:820765/1‑149 (MQ=255) tATGGCGGGTGCGCAGCCGTCGGGGCGAAAAGTTGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGGCGTTGTGGCAGGCGGCAACGCAAAATCTCTTGCCGACAGATTCGCTGTTACCGCCGCCCATTGATGGCCTGATGAAt < 2:95237/149‑1 (MQ=255) gCGGGCGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGGCGTTGTGGCAGGCGGCAACGCAAAATCTCTTGCCGACAGATTCGCTGTTACCGCCGCCCATTGATGGCCTGATGAATGGcc > 1:374559/1‑149 (MQ=255) | TGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAATACAGGAGGCGATCTTCTGCGGATTATGGCGGGTGCGCAGACGTCGCGGCGAAAAGTTGCTGGAGGACAAACTGGAGGCTGGCTGCGCGCCGCTGGCGTTGTGGCAGGCGGCAACGCAAAATCTCTTGCCGACAGATTCGCTGTTACCGCCGCCCATTGATGGCCTGATGAATGGCC > CP009273/1032413‑1032706 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |