Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,594,167 | 0 | C | T | 60.4% | ‑2.7 / 43.9 | 48 | S64F (TCT→TTT) | bglJ | DNA‑binding transcriptional activator BglJ |
Reads supporting (aligned to +/- strand): ref base C (5/14); new base T (29/0); total (34/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.41e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-08 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGC > NZ_CP009273/4594078‑4594254 | cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATtttttctttttc < 1:598413/90‑1 (MQ=255) cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTg < 2:251035/90‑1 (MQ=255) tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGcc < 1:40478/90‑1 (MQ=255) gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTGCCATgag > 2:285758/1‑90 (MQ=255) cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTTCCCTgaga > 1:359855/1‑90 (MQ=255) cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTTCCATtaga > 1:546578/1‑90 (MQ=255) cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATgaga > 1:591775/1‑90 (MQ=255) cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑GCATTTTTTTTTTTTTTTCCATgaga > 1:621520/1‑90 (MQ=255) ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTTCCCTGAGAAg > 2:959844/1‑90 (MQ=255) ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTCTGcc > 1:512243/1‑82 (MQ=255) ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGcc < 2:512243/82‑1 (MQ=255) ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAg < 2:591775/90‑1 (MQ=255) cATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGt < 1:353984/90‑1 (MQ=255) tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTTGCATGAGAAGTGAgcgc > 1:569570/1‑90 (MQ=255) tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTGCCCTGAGGAGGGAgcgc > 2:838020/1‑90 (MQ=255) tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTGCCCTGAGAAGTGAgcgc > 2:400034/1‑90 (MQ=255) tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTGCCATGAGAAGGGAgcgc > 2:222403/1‑90 (MQ=255) ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTTCCATGAGAAGGGAgcgcc > 2:481510/1‑89 (MQ=255) ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCa > 2:262833/1‑90 (MQ=255) cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagaga < 1:639292/90‑1 (MQ=255) aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCGCATTTTTT‑TTTTTTTGCCCTAAGGAGGGAGAGCcgagag > 2:569773/1‑90 (MQ=255) aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGGGCagaggg > 1:427947/1‑88 (MQ=255) aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCagagag > 1:260400/1‑90 (MQ=255) aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCagagag > 2:934253/1‑90 (MQ=255) aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑GCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGGGCaggggg > 1:865860/1‑86 (MQ=255) gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGGGCAGAGAgg > 1:270202/1‑90 (MQ=255) gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGAGCAGAGAgg > 2:765132/1‑90 (MQ=255) gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTGCCATGAGAAGTGGGGGCAGAGAgg > 2:383725/1‑90 (MQ=255) gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTGCCATGAGAAGGGAGGGCAGAGAgg > 2:409391/1‑90 (MQ=255) gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCGGAGAgg > 1:881721/1‑90 (MQ=255) gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTCTGCCATGAGAAGGGAGCGCAGAGAgg > 2:337067/1‑90 (MQ=255) tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGGGCAGAGAggg > 2:71725/1‑90 (MQ=255) ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTTCCATGAGAAGGGAGGGCAGAGAGGGGtt > 2:721646/1‑90 (MQ=255) ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGGGGAGAGAGGGGtt > 2:884495/1‑90 (MQ=255) gTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATc < 2:919663/90‑1 (MQ=255) aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑GCATTTTTTTTTTTTTTTCCCTGGGAAGTGAGGGCAGAGAGGGGTTATCCTg > 1:548550/1‑90 (MQ=255) ggATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgac < 1:409391/90‑1 (MQ=255) gATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgact < 1:583864/90‑1 (MQ=255) aTGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTGCCATGAGAAGGGAGGGCAGAGAGGGGTTATTTTTCctgactg > 2:747628/1‑90 (MQ=255) tGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCGGAGGGGGGTTATCTTGTctgactga > 1:799202/1‑90 (MQ=255) ccATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTTCCTGACTGaaa > 2:916541/1‑89 (MQ=255) ccATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGGAGTGGGGGCAGAGAGGGAGTATCTTTCCTTACTGaaa > 2:615082/1‑89 (MQ=255) cATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACt < 2:244467/90‑1 (MQ=255) tctcGTTTGCGGCGG‑TCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGGTTATTTTTCCTGACTGAACTGGCGATTTAGt > 1:353827/1‑90 (MQ=255) tcGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGttt > 2:901908/1‑90 (MQ=255) cgg‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAccc < 1:672208/90‑1 (MQ=255) gg‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCg < 1:912074/90‑1 (MQ=255) gg‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCg < 2:865860/90‑1 (MQ=255) aTTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGt < 1:959844/90‑1 (MQ=255) ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGt < 1:945901/90‑1 (MQ=255) ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGt < 1:522582/90‑1 (MQ=255) ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGt < 1:270193/90‑1 (MQ=255) tttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTg < 2:513060/90‑1 (MQ=255) ttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCTTACCCGGCGTTTAGTTATTGc < 1:481510/90‑1 (MQ=255) ttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGc < 1:883719/90‑1 (MQ=255) ttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGc < 1:838020/90‑1 (MQ=255) | CCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGC > NZ_CP009273/4594078‑4594254 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |