Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 6,040,850 | T→A | 21.8% | intergenic (+19/+15) | PP_5292 → / ← PP_5293 | catabolite repression control protein/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 6,040,850 | 0 | T | A | 21.8% | 117.0 / 29.8 | 55 | intergenic (+19/+15) | PP_5292/PP_5293 | catabolite repression control protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (20/23); new base A (5/7); total (25/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01 |
GTTCGACTACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGCCCGATAAACACGAACAT > NC_002947/6040705‑6040999 | gTTCGACTACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCa < 2:399099/150‑1 (MQ=255) aCCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCAATGCAGCCCCAAt < 1:1199359/150‑1 (MQ=255) ccAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCAtt > 1:695286/1‑138 (MQ=255) ccAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCAtt < 2:695286/138‑1 (MQ=255) ccAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATg > 1:495311/1‑150 (MQ=255) ccAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCAATGCAGCCCCAATg < 1:229062/150‑1 (MQ=255) tCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTg > 1:265028/1‑150 (MQ=255) cccGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCg > 1:1081853/1‑150 (MQ=255) cGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGgc < 1:291118/150‑1 (MQ=255) gccgcTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCAATGCAGccc < 2:483161/122‑1 (MQ=255) gccgcTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCAATGCAGccc > 1:483161/1‑122 (MQ=255) gcggcgcAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGt < 1:1102946/147‑1 (MQ=255) cGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGGTTATCCCCGCGTCCCCATCTGATCCTCGACGATTACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGAAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCgtgt < 2:265028/150‑1 (MQ=255) gcAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAAt < 2:598065/114‑1 (MQ=255) gcAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAAt > 1:598065/1‑114 (MQ=255) gcAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGa > 1:1142122/1‑150 (MQ=255) gcAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCAAGAATGCGATGCTGTTGACATACTATCACTGGACGCTTACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGa > 1:22344/1‑150 (MQ=255) cAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGAt > 2:80156/1‑150 (MQ=255) aCGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATaa < 2:1081853/150‑1 (MQ=255) gCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTg > 2:257204/1‑150 (MQ=255) ccTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGc < 2:1179323/150‑1 (MQ=255) cTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCAATGCAGCCCCAAt < 1:284144/103‑1 (MQ=255) cTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCAATGCAGCCCCAAt > 2:284144/1‑103 (MQ=255) tGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGccc > 2:525895/1‑150 (MQ=255) ctcCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATgg > 1:568978/1‑85 (MQ=255) ctcCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATgg < 2:568978/85‑1 (MQ=255) tcCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACcgccgc > 1:514466/1‑150 (MQ=255) cAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCgat < 2:514466/150‑1 (MQ=255) gCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCACCGGGTAGCGATAAGACTTGCCCTCATTCGCCCGCACCGCCGCgatga > 2:436246/1‑150 (MQ=255) ccGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCaa < 1:506639/150‑1 (MQ=255) cGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAat > 2:289797/1‑150 (MQ=255) aTCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATatcat < 1:257204/150‑1 (MQ=255) aTCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATatcat < 2:1142122/150‑1 (MQ=255) tcgtcgACTATGACTTGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATatcatc < 2:1123839/150‑1 (MQ=255) gACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGt < 2:715284/150‑1 (MQ=255) tGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAg > 1:644237/1‑150 (MQ=255) aCGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCAt < 1:775778/150‑1 (MQ=255) aCGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCAt < 2:940899/150‑1 (MQ=255) cTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCa < 2:495311/150‑1 (MQ=255) cATTGGGGCTGCAATGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGc < 2:884799/72‑1 (MQ=255) cATTGGGGCTGCAATGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGc > 1:884799/1‑72 (MQ=255) aTTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCAt < 2:898388/129‑1 (MQ=255) aTTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCAt > 1:898388/1‑129 (MQ=255) ttGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGccc < 2:282554/146‑1 (MQ=255) ttGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGccc > 1:282554/1‑146 (MQ=255) ttGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGCCCGATa > 2:845518/1‑150 (MQ=255) ttGGGGCTGCAATGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATg > 1:1014629/1‑129 (MQ=255) ttGGGGCTGCAATGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATg < 2:1014629/129‑1 (MQ=255) ggggCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCCTTTGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGCCCGATaaa < 2:644237/150‑1 (MQ=255) ggCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGATCATGGCAAACAGCACCAGCCCGATAAaca < 1:525895/150‑1 (MQ=255) gatGCAATGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATgg < 1:604874/123‑1 (MQ=255) gCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGCCCGATAAacac > 2:1202841/1‑150 (MQ=255) gCTGCAATGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATgg > 2:604874/1‑125 (MQ=255) tGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACt > 2:762591/1‑107 (MQ=255) tGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACt < 1:762591/107‑1 (MQ=255) tGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGCCCGATAAACACGa > 1:714980/1‑150 (MQ=255) gCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGCCCGATAAACACGaa < 1:289797/150‑1 (MQ=255) ttGCAGCACCAATCGTTATTTGATCGGCCGCCAGATCATAGGCTAGCGAGAAGGCTTCCATTGATTCGCGCGCAGCGCCGCGGTGATGAGCAATCTCATCCCGGACCCCGCGAGGCAGGCACAGCGCACCAGCCCGCTCAACACGAACAt > 2:96063/1‑150 (MQ=255) | GTTCGACTACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGCCCGATAAACACGAACAT > NC_002947/6040705‑6040999 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |