| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 4,575,319 | C→T | intergenic (+1388/+1593) | yjiT → / ← mcrC | putative uncharacterized protein YjiT/Type IV methyl‑directed restriction enzyme EcoKMcrBC subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,575,319 | 0 | C | T | 90.0% | 50.6 / ‑5.0 | 20 | intergenic (+1388/+1593) | yjiT/mcrC | putative uncharacterized protein YjiT/Type IV methyl‑directed restriction enzyme EcoKMcrBC subunit |
| Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (8/10); minor base G (0/2); total (8/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.95e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01 | |||||||||||
GAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGAAAAAGGATTAACGCTGGATGACTGGCGCGATTTTCTCAAGGTTACGTTGGATCATTACGTTCGCGAGAGCAACTTCACACAGCTGGATGATGAGCTGAAAAACTGGATTGGTAGCCGTTTTTCATCAAAATTCGTCCGTAACCCGGAATCAAAAGATCCTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCA > NC_000913/4575083‑4575533 | gAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGAAAAAGGATTAACGCTGGATGACTGGCGCGATTTTCTCAAGGTTACGTTGGATCATTACGTTCGCGAGAGTAACTTCACACAg < 1:200689/249‑1 (MQ=255) cTTAAGCTTTCTGACAAGTTGCTGTTACGGGAATTCCTGCGTCGCCCTACACGGACTAGCCGCCTGGGCCCACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGACCATACATCAGAGCCCCCTGCACTGGCAGGACAAAGGATTACCGCTGGTAGACCGGCGCGATTTTCTCAAGGCTACGTTGGATCATTATGTTCGCGGGAGTAACTTACCACAGCTGGCTGGTGAGCTGAAACACTGGATTgg > 2:119870/1‑251 (MQ=255) cTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGAAAAAGGATTAACGCTGGATGACTGGCGCGATTTTCTCAAGGTTACGTTGGATCATTACGTTCGCGAGAGTAACTTCACACAGCTGGATGATGAGCTGAAAAACTGGATTGGTAGCCg > 2:147325/1‑251 (MQ=255) tCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGAAAAAGGATTAACGCTGGATGACTGGCGCGATTTTCTCAAGGTTACGTTGGATCATTACGTTCGCGAGAGTAACTTCACACAGCTGGATGATGAGCTGAAAAACTGGATTGGTAGCCGTTTTTCATCCAAATTCGTCCGTAACCCg > 1:192359/1‑251 (MQ=255) tCATGCGTCGCCCTAAACGGAATAACAGCCTGGAAACACAGGTGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATTAGAGCCCCTTGCACTGGCAGGAAAAAGGATTAACGCTGGATGACTGGCGCGATTTTCTCAAGGTTACGTTGGATCATTACGTGCGCGAGAGGAACTTCACACAGCTGGATGATGAGCTGAAAATCTGGATTGGTAGCCGTTTTGCATCAAAATTCGTCCGTAAccc < 2:222173/250‑1 (MQ=255) cGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGAAAAAGGATTAACGCTGGATGACTGGCGCGATTTTCTCAAGGTTACGTTGGATCATTACGTTCGCGAGAGTAACTTCACACAGCTGGATGATGAGCTGAAAAACTGGATTGGTAGCCGTTTTTCATCAAAATTCGTCCGTAACCCGGAATc < 1:147325/251‑1 (MQ=255) cGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCCGGGCCTGGAGAAAATACATCAGAGCCACTTGCAATGGCAGGAAAACGGATTAAAGCTGGCTGACTGGCGCGATTTTCTCAAGGTTACGTTGGATCATTACGTTCGAGAGAGTAACTTCACACAGCTGGATGATGAGCTGAAAAACTGGATTGGTAGCCGTTTTTCATCAAAATTCGTCCGTAACCCGGAATCAACAGATCCTg > 1:120854/1‑250 (MQ=255) ggACTAACAGCCAGTAAACACAGGGGGTGGGTAAGGTGGGTTGCCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGAAAAAGGATTAACGCGGGATGACTGGCGCGATTTTCTCAAGGTGACGTTGGATGATTACGTTCGCGAGAGGAACTTCGCACAGCTGGATGATGAGCTGAAAAACTGGATTGGTAGCCGGCTTTCAGCAAAATTCGTCCGTAGCCCGAAATCAAAAGaggctgaa < 2:122554/251‑8 (MQ=255) aGCCGGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGAAAAAGGATTAACGCTGGATGACTGGCGCGATTTTCTCAAGGTTACGTTGGATCATTACGTTCGCGAGAGTAACTTCACACAGCTGGATGATGAGCTGAAAAACTGGATTGGTAGCCGTTTTTCATCAAAATTCGTCCGTAACCCGGAATCAAAAGATCCTGAAGATAATCa < 1:141599/251‑1 (MQ=255) gCTGGGTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGAAAAAGGATTAACGCTGGATGACTGGCGCGATTTTCTCAAGGTTACGTTGGATCAGTACGTTCGCGAGAGTAACTTCACACAGCTGGATGATGAGCTGAAAAACTGGATTGGTAGCCGTTTTTCATCAAAATTCGTCCGTAACCCGGAATCAAAAGATCCTGAAGATAATCAGAACAGACGCTGGCCTc < 1:178036/251‑1 (MQ=255) gttggttACCAGGGGCTGAAGAAAATACATAAGAACCCCTTGCACTGGCAGGAAAAAGGATTAACGCTGGATGACTGGCGCGATTTTCTCAAGGTTACGTTGGATCATTACGTTCGCGAGAGTAACTTCACACAGCTGGATGATGAGCTGAAAAACTGGATTGGTAGCCGTTTTTCATCAAAATTCGTCCGTAACCCGGAATCAAAAGATCCTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTa > 1:71490/1‑250 (MQ=255) ggAGAAAATACATAAGAGCCCCTTGCACTGGCAGGAAAAAGGATTAACGCTGGATGACTGGCGCGATTTTCTCAAGGTTACGTTGGATCATTACGTTCGCGAGAGTAACTTCACACAGCTGGATGATGAGCTGAAAAACTGGATTGGTAGCCGTTTTTCATCAAAATTCGTCCGTAACCCGTAATCAAAAGATCCTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATc < 1:123425/251‑1 (MQ=255) aTACATAAGAACCCCTTGCACTGGCAGGAAAAAGTATTAACGCTGGATGACTGGCGCGATTTTCTCAAGGTTACGTTGGATCATTACGTTCGCGAGAGTAACTTCACACAGCTGGATGATGAGCGGAAAAACTGGATTGGAAGCCGTTTTTCATCAAAATTCGTCCGTAACCCGGATTCAAATGATCCTGAAGATAATCAGAAGAGACGGTGGCCTCCAATTCGTAATGGCAACGTATCCCATCGTTTAgg > 2:25646/1‑250 (MQ=255) tgctatGGCAGGAAAAAGGAGTAACGCTGGAGGACTGGCGCGATTTTCGCAAGGTTACGTTTGATCATTACGTTCGCGAGAGTAACTTCACACAGCTGGATGATGAGCTGAAAAACTGGATTGGTAGCCGTTTTTCATCAAAATTCGTCCGTAACCCGGAATCAAAAg < 2:115080/163‑1 (MQ=255) tGCACTGGCAGGAAAAAGGATTCACGCTGGATGACTGGCGCGATTTTCTCAAGGTTACGTTGGATCATTACGTTCGCGAGAGTAACTTCACACCGCTGGATGATGAGCTGAAAAACTGGATTGGTAGCCGTTTTTCATCAAAATTCGTCCGTAACCCGGAATCAAAAg > 1:115080/1‑168 (MQ=255) ggCAGGAAAAAGGATTAACGCTGGATGACTGGCGCGATTTTCTCAAGGTTACGTTGGATCATTACGTTCGCGAGAGTAACTTCACACAGCTGGATGATGAGCTGAAAAACTGGATTGGTAGCCGTTTTTCATCAAAATTCGTCCGTAACCCGGAATCAAAAGATCCTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGc < 1:68140/250‑1 (MQ=255) cgcgATTTTCTCAAGGTTACGTTGGATCATTACGTTCGCGAGAGTAACTTCACACAGCTGGATGATGAGCTGAAAAACTGGATTGGTAGCCGTTTTTCATCAAAATTCGTCCGTAACCCGGAATCAAAAGATCCTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGa < 1:36232/251‑1 (MQ=255) cacAAGGTTACGTTGGATCATTACGTTCGCGAGAGTAACTTCACACAGCTGGATGATGAGCTGAAAAACTGGATTGGTAGCCGTTTTTCATCAAAATTCGTCCGTAACCCGGAATCAAAAGATCCTGAAGATTATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAAGTATTGATATTATCa < 2:123416/248‑1 (MQ=255) cttACGTTCGCGAGAGTAACTTCACACAGCTGGATGATGAGCTGAAAACCTGGATTGGTAGCCGTTTTTCATCACAATTCGTCCGTAACCCGGAATCAAAAGATCCTGCAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGACGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTc > 1:101655/2‑206 (MQ=255) aTTACGTTCGCGAGAGTAACTTCACACAGCTGGATGATGAGCGTAAAAACTGGAGTGGTAGCCGTTTTTCAGCAAAATTCGTCCGTAACCCGGAATCAAAAGATCCTGAAGATAATCAGAAGAGACGCTGGCCTCAAATTCGTAATTGCAACGTAGCCCATCGTTGAGCGACGTTGCTGATGCTGGGGGCGGGATTCaaaagagtc < 2:101655/206‑6 (MQ=255) | GAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCAGGAAAAAGGATTAACGCTGGATGACTGGCGCGATTTTCTCAAGGTTACGTTGGATCATTACGTTCGCGAGAGCAACTTCACACAGCTGGATGATGAGCTGAAAAACTGGATTGGTAGCCGTTTTTCATCAAAATTCGTCCGTAACCCGGAATCAAAAGATCCTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCA > NC_000913/4575083‑4575533 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |