Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 2,994,461 | (T)8→7 | pseudogene (443/444 nt) | ygeK ← | putative DNA‑binding transcriptional regulator YgeK |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,994,454 | 0 | T | . | 75.0% | 16.6 / 10.8 | 12 | intergenic (+13/+6) | pbl/ygeK | transglycosylase domain‑containing protein Pbl/putative DNA‑binding transcriptional regulator YgeK |
Reads supporting (aligned to +/- strand): ref base T (1/2); new base . (6/3); total (7/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.23e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAAT > NC_000913/2994380‑2994526 | aaaTATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCttttttt > 1:505530/1‑81 (MQ=255) atatatatTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCtttttttat < 1:505503/81‑3 (MQ=255) tatataCCTGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTATATAGTGCaa < 2:505504/80‑1 (MQ=255) tatataCCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCa > 1:505522/1‑80 (MQ=255) atataCCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTATATAGTGCaaa > 2:505521/1‑80 (MQ=255) ccAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTATATAGTGCAAACACCCa < 2:505505/81‑1 (MQ=255) ggCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTATATAGTGCAAACACCCATAc < 1:505506/81‑1 (MQ=255) aCTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGt > 2:505526/1‑81 (MQ=255) tCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCa < 1:505507/81‑1 (MQ=255) aaGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTc > 1:505523/1‑80 (MQ=255) ttCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATAc > 2:505533/1‑80 (MQ=255) gAGCCTAAAGCCTC‑TTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACAt > 2:505536/1‑81 (MQ=255) ccTAAAGCCTC‑TTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCt > 1:505534/1‑80 (MQ=255) aaGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAAt < 1:505508/80‑1 (MQ=255) | AAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAAT > NC_000913/2994380‑2994526 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |