Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 1,496,850 | Δ1 bp | intergenic (+219/‑6) | ydcJ → / → opgD | DUF1338 domain‑containing protein YdcJ/glucan biosynthesis protein D |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,496,848 | 0 | T | . | 93.3% | 49.4 / ‑1.8 | 15 | intergenic (+217/‑8) | ydcJ/opgD | DUF1338 domain‑containing protein YdcJ/glucan biosynthesis protein D |
Reads supporting (aligned to +/- strand): ref base T (0/0); major base . (9/5); minor base A (1/0); total (10/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01 |
TCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGCTATGGCCGCCGTGTGCGGTACCAGCGGCAT > NC_000913/1496777‑1496920 | tCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCAC‑TTCAGGTATg < 2:282187/81‑1 (MQ=255) aataaGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCAC‑TTCAGGTATGGATCg > 1:282206/1‑81 (MQ=255) ataaGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCAC‑TTCAGGTATGGATCg > 2:282207/1‑80 (MQ=255) gCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCAC‑TTCAGGTATGGATCGTAGACGATttat < 2:282188/78‑1 (MQ=255) gAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCAC‑TTCAGGTATGGATCGTAGACGATTTATTAAAGGtt > 1:282211/1‑80 (MQ=255) aCCAGCGCGGGACCATTCACAACACCAGAAGGACTCAC‑TTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGc > 1:282218/1‑80 (MQ=255) cAGCGCGGGACCATTCACAACACCAGAAGGACTCAC‑TTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGCTAt > 2:282212/1‑81 (MQ=255) gcgcGGGACCATTCACAACACCAGAAGGACTCAC‑TTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGCTATgg < 2:282189/81‑1 (MQ=255) gcgGGACCATTCACAACACCAGAAGGACTCAC‑TTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGCTATGgcc < 2:282190/81‑1 (MQ=255) gggACCATTCACAACACCAGAAGGACTCAC‑TTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGCTATGgccgc > 1:282198/1‑81 (MQ=255) gACCATTCACAACACCAGAAGGACTCA‑ATTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGCTATGGccgccg > 2:282205/1‑81 (MQ=255) aacaCCAGAAGGACTCAC‑TTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGCTATGGCCGCCGTGTGCGGTAc < 1:282191/81‑1 (MQ=255) aacaCCAGAAGGACTCAC‑TTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGCTATGGCCGCCGTGTGCGGTAc > 2:282208/1‑81 (MQ=255) aacaCCAGAAGGACTCAC‑TTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGCTATGGCCGCCGTGTGCGGTAc > 2:282216/1‑81 (MQ=255) aGGACTCAC‑TTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGCTATGGCCGCCGTGTGCGGTACCAGCGGCAt > 2:282215/1‑81 (MQ=255) | TCATAAATAAGCACTTTCTTGCCGCTGAAAACGACCAGCGCGGGACCATTCACAACACCAGAAGGACTCACTTTCAGGTATGGATCGTAGACGATTTATTAAAGGTTCAATGGCTATGGCCGCCGTGTGCGGTACCAGCGGCAT > NC_000913/1496777‑1496920 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |