Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,604,344 | 0 | C | T | 41.2% | 46.6 / 30.0 | 51 | S62F (TCT→TTT) | bglJ | DNA‑binding transcriptional regulator BglJ |
Reads supporting (aligned to +/- strand): ref base C (4/26); new base T (21/0); total (25/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.11e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCG‑ATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCA > NC_000913/4604209‑4604483 | gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑Ttttttctttt > 1:1082327/1‑140 (MQ=255) cGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑tttttttttttt > 2:5133/1‑140 (MQ=255) cGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑Ttttttcttttt < 2:1082327/140‑1 (MQ=255) cGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTAAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑Ttttttcttttt < 1:759573/140‑1 (MQ=255) gTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATA‑CAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTttttctttttct < 2:1127658/140‑1 (MQ=255) tCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTttttctttttct < 2:698093/140‑1 (MQ=255) aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGc < 2:155658/140‑1 (MQ=255) aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGc < 2:32066/140‑1 (MQ=255) gtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATgaga < 2:475249/140‑1 (MQ=255) tCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATTCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGt < 1:207596/140‑1 (MQ=255) tCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGt < 2:927598/140‑1 (MQ=255) cGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTTCCCTTAGaaagg > 2:1135720/1‑137 (MQ=255) gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTTCCATTAGAagagag > 2:207596/1‑140 (MQ=255) gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAg < 1:790849/140‑1 (MQ=255) aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTTTTCCATGAGAAGTGAgcg > 2:599181/1‑140 (MQ=255) aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAgcg < 2:442011/140‑1 (MQ=255) aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGGCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCA‑TTTTTTCTTTTTTTTCCGTAAGAAGTGaagg > 2:489922/1‑137 (MQ=255) gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑ATTTTTTTTTTTTTGTCCTGAGAAGGGAgggcgg > 2:234731/1‑138 (MQ=255) gagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagag < 2:801732/140‑1 (MQ=255) gagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagag < 2:7182/140‑1 (MQ=255) ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGAggg > 2:539702/1‑140 (MQ=255) aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTGCCCTGAGAAGGGGGCGCCGAGAGGGAt > 1:819816/1‑140 (MQ=255) tCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTa < 1:530263/140‑1 (MQ=255) aaaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTGCCTTG‑GAAGTGGGCGCCGCGAGGGGTTAtttt > 1:581898/1‑140 (MQ=255) aaGTTTGCGGGTAACCCTTATAAGCTCCATTCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTg < 2:1077076/140‑1 (MQ=255) tttGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCct < 1:357610/140‑1 (MQ=255) tttGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCct < 1:1098663/140‑1 (MQ=255) tAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAc < 2:539458/140‑1 (MQ=255) cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTTGCG‑ATTAAGtt > 2:974335/1‑140 (MQ=255) tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTCTTCCCTGAGAAGTGAGGGCAGAGGGGGGTTTTCTTTCCTGACTGAACTGGCG‑GTTAAGttt > 2:716027/1‑140 (MQ=255) tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCACGAGAAGTGAGCGCAGAGAGGGAATATCTTGGCTGGCTGTACTGGCG‑ATTAAGGTTTcgc > 2:545334/1‑140 (MQ=255) aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCG‑ATTAAGTTTCcgcg < 1:441522/140‑1 (MQ=255) ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCG‑ATTAAGTTTCCGCGTAcc > 1:713008/1‑140 (MQ=255) taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCAGAGAGGGGTTATTTTTCCTGACTGAACTGGCG‑GTTAAAGTTTCGCGTAccc > 2:94925/1‑140 (MQ=255) taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCG‑ATTAAGTTTCCGCGTAccc < 1:1186745/140‑1 (MQ=255) gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTGCCATGAGAAGTGAGGGCAGAGGGGGGTTATCTTGCCTGACTGAACTGGCG‑ATTAAGTTTTCGCGTACCCcgcgg > 2:834124/1‑139 (MQ=255) gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCG‑ATTAAGTTTCCGCGTACCCGGCGt < 1:1053685/140‑1 (MQ=255) gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCT‑TTTTTTTTTTTTTTGCCATGATGAATGATACCAGAGAGGGGTTATCTCTTCTTACTGGAGTGTCGCATTAATTTTCCCCGATCCCGGCg > 2:1017372/1‑140 (MQ=255) tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTTCCATGAGAAGTGGGGGCAGAGAGGGATTATCTTTCCTGACTGAACTGGCG‑AATAAGTTTTTGCGTACCCGGCGtt > 2:605992/1‑140 (MQ=255) tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTGCCCTGTGAAGGGGGGGCAGAGAGGGGTTAACTTTCCGTACTGAACCGGGG‑ATTTAGTTTTCGCGTACCCcgcggt > 2:522914/1‑138 (MQ=255) cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCC‑TTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCAGAGAGGGGGTATCTTTCCTGACTGAACTGGCG‑GTTAAATTTTCGCGGACCCGGCGgtt > 2:203420/1‑140 (MQ=255) cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCC‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCG‑ATTAAGTTTCCGCGTACCCGGCGttt > 2:612407/1‑140 (MQ=255) cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCA‑ATTTTTTTTTTTTTACCATCAGAAAGTAACGGCGAGAGGGACTTTCTTACCCGGCAGAACTTGCG‑AGTTAATTTCCCCGTTCCCCGccgtt > 2:257409/1‑136 (MQ=255) tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCG‑ATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGc > 1:971477/1‑140 (MQ=255) ttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCG‑ATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCg < 1:1024886/140‑1 (MQ=255) gATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCG‑ATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatga > 2:1058775/1‑140 (MQ=255) aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTTCCATGAGAAGTGAGCGGAGAGAGGGGTTATTTTGCCTGACTGCACTGGCG‑ATTAAGTTTTCGCGTACCCGGCGGTTAGTTATTGCGggtgag > 2:786151/1‑139 (MQ=255) gCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTGCCATGAGGAGGGAGCGCAGAGAGGGGTTTTTTTTCCTGACTTAAATGGGG‑ATTTAGTTTTCGCGGACCCGGCGGTTTGTTATTGCGggtggtga > 2:559187/1‑140 (MQ=255) gCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTCCGCCTTTAGGAATGAGGGGAGGGGGGGGTTATTTTTCCGGCCTTAAAAAAGG‑AATAAGCTTGCCCGCACCCGGCGTTTAAATGTTGCgggttgtga > 2:116106/1‑133 (MQ=255) cATGTCGCGGATCTCGTTTGCGGCGGTCC‑TTTTTTTTTTTTTTTCCCTGAGAAGTGGGGGCAGAGAGGGGTTATTTTTCCTTACCGAACTGGCG‑GATAAGTTTTCGCCGACCCCGCGTTTTGTTTTTGCGgatgatgatt > 2:855251/1‑139 (MQ=255) aTGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCG‑ATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGatat < 2:363556/140‑1 (MQ=255) gcgGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCG‑ATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAg < 2:679741/140‑1 (MQ=255) gATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCG‑ATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTc < 2:581898/140‑1 (MQ=255) tctcGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCG‑ATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCgg < 1:94925/140‑1 (MQ=255) tctcGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCG‑ATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCgg < 1:1017372/140‑1 (MQ=255) gTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGACCTGGCG‑ATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGa < 1:116106/140‑1 (MQ=255) gTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCG‑ATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGa < 1:215610/140‑1 (MQ=255) gTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCG‑ATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCAtt < 1:257409/140‑1 (MQ=255) ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCG‑ATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt < 1:153451/140‑1 (MQ=255) ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCG‑ATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt < 1:612407/140‑1 (MQ=255) ctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCG‑ATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCa < 1:1135720/140‑1 (MQ=255) | GCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCG‑ATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCA > NC_000913/4604209‑4604483 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |