| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,604,350 | 0 | C | T | 46.2% | 18.9 / 14.0 | 26 | S64F (TCT→TTT) | bglJ | DNA‑binding transcriptional regulator BglJ |
| Reads supporting (aligned to +/- strand): ref base C (4/10); new base T (12/0); total (16/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.01e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGG‑TTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCG‑ACCGACCGACTGAACAATCAGTGGTAC > NC_000913/4604109‑4604597 | gcAGGGTGCAAGAGCACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGACAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCC‑Ttttttttttttcatccatc > 2:450566/1‑242 (MQ=255) gagaACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGGATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTAATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAATTATTTTTTTTTTTGCCAT‑AGAAGTgtgtgcc > 2:397249/1‑245 (MQ=255) agaACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCag < 1:253077/251‑1 (MQ=255) tGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGGATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATTAGAAGTGAGCGCAGAGAGGGAtt < 2:165900/249‑1 (MQ=255) aTGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGCGTTTATTCAGCAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTAA‑TTTTTTCTTTTTCTTCCCTGTGAAGTGAGTGCagggggggcttttgtttcctcgctg > 2:181689/1‑228 (MQ=255) tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTGCCATGAGGAGGGAGGGCCGGGGGGGAGTATTTTGCCTGAATGAACTggg > 1:305089/1‑250 (MQ=255) cAGAATGGAACACAGCCCAATTGAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTa < 1:2670/251‑1 (MQ=255) acaGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTA‑TTTTTTTTTTTTTTTCCATGTGGAGTGGGCGGCGAGAGGGGTTATCTTTTCTGTCTGAACTGGGGATTTAGTTTTcgcgg > 2:13868/1‑250 (MQ=255) tAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTGCCATGAGAAGTGGGGGCAGAGAGGGGTTATCTTTCCTGACTGAACTGGCGATTAAGTTTCCCCGTACCCGGCGtt > 1:212274/1‑251 (MQ=255) tCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCC‑TTTTTTTTTTTTTTGCCCTGAGAAGGGAGGGCAGGGAGGGGTTATCTTTCCTGGCTGAACTTGCGGTTAAAATTCCGCGGCCCCGGCGttttggttttttgggg > 2:247745/1‑238 (MQ=255) aTAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGatat < 2:257934/251‑1 (MQ=255) cATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTTTGCCATGAGAAGTGAGGGCCGGGAGGGATTATCTTGCCTGACTGAACTGGCGGTTAAAGTTCCGCGTACCCGGCGGTTTGGTATTGCGGGTGATGATATTTAAGCTCGGCTg > 1:301780/1‑250 (MQ=255) agtagtATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCTGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTgg < 1:64849/251‑1 (MQ=255) gtagtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCAGGGAGGGGTTATCTTTCCTGTATGTACTTGGGCTTTAGTTTCCCCGGACCCGGCGTTTTGGTTTTGCGGGTGATGTTATTGAAGCTCGGCTGGTTggg > 2:44178/1‑250 (MQ=255) aTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTGCCCTGAGGAGTGAGCGGCGAGAGGGGTTATCTTTCCTGACTGAACTGGCGGTTAAATTTTCGCGTACCCCGCGGTTTGTTTTTGTGGGTGGTGATATTGAAGCTCGGCTGGTTGG‑GTCAtt > 2:351333/1‑251 (MQ=255) tGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGG‑TTCATTGTCGCCa < 1:351333/250‑1 (MQ=255) ttCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTTTGCCATGAGAAGCGGGCGCAGAGAGGGGTTATTTTGCCTGACTGAACTGGGGATTAAGGTTCCGCGTACCCCGCGGTTTGGTATTTCGGGTGATGATATTGAAGCTCGGCTGGTTGGGTTC‑TTTTCGCCATC‑ACCGCTggg > 1:117180/1‑250 (MQ=255) ccATGCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGG‑TTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGgag < 1:141108/250‑1 (MQ=255) ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTGCCACGAGAAGTGAGGGCCGGGAGGGGTTTTCTTTCCTGACTGAAATGGGGGTTAAGGTTTCGCGGATCCGGCGGTTAGGTATTTCGGGTGGTAATATTTAAGCCCGGGTGATTTG‑TTCATTGTTGGCATC‑ACCCCTGGGCGGGGGGTTTATTAAAACGTCAACGCGGGGGAttttttat > 2:400729/1‑248 (MQ=255) ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGG‑TTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTtcat < 1:13868/251‑1 (MQ=255) ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCC‑CCCTTTCCCTTTCTCCCCTGCGCAGGGAGCGCCGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTTAGGTTCCGCGTACCCGGCGTTTTGTTTTTGCGGGTGATGGTATTTAAGCTCGGATGATTGG‑TTCCTTGTTGCCATC‑ACCGCTGGACGGGGTATTTAGTAAAGc > 1:78696/1‑219 (MQ=255) ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTGGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTTCCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGG‑TTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGc < 2:78696/219‑1 (MQ=255) ttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGGTTATCTTGCCTGACTGCACTGGCGGTTAAATTTTCGCGTACCCCGCGGTTTGTTTTTGCGGGTGGTGGTATTTAAGCTCGGCTGGTTGG‑TTTATTGTCGCCATCCACC‑CCGGGCGGGGTATTAAGTAAAACGTCCACACTGGAGATTTTTTATCAGGAACCCtttttttc > 2:444350/1‑247 (MQ=255) cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGG‑TTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCAtt < 1:42562/251‑1 (MQ=255) ggATGCCATGTCGCGGATCTCGTTTGCGGCGGGAG‑TTTTTTTTTTTTTTGCCCTGGGAAGGGAGCGCAGAGAGGGGTTGTCTTGGCTGACTGTAAAGGGGATTTAGTTTTCCCGTGCACGGCGGTTTGTTTTTACGGGTGAAGATATTTAAACTCGCCTGATTTG‑TTTATTTTCGCCCTC‑ACCCCTGGGCGGTGTGTTTATGGAAGAGGCAACGCTGGGGATTTTTTTTCCGGGACCCTTTTTGTCAttt > 2:430804/1‑249 (MQ=255) cA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGG‑TTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGAGTACGTCAGGCG‑ACCGACCGACTGAAc > 2:244857/1‑250 (MQ=255) a‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGCCGTTTGGCTATTGCGGTTGATGATATTGAAGCTCGGCTGATTGG‑TTCATTGTCGCCATC‑ACCGGTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTAAGTCAGGCGTACC‑ACCGACTGaata > 2:375518/1‑248 (MQ=255) ttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGG‑TTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCG‑ACCGACCGACTGAACAATCAGTgg < 1:412331/249‑1 (MQ=255) ttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGG‑TTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCG‑ACCGACCGACTGAACAATCAGTGGTAc < 1:181689/250‑1 (MQ=255) | GCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGG‑TTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCG‑ACCGACCGACTGAACAATCAGTGGTAC > NC_000913/4604109‑4604597 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |