Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 35,009 | A→G | 15.9% | intergenic (+27/+31) | czrSA → / ← PP_0031 | metal‑sensor histidine kinase/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 35,009 | 0 | A | G | 15.9% | 68.8 / 6.4 | 32 | intergenic (+27/+31) | czrSA/PP_0031 | metal‑sensor histidine kinase/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (13/13); new base G (2/3); total (16/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01 |
GTCGATCATGAGTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAGCATTCACTT > NC_002947/34862‑35155 | gTCGATCATGACTCTGAACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCTGAAGCCCAGCAACTGCAGCGGGCAATAACCGCCGCTTAACTTGCGCAGACACGACTCGGTCCGTCCAt > 2:306853/1‑150 (MQ=255) tCATGAGTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACg > 1:671553/1‑150 (MQ=255) tGAGTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCGATTGACGGGc < 1:682867/150‑1 (MQ=255) tGAGTCTGCACCAGGGATCGGCGAGCGTGGACGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAACGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTTGATTGACGGGc < 2:293286/150‑1 (MQ=255) gCACCAGGGATCGCCGAGCGTGGAGGGGCCCGAAGAAAACGTCACCTGGTTCCACCGTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCCGTCACCATCACGCCCGTCAATCAACGCGCCttttat > 2:53202/1‑148 (MQ=255) gggATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCt < 1:495745/150‑1 (MQ=255) gCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGc < 2:690412/130‑1 (MQ=255) gCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGc > 1:690412/1‑130 (MQ=255) cGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCCGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCAt < 1:88378/150‑1 (MQ=255) ggggACCGAAGCAAACTTCACCTGGTTCCACCTTCACTTTCCGGAACCCTCGCTCGTGAACCCTGCCGTTGCGGCCCATTCAATTGCTCGGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGaa < 1:53202/147‑1 (MQ=255) cGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCGATTGACGGGc > 2:465459/1‑113 (MQ=255) cGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCGATTGACGGGc < 1:465459/113‑1 (MQ=255) aagaaAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACcaca > 2:481929/1‑150 (MQ=255) gaaAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAgg < 1:481929/150‑1 (MQ=255) ttCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCAt > 1:447885/1‑150 (MQ=255) cTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCg > 2:622146/1‑150 (MQ=255) ggTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTAAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGAc < 1:2946/150‑1 (MQ=255) aCCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGa > 1:452612/1‑150 (MQ=255) ccTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGaa > 2:137647/1‑150 (MQ=255) cACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCAc < 1:109189/150‑1 (MQ=255) cccAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAAtc < 1:137647/150‑1 (MQ=255) ccAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCTTTTTGTTTGACTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAAtct > 1:52753/1‑150 (MQ=255) cAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAAtctt < 2:447885/150‑1 (MQ=255) aTGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTc < 2:671553/150‑1 (MQ=255) cgccgcTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAAt < 1:317729/64‑1 (MQ=255) cgccgcTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAAt > 2:317729/1‑64 (MQ=255) gcTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCca < 1:651902/150‑1 (MQ=255) acGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTACCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACt > 1:176739/1‑150 (MQ=255) ggCCCGTCGATTGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAGCAtt > 1:413180/1‑150 (MQ=255) ggCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAGCAtt > 2:647280/1‑150 (MQ=255) ccGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAGCATTCAc > 2:596659/1‑150 (MQ=255) gTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAGCATTCACtt < 1:622146/150‑1 (MQ=255) | GTCGATCATGAGTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAGCATTCACTT > NC_002947/34862‑35155 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |