Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,642,313 | G→C | 13.2% | intergenic (‑99/+53) | PP_1442 ← / ← lon‑I | hypothetical protein/DNA‑binding ATP‑dependent protease |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,642,313 | 0 | G | C | 13.2% | 154.5 / 12.5 | 53 | intergenic (‑99/+53) | PP_1442/lon‑I | hypothetical protein/DNA‑binding ATP‑dependent protease |
Reads supporting (aligned to +/- strand): ref base G (23/23); new base C (2/5); total (25/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.26e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01 |
CAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCG > NC_002947/1642166‑1642459 | cAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGATGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAggg < 1:191935/150‑1 (MQ=255) gAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGt > 2:387064/1‑150 (MQ=255) cAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGt < 2:351369/150‑1 (MQ=255) gggAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCg < 1:358418/150‑1 (MQ=255) ggAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCACGGGTACGC‑TACCCCTGCAGGCGCgg < 2:351045/150‑1 (MQ=38) gAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCggg < 2:322409/150‑1 (MQ=255) gAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTg > 2:52911/1‑150 (MQ=255) aGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGc < 2:603998/150‑1 (MQ=255) gaggagAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAggg > 2:353624/1‑114 (MQ=255) gaggagAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAggg < 1:353624/114‑1 (MQ=255) gagGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAggg < 1:241247/150‑1 (MQ=255) agGTAATGGTGGCGCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGGGCAAGCCCGCACCTGCACGGGTA‑CGCTACCCCTGCAGGCGCGGGCTTACTCGCGaaaggacc < 1:93337/150‑7 (MQ=17) tCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGt < 1:613479/150‑1 (MQ=255) gtggCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAgg < 2:373801/150‑1 (MQ=255) cTCCGAAGGCTGGCGATGGCCTGAGGATATCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAg > 1:314456/1‑150 (MQ=255) cTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGAGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAg > 1:292760/1‑150 (MQ=255) cGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAAc > 1:467473/1‑150 (MQ=255) gAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACa > 2:615288/1‑150 (MQ=255) tGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACtt > 1:22112/1‑150 (MQ=255) tGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACtt > 1:163945/1‑150 (MQ=255) gATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGc < 1:52911/150‑1 (MQ=255) aTACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCaaa < 1:562949/150‑1 (MQ=255) gAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGaaa > 1:500876/1‑150 (MQ=255) gAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCg > 2:611965/1‑97 (MQ=255) gAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCg < 1:611965/97‑1 (MQ=255) gAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGc > 2:375941/1‑150 (MQ=255) aaaTTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCa < 1:387064/150‑1 (MQ=255) aaaTTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCa < 2:22112/150‑1 (MQ=255) aaTTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCAc < 2:500876/150‑1 (MQ=255) gCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCgcgc > 1:253912/1‑150 (MQ=255) ccGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCa < 1:185899/150‑1 (MQ=255) cGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAaga > 2:81387/1‑72 (MQ=255) cGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAaga < 1:81387/72‑1 (MQ=255) gTCCTTTCGCGAGTCAGACCTCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGt > 2:458560/1‑150 (MQ=255) cTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTc > 2:210491/1‑95 (MQ=255) cTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTc < 1:210491/95‑1 (MQ=255) cTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGt > 2:542222/1‑150 (MQ=255) cTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGt > 1:247324/1‑150 (MQ=255) tttCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTc < 1:615288/150‑1 (MQ=255) tttCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTc > 1:527732/1‑150 (MQ=255) tttCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTc > 2:164800/1‑150 (MQ=255) ttCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTcc < 2:467473/150‑1 (MQ=255) cgAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCa < 2:543664/150‑1 (MQ=255) gcaaGCCCGCACCTGCACGGGTACGC‑TACCCCTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCa < 1:34711/145‑1 (MQ=38) gcaaGCCCGCACCTGCACGGGTA‑CGCTACCCCTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCa > 2:34711/3‑147 (MQ=38) gCCCGCGCCTGCAGGGGTAGCG‑TACCCGAGCAGGTGCGGGCTTGCCCGCGAAGAAGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGttcttc > 1:22197/1‑150 (MQ=255) cccGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCg < 1:358348/150‑1 (MQ=255) ccGCGCCTGTAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTcgc > 1:74019/1‑126 (MQ=255) ccGCGCCTGTAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTcgc < 2:74019/126‑1 (MQ=255) ccGCACCTGCACGGGTACGC‑TACCCCTGCAGGCGCGGGCTTACTCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGttcttc > 2:215455/1‑148 (MQ=25) ccGCACCTGCACGGGTACGC‑TACCCCTGCAGGCGCGGGCTTACTCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGttcttc < 1:215455/148‑1 (MQ=25) gcgcCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAg > 1:214064/1‑150 (MQ=255) cTGCACGGGTACGC‑TACCCCTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGttcttc < 2:429894/142‑1 (MQ=39) tGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTcc > 2:112503/1‑130 (MQ=255) tGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTcc < 1:112503/130‑1 (MQ=255) gCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTcgccg < 2:163945/150‑1 (MQ=255) | CAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCG > NC_002947/1642166‑1642459 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |