Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 1,759,680 | (C)5→6 | coding (844/1221 nt) | sufS ← | cysteine desulfurase, stimulated by SufE, selenocysteine lyase, PLP‑dependent |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,759,675 | 1 | . | C | 92.9% | 136.9 / 3.0 | 42 | coding (849/1221 nt) | sufS | cysteine desulfurase, stimulated by SufE, selenocysteine lyase, PLP‑dependent |
Reads supporting (aligned to +/- strand): ref base . (1/2); new base C (23/16); total (24/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.67e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01 |
TTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑CCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTT > NC_000913/1759537‑1759808 | ttGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAatgatg < 1:1445824/139‑1 (MQ=255) gggCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑cccc > 2:545077/1‑139 (MQ=255) agTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTggg > 2:226419/1‑139 (MQ=255) gTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGgt > 2:1587020/1‑139 (MQ=255) gaTCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTAc > 2:617116/1‑139 (MQ=255) gaTCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTAc < 2:1389666/139‑1 (MQ=255) tCCGGTACAGATTCCAGCTGCGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCg > 2:765159/1‑139 (MQ=255) ggTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCt > 1:1108495/1‑139 (MQ=255) aTTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑ccccccggt > 2:276121/1‑112 (MQ=255) aTTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑ccccccggt < 1:276121/116‑5 (MQ=255) tCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑ccccc > 2:195526/1‑110 (MQ=255) tCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑ccccc < 1:195526/110‑1 (MQ=255) cgcATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGt > 2:767098/1‑139 (MQ=255) gcATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTc < 2:436832/139‑1 (MQ=255) cAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGcc > 1:1268332/1‑139 (MQ=255) aTTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTc > 1:641369/1‑139 (MQ=255) gCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCa < 2:785935/96‑1 (MQ=255) gCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCa > 1:785935/1‑96 (MQ=255) gCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTc < 2:1356846/109‑1 (MQ=255) gCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTc > 1:1356846/1‑109 (MQ=255) gCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTg > 2:1105374/1‑139 (MQ=255) tATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGc < 2:1268332/139‑1 (MQ=255) aTGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCg > 1:1596715/1‑139 (MQ=255) gCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAAc < 2:544946/139‑1 (MQ=255) ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑ccccccg > 2:1480701/1‑51 (MQ=255) ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑ccccccg < 1:1480701/53‑3 (MQ=255) ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑cc < 1:1374165/48‑1 (MQ=255) ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑cc > 2:1374165/1‑48 (MQ=255) cccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCg < 1:1587020/139‑1 (MQ=255) ccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt > 2:535284/1‑87 (MQ=255) ccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt < 1:535284/87‑1 (MQ=255) ccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGc > 1:1336607/1‑139 (MQ=255) gcgcCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt > 2:1288264/1‑84 (MQ=255) gcgcCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt < 1:1288264/84‑1 (MQ=255) gcgcCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGcccc > 2:1230463/1‑139 (MQ=255) tCCAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACggcggc < 1:617116/139‑1 (MQ=255) ccAGCGCCGCGCCAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCa < 1:528754/139‑1 (MQ=255) gcgccgcgccAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTcc > 2:964795/1‑126 (MQ=255) gcgccgcgccAAGACCAATGATG‑CCCCCAGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTcc < 1:964795/126‑1 (MQ=255) gccgcgccAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTcc > 2:1404747/1‑139 (MQ=255) cgcgccAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCtt < 1:988345/63‑1 (MQ=255) cgcgccAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCtt > 2:988345/1‑63 (MQ=255) cgcgccAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt > 1:1155373/1‑62 (MQ=255) cgcgccAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt > 2:304025/1‑62 (MQ=255) cgcgccAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt < 2:1155373/62‑1 (MQ=255) cgcgccAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt < 1:304025/62‑1 (MQ=255) cgccAAGACCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGca > 2:1124375/1‑139 (MQ=255) cgccAAGACCAATGATG‑CCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGc < 2:763557/97‑1 (MQ=255) cgccAAGACCAATGATG‑CCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGc > 1:763557/1‑97 (MQ=255) aCCAATGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGc < 2:1596715/139‑1 (MQ=255) aaTGATGCCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGt > 1:990217/1‑61 (MQ=255) aaTGATGCCCCCC‑GGTATTCGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGt < 2:990217/61‑1 (MQ=255) atgatgcCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGcttct > 2:644760/1‑139 (MQ=255) tgatgcCCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGcttctt < 1:1404747/139‑1 (MQ=255) | TTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑CCCCC‑GGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTT > NC_000913/1759537‑1759808 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |