Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,897,782 | G→C | 8.9% | intergenic (‑182/+52) | PP_5172 ← / ← PP_5173 | hypothetical protein/RND family transporter |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,897,782 | 0 | G | C | 8.9% | 140.9 / 5.5 | 55 | intergenic (‑182/+52) | PP_5172/PP_5173 | hypothetical protein/RND family transporter |
Reads supporting (aligned to +/- strand): ref base G (23/25); new base C (3/2); total (28/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.69e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.54e-01 |
TTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCG‑ATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATCCC > NC_002947/5897638‑5897926 | ttGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATAATGTTTAAGAATTTTTATATAACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTAACCCCGAATACTGTGGTCG‑ATTCAc > 5:19301/1‑148 (MQ=255) tGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCG‑ATGCACCa > 7:313981/1‑149 (MQ=255) gTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCG‑ATGCACCATc > 3:26603/1‑144 (MQ=255) gTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCG‑ATGCACCATc < 4:26603/144‑1 (MQ=255) ttttGCGCCTGGCTTCATGGCAATAGTAGTGAGGCTGGTTTATTATTAAAAATACTGTTTAAGAATTTTTATATAAACATTTAA‑GAGAACACTCCTCTGGTAGGAGAGGGTTTAACCGCGAATACGGTTGTCG‑ATGCACCATCGCATTc < 2:187270/149‑1 (MQ=255) tGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTTCCCCCGGAAAAGGTGGGCG‑GTGACCCCTCGCGTTAGCCGGTaaa > 6:9840/1‑149 (MQ=255) gCTTGAGTGCAATAGTAGTGAGGCGGTCTTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCG‑ATGCACCATCGCATTCGCGGGTAAAc < 5:253004/148‑1 (MQ=255) gCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGCATGCACCATCGCATTCGCGGGTAAAc < 1:179353/149‑1 (MQ=255) gCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCG‑ATGCACCATCGCATTCGCGGGTAAAc < 8:313981/148‑1 (MQ=255) gCTTGAGGGCAATAGTAGTGAGGCGGGATTAATCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCG‑ATGCACCATCGCATTCGCGGGTAAAcc < 3:86647/149‑1 (MQ=255) cTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCG‑ATGCACCATCGCATTCGCGGGTAAAccc < 8:42413/149‑1 (MQ=255) ttGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTTACCGCGGATAAGGTTGTGG‑ATTGACCCAAGCATTACCGGGTTAAACCg > 3:339138/1‑149 (MQ=255) aGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCG‑ATGCACCATCGCATTCGCGGGTAAACCCGCt < 3:145308/148‑1 (MQ=255) ggCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCG‑ATg > 6:349294/1‑118 (MQ=255) ggCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGCTTAAGAATTTTTATATAACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCG‑ATg < 5:349294/118‑1 (MQ=255) aaTAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACAATTTAA‑GAGAACACGCCTCTGGTAGGAGAGGGTTTACCCTCGAATGCGATGGTGC‑ATCGACCACCGTATTcg > 3:71787/1‑129 (MQ=21) gtagtGAGGCGGGATTATTCTTAAAAATACTTTTTAATAATTTTAATATAACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCG‑ATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAggg < 6:137754/149‑1 (MQ=255) tgtctcgCGGTATAATTCTAAAAAATA‑TTTTTAACAAATTTTAAAAAAACAATTAACCACACCCCGCCTTGGTTTGCTGCGGGTTTACCCGCTAATAGGGTGGTCT‑ATCCACCATCGCATTCGGGGTAAAACCCGCTCCTACCAgggg < 8:47879/142‑1 (MQ=255) aGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTTAACGCGAATAAGGTGGTCG‑ATGGACCAACGCATTAGCGGGTTAACCCGCGCCTAACAGGGGCCaaa > 6:42041/1‑147 (MQ=255) cGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCG‑ATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATgg < 5:132692/149‑1 (MQ=255) gggATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTAACCGCGAATAAGGTGGTCG‑ATGCACCATCGCATTAGCGGGTAAACCCGCTACTACCAGGGGCCATATggg > 3:242707/1‑149 (MQ=255) ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCG‑ATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCg > 4:115168/1‑149 (MQ=255) ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCG‑ATGCACCATAGCATTAGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCg > 1:138483/1‑149 (MQ=255) ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATAAGGTGGTCG‑ATGCACCATCGCATTAGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCg > 7:15691/1‑149 (MQ=255) tttAAGAATTTTTATATAACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCG‑ATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGc > 3:85345/1‑149 (MQ=255) tAAGAATTTTTATATAACCATTTAA‑GAGAAAACGCCTCTGGTAGGAGCGGGTTTAACCTACAATACGGTGCTCG‑ATGCACCATCGCATTCGCGGCTAAAACCGCTCATACCATGGTCCATATGCGCCTTCAGGCGTGTTTATGCCGCtg > 2:46594/1‑148 (MQ=255) ttttATATAACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCG‑ATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTcgcg < 7:313867/145‑1 (MQ=255) ttttATATAACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCG‑ATGCACCATCGCATTCGCGGGTAAACCCGCACCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTcgcg > 8:313867/1‑145 (MQ=255) tatataACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCG‑ATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGc < 2:91683/107‑1 (MQ=255) tatataACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATAAGGTGGTCG‑ATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGc > 1:91683/1‑107 (MQ=255) atataACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCG‑ATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCTGGCGTGTTTCTGCGGCTCGCGGATCTTg < 1:193188/148‑1 (MQ=255) atataACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCG‑ATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTg < 8:27435/148‑1 (MQ=255) atataACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCG‑ATGCACCATCGCATTAGCGGGTAAACCCGCTCCAACCAAGGGCCATATGGGCCTTAAGGCGTGTTTATTCGGCTACCGGATCTTg > 8:301423/1‑148 (MQ=37) tataACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCG‑ATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTa < 3:115168/149‑1 (MQ=255) tataACATTTTAA‑GAGACCAGGCCTGTGGTAGGAGCGTGTTTACCCGCTAATACGGTGGGCT‑ATGCACCATCGCATTCGGGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGt < 7:301423/148‑1 (MQ=37) aaCCATTTAA‑GAGAACACGCCTCTGTTAGGAGCGGGTTTACCCGCGAATACGGTGGTCG‑ATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCgtgt < 6:285801/123‑1 (MQ=255) aaCCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCG‑ATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCgtgt > 5:285801/1‑123 (MQ=255) aCCATTTAA‑GCCAACACGCCTCTGGTCGGAGCGGGTTTACCCGCTAATACGGTGGCCG‑ATGCACCATCGCATGCGCGGTTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAg < 8:15691/149‑1 (MQ=255) aCCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCG‑ATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAg > 3:90253/1‑149 (MQ=255) tttAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATAAGGTGGTAG‑ATGCACAATCGCAATCGCGGGTAAACCCGCTCCTACCAGGGGCCCTAAGGGCCTTCACGCGTGTTTCTGCGGCTCGCGGAGCTTGTACCAGCCCa > 6:271208/1‑149 (MQ=37) aa‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCG‑ATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACg > 5:298536/1‑148 (MQ=255) gagaACACGCCTCTGGTAGGAGCGGGTTTTACCGCGAATACGGGGGGGG‑ATGCACCATGGCATTAGCGGGTTAACCCGCGGCTACCAGGGGGCATATGGGCCTTAAGGCGTGTTTGTGCGGGTGGCGGCTCTTGTAACAGGCCACGTa > 8:46774/1‑148 (MQ=255) aaCACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCG‑ATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAgc < 4:85345/149‑1 (MQ=255) acGCCTCTGGTAGGAGCGGGTTTACGCGAGAATCCGGTGGTCG‑ATGCCCCATCGCATGCGCGGGTAAACCCGCTCCTACCAGGGGCCAGATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGcc < 5:242942/149‑1 (MQ=38) gTAGGTGCGGTTTTACCCGCTAATCCGGGGGTCG‑TTGCACCATCGAATTGGGGGGAAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGaa < 5:9840/148‑1 (MQ=255) ggTTTACCCGCGAATACGGTGGTCG‑ATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGt > 4:314984/1‑148 (MQ=255) ccGCGAATGCGATGGTGC‑ATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAgc > 5:214625/1‑121 (MQ=21) ccGCGAATGCGATGGTGC‑ATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAgc < 6:214625/121‑1 (MQ=21) ccGCGAATACGGTGGTCG‑ATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAgg > 4:221021/1‑149 (MQ=255) ccGCGAATACGGTGGTCG‑ATGCACCATAGCATTCGCGGGTAAACCCGCACCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCACGCGGATCATGTACCAGGCCACGTCCAGCGCCGGCAGGAACAGCAGCGTCAGCAAgg > 2:148436/1‑149 (MQ=255) cgcgAATGCGATGGTGC‑ATCGACCACCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt > 7:126673/1‑149 (MQ=37) gcgAATACGGTGGTCG‑ATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAgg < 2:157196/128‑1 (MQ=255) gcgAATACGGTGGTCG‑ATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAgg > 1:157196/1‑128 (MQ=255) gcgAATACGGTGGGCG‑ATGCACCATGGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTAGAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTc < 1:46594/149‑1 (MQ=255) ggtgggCG‑ATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACgat < 4:242707/149‑1 (MQ=255) gtCG‑ATGCACCATCGCATTCGCGTGTAAACCCGCTCCTACCAGGGGCCATAGGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATc < 5:271208/149‑1 (MQ=255) cG‑ATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGTCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATccc > 8:16750/1‑149 (MQ=255) | TTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAA‑GAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCG‑ATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATCCC > NC_002947/5897638‑5897926 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |