Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,489,374 | C→T | 13.3% | intergenic (+166/‑37) | ybgI → / → cysD | metal‑binding hydrolase/oxidase/sulfate adenylyltransferase subunit 2 |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,489,374 | 0 | C | T | 13.3% | 108.5 / 8.9 | 45 | intergenic (+166/‑37) | ybgI/cysD | metal‑binding hydrolase/oxidase/sulfate adenylyltransferase subunit 2 |
Reads supporting (aligned to +/- strand): ref base C (19/20); new base T (2/4); total (21/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.70e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01 |
TTCGCGGGCCTCATAGGAGCGGGTTCACCCGCGAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGAGTTCGACAACCCGGTGATGCTGTACTCGATCG > NC_002947/1489226‑1489516 | ttCGCGGGCCTCATAGGAGCGGGTTCACCCGCGAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCg > 1:159117/1‑150 (MQ=255) cgGGCCTCATAGGAGCGGGTTCACCCGCGAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAgg < 1:905816/150‑1 (MQ=255) ggCCTCATAGGAGCGGGTTCACCCGCGAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGcc > 2:287253/1‑150 (MQ=255) tCATAGGAGCGGGTTCACCCGCGAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATACAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCTGCGGGCCGtgt < 1:1323330/150‑3 (MQ=255) aTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCTGCGGGCCGtgttc < 1:1139367/150‑5 (MQ=255) aTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCTGCGGGCCGtgttc < 2:1393187/150‑5 (MQ=255) gggTTCACCCGCGAAGAGGCCGGCCTTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGt < 2:5841/150‑1 (MQ=255) ccGCGAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTGGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTg < 1:859037/150‑1 (MQ=255) ccGCGAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTGGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTg > 2:859037/1‑150 (MQ=255) ccGCGAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTg > 1:283533/1‑150 (MQ=255) ccGCGAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTg > 2:451744/1‑150 (MQ=255) gAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTa < 1:1021096/150‑1 (MQ=255) gagGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGcc > 2:363462/1‑150 (MQ=255) agGCCGGCCCTGCCAGCCCCAACCCCAAGTCATATCGTTAGACTTTTTCGCTCTATCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCa < 1:653950/150‑1 (MQ=255) agGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCTGCGGGCCGt > 1:337145/1‑122 (MQ=255) agGCCGGCCCTGCCAGCACAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCTGCGGGCCGt < 2:337145/122‑1 (MQ=255) ggccCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTc < 1:1024474/150‑1 (MQ=255) gCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACaaa < 1:134821/150‑1 (MQ=255) ccAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAAc < 1:305470/150‑1 (MQ=255) cccccAGTCCTATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACtt < 2:164155/150‑1 (MQ=255) cccccAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACtt < 2:159117/150‑1 (MQ=255) ccccAGTCATATCGTTAGACTTTTTCTATTTATCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTTGACAAACTGACGCACTTg > 2:721660/1‑150 (MQ=255) tatCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTg < 2:883130/150‑1 (MQ=255) atCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGt > 2:74503/1‑108 (MQ=255) atCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGt < 1:74503/108‑1 (MQ=255) aCTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGgag < 1:451744/150‑1 (MQ=255) tCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCcaca < 2:283533/150‑1 (MQ=255) tcctAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAg < 1:287253/150‑1 (MQ=255) tcctAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAg > 2:352372/1‑150 (MQ=255) tAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTc < 1:1082580/90‑1 (MQ=255) tAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTc > 2:1082580/1‑90 (MQ=255) gctgctGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTgg > 1:77314/1‑106 (MQ=255) gctgctGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTgg < 2:77314/106‑1 (MQ=255) tgtgATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAACATCCACATCATCCGCGAGGTGGCCGCCGAGTTGGACaa > 2:566947/1‑150 (MQ=255) tgATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGAGTTCGACAAcc > 1:160603/1‑150 (MQ=255) tAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGAGTTCGACAACCCgg < 1:1307151/150‑1 (MQ=255) ggCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGAGTTCGACAACCCGGTGATGc < 1:363462/150‑1 (MQ=255) gctcgctcGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGAGTTCGACAACCCGGTGATGCt > 1:232459/1‑150 (MQ=255) tcgctcgAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGAGTTCGACAACCCGGTGATGCTGt > 2:1016449/1‑150 (MQ=255) tcgctcgAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGAGTTCGACAACCCGGTGATGCTGt > 2:612182/1‑150 (MQ=255) tcgctcgAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGAGTTCGACAACCCGGTGATGCTGt > 1:1280693/1‑150 (MQ=255) gctcgAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGAGTTCGACAACCCGGTGATGCTGTAc > 1:400191/1‑150 (MQ=255) acacGGCCCGCAGGCCGTCCATAAGATCGTTTTTCATGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATcc > 1:696373/1‑103 (MQ=255) acacGGCCCGCAGGCCGTCCATAAGATCGTTTTTCATGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATcc < 2:696373/103‑1 (MQ=255) gacGGCCTGCGGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGAGTTCGACAACCCGGTGATGCTGCACTCGATCg > 1:830711/2‑150 (MQ=255) | TTCGCGGGCCTCATAGGAGCGGGTTCACCCGCGAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGAGTTCGACAACCCGGTGATGCTGTACTCGATCG > NC_002947/1489226‑1489516 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |