Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,102,760 | A→G | 8.2% | M1M (ATG→GTG) † | phhB → | pterin‑4‑alpha‑carbinolamine dehydratase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,102,760 | 0 | A | G | 8.2% | 160.6 / 3.2 | 57 | M1M (ATG→GTG) † | phhB | pterin‑4‑alpha‑carbinolamine dehydratase |
Reads supporting (aligned to +/- strand): ref base A (28/24); new base G (3/2); total (31/26) | |||||||||||
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 |
ACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGAAGTACGTGACGGCCACATGGAGCTGGAGCGCGTGTT > NC_002947/5102611‑5102905 | aCATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCa < 1:1078588/150‑1 (MQ=255) aTGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCGTg > 1:91774/1‑150 (MQ=255) cTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTg > 1:520821/1‑150 (MQ=255) gTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGACAACACCTCATGAATGCCTTGaa < 1:583885/150‑1 (MQ=255) aCCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAg > 1:1042539/1‑150 (MQ=255) ccACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTg < 2:250685/150‑1 (MQ=255) aCGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCg > 1:742255/1‑150 (MQ=255) gAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGaa > 2:1055685/1‑150 (MQ=255) atcatGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCCCCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGc > 1:412744/1‑150 (MQ=255) atcatGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCt > 1:371227/1‑124 (MQ=255) atcatGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCt < 2:371227/124‑1 (MQ=255) tGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCcgcg > 1:611843/1‑150 (MQ=255) cATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGaa < 1:268900/111‑1 (MQ=255) cATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGaa > 2:268900/1‑111 (MQ=255) aTAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGa > 1:505776/1‑150 (MQ=255) aaaGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGaaa < 2:520821/150‑1 (MQ=255) ccATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGGCTGCCGCGCCGATGCGCCGAAAGTCa > 2:4725/1‑150 (MQ=255) aTGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGaaa > 2:1058724/1‑144 (MQ=255) aTGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGaaa < 1:1058724/144‑1 (MQ=255) gTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGc < 2:909085/150‑1 (MQ=255) tgcatgcaCCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTgg < 1:735047/150‑1 (MQ=255) atgcaCCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCg > 2:684710/1‑150 (MQ=255) aCCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCt > 1:386800/1‑150 (MQ=255) aCCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCt > 1:564526/1‑150 (MQ=255) aCCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCt > 1:564519/1‑150 (MQ=255) aaGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATc > 2:1099130/1‑150 (MQ=255) aCCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTTAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAAt > 2:513360/1‑150 (MQ=255) aaGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTTAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCg < 1:513360/150‑1 (MQ=255) aaGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCg < 1:36036/150‑1 (MQ=255) ggTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGa < 2:611843/150‑1 (MQ=255) gCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGTCGAGCTGATCCGCGAAATCCCGGACTgg < 1:1055685/150‑1 (MQ=255) gCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTgg < 2:412744/150‑1 (MQ=255) cTGCCTGAGCGACCTTGCTGATAACAACTCGAACCGGAAAACACCTCGTGAATGCCTTGAACCAAGCCCATTGCGAAGCCTTCCGCGCCGATGCGCCGAAAGTCa > 2:533612/1‑105 (MQ=255) cTGCCTGAGCGACCTTGCTGATAACAACTCGAACCGGAAAACACCTCGTGAATGCCTTGAACCAAGCCCATTGCGAAGCCTTCCGCGCCGATGCGCCGAAAGTCa > 2:533602/1‑105 (MQ=255) cTGCCTGAGCGACCTTGCTGATAACAACTCGAACCGGAAAACACCTCGTGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCa < 1:533602/105‑1 (MQ=255) cTGCCTGAGCGACCTTGCTGATAACAACTCGAACCGGAAAACACCTCGTGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCa < 1:533612/105‑1 (MQ=255) tGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGaa > 1:982994/1‑150 (MQ=255) gCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGaa < 2:1042539/150‑1 (MQ=255) cGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGATCGGTCCCGGTTTCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCNAAg < 2:705021/150‑1 (MQ=255) ccTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGAAGTAc > 1:783244/1‑150 (MQ=255) cTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGAAGTACg < 1:509506/150‑1 (MQ=255) ttGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGaa < 1:344428/144‑1 (MQ=255) ttGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGaa > 2:344428/1‑144 (MQ=255) ttGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGAAGTACGt < 2:742255/150‑1 (MQ=255) ttGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGAAGTACGt > 2:437067/1‑150 (MQ=255) ttGCCGATAACAACTCGAACCGGAAAACAACCCATAAAAGCTTATCATGGATATCGTTGCGAAGCATGCCTGTCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGCACATCGAAGTActt > 1:193368/1‑148 (MQ=255) cGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGAAGTACGTGACg > 2:895042/1‑150 (MQ=255) gATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGAAGTACGTGACgg > 1:104655/1‑150 (MQ=255) aacTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTgg < 2:752923/95‑1 (MQ=255) aacTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTgg > 1:752923/1‑95 (MQ=255) acTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGAAGTAc < 1:640337/136‑1 (MQ=255) acTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGAAGTAc > 2:640337/1‑136 (MQ=255) tCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGAAGTACGTGACGGCCACAtgga < 2:982994/150‑1 (MQ=255) tCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGAAGTACGTGACGGCCACAtgga > 2:293085/1‑150 (MQ=255) aaaCACCTCAGGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGAAGTACGTGACGGCCACATGGAGCTGGAgcgc < 1:461763/150‑1 (MQ=255) aaCACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGAAGTACGTGACGGCCACATGGAGCTGGAGcgcg > 2:760492/1‑150 (MQ=255) acCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGAAGTACGTGACGGCCACATGGAGCTGGAGTGCgtgt < 1:895042/150‑1 (MQ=255) ccTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAACCACGGACTGGAACATCGAAGTACGTGACGGCCACATGGAGCTGGAGCGCGTGtt < 1:202383/150‑1 (MQ=255) | ACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAAGCCTGCCGCGCCGATGCGCCGAAAGTCACCGACGAAGAGCTGGCCGAGCTGATCCGCGAAATCCCGGACTGGAACATCGAAGTACGTGACGGCCACATGGAGCTGGAGCGCGTGTT > NC_002947/5102611‑5102905 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |