| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 863,898 | T→A | 8.5% | intergenic (+15/+21) | hemH → / ← uraA | ferrochelatase/uracil permease |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 863,898 | 0 | T | A | 8.5% | 151.4 / 5.5 | 59 | intergenic (+15/+21) | hemH/uraA | ferrochelatase/uracil permease |
| Reads supporting (aligned to +/- strand): ref base T (30/23); new base A (2/3); total (33/26) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.48e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-01 | |||||||||||
AGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTGCC > NC_002947/863759‑864044 | agaTAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCAAAGCGGccc > 8:33351/1‑149 (MQ=255) tAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCAAAGCGGCCCCtt < 3:394233/149‑1 (MQ=255) tAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCAAAGCGGCCCCtt < 8:110608/149‑1 (MQ=255) tAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCAAAGCGGCCCCtt < 7:33351/149‑1 (MQ=255) aGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTc > 2:220738/1‑149 (MQ=255) gCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGGTTTGCGGCCCCTTcc > 5:12420/1‑148 (MQ=255) cATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGa > 7:342529/1‑149 (MQ=255) cGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCg < 5:233137/148‑1 (MQ=255) cGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACAt > 1:15980/1‑149 (MQ=255) gAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATc > 3:176996/1‑149 (MQ=255) gAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATc > 4:187415/1‑149 (MQ=255) ctgTTTGTCGAGGCGGGCGGCCAGGCGCTGGTGCTGGTGCCTTGCCTGAATGATCCCCCGTAGTGGGTGAGGGTGTTGGCTGATATGTGTCAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAg < 7:328164/147‑1 (MQ=255) gTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGct > 3:50067/1‑149 (MQ=255) gTCGAGGCGGGCGGCCTGTAGCTGGTGCTGGGGCCGTGCCTGAATGATCACCCGGAGTGGGGGAGGGGGTGGGCTGATATGTGTTAAAAGGCCTGAAGCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGcttct < 2:222411/148‑1 (MQ=255) tCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGcttctt < 4:124164/148‑1 (MQ=255) gAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTtgtg > 1:54902/1‑149 (MQ=255) cgggcggCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTc < 3:187415/149‑1 (MQ=255) gggcggCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCACCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGttcttc > 6:94301/1‑148 (MQ=255) ggcggcCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCa < 1:220738/149‑1 (MQ=255) gcCTGGAGCTGGTGCTGGTGCCGTGCCTGATTGATCACGGGGTGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGTAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTg < 5:94301/149‑1 (MQ=255) ccAGGAGCTGGTGCTCGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGCCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGt > 2:69751/1‑149 (MQ=255) gctggtgctggtGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATAc < 5:239808/117‑1 (MQ=255) gctggtgctggtGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATAc > 6:239808/1‑117 (MQ=255) ctgtggctggtgcCGGGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTAc < 6:12420/143‑1 (MQ=255) ctggtgctggtgCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTAcc > 5:83840/1‑149 (MQ=255) tgctggtgcCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCGTTGCGGCCCCTTCCGACCGGTCATACATCAAGCTTCTTGTGCTTCCAGCTTTCATTACCAgg > 7:345222/1‑148 (MQ=255) ctggtgcCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGACTCGAAGGGGCCGCTTTGCGGCCCCGTGCGCCCGGGAATAAATCCAGCTTATTGTGCTTGCCGCTGTCATTAccccg > 8:334301/1‑144 (MQ=255) ggtgcCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGaa > 5:373101/1‑149 (MQ=255) gtgccgtgccTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTACGAACGGTCATACATCCAGATTCTTGTGCTTCAAGATGTCATTAACAGgaagaa > 3:28537/1‑148 (MQ=255) tgccgtgccTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATc < 6:373101/149‑1 (MQ=255) ccgtgccTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAg < 8:342529/149‑1 (MQ=255) tGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCg < 6:324806/132‑1 (MQ=255) tGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCg > 5:324806/1‑132 (MQ=255) aaTGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAgg < 1:69751/147‑1 (MQ=255) tGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCaa > 5:197041/1‑149 (MQ=255) tCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTGCAGCTGTCATTTCCCGGCAGAAGCAGGTTCAGGGCAATAg > 4:360788/1‑149 (MQ=255) tCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAg > 1:57673/1‑149 (MQ=255) tCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAg < 2:15980/149‑1 (MQ=255) aCCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGGCCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGc > 1:51534/1‑148 (MQ=255) ggTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGc > 4:311161/1‑94 (MQ=255) ggTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGc < 3:311161/94‑1 (MQ=255) ggCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATg < 6:197041/148‑1 (MQ=255) cTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGccc > 1:51186/1‑149 (MQ=255) atatGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCATCTGTCATTACCCGGCAGAATCAGTTTCAGTGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTc < 2:57673/149‑1 (MQ=255) atatGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCtt > 5:366066/1‑148 (MQ=255) atatGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCtt > 5:30697/1‑148 (MQ=255) tatGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGAGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTTTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCa < 4:50067/149‑1 (MQ=255) tgtgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGcc < 4:176996/149‑1 (MQ=255) gtgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCATCCAGCTGTCATTACCCGGCAGAATCAGGTTAAGTGCAATAGCCACAATGGCGCACAGCGCGATTCCAGTCAGGccc > 2:109615/1‑149 (MQ=255) tgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTTCCCGGCAGAATc > 8:86834/1‑97 (MQ=255) tgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATc < 7:86834/97‑1 (MQ=255) tgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGccc < 6:366066/148‑1 (MQ=255) tgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGccc > 5:146273/1‑148 (MQ=255) cTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCg > 6:393823/1‑149 (MQ=255) aaGGGGCCGCAAAGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCa > 5:8108/1‑149 (MQ=255) ggggCCGCTTTGCGGCCGCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCAc < 8:345222/148‑1 (MQ=255) gCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCt < 2:54902/149‑1 (MQ=255) cGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTGc > 8:307617/1‑149 (MQ=255) gCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTGcc > 3:195044/1‑149 (MQ=255) | AGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTGCC > NC_002947/863759‑864044 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |