Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,604,344 | 0 | C | T | 39.2% | 56.3 / 34.8 | 51 | S62F (TCT→TTT) | bglJ | DNA‑binding transcriptional regulator BglJ |
Reads supporting (aligned to +/- strand): ref base C (5/26); new base T (19/1); total (24/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.64e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGA‑GGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCA > NC_000913/4604205‑4604483 | aaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTtttttc > 1:75451/1‑140 (MQ=255) gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttttttt > 1:481211/1‑140 (MQ=255) tCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct < 1:142760/140‑1 (MQ=255) aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGc < 1:684803/140‑1 (MQ=255) aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGc < 2:730127/140‑1 (MQ=255) aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGc < 1:611932/140‑1 (MQ=255) taTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGaa > 1:8287/1‑140 (MQ=255) tCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGt < 2:267899/140‑1 (MQ=255) tCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGt < 2:195350/140‑1 (MQ=255) aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcg < 2:98874/140‑1 (MQ=255) aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcg < 2:540941/140‑1 (MQ=255) aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcg < 1:31553/140‑1 (MQ=255) tGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCa < 2:75451/140‑1 (MQ=255) agTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCATCAGAAATGGGCGCagagg > 2:31553/1‑139 (MQ=255) gTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagaga‑g < 2:596846/140‑1 (MQ=255) tATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGA‑GGGa < 2:481211/140‑1 (MQ=255) ccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGA‑GGGATTATCTTGCCTGACTGAACTGGc < 2:8287/140‑1 (MQ=255) taAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGA‑GGGATTATCTTGCCTGACTGAACTGGCGATTa > 1:711246/1‑140 (MQ=255) tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGATCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGA‑GGGATTATCTTGCCTGACTGAACTGGCGATTAAGttt > 2:600555/1‑140 (MQ=255) tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGGGAAGTGAGGGCCGAGA‑GGGGTTATTTTGCCTGACTGAACTGGCGATTAAGTTTCcgc > 2:405981/1‑140 (MQ=255) tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGA‑GGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAc < 1:187864/140‑1 (MQ=255) ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGA‑GGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAcc < 1:141519/140‑1 (MQ=255) taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCCTTTTTTTTTTTTTTTGGCATGTGAAGTGAGAGGAGTGA‑GGGGTTGTTTTGGCTGGGTTAACTAGCGACTAAGTTTCtgcctgtcc > 2:443452/1‑134 (MQ=255) cagaGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGA‑GGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCg < 1:75519/137‑1 (MQ=255) ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGA‑GGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCgg < 1:112016/140‑1 (MQ=255) gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCCTTAGAAAAGAGCGCCGAGA‑GGGGTTATTTTTTCCGACTGAACCGGGGGTTTAAATTCCGCGGCCCCcgcgg > 2:562953/1‑139 (MQ=255) gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGA‑GGGGTTATCTTGCCTGACTGAACTGGCGATTAAATTTCCGCGTACCCcgcgg > 2:305821/1‑139 (MQ=255) gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGAAGTGAGCGCAGAGA‑GGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGt > 2:280075/1‑140 (MQ=255) tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTTCACTTAGAAGGGAACGCAGAGC‑GGGATTTTCTTTGCTGACTGAAATTTGCATTAAATTTCCGTCTGCCCGGCgct > 2:112500/1‑138 (MQ=255) cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCAGAGA‑GGGGTTATCTTGCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCCGCGgtt > 2:154437/1‑140 (MQ=255) aGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTTCATGAGGAGGGAGCGCAGAGG‑GGGGTTTTTTTGCCTGACTGAAATGGCGGTTTAATTTTCGCCGACCCCGCGgttt > 2:187864/1‑139 (MQ=255) aGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCAGAGG‑GGGGGTATCTTTCCTGACTGAAATGGCGGTTAAGTTTCCCCGGACCCGGCGgttt > 2:397696/1‑139 (MQ=255) aGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGGGCCGAGA‑GGGGTTATCTTTCCTGACTGGACTGGCGATTAAGTTTTCGCGGACCCGGCGgttt > 2:711853/1‑139 (MQ=255) ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTTTTTCATGTGGAATGAGCGGTGAGAGGGGATAATC‑TGGCTGACTGAACATGCGTATAAATTTCCCCCAACCCGGCGGTTAg > 2:134378/1‑140 (MQ=255) tCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGGGCCGAGG‑GGGGTTATCTTTCCTGGCTGGAATGGCGGTTAAAATTTCTCGTACCCGGCGCTTTGtttt > 2:169978/1‑138 (MQ=255) aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCTTTGGAAGTGAGCGCGGAGG‑GGGGTTTTCTTGCCTGCCTGAACTGGCGGTTAAGTTTTCCCGCAACCCGCGTTTAGGTAttt > 2:141711/1‑139 (MQ=255) aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGGCGAGA‑GGGGTTATCTTTCCTGACTGAACTGGGGATTAAGTTTTCGCGTACCCCGCGTTTTGTTaatt > 2:75519/1‑137 (MQ=255) aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGA‑GGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTg < 1:711853/140‑1 (MQ=255) aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTCTTCTGCCATGAGAAGTGAGCGCAGAGA‑GGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCCCGTACCCGGCGTTTAGTTATTg < 2:371291/140‑1 (MQ=255) tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCAGAGG‑GGGGTTATCTTTCCTGACTGAACTGGCGATTAAGTTTTCGCGGACCCGGCGTTTTGTTATTgt > 2:239217/1‑139 (MQ=255) tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCCGAGA‑GGGGGTATCTTGCCTGACTGAACTGGCGGTTTAGTTTCCCCGTACCCCGCGGTTTGTTATTgg > 2:72790/1‑139 (MQ=255) ttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGA‑GGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCg < 1:305821/140‑1 (MQ=255) ttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGA‑GGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCg < 1:405981/140‑1 (MQ=255) cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGA‑GGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGga < 2:711246/140‑1 (MQ=255) cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGA‑GGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGga < 2:706310/140‑1 (MQ=255) gATGCCATGTCGCGGATCTCGTTTGCGGCGATCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGA‑GGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatga < 1:600555/140‑1 (MQ=255) aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGA‑GGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgat < 2:219054/140‑1 (MQ=255) aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGA‑GGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgat < 1:169978/140‑1 (MQ=255) gATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCAGGGA‑GGGGGTATTTTGCCTGACTGAAATGGGGGGTAAGTTTTCGCGGACCCGGCGTTTTGTTATTTTGGGTGATGATATTTAAGCTc > 2:623604/1‑140 (MQ=255) tctcGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGA‑GGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCgg < 1:280075/140‑1 (MQ=255) tcGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGA‑GGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCt < 1:72790/140‑1 (MQ=255) cggcggTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGG‑GGGGGTTTTTTTTCTGACTGAACTGGCGGTTTAGTTTTCGCGTACCCCGCGGTTTGTTATTTTGGATGATGATATTGAAGATCGGCTGATTTGt > 2:397559/1‑140 (MQ=255) cggcggTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGA‑GGGATTATATTTCCTGACTGAACTGGCGGATAAATTTCCGCGGACCCGGCGTTTAGTTATTTTGGATGGTGAAATTTAAGCTCCGCTTTTTGGt > 2:81860/1‑140 (MQ=255) cggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGA‑GGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCa < 1:562953/140‑1 (MQ=255) gataatttttcttttttTTCCATGAGAAGTTAGCGCAGCGA‑GGGATTATCTTTTCTGACTGAACCGGCGATTAAGATTTCGCGCTCCCCGCGGTTTGTTATTGCGGATGATGATATTTAAGCTCGGCTGATTGGTTCAtt > 2:527103/6‑140 (MQ=255) gTCATTTTTTTTTTTTTTTCCCTGAGGAATGAACGCAGAGA‑GGGATTATCTTGCCTGTCTTAACCGGCGAATAAGTTTTCGCGTACCCCGCGGTTTGTTATTGCGGGTGATGATATTGAAGGTCGGCTGATTGGTTCAtt > 2:667637/1‑140 (MQ=255) gTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCAGGGG‑GGGGTTATCTTTCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCGGCGGTTTGTTATTGCGGATGATGATATTGAAGCTCGGGTGATTTGTTCAtt > 2:448745/1‑140 (MQ=255) gTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGA‑GGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCAtt > 1:28192/1‑140 (MQ=255) ttctttttctGCCATGAGAAGTGAGCGCAGAGA‑GGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt < 1:255325/140‑1 (MQ=255) ttctttttctGCCATGAGAAGTGAGCGCAGAGA‑GGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt < 1:112500/140‑1 (MQ=255) ttctttttctGCCACGAGAAGTGAGCGCAGAGA‑GGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt < 2:561875/140‑1 (MQ=255) ctttttctGCCATGAGAAGTGAGCGCAGAGA‑GGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCa < 1:213738/140‑1 (MQ=255) | AAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGA‑GGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCA > NC_000913/4604205‑4604483 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |