Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,604,344 | 0 | C | T | 40.9% | 48.1 / 24.0 | 44 | S62F (TCT→TTT) | bglJ | DNA‑binding transcriptional regulator BglJ |
Reads supporting (aligned to +/- strand): ref base C (3/23); new base T (18/0); total (21/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.29e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAT > NC_000913/4604207‑4604481 | aTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttttt > 1:81591/1‑140 (MQ=255) aTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttttt > 1:522520/1‑140 (MQ=255) cGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttcttttt < 1:529531/140‑1 (MQ=255) tCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct < 1:860701/140‑1 (MQ=255) cATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAttttttttttttttt > 2:971917/1‑137 (MQ=255) aTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAg < 1:291715/140‑1 (MQ=255) tCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAgg > 1:119693/1‑139 (MQ=255) tCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGt < 2:647673/140‑1 (MQ=255) gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCATGAGAAGTGAg > 2:265291/1‑140 (MQ=255) gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAg < 1:265291/140‑1 (MQ=255) aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAgcg > 2:575001/1‑140 (MQ=255) aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcg < 1:762674/140‑1 (MQ=255) tGAGAGTTTATTCAGAAAGTTTGCGGGTTACGCTTATAAGCGCCATTCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTTAGCGCa < 1:79218/140‑1 (MQ=255) gagaGTTTATTCAGAAAGTTTGC‑‑GTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagag > 1:609015/1‑140 (MQ=255) gTTTATTCAGAAAGTTTGC‑‑GTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAggg < 2:609015/140‑1 (MQ=255) tttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTCTTCCCTTAGAACTGAGCGGCGAGAgg > 2:706580/1‑140 (MQ=255) ttATTCAGAAAGTTTGCGGGTAACCTTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAggg < 1:280599/140‑1 (MQ=255) aGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATc < 2:775944/140‑1 (MQ=255) tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCAGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑Gttt > 2:923918/1‑140 (MQ=255) ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGt < 2:383879/140‑1 (MQ=255) aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCg < 2:81591/140‑1 (MQ=255) aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCg < 1:410319/140‑1 (MQ=255) ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGGGCGGAGGGAGGGGTTATCTTGCCTGACTGAACTGGCGGTTAA‑GTTTTCCCGTAccccg > 2:665031/1‑138 (MQ=255) gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCAGAGAGGGGTTATCTTTCCTGACTGAACTGGCGATTAA‑GTTTTCGCGTACCCGGCgg > 2:457894/1‑139 (MQ=255) gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCAGAGGGGGGTTTTTTTGCCTGACTGAACTGGCGGTTAA‑GTTTCCGCGTACCCcgcgg > 2:962116/1‑139 (MQ=255) gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGt < 1:125197/140‑1 (MQ=255) tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGAGAGGGCAGAGAGGGGTTATTTTGCCTGACTCAAATGGGGAATGA‑GTTTTCGCGGACCCGGCggt > 2:233800/1‑138 (MQ=255) tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTTCATGTGAAGTGAGGGGGGGGGGGGGGTATTTTTTCTGACTGAAATTGGGGTTTA‑GATTTCGCGGATCCGGGGtt > 2:125197/1‑140 (MQ=255) tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCATTAGAAGTGAGCGCAGGGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGGACCCGGCGtt > 2:607044/1‑140 (MQ=255) cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCAGAGAGGGGTTATCTTTCCTGACTGAACTGGCGATTAA‑GTTTTCGCGGACCCGGCGttt > 2:604065/1‑140 (MQ=255) cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATTAGAAGTGAGCGCAGAGAGGGGGTATCTTGCCTGACTGAACTGGCGGTTAA‑GTTTTCGCGTACCCGGCGgtt > 2:157493/1‑140 (MQ=255) cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCTGCGttt < 1:477414/140‑1 (MQ=255) ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAg < 1:604065/140‑1 (MQ=255) tCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTAt < 1:665031/140‑1 (MQ=255) aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTg < 1:828456/140‑1 (MQ=255) aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTg < 1:329037/140‑1 (MQ=255) aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTg < 1:535809/140‑1 (MQ=255) aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTg < 1:469385/140‑1 (MQ=255) tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTTCCTGACTGAACTGGCGATTAA‑GTTTTCGCGTACCCGGCGGTTTGTTAtttt > 2:88202/1‑138 (MQ=255) tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCATGAGAAGGGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAA‑ATTTTCGCGTACCCGGCGGTTAGTTTTTgt > 2:329037/1‑139 (MQ=255) aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCGTTGGAAGTGGGGGCCGAGGGGGATTGTCTTTCCTGCCTGAACGGGCGCTTAAGGTTTCCG‑GCACCCCGCGGTTAGGTATTGCgggg > 2:462187/1‑138 (MQ=255) aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGGTTATCTTTCCTGACTGAACTGGCGGTTAA‑GTTTTCGCGTACCCCGCGTTTAGTTATTTCGgatggt > 2:535809/1‑138 (MQ=255) gCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgatga < 2:364336/140‑1 (MQ=255) cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGGGCAGAGAGGGGTTATCTTTTCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCCGCGTTTTGTTATTTCGGGTGATGata > 2:700569/1‑140 (MQ=255) gTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGa < 2:885929/140‑1 (MQ=255) tttGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTGA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGAt < 1:233800/140‑1 (MQ=255) tGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGGGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTTCGCGTACCCGGCGTTTTGTTATTTCGGGTGATGATATTGAAGCTCGGCTGGTTg > 2:451680/1‑140 (MQ=255) cggcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTTGTTATTGCGGGTGATGATATTTAAGCTCGGCTGATTGGt > 1:573782/1‑140 (MQ=255) gTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGGGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTTCGGGTGATGATATTTAAGCTCGGCTGATTGGTTCAtt > 2:190792/1‑140 (MQ=255) tCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTg < 1:190792/140‑1 (MQ=255) ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt < 1:457894/140‑1 (MQ=255) | ATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAT > NC_000913/4604207‑4604481 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |