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.0% | 38.1 / 24.0 | 41 | S62F (TCT→TTT) | bglJ | DNA‑binding transcriptional regulator BglJ |
Reads supporting (aligned to +/- strand): ref base C (6/19); new base T (16/0); total (22/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.33e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATT > NC_000913/4604207‑4604473 | aTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCCTtttttttt > 2:974144/1‑140 (MQ=255) gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTtttttttttt > 1:188946/1‑140 (MQ=255) tCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct < 1:910539/140‑1 (MQ=255) aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGc < 1:645627/140‑1 (MQ=255) aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGc < 2:448687/140‑1 (MQ=255) aTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAg < 1:663722/140‑1 (MQ=255) tCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGt < 1:339379/140‑1 (MQ=255) gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAg > 1:601019/1‑140 (MQ=255) gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCCTTTTTTTTTTTTTTTTTCATGAAAAGTGAg > 2:249512/1‑140 (MQ=255) tATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgc < 1:500865/140‑1 (MQ=255) aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGGAGTGAgcg > 2:275075/1‑140 (MQ=255) aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcg < 1:664204/140‑1 (MQ=255) gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGGgcggcg > 2:772805/1‑137 (MQ=255) agTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGAGAGCGCagagg > 2:833863/1‑139 (MQ=255) tttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGAgg > 1:374364/1‑140 (MQ=255) ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAggg > 1:256376/1‑140 (MQ=255) aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGAt < 1:210033/140‑1 (MQ=255) tCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCTATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTa < 2:596819/140‑1 (MQ=255) gAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGGGCGCAGAGAGGGGGTATCt > 2:103116/1‑140 (MQ=255) aaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGAGGGGTTATCttt > 2:66580/1‑139 (MQ=255) gcccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTg < 1:163371/139‑1 (MQ=255) ccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGc < 2:886086/140‑1 (MQ=255) cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCAGAGAGGGGTTATTTTTCCTGACTGAACTGGCGGTTAAGtt > 2:497849/1‑140 (MQ=255) tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTTCCTGACTGAACTGGCGATTTAGttt > 2:163371/1‑140 (MQ=255) tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCAGAGAGGGGTTATCTTTCCTGACTGAACTGGCGGTTAAGTTTTcgc > 2:927585/1‑140 (MQ=255) taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAccc < 1:930775/140‑1 (MQ=255) aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCg < 1:974144/140‑1 (MQ=255) gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGGAGTGAGCGCAGAGAGGGGTTATTTTTCCTTACTGAACTGGCGATTTAGTTTCCGCGTACCCCGCGt > 2:598050/1‑140 (MQ=255) gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGt > 1:821571/1‑140 (MQ=255) tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGGAGTGAGCGCAGAGAGGGATTATTTTTCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCGGCggt > 2:314633/1‑138 (MQ=255) tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGGTTATCTTTCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCGGCggt > 2:509351/1‑138 (MQ=255) tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCCATTTTTTTTTTTTTTGCCCTGCGAAGTGGGCGCAGAGAGGGGGTATCTTGGCTGACTGAACTGGCGGTTAAGGTTTCGCGTCCCCGGCGtt > 2:755261/1‑140 (MQ=255) cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGGGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGGTTAAGTTTTCGCGTACCCGGCGttt > 2:73994/1‑140 (MQ=255) cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTTAGAAGTGAGCGCAGAGAGGGGGTATCTTTCCTGACTGAACTGGCGATTAAATTTTCGCGGACCCCGCGgtt > 2:589729/1‑140 (MQ=255) cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGttt < 1:497849/140‑1 (MQ=255) cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGga < 1:788413/140‑1 (MQ=255) tGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgatg < 1:275075/140‑1 (MQ=255) ccATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCCTGAGAAGTGAGCGCCGAGAGGGGTTATCTTTCCTGACTGGACTGGCGATTAAGTTTTCCCGTACCCGGCGTTTTGTTATTTTGggtgaggat > 2:788413/1‑140 (MQ=255) aTGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGatat < 2:601019/140‑1 (MQ=255) gcgGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAg < 1:103116/140‑1 (MQ=255) gATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGGGGGGTTTTTTTTCCTGACTGAACTGGCGGTTTAGTTTTCGCGTACCCGGCGGTTTGTTTTTTCGGATGATGGTATTTAAGCTc > 2:264327/1‑140 (MQ=255) tctcGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCgg < 1:598050/140‑1 (MQ=255) tcGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTTACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCt < 1:66580/140‑1 (MQ=255) cggcggTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGGGATTAAGTTTCCGCGGACCCGGCGGTTTGTTATTGCGGGTGATGATATTGAAGCTCGGCTGATTGGt > 2:556102/1‑140 (MQ=255) gTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGTTTTCGCGGACCCGGCGGTTTGTTATTGCGGGTGATGATATTGAAGCTCGGCTGATTGGTTCAtt > 2:519519/1‑140 (MQ=255) | ATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATT > NC_000913/4604207‑4604473 |
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: ‑ |