| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 2,873,755 | C→T | A80T (GCC→ACC) | cysC ← | adenosine 5'‑phosphosulfate kinase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,873,755 | 0 | C | T | 97.3% | 112.8 / NA | 37 | A80T (GCC→ACC) | cysC | adenosine 5'‑phosphosulfate kinase |
| Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (13/23); minor base A (1/0); total (14/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.78e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01 | |||||||||||
TCCTACGCGTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCC > NC_000913/2873624‑2873890 | tcCTACGCGTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTa < 2:670427/139‑1 (MQ=255) cgTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCga < 1:267303/139‑1 (MQ=255) cgTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCga < 2:642703/139‑1 (MQ=255) tCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTTAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGAt < 2:533470/139‑1 (MQ=255) tCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGAt > 2:251505/1‑139 (MQ=255) gcgAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCg < 1:34248/139‑1 (MQ=255) aCCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGc < 1:234855/139‑1 (MQ=255) tCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATaa > 2:34454/1‑139 (MQ=255) cTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTAACCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCAACGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAAt < 2:550708/139‑1 (MQ=255) aTAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATc > 2:558886/1‑137 (MQ=255) aTAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATc < 1:558886/137‑1 (MQ=255) tAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAATACCGAGATCGCTGCAGAATCCGTGGCGAACATTGTCGCCATCcag < 2:536814/139‑1 (MQ=255) gTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTTTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATcc < 2:659107/107‑1 (MQ=255) gTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATcc > 1:659107/1‑107 (MQ=255) tCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGt > 1:914255/1‑139 (MQ=255) aGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTAc < 1:585198/139‑1 (MQ=255) caccacCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCg < 2:26270/91‑1 (MQ=255) caccacCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCg > 1:26270/1‑91 (MQ=255) accaGACCGGCTTCAACCATCAAATTCGCGACTTCACCGACGCGTCGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGc < 2:898040/139‑1 (MQ=255) ccGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGtt > 1:338098/1‑139 (MQ=255) ccGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACAGACGCCGAGtt > 1:602984/1‑139 (MQ=255) tCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTa < 2:658292/139‑1 (MQ=255) cAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCcagcag > 1:599247/1‑105 (MQ=255) cAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCcagcag < 2:599247/105‑1 (MQ=255) tCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAg < 1:34454/139‑1 (MQ=255) cgacgGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCg > 2:385857/1‑62 (MQ=255) cgacgGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCg < 1:385857/62‑1 (MQ=255) cgacgGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGAc < 2:907329/134‑1 (MQ=255) cgacgGATATTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGAc > 1:907329/1‑134 (MQ=255) atatTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTAc < 2:602984/139‑1 (MQ=255) atatTCTCTTTACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTAc < 1:900826/139‑1 (MQ=255) ttACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCcag < 1:31282/65‑1 (MQ=255) ttACGATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCcag > 2:31282/1‑65 (MQ=255) gATCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGagag < 2:338098/139‑1 (MQ=255) tCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGc > 1:699607/1‑139 (MQ=255) tCGGTATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGc < 2:669130/139‑1 (MQ=255) cGGAATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTACTTCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGcc > 2:599355/1‑139 (MQ=255) | TCCTACGCGTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCC > NC_000913/2873624‑2873890 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |