| 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 | 98.5% | 210.0 / ‑5.2 | 66 | *151* (TAG→TAA) | moaE | molybdopterin synthase, large subunit |
| Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (31/34); total (31/35) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTTATATCTT > NC_000913/819612‑819871 | tGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAAttt > 2:1105022/1‑139 (MQ=255) cGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGCAGAGCGATCAGCAGGCGGCAAAACGCTGGTAAttttt > 1:1027691/1‑139 (MQ=255) tttATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGc < 1:268199/139‑1 (MQ=255) tctcAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCgg > 1:714233/1‑139 (MQ=255) tctcAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCgg > 1:310122/1‑139 (MQ=255) tctcAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCgg > 1:1428947/1‑139 (MQ=255) aaaCCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAg < 1:1055668/109‑1 (MQ=255) aaaCCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAg > 2:1055668/1‑109 (MQ=255) aaCCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAAttt > 2:182458/1‑139 (MQ=255) cgcACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGc < 2:1582284/139‑1 (MQ=255) gcACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCt < 2:1118389/139‑1 (MQ=255) cACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTc < 2:1428947/139‑1 (MQ=255) ttCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCt < 2:640193/133‑1 (MQ=255) ttCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCt > 1:640193/1‑133 (MQ=255) ttCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTAc > 2:1396944/1‑139 (MQ=255) ttCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTAc > 2:189527/1‑139 (MQ=255) gcgcGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGttt > 2:364123/1‑139 (MQ=255) gcgcGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGttt > 2:1478441/1‑139 (MQ=255) cgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTcc < 1:1551737/107‑1 (MQ=255) cgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTcc > 2:1551737/1‑107 (MQ=255) cgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCaa < 1:1533199/112‑1 (MQ=255) cgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCaa > 2:1533199/1‑112 (MQ=255) cgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCt > 2:663128/1‑122 (MQ=255) cgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCt < 1:663128/122‑1 (MQ=255) cgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACg > 2:88035/1‑139 (MQ=255) cgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACg > 1:696940/1‑139 (MQ=255) cgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACg > 2:1612107/1‑139 (MQ=255) aGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTg < 2:168129/139‑1 (MQ=255) gAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGca < 1:1478441/139‑1 (MQ=255) gACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCt < 2:175995/139‑1 (MQ=255) aCCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCtc > 2:865525/1‑139 (MQ=255) ccGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGTACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCtca < 1:1105022/139‑1 (MQ=255) ccGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCtca < 2:1027691/139‑1 (MQ=255) ccGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCtca < 1:189527/139‑1 (MQ=255) ttGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAg > 1:509179/1‑55 (MQ=255) ttGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAg < 2:509179/55‑1 (MQ=255) aaGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAg < 1:657008/139‑1 (MQ=255) aaGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAg < 2:275942/139‑1 (MQ=255) aaGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAg < 1:182458/139‑1 (MQ=255) aGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGGTGCGTCCGGCAATTTCACCGCTCAGTACATAAGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGt < 2:979748/139‑1 (MQ=255) aGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGt < 1:865525/139‑1 (MQ=255) aGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGt < 1:1612107/139‑1 (MQ=255) aGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGt < 2:696940/139‑1 (MQ=255) cGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCt > 2:1242584/1‑85 (MQ=255) cGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCt > 1:616851/1‑85 (MQ=255) cGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCt < 2:616851/85‑1 (MQ=255) cGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCt < 1:1242584/85‑1 (MQ=255) ggAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGcc < 2:102951/47‑1 (MQ=255) ggAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGcc > 1:102951/1‑47 (MQ=255) gCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGat < 1:1396944/139‑1 (MQ=255) tCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGata < 2:81702/136‑1 (MQ=255) tCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGata > 1:81702/1‑136 (MQ=255) agcagGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGcaca < 1:201690/104‑1 (MQ=255) agcagGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGcaca > 2:201690/1‑104 (MQ=255) agcagGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCt > 2:1405382/1‑108 (MQ=255) agcagGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCt < 1:1405382/108‑1 (MQ=255) agGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGGAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTtat < 1:88035/139‑1 (MQ=255) agGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTtat < 1:1335716/139‑1 (MQ=255) ggcggcAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCa > 1:701232/1‑116 (MQ=255) ggcggcAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCa < 2:701232/116‑1 (MQ=255) ggcggcAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTtata > 2:549928/1‑139 (MQ=255) gcAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTTATATCtt > 1:1515805/1‑139 (MQ=255) ggTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTtat < 2:310498/124‑1 (MQ=255) ggTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTtat > 1:310498/1‑124 (MQ=255) gTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCt > 1:1182040/1‑53 (MQ=255) gTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCt < 2:1182040/53‑1 (MQ=255) | TGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTTATATCTT > NC_000913/819612‑819871 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |