Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,619,711 | C→T | 7.3% | intergenic (+934/‑1326) | PP_5652 → / → PP_5654 | hypothetical protein/phosphotriesterase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,619,711 | 0 | C | T | 7.3% | 152.7 / 2.7 | 55 | intergenic (+934/‑1326) | PP_5652/PP_5654 | hypothetical protein/phosphotriesterase |
Reads supporting (aligned to +/- strand): ref base C (22/28); new base T (2/2); total (24/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01 |
AATTATTAGAGTATATAAAAAAAATTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGAAAAAATATTCTCTGAAACCAAGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCGGTATGAAGGCTCACATTGGTGTGGATGCGGAGTCGGTCTTGGTGCACTGCGTGGTAGGCACGTCGG > NC_002947/4619562‑4619859 | aaTTATTAGAGTATATAAAAAAAATTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGAAAAAATATTCTCTGAAACCAAGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCAc < 2:568293/150‑1 (MQ=255) taAAAAAAATTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGAAAAAATATTCTCTGAAACCAAGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGAcc < 1:267825/150‑1 (MQ=255) taAAAAAAATTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGAAAAAATATTCTCTGAAACCAAGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGAcc < 1:598369/150‑1 (MQ=255) aaaTTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGAAAAAATATTCTCTGAAACCAAGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAAc < 1:550014/150‑1 (MQ=255) cgtgattACCCGCTCCCCAACAATACACGAGAGAAATTGAAAAAATATTCTCTGAAACCAAGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCaa > 2:1604122/7‑150 (MQ=255) cccAACAATACACGAGAGAAATTGAAAAAATATTCTCTGAAACCAAGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCAGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGcc > 1:994226/1‑150 (MQ=255) cacGAGAGAAATTGAAAAAATATTCTCTGAAACCAAGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGc < 1:155456/150‑1 (MQ=255) acGAGAGAAATTGAAAAAATATTCTCTGAAACCAAGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCa > 2:102976/1‑150 (MQ=255) cGAGAGAAATTGAAAAAATATTCTCTGAAACCAAGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCAc < 1:1212480/150‑1 (MQ=255) cGAGAGAAATTGAAAAAATATTCTCTGAAACCAAGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCAc < 2:1923851/150‑1 (MQ=255) ttGAAAAAATATTCTCTGAAACCAAGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAg > 1:308477/1‑150 (MQ=255) aaaaaaTATTCTCTGAAACCAAGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACCTAAATGCACCAAGCTAAGaaa > 2:642556/1‑150 (MQ=255) aaaaTATTCTCTGAAACCAAGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAgg < 1:1750817/150‑1 (MQ=255) aaaaTATTCTCTGAAACCAAGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACAAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAgg < 1:520599/150‑1 (MQ=255) ttCTCTGAAACCAAGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCg < 1:1604122/150‑1 (MQ=255) tctGAAACCAAGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCggtg < 1:1246725/150‑1 (MQ=255) tctGAAACCAAGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACACAAATGCACCAAGCCAAGAAAGGCAATCggta > 2:396553/1‑149 (MQ=255) aaaCCAAGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCACGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTAc > 2:1814256/1‑150 (MQ=255) ccAAGCCCTCGTGAGTAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTc < 1:1405529/150‑1 (MQ=255) ccAAGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCAGCGGTACCTc > 2:663043/1‑150 (MQ=255) aaGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCATGCGCCGAGTTCGACCGAGAACCAAGACAGCAAgcgc > 2:1396597/1‑103 (MQ=255) aaGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCATGCGCCGAGTTCGACCGAGAACCAAGACAGCAAgcgc > 2:1396347/1‑103 (MQ=255) aaGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCATGCGCCGAGTTCGACCGAGAACCAAGACAGCAAgcgc < 1:1396597/103‑1 (MQ=255) aaGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCATGCGCCGAGTTCGACCGAGAACCAAGACAGCAAgcgc < 1:1396347/103‑1 (MQ=255) tGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCGGTATGAAGGc < 1:1721019/150‑1 (MQ=255) tGAGCAAAAGCGTCATTCGACGCTGGATGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCGGTATGAAGGc < 2:408757/150‑1 (MQ=255) aTGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCGGTATGAATGCTCACAttt < 2:840662/150‑2 (MQ=255) tCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCGGTATGAAGGCTCACATTGgtgt > 2:1008584/1‑150 (MQ=255) aTTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCGGTATGAAGGCTCACATTGGTGTgg > 2:352213/1‑150 (MQ=255) ttCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGAGCACGATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCGGTATGAAGGCTCACATTGGTGTGGa > 1:595192/1‑150 (MQ=255) ggCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCGGTATGAAGGCTCACATTGgt < 1:722492/142‑1 (MQ=255) ggCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCGGTATGAAGGCTCACATTGgt > 2:722492/1‑142 (MQ=255) ggCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCGGTATGAAGGCTCACATTGGTGTGGATGc < 2:1391327/150‑1 (MQ=255) ggaggaAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCGGTATGAAGGCTCACATTGGTGTGGATGCGGAGTc > 1:1845956/1‑150 (MQ=255) ggaggaAACGCTGTACGAGGCCACGCAGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCGGTATGAAGGCTCACATTGGTGTGGATGCGGAGTc < 2:994226/150‑1 (MQ=255) gaAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCGGTATGAAGGCTCACATTGGTGTGGATGCGGAGTCGGTc > 2:514823/1‑150 (MQ=255) aaCGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCTGTATGAAGGCTCACATTGGTGTGGATGCGGAGTCGGTCtt < 2:121865/150‑1 (MQ=255) aCGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCGGTATGAAGGCTCACATTGGTGTGGATGCGGAGTCGGTCTTg < 1:102976/150‑1 (MQ=255) aCGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCGGTATGAAGGCTCACATTGGTGTGGATGCGGAGTCGGTCTTg < 2:308477/150‑1 (MQ=255) gCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCACTCGGTGGTACTTCGGTATGAAGGCTCACATTGGTGTGGATGCGGAGTCGGTCTCggg > 1:739896/1‑149 (MQ=255) gTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCGGTATGAAGGCTCAAATTGGTGTGGATGCGGAGTCGGTCTTGGTGCa < 1:252032/150‑1 (MQ=255) gAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCGGTATGAAGGCTca < 2:208169/113‑1 (MQ=255) gAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCGGTATGAAGGCTca > 1:208169/1‑113 (MQ=255) gAGGCCACGCGGATCAACGCGCCGAGTTCGACTGAGAACCAAGACAGCAAGCGCGACACATAAATGCACCAAGCTAAAAAAGGCAATCGGTGGTACTTCGGTATGAAGGCTAACATTGGTGTGAATGCGGAATCGGTCTTGGTGCACTGc < 2:863027/150‑1 (MQ=255) gCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCGGTATGAAGGCTCACATTGgtgt < 2:1637357/119‑1 (MQ=255) gCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCGGTATGAAGGCTCACATTGgtgt > 1:1637357/1‑119 (MQ=255) gCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCGGTATGAAGGCTCACATTGGTGTGGATGCGGAGTCGGTCTTGGTGCACTGCGTg < 1:1065763/150‑1 (MQ=255) ccACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCGGTATGAAGGCTCACATTGGTGTGGATGCGGAGTCGGTCTTGGTGCACTGCGTgg > 1:985175/1‑150 (MQ=255) cACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCGGTATGAAGGCTCACATTGGTGTGGATGCGGAGTCGGTCTTGGTGCACTGCGTGGt > 2:1439790/1‑150 (MQ=255) cACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCGGTATGAAGGCTCACATTGGTGTGGATGCGGAGTCGGTCTTGGTGCACTGCGTGGt > 1:1181007/1‑150 (MQ=255) aCGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCgg < 2:1393488/95‑1 (MQ=255) aCGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCgg > 1:1393488/1‑95 (MQ=255) cgcgGGTCCAAGCGGCGAGTCCGACCGTGAACCAAGGCAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTTGGTATGGCGGGTCACATTGGTGTGGATGCGGAGTCGGTCTTGGGGCACTGGGTGGGAg < 2:739896/150‑1 (MQ=255) gcgGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAACAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCGGTATGAAGGCTCACATTGGTGTGGATGCGGAGTCGGTCTTGGTGCACTGCGTGGTAgg > 2:899486/1‑150 (MQ=255) ccACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCGGTATGAAGGCTCACATTGGTGTGGATGCGGAGTCGGTCTTGGTGCACTGCGTGGTAGGCACGTc < 1:1428956/150‑1 (MQ=255) ccACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCGGTATGAAGGCTCACATTGGTGTGGATGCGGAGTCGGTCTTGGTGCACTGCGTGGTAGGCACGTc < 1:1439419/150‑1 (MQ=255) ccgcgCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCGCCAAGCTAAGAAAGGCAATCGGTGGTACTTCGGTATGAAGGCTCACATTGGTGTGGATGCGGAGACGGTCTTGGTGCACTGCGTGGTAGGCACGTAgg > 2:1151256/2‑150 (MQ=255) | AATTATTAGAGTATATAAAAAAAATTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGAAAAAATATTCTCTGAAACCAAGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTTCGACCGAGAACCAAGACAGCAAGCGCGCCACATAAATGCACCAAGCTAAGAAAGGCAATCGGTGGTACTTCGGTATGAAGGCTCACATTGGTGTGGATGCGGAGTCGGTCTTGGTGCACTGCGTGGTAGGCACGTCGG > NC_002947/4619562‑4619859 |
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: ‑ |