| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 1,369,921 | A→G | E78G (GAG→GGG) | pspE → | thiosulfate:cyanide sulfurtransferase (rhodanese) |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,369,921 | 0 | A | G | 93.5% | 139.1 / 1.0 | 46 | E78G (GAG→GGG) | pspE | thiosulfate:cyanide sulfurtransferase (rhodanese) |
| Reads supporting (aligned to +/- strand): ref base A (2/1); new base G (22/21); total (24/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01 | |||||||||||
AGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACC > NC_000913/1369825‑1370014 | aggaagTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGGGAt > 1:739007/1‑100 (MQ=255) aagTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGGGATCCt > 1:1097962/1‑100 (MQ=255) gAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGGGATCCTTAGc > 2:1290062/1‑100 (MQ=255) aaGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGGGATCCTTAGCga > 1:727024/1‑100 (MQ=255) agagCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGGGATCCTTAGCgag < 2:45292/100‑1 (MQ=255) gagCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGGGATCCTTAGCg > 2:1121302/1‑97 (MQ=255) agCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGGGATCCTTAGCGAGAt > 1:521213/1‑100 (MQ=255) gcgcATTGCCACCGCCGTGCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGGGATCCTTAGCGAGATg < 2:454657/100‑1 (MQ=255) gcATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGGGATCCTTAGCGAGATggg < 2:1042604/100‑1 (MQ=255) ttGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGAGCAAGCAAAAGGGATCCTTAGCTAGATGTGata < 2:892076/100‑1 (MQ=255) cACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGGGATCCTTAGCGAGATGGGATATAcc > 1:659029/1‑100 (MQ=255) tCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGGGATCCTTAGCGAGATGGGATATACCCACGTTgag < 1:1290073/100‑1 (MQ=255) cGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGGGATCCTTAGCGAGATGGGATATACCCACGTTGAGaa < 1:347293/100‑1 (MQ=255) cGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGGGATCCTTAGCGAGATGGGATATACACACGTTgag > 2:631147/1‑98 (MQ=255) aaaCGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGGGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGt < 1:929862/100‑1 (MQ=255) aaaCGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGGGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGt > 1:1014412/1‑100 (MQ=255) aaCGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGGGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTg < 1:654305/100‑1 (MQ=255) gACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGGGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGcc < 2:971985/100‑1 (MQ=255) cacCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGGGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTg > 1:857351/1‑100 (MQ=255) tGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGGGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGa < 1:17103/100‑1 (MQ=255) gAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGAc > 2:103888/1‑100 (MQ=255) aGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGGGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATc < 2:203542/100‑1 (MQ=255) tgtATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGGGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGc < 1:585436/100‑1 (MQ=255) ttGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGGGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATg > 2:1193108/1‑100 (MQ=255) tGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGGGATCCTTAGCTAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGc < 2:1189828/100‑1 (MQ=255) tGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGGGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGc < 1:985870/100‑1 (MQ=255) tGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGc < 1:1092851/100‑1 (MQ=255) gCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGGGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGcc > 2:319643/1‑100 (MQ=255) cAATGCCGGCCGCCAGTCAGGGCAAGCAAAAGGGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCACTGCCg > 2:1190027/1‑100 (MQ=255) cAATGCCGGACGCCAGTCAGGGCAAGCAAAAGGGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCg > 2:405752/1‑100 (MQ=255) aTGCCGGACGCCAGTCAGGGCAAGCAAAAGGGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGaa < 2:677043/100‑1 (MQ=255) ggACGCCAGTCAGGGCAAGCAAAAGGGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCa < 2:457370/100‑1 (MQ=255) ggACGCCAGTCAGGGCAAGCAAAAGGGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCa < 2:951202/100‑1 (MQ=255) ggACGCCAGTCAGGGCAAGCAAAAGGGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCa > 1:1275731/1‑100 (MQ=255) cGCCAGTCAGGGCAAGCAAAAGGGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAg < 1:709878/100‑1 (MQ=255) gCCAGTCAGGGCAAGCAAAAGGGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAgg > 1:899338/1‑100 (MQ=255) ngTCAGGGCAAGCAAAAGGGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTa > 2:819755/2‑100 (MQ=255) aGTCAGGGCAAGCAAAAGGGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTa > 2:223682/1‑100 (MQ=255) aGGGCAAGCAAAAGGGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAg > 2:292726/1‑100 (MQ=255) gggCAAGCAAAAGGGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGa > 2:372052/1‑100 (MQ=255) gcaagcaaAAGGGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAAt < 1:1297414/100‑1 (MQ=255) aagcaaAAGGGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTc < 2:857345/100‑1 (MQ=255) agcaaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCa > 1:801346/1‑100 (MQ=255) gcaaAAGGGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCAc > 1:449803/1‑100 (MQ=255) caaAAGGGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTACAAGAATTCAcc > 1:575809/1‑100 (MQ=255) caaAAGGGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCAcc < 2:895178/100‑1 (MQ=255) | AGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACC > NC_000913/1369825‑1370014 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |