Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,604,344 | 0 | C | T | 43.1% | 49.7 / 33.9 | 51 | S62F (TCT→TTT) | bglJ | DNA‑binding transcriptional regulator BglJ |
Reads supporting (aligned to +/- strand): ref base C (5/24); new base T (22/0); total (27/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.69e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAT > NC_000913/4604207‑4604481 | aTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctt > 1:358451/1‑140 (MQ=255) tCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct < 1:1055300/140‑1 (MQ=255) caggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGt < 1:562691/138‑1 (MQ=255) gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAg > 2:1122734/1‑140 (MQ=255) aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATTAGAAGGGAgcg > 2:782960/1‑140 (MQ=255) aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAgcg > 1:548578/1‑140 (MQ=255) aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcg < 2:358451/140‑1 (MQ=255) tGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCa > 1:924657/1‑140 (MQ=255) tGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCa < 2:605113/140‑1 (MQ=255) tGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCa < 1:3153/140‑1 (MQ=255) tGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCa < 1:1122734/140‑1 (MQ=255) gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCag > 2:3153/1‑140 (MQ=255) gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCag > 2:1234567/1‑140 (MQ=255) gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCag > 1:939265/1‑140 (MQ=255) agagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACTTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCaga < 2:242511/140‑1 (MQ=255) agagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACTTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCaga > 1:242511/1‑140 (MQ=255) ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCagagactg > 2:562691/1‑137 (MQ=255) ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCATGAGAAGAGAGCGCAGAGAggg > 2:184637/1‑140 (MQ=255) ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAggg < 1:1103914/140‑1 (MQ=255) ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAggg < 1:713838/140‑1 (MQ=255) ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAggg < 1:782960/140‑1 (MQ=255) ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAggg > 1:377052/1‑140 (MQ=255) gAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCAGGGAGGGGGTATCt > 1:557472/1‑140 (MQ=255) aaaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCATGAGAAGTGAGCGCAGAGAGGGGTTATCtt > 1:467436/1‑140 (MQ=255) aaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTg < 1:665608/140‑1 (MQ=255) aCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTg > 1:959974/1‑140 (MQ=255) taAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTa < 1:555625/140‑1 (MQ=255) cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGtt < 2:377052/140‑1 (MQ=255) aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCAGAGAGGGGTTATCTTTCCTGACTGAACTGGC‑GGTTAAGTTTTcgcg > 2:909482/1‑140 (MQ=255) agacagaGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTAc < 1:71618/134‑1 (MQ=255) ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGGGGGGTTATTTTGCCTTACTGAACTGGC‑GATTAAGTTTTCGCGTAcc > 2:383081/1‑140 (MQ=255) cccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCg < 1:449908/139‑1 (MQ=255) ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGGAGGGAGCGCAGGGAGGGGGTATTTTTCCTGACTTAACTGGC‑GGTTTAGTTTCCGCGTACCCgg > 2:40111/1‑140 (MQ=255) gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGACATTTTTTCTTTTTCTGCCCTGGGGAGCGAGCGCCGTGGGGGGTTATCTTGCCTGACTTGAGTGGG‑GATTAAGTTTCCGCGTACCCGGccg > 2:449908/1‑138 (MQ=255) tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCGTTTTTTTCTTTTTTTTTCCAAATAAATGGACGCAAGGGGGGGGTATCTTGTCTTACTGAAATGGC‑GGATAAGGTTTCGGCTTCCCGGCGtt > 2:171509/1‑140 (MQ=255) cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGttt > 2:34355/1‑140 (MQ=255) aGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCCGACTGAACTGGC‑GATTAAAATTCCGCGTACCCGGCGtttt > 2:760578/1‑139 (MQ=255) ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAg < 1:171509/140‑1 (MQ=255) gAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTTAGAAGGGAGCGCAGAGAGGGGTTATCTTTCCTGACTGAACTGGC‑GATTAAGTTTCCGCGGACCCGGCGGTTTGt > 2:140368/1‑140 (MQ=255) aGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGGAAGGGGCGCAGAGAGGGGATTTATTTCCTTAGTTAAATTGCGGGTTAAGTTT‑CGCGTACCCCGCGGTTTGtt > 2:596893/1‑140 (MQ=255) cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTAtt < 1:1216339/140‑1 (MQ=255) aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTCTTCCCTGGGAAGGGAGGGGAGAGGGGGGGTTTTTTTCCTTGCTGGACTGGG‑GGTTAAGTTTTCGCGGACCCCGCGGTTTGTTAttt > 2:174323/1‑139 (MQ=255) ttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCg < 1:184637/140‑1 (MQ=255) aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCAGAGGGGGGTTATCTTTCCTGACTTAACTGGC‑GGTTAAGTTTCCGCGTACCCGGCGGTTAGTTTTTGCgggt > 2:458747/1‑138 (MQ=255) gATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatga < 2:924657/140‑1 (MQ=255) gATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatga < 1:383081/140‑1 (MQ=255) tGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTCTTTTTCTGCCTTGAGAGGCTAGTGTAGAGAGGGGATATTTTCGCTGACTGGACTGGC‑GATTAAGTTTCCGGGGACCCGGCGTTTAATTATTGCGggtgatg > 2:1068002/1‑140 (MQ=255) aTGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGatat < 2:939265/140‑1 (MQ=255) aTGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGatat < 2:557472/140‑1 (MQ=255) tGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCCTGAGAAGGGAGCGCAGAGGGGGGTTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTTGTTATTTCGGGTGATGATAtt > 2:1163034/1‑140 (MQ=255) ggATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCt < 1:140368/140‑1 (MQ=255) tctcGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCgg < 1:760578/140‑1 (MQ=255) cggcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGt < 1:1234567/140‑1 (MQ=255) gTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATTTTTCCTGACTGAACTGGC‑GATTTAGTTTTCTCGTACCCGGCGGTTTGTTATTTCGGATGATGATATTTAAGCTCGGCTGGTTGGTTCAtt > 2:1068564/1‑140 (MQ=255) ttttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTTAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGcc < 1:458747/140‑1 (MQ=255) ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCTCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt < 1:1068564/140‑1 (MQ=255) ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt < 2:548578/140‑1 (MQ=255) ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt < 2:739179/140‑1 (MQ=255) | ATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAT > NC_000913/4604207‑4604481 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |