Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,415,081 | G→A | 12.5% | intergenic (+146/+22) | PP_5644 → / ← PP_3899 | membrane protein/hicB family protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,415,081 | 0 | G | A | 12.5% | 112.6 / 8.9 | 48 | intergenic (+146/+22) | PP_5644/PP_3899 | membrane protein/hicB family protein |
Reads supporting (aligned to +/- strand): ref base G (21/21); new base A (2/4); total (23/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.68e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.66e-01 |
GATGGCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATG‑CATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCTCACATTGAA > NC_002947/4414935‑4415225 | gATGGCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATG‑CATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGc < 1:597119/150‑1 (MQ=255) ggCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATG‑CATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCACTGCAATGCggg < 1:840187/150‑1 (MQ=255) ggCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATG‑CATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCACTGCAATGCggg < 1:711362/150‑1 (MQ=255) gCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATG‑CATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCttt > 2:433761/1‑150 (MQ=255) cACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATG‑CATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGAGCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTCCtttt > 1:52426/1‑149 (MQ=255) ccctcTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATG‑CATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGggc > 1:64424/2‑150 (MQ=255) cctcTTGAAGCCGCAATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATG‑CATTCACGCTTTTTTCACGCCCGCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGggcg < 2:622482/150‑1 (MQ=255) ttGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATG‑CATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTa > 1:201268/1‑150 (MQ=255) aGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATA‑CTGACATGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTAcgcg < 1:480325/150‑1 (MQ=255) ccgccATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATG‑CATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTATCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGtg > 1:168665/1‑150 (MQ=255) cATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATG‑CATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGggc < 2:795065/135‑1 (MQ=255) cATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATG‑CATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGggc > 1:795065/1‑135 (MQ=255) aaGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATG‑CATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCa > 2:615956/1‑150 (MQ=255) aGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGACATT‑ACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGTCGGCTACGCGTGATGACCCAg < 1:224777/150‑1 (MQ=255) aGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATG‑CATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAg > 1:446685/1‑150 (MQ=255) cTCGGCGGGCTTCACCATG‑CATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCttt < 1:11967/150‑1 (MQ=255) acaaagacTTCACCATG‑CATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGc < 1:759951/143‑1 (MQ=255) gcggGCTTCACCATG‑CATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCgt < 1:259194/150‑1 (MQ=255) ggCTTCACCATG‑CATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTc < 2:446685/150‑1 (MQ=255) ttCACCATG‑CATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGt < 1:572304/112‑1 (MQ=255) ttCACCATG‑CATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGt > 2:572304/1‑112 (MQ=255) g‑CATTCACGCTTTTTTCACGCCCTTCCCTGCTCAGTGCAGTCTCAACCTTTTCCCTACGTGAGTCCGCAGTGCAGTGCGGGCTTTTCTTTGGGCGGCCACGCGTGATTACCCAGGTACTCACTAAGCGCTTTGCGGGTCAGATCATTCaa < 2:52426/150‑1 (MQ=255) g‑CATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCCGATCATTCaa < 2:168665/150‑1 (MQ=255) cATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAg < 1:334719/150‑1 (MQ=255) aTTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTgg < 1:417604/90‑1 (MQ=255) aTTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTgg > 2:417604/1‑90 (MQ=255) gCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCttt > 2:34134/1‑124 (MQ=255) gCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCttt < 1:34134/124‑1 (MQ=255) gCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTtggtg > 1:273078/1‑150 (MQ=255) cccTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGcc > 1:694700/1‑150 (MQ=255) cccTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGcc > 2:395098/1‑150 (MQ=255) tGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCg > 2:535558/1‑150 (MQ=255) ctggCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAg > 2:525355/3‑77 (MQ=255) ctggCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAg < 1:525355/75‑1 (MQ=255) ggCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAg < 2:64424/150‑1 (MQ=255) aTCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCa > 2:366463/1‑150 (MQ=255) ccGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAg > 1:685844/1‑105 (MQ=255) ccGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAg < 2:685844/105‑1 (MQ=255) cTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCg > 1:631955/1‑150 (MQ=255) gaaaaGCCCGCACTGCAATGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGaaa > 1:786721/5‑98 (MQ=39) gaaaaGCCCGCACTGCAATGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGaaa < 2:786721/94‑1 (MQ=39) aaGCCCGCACTGCAATGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCaa > 1:760485/2‑90 (MQ=255) aaGCCCGCACTGCAATGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCaa < 2:760485/89‑1 (MQ=255) aTTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCAtt > 1:518279/1‑79 (MQ=255) aTTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCAtt < 2:518279/79‑1 (MQ=255) aTTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCTCACATTGa > 1:797966/1‑150 (MQ=255) aTTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCTCACATTGa < 1:433761/150‑1 (MQ=255) ttGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCTCACATTGaa < 1:535558/150‑1 (MQ=255) | GATGGCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATG‑CATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCTCACATTGAA > NC_002947/4414935‑4415225 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |