| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,604,350 | 0 | C | T | 35.5% | 38.3 / 11.2 | 32 | S64F (TCT→TTT) | bglJ | DNA‑binding transcriptional regulator BglJ |
| Reads supporting (aligned to +/- strand): ref base C (6/14); new base T (11/0); total (18/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.50e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTC‑ATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGG > NC_000913/4604119‑4604594 | agagAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTC‑ATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAgcgc < 1:335462/251‑1 (MQ=255) nagaACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTC‑ATTTTTTTTTTTTTTGCCA‑Tgggaggttgggagg > 2:244979/2‑238 (MQ=255) gagaACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTC‑ATTTTTTCTTTTTCTGCCA‑Tgag < 1:244979/240‑1 (MQ=255) acttGCCTTCTGTTTATGGAGATGCCGCAGAATGGCACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTC‑ATTTTTTTTTTTTTTTCCA‑TGTGAAGGGAGCGCCGCGAGGGAt > 1:114988/1‑251 (MQ=255) acttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTC‑ATTTTTTTTTTTTCTGCCC‑TGAGGAGGGAGCGCAGAGAggggt > 1:500/1‑249 (MQ=255) tATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTC‑ATTTTTTTTTTTTTTTCCA‑TGAGAAGGGGGCGCCGAGAGGGGTTATCTTGCctgact > 1:15071/1‑251 (MQ=255) agaTGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTC‑ATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACt < 2:124563/251‑1 (MQ=255) gaTGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATCGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTC‑TTTTTTTTTTTTTTTTTCC‑TGCGAAGTGGGCGCAGAGAGGGGTTCTCTTTCctttctgaaacg > 2:465847/1‑241 (MQ=255) ccGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTC‑ATTTTTTTTTTTTCTGCCC‑TGAGAAGTGAGCGCAGAGAGGGGTTATCTTTCCTGACTGAACTGGCgg > 1:249768/1‑250 (MQ=255) ggAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTC‑ATTTTTTTTTTTTTTTCCC‑TGCGAAGTGAGCGCGGAGAGGGGTTATTTTGCCTGCCTGAACTgggggtttaatttt > 2:541415/1‑239 (MQ=255) cGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTC‑ATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAG‑TTTCCGCGTACCCg < 1:33813/251‑1 (MQ=255) aagaagAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTC‑ATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAG‑TTTCCGCGTACCCGGCGttt < 1:35653/251‑1 (MQ=255) gaagaGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTC‑ATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAG‑TTTCCGCGTACCCGGCGTTTAg < 2:187331/251‑1 (MQ=255) agagAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTC‑ATTTTTTCTTTTTTTGCCA‑TGAGAAGTGAGCGCAGAGAGGGGTTTTTTTTCCTTAATGAACTGGGGATTAAG‑TTTTCGCGTACCCGGCggtttgtt > 1:146734/1‑244 (MQ=255) gagaAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTC‑ATTTTTTTTTTTTCTGCCA‑TGAGGAGGGGGCGGCGGGAGGGGTTATCTTTCCTGGCTGGACCGGGGGTTAAG‑TTTTCGCGGACCCGGCggttcgttt > 1:213368/1‑243 (MQ=255) aGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTC‑ATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAG‑TTTCCGCGTACCCGGCGTTTa < 1:494289/229‑1 (MQ=255) aGAAAAATGCGCCATGAGTAGTACCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCCTTTCAGGCTGCCATGTCGCGGATCTCGTTTGCGGCGGTC‑ATTTTTTTTTTTTTTTTCC‑TGAGAAGTGAGTGCAGAGGGGGATTTTCTTTCCTGGCTGGACCTGCGATTTCG‑TTTTCGCGCACCCGGCGGTTaatagcggcagcgggtcggcg > 2:494289/1‑229 (MQ=255) aaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTC‑ATTTTTTTTTTTTTTGCCC‑TGAGAAATGCGCGCAGAGAGGGATTATTTTGCCTTACTGGAATGGCGATTAAG‑TTTTCGCCTACCCGGCGGTTAGTTATTTCgggttgttttatttaagcc > 1:503428/1‑234 (MQ=255) gAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTC‑ATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTg < 2:114988/251‑1 (MQ=255) agtaTCGGTATTGAGAGATTATTAAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGC‑ATTTTTTTTTTTTTTGCCA‑TGAGAAGTGAGCGCCGGGGGGGGTTATCTTGCCTGATTGAAATGGCGGTTAAA‑TTTTCGCGTAACCGGCGTTTTGTTATTGCGGGTGGAGATATTTAAGATCGGATGATTGGttt > 2:212505/1‑250 (MQ=255) cGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTC‑ATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTg < 1:541415/250‑1 (MQ=255) aaaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTC‑ATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGgtg < 1:202677/251‑1 (MQ=255) gTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTC‑ATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTAt < 1:147861/251‑1 (MQ=255) ttGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTC‑ATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTAtt < 1:210063/250‑1 (MQ=255) tGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTGTTTTTTATTTTT‑TGCCC‑TGAGAAGGGGGAGCAGAGAGGGATTATTTTTGCTGACTGAACTGGCGATTAAG‑TTTTCGCGTACCCGGCGTTTTTTTATTGCGGATGATGTGATTTAAACTCGGCTGATTGGTTCATTGTCTCCATAACCGcctggaggggggttt > 1:539604/1‑236 (MQ=255) tataAGCTCCATACCTATTCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTC‑ATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAG‑TTTCCGCGTGCCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTTGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGt < 2:289738/240‑1 (MQ=255) tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTC‑ATTTTTTTTTTTTCTGCCC‑TGAGAAAGGGGCGCACAGAGGGATTATCTTTCCTGGCTGAACTGGCGGATAAG‑TTTTCGCGAAACCGGCGGTTAGTTATTGCGGATAAGGGTATTTAAGCTCGGCTGATTGGTTCATTGTCGCCCTCACCGCTGGACGGGGGATTAAgaaggacggaagg > 1:289738/1‑239 (MQ=255) ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTC‑ATTTTTTTTTTTTTTGCCA‑TGAGAAGTGAGGGCAGAGAGGGATTATCTTGGCTGACTGAACTGGCGATTAAG‑TTTCCGCGGACCCGGCGTTTAGTTATTGCGGGTGATGATATTGAAGCTCGGCTGGTTGGTTCATTGTCGCCATCACCGCTGGGCGGGGTATTTAGTAAAGCGTCAACGCTGGAGAatttt > 1:477860/1‑251 (MQ=255) ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTC‑CTTTTTTTTTTTTCTTCCACTG‑GGAGTGGGGGCCGCGAGGGGATATCTTGTCTTACTGAAATTGCGCTTAAA‑TTTCCGCGTCCCCCGCGTTTTGTTTTTTCTGAGGATGATAATGAAGCTCGGCTGATTTGTTCATTTTTTCCCTCACCGCTGGGCGGTGGTTTTAGTTAAACGTGCACGCTGGCGGTTTTTCATCAGGAActct > 2:406028/1‑251 (MQ=255) aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTT‑TTTTTTTTTTTTTTTGACA‑TGAGAAGGGAGCGCAGAGGGGGAATATTTTGCCTGACCGAACTGGCGAATAAGGTTTCTG‑GTACCCGGCGGTTTGTTTTTGCGGATGATGAAATTTAAACTCGGCTGGTTTGTTCGTTTTTGGCATCACCTCCGGGCAGGGTATTAAGTAAAACGGCCACGCTGGCGGtttttttttagggaatcttttttttcatt > 1:472848/1‑227 (MQ=255) gcggTC‑ATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGCCCGACCGACt < 2:498488/251‑1 (MQ=255) gTC‑ATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGa > 2:227431/1‑250 (MQ=255) tttttctttttcTGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAAt < 1:213301/250‑1 (MQ=255) tttctttttctGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATc > 2:414708/1‑249 (MQ=255) ctttttctGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTgg > 1:361511/1‑251 (MQ=255) | AGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTC‑ATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGG > NC_000913/4604119‑4604594 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |