Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,130,198 | A→T | 19.1% | intergenic (‑247/+57) | gcvH‑I ← / ← PP_0990 | glycine cleavage system protein H/DoxX family protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,130,198 | 0 | A | T | 19.1% | 105.2 / 19.8 | 47 | intergenic (‑247/+57) | gcvH‑I/PP_0990 | glycine cleavage system protein H/DoxX family protein |
Reads supporting (aligned to +/- strand): ref base A (20/18); new base T (5/4); total (25/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01 |
CAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCCAGGTCAGGTGGGTGG > NC_002947/1130057‑1130346 | cAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGa > 1:579959/1‑150 (MQ=255) aaTAGATAACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGaa > 2:103496/1‑150 (MQ=255) tAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGaaaa < 2:1108535/150‑1 (MQ=255) tAAATTTAGGGAGAAGACATATTTATAGTGTGGGATGATATCGTAACAGATTCCTCTGAGAGGAAAAGGGTTGGAGCATCAAAGCCTTGTGCACCGACGGCGAAGAGGCGGTGTATCGATTTCGCTACGACGTAGCGAAAACGATACAtg > 2:321175/1‑150 (MQ=255) aGCCGTATTTATTGTGTGGAATTATATCTTTACAGGTGCCTCTGAGAGGAAAAGGGTTTCAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTg < 1:241076/150‑1 (MQ=255) gTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCCGTGTATCGTTTTCGCTACGTCGTAGCGaaaa < 1:77756/123‑1 (MQ=255) gTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCCGTGTATCGTTTTCGCTACGTCGTAGCGaaaa > 2:77756/1‑123 (MQ=255) tatttattGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGc < 1:1060405/137‑1 (MQ=255) tatttattGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGc > 2:1060405/1‑137 (MQ=255) ggAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGcc < 2:834456/150‑1 (MQ=255) aaTGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCtt > 1:549524/1‑150 (MQ=255) gATAACGTTACAGGTGCCGCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACAGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGAGAAAACGATACATGTGGTGCCTGTACTGGCGAGGTGCTAATGGCCTTCgg < 2:263572/150‑1 (MQ=255) ttACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGTCGTAGCGaaaa > 2:5692/1‑98 (MQ=255) ttACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGTCGTAGCGaaaa > 2:247060/1‑98 (MQ=255) ttACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGTCGTAGCGaaaa < 1:247060/98‑1 (MQ=255) ttACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGTCGTAGCGaaaa < 1:5692/98‑1 (MQ=255) gCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTcag > 2:628969/1‑150 (MQ=255) gAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCa > 2:736504/1‑150 (MQ=255) aaaGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGc < 2:918407/150‑1 (MQ=255) tttGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAAc > 2:588600/1‑150 (MQ=255) gCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGcc < 2:549524/150‑1 (MQ=255) tCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGgc > 2:1063572/1‑150 (MQ=255) aaCCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGcc > 1:999980/1‑150 (MQ=255) cccTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTc < 1:946979/150‑1 (MQ=255) ccTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGAc < 1:435082/78‑1 (MQ=255) ccTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGAc > 2:435082/1‑78 (MQ=255) ccTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCg > 2:223800/1‑150 (MQ=255) tttGTGCACCGAGGGCGAGAAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCgg < 1:56787/149‑1 (MQ=255) ttGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCCGAAGATGATCCAGCGAAAACGCGCCGGCGCCTCggg > 1:113548/1‑150 (MQ=255) tgtgCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGc < 2:579959/150‑1 (MQ=255) gCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGATTGGTCAGTTGCCAATGACGATCCCTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAAc < 1:981490/150‑1 (MQ=255) aCCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCtt > 1:593526/1‑91 (MQ=255) aCCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCtt < 2:593526/91‑1 (MQ=255) aCCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATAAATGTGGTGACTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCCCTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTGGGGGAACCa > 2:36648/1‑150 (MQ=255) ccGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAg < 2:447711/150‑1 (MQ=255) gCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGcc < 2:999980/150‑1 (MQ=255) tgtgTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCt < 1:628969/150‑1 (MQ=255) tATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCAAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGccc > 2:662018/1‑150 (MQ=255) aTCGTTTTCGCTACGACGTAGCGTAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCAtc > 2:554148/1‑77 (MQ=255) aTCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCAtc < 1:554148/77‑1 (MQ=255) aTCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCca > 2:754565/1‑150 (MQ=255) gTTTTCGCTACGTCGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCGCTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGagca < 1:980150/131‑1 (MQ=255) gTTTTCGCTACGTCGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCGCTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGagca > 2:980150/1‑131 (MQ=255) ttttCGCTACGTCGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCACGCACCGCGGGCAACCAGGAGCAGCAAGAGCCCTGCCcaggtc > 1:851913/1‑150 (MQ=255) tACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCCAGGTCAggtggg > 2:49735/1‑150 (MQ=255) tACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCCAGGTCAggtggg > 1:32652/1‑150 (MQ=255) tACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCCAGGTCAggtggg < 1:1063572/150‑1 (MQ=255) acgacgTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCCAGGTCAggtgggt < 1:698365/150‑1 (MQ=255) gacgTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCCAGGTCAGgtgggtgg > 2:638523/1‑150 (MQ=255) | CAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCCAGGTCAGGTGGGTGG > NC_002947/1130057‑1130346 |
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: ‑ |