Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 2,653,105 | G→C | 8.7% | W204C (TGG→TGC) ‡ | sseA → | 3‑mercaptopyruvate sulfurtransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,653,105 | 0 | G | C | 8.7% | 117.4 / 6.0 | 46 | W204C (TGG→TGC) ‡ | sseA | 3‑mercaptopyruvate sulfurtransferase |
Reads supporting (aligned to +/- strand): ref base G (28/14); new base C (2/2); total (30/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.02e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.89e-02 |
TATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCAC > NC_000913/2652975‑2653244 | tattGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTCCACGGAAcattca < 1:631070/143‑6 (MQ=255) tattCGCAACCCATGAAACTACAGCGCACATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGCTCATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTCCACGGAAcattca < 1:14297/143‑6 (MQ=255) ttGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGGTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCGGGAGCACTGAATGTTACGTGGACGGAACTgggacg > 2:268407/1‑139 (MQ=255) ccATGAAAATACTGCGCAAAGTAGTGATGCGCGCCCGGCTGCACGTCTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCg < 2:126595/143‑1 (MQ=255) tGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAAc > 2:624600/1‑143 (MQ=255) gAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACt > 2:799032/1‑143 (MQ=255) gAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACt > 1:612989/1‑143 (MQ=255) tACGGCGCAAATTATGGATGCCCGCCCGGCTGCACGTTTTAGCGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTaaaaa < 2:507443/143‑1 (MQ=255) gtgctcgccgcccgGCTGGACGTTTTGACGCAGGAGTTCATGAACCTCGCCCAGGCGTACGTCGCGGACATATTCCCGGAGCATTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAACCGACCGATGAACTg < 2:707775/136‑1 (MQ=255) aTGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGAt > 2:616589/1‑143 (MQ=255) aTGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGAGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGAt > 1:238631/1‑143 (MQ=255) tGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATg < 1:439307/143‑1 (MQ=255) tGCCCGCCCGGCTGCACGTTTCAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATg < 1:443436/143‑1 (MQ=255) gCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATAttttt < 2:738738/143‑1 (MQ=255) gCTGCACGTTGTATCGCAGAAGTGGCGGGACCTCGCCCAGGCTTACGTCGCGGACATATTCCCGGAGCACTGAATGTCCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATAttttt < 2:801876/143‑1 (MQ=255) aCGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGCAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAACACGACCGCTGAACTGGCTGCGATATTTTTTGGTc > 2:510227/1‑143 (MQ=255) ttAACGCAGAAGCTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCgcggc < 2:209868/143‑1 (MQ=255) aaCGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGGACTGGATGCGATATTTTTTGGTCGCGGCGt > 1:686841/1‑143 (MQ=255) aaCGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTTGTCGCGGCGt > 1:178619/1‑143 (MQ=255) aGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGGCGGAACTGGTGCGCGAAGGCGAACCAAACACGCTCGAGGAACTGGATGCGATAGTGTTTGCTCGCGGCGTCAGCt > 1:670184/1‑143 (MQ=255) aGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTCGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGCTGCGATATTTTTTGGTCGCGGCGTCAGCt > 1:240024/1‑143 (MQ=255) aaGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGCACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATGTTTTTTGGTCGCGGCGTCGGCTAc > 2:75907/1‑143 (MQ=255) gTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTAAGGGGTCGCGGCGTCAGCTACGa > 2:522646/1‑143 (MQ=255) atgaACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACaaa > 1:417598/1‑143 (MQ=255) tgaACCTCGCCCAGCTTTACGTCGCGGACATATTCCCGGAGCACTGCATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACAGATGAACTGGATGCGACATTTTTTGGTCGCGGCGTCCGCTACGACAAAc > 2:349005/1‑143 (MQ=255) aCCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGGCAGCTACGACAAGCCag > 1:424994/1‑142 (MQ=255) gcgcAGGTTTACGTCGCCGATATATTCCCGGAGCACTGACGGTTTCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATc < 2:122600/140‑1 (MQ=255) gCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTTAATGTTCCGTGGACGGAACTGGTGCGAGACCGGGAACTAAAAACGACCGATGAACTGGATGCAATATTGTGTCGTCGCGCCGTCTGCTACGACAAACCAATTTTc > 2:674613/1‑143 (MQ=255) nggTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGAGAAGGCGAACTAAAAACGACCGATGAACTGGACGCGCTATCCTTTGGTCGCGGCGTCAGGTACGCAAAACCAGTGATCggcg > 1:67809/2‑140 (MQ=255) gTTTACGTCGCGGACATAGTCCCGGAGCACAGAATGTTCCGTGGACGGAACTGGTGAGCGAAGGCGAACTAAAAACGACAGATGAACTGGACGCGATATTATTTGGTCGCGGCGGCAGCTACGACCAACCCATTAGCGGCacc > 2:42080/1‑141 (MQ=255) tttACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCt > 1:691990/1‑143 (MQ=255) cGTCGCGGACATATGCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCgg < 2:69423/143‑1 (MQ=255) gcgGACATATTCCCGGAGCACTGAATGGTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGctct < 1:560907/143‑1 (MQ=255) ttCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCGCTGGTGTAGCg > 1:252863/1‑143 (MQ=255) cccGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGc < 2:417598/143‑1 (MQ=255) ccGGAGCACTGAATGTTCCGTGGACGGAACTCGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCa < 1:624600/143‑1 (MQ=255) cACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGACCGCCGta > 2:381085/1‑142 (MQ=255) cTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGt < 1:396326/143‑1 (MQ=255) cTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGt > 2:83166/1‑143 (MQ=255) cTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGAGGAACTGGACGCGACAGTTTTTGGTCGCGGCGTCAGCTACGAGCAACCAATTATCGTCAGCTGCGCCGCGGGTGAAACGCCCGCCGCGGt > 2:343238/1‑143 (MQ=255) tGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTGACCGCAGCCGTGGtt > 1:59905/1‑143 (MQ=255) ccagTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGtttt > 2:415556/4‑143 (MQ=255) ccagTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGtttt > 2:122317/4‑143 (MQ=255) ttCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGCCTCTGGTGTCACGGCAGCCGTGGTTTTGATa > 2:343134/1‑143 (MQ=255) cGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCa > 1:599744/1‑143 (MQ=255) gTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTGTGTTAGCAc > 1:317842/1‑143 (MQ=255) | TATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCAC > NC_000913/2652975‑2653244 |
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: ‑ |