Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,698,852 | Δ1 bp | 24.1% | coding (558/1014 nt) | cheBC → | chemotaxis response regulator protein‑glutamate methylesterase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,698,852 | 0 | G | . | 24.1% | 69.8 / 21.3 | 29 | coding (558/1014 nt) | cheBC | chemotaxis response regulator protein‑glutamate methylesterase |
Reads supporting (aligned to +/- strand): ref base G (9/13); new base . (4/3); total (13/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.66e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01 |
GGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATG‑TGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGCCT > NC_002947/1698714‑1698992 | ggTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATG‑TGGACCa > 2:331975/1‑150 (MQ=255) aCCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATG‑TGGACCAGGTGTTTGc < 1:267092/150‑1 (MQ=255) ggcgcgcGGCTCCCCAGCGACGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGTGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATG‑TGGACCAGGTGTTTGCCGCAg < 2:236976/149‑1 (MQ=255) ggCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATG‑TGGACCAGGTGTTTGCCGCAGGCATGGc < 1:429220/150‑1 (MQ=255) cccAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATG‑TGGACCAGGTGTTTGCCGCAGGCATGGCCGAAt < 2:501539/150‑1 (MQ=255) ggTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATG‑TGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAgcagcg > 2:530646/1‑150 (MQ=255) cAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATG‑TGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGcc < 2:296079/150‑1 (MQ=255) aaTCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATG‑TGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCt > 2:327263/1‑150 (MQ=255) aaTCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCA‑CATGCTGGACCa > 2:456316/1‑102 (MQ=255) aaTCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCA‑CATGCTGGACCa < 1:456316/102‑1 (MQ=255) gggTCTTCGGCGGGTGGCCCGGCAGCAGTCGAAGTGCTGCTCAAGCGTTTGCCGGCAGCGTTTACGGCGGCCATCGTGCTGGTCCAGCATG‑TGGACCAGGTGTTTGCCTCAGCCATGGCCTAATGGGTCACCAGCGCAGCAGGCCTGCCg > 1:361494/1‑150 (MQ=255) gTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCA‑CATGCTGGACCAgc < 2:492098/98‑2 (MQ=255) gTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCA‑CATGCTGGACCAgc > 1:492098/1‑97 (MQ=255) ggTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATG‑TGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTTCGCCttcac > 2:396414/1‑146 (MQ=255) cccGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATG‑TGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGa > 1:353824/1‑150 (MQ=255) ggCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATG‑TGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAgg < 1:331975/150‑1 (MQ=255) gcagcaCTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATG‑TGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGc < 1:579810/150‑1 (MQ=255) cAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATG‑TGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCa < 1:530646/150‑1 (MQ=255) aaGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATG‑TGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAg > 2:116758/1‑150 (MQ=255) aaGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATG‑TGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAg > 1:519376/1‑150 (MQ=255) aaGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATG‑TGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAg < 1:441625/150‑1 (MQ=255) ggCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATG‑TGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGc < 1:116758/150‑1 (MQ=255) gTTTCCGGCGGCCATCGTGCTGGTCCAGCATG‑TGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTac < 1:327263/150‑1 (MQ=255) ggcCATCGTGCTGGTCCAGCATG‑TGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAAccacca < 2:502799/150‑1 (MQ=255) caccTGGTCCA‑CATGCTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTac < 2:334753/131‑1 (MQ=255) caccTGGTCCA‑CATGCTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTac > 1:334753/4‑134 (MQ=255) ccTGGTCCA‑CATGCTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGCCt > 2:119463/2‑150 (MQ=255) cAGCATG‑TGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAaccac > 2:301477/1‑132 (MQ=255) cAGCATG‑TGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAaccac < 1:301477/132‑1 (MQ=255) | GGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATG‑TGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGCCT > NC_002947/1698714‑1698992 |
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: ‑ |