| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 819,747 | G→A | *151* (TAG→TAA) | moaE → | molybdopterin synthase, large subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 819,747 | 0 | G | A | 97.4% | 113.5 / ‑3.9 | 38 | *151* (TAG→TAA) | moaE | molybdopterin synthase, large subunit |
| Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (17/20); total (18/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.74e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01 | |||||||||||
CGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTTATATCTTCAGGAGATCGT > NC_000913/819614‑819882 | cGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGttttt > 1:419394/1‑139 (MQ=255) ttGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGt > 1:737248/1‑139 (MQ=255) tATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTcc > 1:629499/1‑139 (MQ=255) tATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTcc > 2:232706/1‑139 (MQ=255) cccGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCa < 2:737248/139‑1 (MQ=255) cgcgCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCAcc < 1:232996/139‑1 (MQ=255) gcgcACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCg < 2:894739/139‑1 (MQ=255) gcgcACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCg < 2:726902/139‑1 (MQ=255) ttCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTAc > 1:758622/1‑139 (MQ=255) gcgcGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAg < 2:748799/87‑1 (MQ=255) gcgcGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAg > 1:748799/1‑87 (MQ=255) gcgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGt > 1:335436/1‑82 (MQ=255) gcgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGt < 2:335436/82‑1 (MQ=255) cgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACg > 2:54447/1‑139 (MQ=255) ccGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTACGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCtca < 2:889546/139‑1 (MQ=255) aGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGt < 1:232706/139‑1 (MQ=255) gggAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGTCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCt < 2:857149/139‑1 (MQ=255) gggAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCt < 2:859568/139‑1 (MQ=255) ggAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGc < 2:44869/117‑1 (MQ=255) ggAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGc > 1:44869/1‑117 (MQ=255) ggAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCa < 2:47936/131‑1 (MQ=255) ggAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCa > 1:47936/1‑131 (MQ=255) agCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGa > 1:823266/1‑139 (MQ=255) gCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGat < 2:823266/139‑1 (MQ=255) cGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCa > 1:582321/1‑79 (MQ=255) cGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCa < 2:582321/79‑1 (MQ=255) cGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGata < 1:54447/139‑1 (MQ=255) cGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGata < 2:758622/139‑1 (MQ=255) ggcggcAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCa < 2:215455/63‑1 (MQ=255) ggcggcAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCa > 1:215455/1‑63 (MQ=255) ggcAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCtcatca > 2:901509/1‑106 (MQ=255) ggcAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCtcatca < 1:901509/106‑1 (MQ=255) ggcAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCa < 1:480136/113‑1 (MQ=255) ggcAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCa > 2:480136/1‑113 (MQ=255) ggcAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTTATATCt > 2:436340/1‑139 (MQ=255) gcAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTTTCCTTAAGATATTCCTTATATCtt > 1:518194/1‑139 (MQ=255) gCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTTATATCTTCAGgaga < 2:518194/139‑1 (MQ=255) gTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTTATATCTTCAGGAGATCGt > 2:617310/1‑139 (MQ=255) | CGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTTATATCTTCAGGAGATCGT > NC_000913/819614‑819882 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |