| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,536,197 | 0 | A | T | 40.0% | 37.9 / 22.1 | 35 | intergenic (‑2144/‑400) | topAI/ythB | toxin of the TopAI‑YjhQ toxin‑antitoxin system, TopA inhibitor/protein YthB |
| Reads supporting (aligned to +/- strand): ref base A (0/21); new base T (14/0); total (14/21) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.31e-10 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GTCTTTGGTGTACCGTGGTGAAGGCGTAAAAAATGAAATAAAGGCAATATATGGTAGTCTGGATCAACGTCAGCTTTAAATGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGC‑TTTTAGGAATGAAAGAAAATGCAAGCA‑ATAATATAATAA‑TGGCCTAA‑TTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCA‑GGACGGT‑GCTGCGAACGCTGTGG‑CCCAAAACGTTCCAGTAAACTGTG‑TAGATGTTCGTAGATTTCCGAGCCATAATCGAAAAGGCGCGGATGTCTATTAAGGCATCC > NC_000913/4535953‑4536427 | gTCTTTGGTGTACCGTGGTGAAGGCGTAAAAAATGAAATAAAGGCAATATATGGTAGTCTGGATCAACGTCAGCTTTAAATGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGC‑TTTTAGGAATg < 2:7109/251‑1 (MQ=255) ttGGTGTACCGTGGTGAAGGCGTAAAAAATGAAATAAAGGCAATATATGGTAGTCTGGATCAACGTCAGCTTTAAATGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGC‑TTTTAGGAATgaa < 1:267258/249‑1 (MQ=255) ggtccGTGGTGAAGGCGTAAAAAATGAAAGAAAGGCAATATATAGTAGGCTGGATCAACGTCAGCTTTAAATGTGATGAGCCGTTCATGATACTGGAGTAGTAGAGGGCATATGTTACCAATGTTAAGGCGACTGTCCGTTTTATTGCGAGTCGCCTGTCCTGATTATAAATCTTATTGGGTACCTAGTCATTTTGCTTTGCATAATGCGCCGTGGATTAATGTCGGGTGC‑TTTTAGGAATGAAAGAAAAt < 2:379393/248‑1 (MQ=255) gAAGGCGTAAAAAATGAAATAAAGGCAATATATGGTAGTCTGGATCAACGTCAGCTTTAAATGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGC‑TTTTAGGAATGAAAGAAAATGCAAGCa‑at < 2:478924/250‑1 (MQ=255) aGGCGTAAAAAATGAAATAAAGGCAATATATGGTAGTCTGGATCAACGTCAGCTTTAAATGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGC‑TTTTTTGGATGAAAGAAAATACAAGCa‑ataa > 1:564723/1‑250 (MQ=255) aaaaaaTGAAATAAAGGCAATATATGGTAGTCTGGATCAACGTCAGCTTTAAATGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGG‑GTTTTTGAATGAAAGAAAAAGAAAGGa‑ataaaaagaa > 2:275749/1‑244 (MQ=255) aaTAAAGGCAATATATGGTAGTCTGGATCAACGTCAGCTTTAAATGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTAC‑TAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGC‑GTTTTTGTATGAAAAAAAATacaaaca‑aaaaatattaaaatgaccaat > 2:36300/1‑231 (MQ=255) tatatGGTAGTCTGGATCAACGTCAGCTTTAAATGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGC‑TTTTTGTAATGAAAAAAAAAAAAAACA‑ATAAaataataa‑tagtctcattttttttta > 1:34566/1‑233 (MQ=255) atGGTAGTCTGGATCAACGTCAGCTTTAAATGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGC‑TTTTAGGAATGAAAGAAAATGCAAGCA‑ATAATATAATAA‑TGGCCTAA‑TTTCTTTGAAcat < 1:275749/250‑1 (MQ=255) tAGTCTGGATCAACGTCAGCTTTAAATGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGG‑TTTTTGGAATATAAGAAAATAAAAGCA‑ATAATATAATAA‑TTGGCTAA‑ATTTTTTTAACATCAttt > 1:87928/1‑250 (MQ=255) aGTCTGGATCAACGTCAGCTTTAAATGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGC‑TTTTAGGAATGAAAGAAAATGCAAGCA‑ATAATATAATAA‑TGGCCTAA‑TTTCTTTGAACATCATTa < 2:429055/250‑1 (MQ=255) cAACGTCAGCTTTAAATGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGG‑TTTTTTGAATGAAAAAAAAAGAAAACA‑ATAAAATAAAAAATGGC‑TAA‑TTTATTTTAACATCAtttaacaaattat > 1:267682/1‑240 (MQ=255) cAACGTCAGCTTTAAATGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGC‑TTTTTGGAATGAAAGAAAAAGCAAAGA‑AAAAAATAATAA‑AAGCCTAA‑Tttttttttacaacatttataatatcat > 1:218910/1‑230 (MQ=255) cAACGTCAGCGTTAAATGTGATGAGCCGTTCATGATACAGGAGTATTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTGATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTTCCTAATCATTTTCCTTTGCATAATGCGCCTTGGATTAATGGCGGGTGC‑TTTTAGGAATGAAAGAAAATGCAAGCA‑ATAATATAATAA‑TTGCCTAA‑TTTCTTTGAACATCATTAAGCATatca < 1:236171/250‑1 (MQ=255) aGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGC‑TTTTAGGAATGAAAGAAAATGCAAGCA‑ATAATATAATAA‑TGGCCTAA‑TTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCg < 2:153798/251‑1 (MQ=255) aTGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGC‑TTTTAGGAATGAAAGAAAATGCAAGCA‑ATAATATAATAA‑TGGCCTAA‑TTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTta < 1:244332/251‑1 (MQ=255) agtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGC‑TTTTTTGAATTAAAAAAAAAGCAAACA‑AAAATATAAAAA‑TAGCCTAA‑TTTTTTTTTAAATCATTAAACATATCATCATCGTCTTTTTTTCCGAATCCCTCCACTTATAATCCCCAGAAc > 1:144244/1‑251 (MQ=255) aaaGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGC‑TTTTAGGAATGAAAGAAAATGCAAGCA‑ATAATATAATAA‑TGGCCTAA‑TTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGc < 1:458944/251‑1 (MQ=255) atatTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGC‑TTTTAGGAATGAAAGAAAATGCAAGCA‑ATAATATAATAA‑TGGCCTAA‑TTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGcc < 2:213399/212‑1 (MQ=255) ccAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGC‑TTTTTGGGATGAAAGAAAAAGAAAGCA‑ATAATATAATAA‑TGGCCTAA‑TTTCTTTTTACATCATTCAACATATCATCATCTTCTTTTTTCTCGAATCCGCCGACTTTTTATCACCCAAACGGCCTATCA‑GGACggggc > 1:58668/1‑248 (MQ=255) ttATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGC‑TTTTAGGAATGAAAGAAAATGCAAGCA‑ATAATATAATAA‑TGGCCTAA‑TTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATAGCA‑GGACGGT‑GCTGCGAACGCTGTGG‑CCCAAAACg < 1:526880/250‑1 (MQ=255) tGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGC‑TTTTAGGAATGAAAGAAAATGCAAGCA‑ATAATATAATAA‑TGGCCTAA‑TTTCTTTGAACATCATTAAGCATATCAGCATCGTCTATTTGCCAGAAGCCGTCGACTTATAATCGCCAGAACGGCAGATCA‑GGACGGT‑GCTGCGAACGCTGTGG‑CCCAAAACGTTcc < 1:28122/250‑1 (MQ=255) ttCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGA‑TTTTTGGAAAGAAAGAAAAAACAAACACATAA‑ATAAAAA‑TTACCTCA‑TTTTTTTTTAAATCATTAAGCATATCATTAAATTCTATTTGTCCGAAACCGTCGACTTATTAACGCCCGAAAGGCAACTCC‑GGACGGC‑GGTGGGAAAGCTGTGG‑GCCCAAACgatatcgtaaatagtgg > 1:566760/1‑234 (MQ=255) aTTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGC‑TTTTAGGAATGAAAGAAAATTCAAGCA‑ATAATATAATAA‑TGGCCTAA‑TTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCATAATCCGTCGACTTATAATCGCCAGAACGGCATATCA‑GGACGGT‑GCTGCGAACGGTGTGG‑CCCAAAACGTTCCTGTAGACTGTG‑TAGATGtt < 2:373274/251‑1 (MQ=255) tataAATCTTATTGGTTACCTAATCATTTTCCGTTGCATAATGCGCCTGGGATTAATGGCGGGTGC‑TTTTAGGAATGAAAGAAAATGCAAGCA‑ATAATATAATAA‑TGGCCTAA‑TTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCA‑GGACGGT‑GATGCGAACGCTGTGG‑CCCAAAACGTTCCAGTAAACtctttgttttttc < 1:433125/250‑13 (MQ=255) tATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGC‑TTTTAGGAATGAAAGAAAATGCAAGCA‑ATAATATAATAA‑TGGCCTAA‑TTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCA‑GGACGGT‑GCTGCGAACGCTGTGG‑CCCAAAACGTTCCAGTAAACTGTG‑TAGATGTTCGTAGAttt < 1:214679/249‑1 (MQ=255) ttACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGC‑TTTTTTGAATGAAAAAAAAAGCAAGAA‑ATAAAAAAATAA‑AAGCCTTA‑TTTTTTTTAAAATAATTAAACATATCATTATAGTCTATTTTCCCGAATCCGTCGACTTTTTATCTCCAGAAACGAATATCA‑GGAAGGG‑GCTGCGGAAGCTGTGGGCCCAAAA‑GTTTCCGTAAAATATG‑TAGATGtttgttggtttttccgccc > 1:284352/1‑234 (MQ=255) ccTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGC‑TTTTTTGAATGAAAAAAAAAGCAAGCA‑AAAATATAATAA‑TGGCCTAAATTTTTTTTAA‑ATAATTTAACATATTATCCTCGTCTATTTTTCAGAATCCCTCGCCTTTTAATCCCCCGAACGGCATATAA‑GGACGGG‑GGTGCGGACGCTGTGGGCCCAAAA‑GTTCCAGTAAAATGTG‑TAGATGttttttgattttccagccata > 1:470443/1‑231 (MQ=255) cTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCGTTTTTGGAA‑GTAAAAAAAAAAAAGCA‑AAAATATAATAA‑TTGACTCA‑TTTTTTTGAAAATTATTTAACATTTCATCATCGTCTTTTTTTTCGAATCCGTCGTCTTTTTATTGCCCGAACGGCACAACA‑GGGAGGT‑GGTGCGGAAGCAGTGGTCCCAAAAC‑TTCCCGTAAAATGTG‑TAGAGGTTGGGTGATTTCTGAGCCATAAt > 1:179272/1‑251 (MQ=255) cTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGC‑TTTTTGGGATGAAAAAAAAAGCAAAAA‑ATTATATAATAAATGCCCCAT‑TTTTTTTTAA‑ATCACTAAACATATCATCATCGTCTATTTTTCCGAATCCCTCGACTTATAATCGCCCGAAAGGCATATCAAGGACGG‑‑GCTGCGGAAGCTGTGGGCCCAAAA‑ATTCCACTAAAATTTG‑TTGAAGTTTGTAGATTTctgtgccctaat > 1:359638/1‑240 (MQ=255) aaTCATTTTCCTTTGCATAATGCGCCTGGTATTAATGGCGGGTGC‑TTTTAGAAATGAAAGAAAATGCAAGAA‑ATAATATAATAA‑TGGCCTAA‑TTTCTTTGAACATCATTAAGCATAGCATCATGGTCTATTTGCCAGAATCCTTCGACTTATAATCGCAAGAAGGGCATATCA‑GG‑CGGTTGCTGCGACCGCTGTGG‑CCCAAAACTTTCCAGTAAACTGTG‑TAGATGTTCGTAGATTTCCGAGCCATAATCg < 1:55548/251‑1 (MQ=255) aaTCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGC‑TTTTAGGAATGAAAGAAAATGCAAGCA‑ATAATATAATAA‑TGGCCTAA‑TTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCA‑GGACGGT‑GCTGCGAACGCTGTGG‑CCCAAAACGTTCCAGTAATCTGTG‑TAGATGTTCGTAGATTTCCGAGCCATAAt < 1:164005/249‑1 (MQ=255) cattggcaaaaggcgcCGGGGATTAATGGGGGGTGC‑TTTTAGGAATGAAAAAAAATGCAAGCA‑AAAATATAATAA‑TGGCCTAA‑TTTCTTTGAACATCATTAAGCATATCAGCATCGTCTATTTGCAAGAATCCGTCAACTTATAACCGCCAGAAGGGCATACCA‑GGAGGGG‑GCTGCGAACGCTGTGC‑CCCAAAACGTCCAAGTAAACGGTG‑TAGGTGTTCGTAGTTTTCCGAGCCATAATCGAAAAGgcgc < 2:227436/239‑1 (MQ=255) gcgcCTGGGATTAATGGCGGTTGC‑TTTTAGGAATGAAAAAAAATGCAAGCA‑AAAATATAATAA‑TGGCCTAA‑TTTCTTTGAACATCATTAAGAATATCATCATCGTCTATTTGCCAGAATCCTTCGACTTATAACCGCCAGAGCGGCATATCA‑GGACGGT‑GCTCCGATCGCGGTGG‑CCAAAAGCGTTCCAGGAAA‑GGTGGTAGATGTTAGTAGATTTCCGAGCCATAATCGAAAAGGCGCGGATGTCTATTa < 2:267682/251‑1 (MQ=255) gATTAAGGGCGGGTGC‑TTTTAGGAATGAAAGAAAATGCAAGCA‑ATAATATAATAA‑TGGCCTAA‑TTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCA‑GGACGGT‑GCTGCGAACGCTGTGG‑CCCAAAACGTTCCAGTAAACTGTG‑TAGATGTTCGTAGATTCCCGAGCCATAATCGAAAAGGCGCGGATGTCTATTAAGGCATcc < 1:38823/251‑1 (MQ=255) | GTCTTTGGTGTACCGTGGTGAAGGCGTAAAAAATGAAATAAAGGCAATATATGGTAGTCTGGATCAACGTCAGCTTTAAATGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGC‑TTTTAGGAATGAAAGAAAATGCAAGCA‑ATAATATAATAA‑TGGCCTAA‑TTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCA‑GGACGGT‑GCTGCGAACGCTGTGG‑CCCAAAACGTTCCAGTAAACTGTG‑TAGATGTTCGTAGATTTCCGAGCCATAATCGAAAAGGCGCGGATGTCTATTAAGGCATCC > NC_000913/4535953‑4536427 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |