| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,862,668 | 0 | A | G | 46.9% | 20.7 / 13.3 | 32 | I312I (ATT→ATC) | glvC | putative PTS enzyme II component GlvC |
| Reads supporting (aligned to +/- strand): ref base A (2/15); new base G (15/0); total (17/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.42e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCAT‑GTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAA‑ACAGCAACGGTGAAA‑TGAA‑CAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAG‑CATGGCGGTT‑AAGG‑TTGCGGGAATCAGC‑AAGCCC‑GCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGA‑TTTTTGAGTTA‑CCGTGCAGGGCAAAA‑CCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAA > NC_000913/3862428‑3862913 | aaTGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCAT‑GTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAA‑ACAGCAACGGTGAAA‑TGAA‑CAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAG‑CATGGCGGTT‑AAGG‑GTGCGGGGATCAGC‑aaac > 1:691646/1‑249 (MQ=255) ggCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCAT‑GTTGCCCACCACACCAAAGATATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAA‑ACAGCAACGGTGAAA‑TGAA‑CAGGAAGGTAAAGTCCAGCGGTTCGGTAATTCCCACCAG‑CATGGCGGTT‑AAGG‑TTGCGGGAATCAGC‑AAGccc < 1:806825/250‑1 (MQ=255) gATTTGGGTTAGCCTCAAATCCGCGGGTTTGCTGAACATCGGGATCCTGTTTTGCGGTAAAACCTAGTCAATCAGGCCTCCCCCCAT‑GGTGCCCCCCACACCAAAGAGGTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAA‑ACAGCAACGGTGAAA‑TGAA‑CAGGAAAGGAAATTCCAGCGGTTCGGTAATTCCCACCAG‑CATGGCGGTT‑AAGG‑TTGCGGGAATCAGC‑AAGCCC‑Gcc < 2:663452/250‑1 (MQ=255) tGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCAT‑GTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAA‑ACAGCAACGGTGAAA‑TGAA‑CAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAG‑CATGGCGGTT‑AAGG‑TTGCGGGAATCAGC‑AAGTCC‑GCCAC‑tt < 1:141753/250‑1 (MQ=255) cgcgtgtttgtaaaccATCGGTACCCATTTTTGCGTAAAAACCTAGACAACCAG‑CCCCCCCCCATTGTGGCCCCCCACACAAAAGAAATACATTAGGGGCGTCATTGAGGCCGCCACCGGGGCGTGTCCCGCAA‑ACAGCA‑GGGTGAAAATGAA‑CGGGAAGGTAAATTCCAGGGGTGCGGAAATTCCCCCTAG‑CATGGGGGTT‑AGGT‑TTGCGTAAATAAGA‑AAGCCC‑GCCAC‑TTTTACCCGATTTTCCGGTg < 1:513923/236‑1 (MQ=255) ttGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCAT‑GTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAA‑ACAGCAACGGTGAAA‑TGAA‑CAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAG‑CATGGCGGTT‑AAGG‑GTGCGGGGATCCGC‑AAGCCC‑CCCAC‑TTTTTCCCGCttttttgggggggcagg > 1:607653/1‑238 (MQ=255) cTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCAT‑GTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAA‑ACAGCAACGGTGAAA‑TGAA‑CAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAG‑CATGGCGGGT‑AAAG‑GTGCGGGGATCAGC‑AAGCAC‑GCCAC‑TTTTTCCCGCttttttggggctggaggg > 2:300500/1‑235 (MQ=255) tGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCAT‑GTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAA‑ACAGCAACGGTGAAA‑TGAA‑CAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAG‑CATGGCGGTT‑AAAG‑GTGCGGGGATCAGC‑CAGCCC‑CCCACCTTTTTCCCGTTTTTcgtggggggcgggatg > 1:73864/1‑236 (MQ=255) gAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCAT‑GTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAA‑ACAGCAACGGTGAAA‑TGAA‑CAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAG‑CATGGCGGTTTAAGG‑T‑GCGGGGAGCAGC‑AAGGCC‑CCCACCTTTTCCCCCATTTTctgtgggggggggg > 2:48957/1‑235 (MQ=255) ccggaaTCCAGTTTTGCGGTAAAACCTATCCAACCACCCCCCCCCCCAT‑GTTCCCCCCCCCCCCAAAGAGATCCATTCCGGCCGACATGGAGGCCGCCAGCACGGCGTGTACCGAAA‑ACAGCAAGGGTGAAA‑GGAA‑CAGGAAGGTAAATCCCGGCGTTTGGGTAATTCCCCCCAG‑CATGGCGGTT‑AAGG‑TTGCGGGAACCAGC‑AAGCCC‑GCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACAt < 1:300500/245‑1 (MQ=255) aTCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCAT‑GTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAA‑ACAGCAACGGTGAAA‑TGAA‑CAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAG‑CATGGCGGGT‑AAGG‑GTGCGGGGATTAGA‑AAACCC‑CCCCC‑TTTTTACCGATTTTCCCGTGCGGCCGGGAAGTACATAGCtt > 1:429649/1‑249 (MQ=255) aaaaCCTAGTCAATCAGACCTCCGCCCAT‑GTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAA‑ACAGCAACGGTGAAA‑TGAA‑CAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAG‑CATGGCGGGT‑AAAG‑TTTCGGGGAACAGC‑AAACCC‑CCCCC‑TTTTATCCGCTTTTCCGGGGGGGGCGGGGAGTACATCGCTAAAGaaaaacccccggc > 1:667204/1‑241 (MQ=255) caataagccccccccccAT‑TTTGCCCCCCACACAAAAGAGATACATCACGCGCGAAATAGAGCCCCCCACCTGGGCGTAAACC‑CAACACAGCAAGGGTAAAA‑AGAA‑CAGAAAGGAAAATCCCAGCGGTTGGGTAATCCCCACCAG‑CATGGGGGTT‑AAGG‑TGGGGGGAAACAAC‑AAGCCC‑GCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGA‑tt < 1:43018/237‑1 (MQ=255) cgaaaTTGAGGCCGCCAGCACGGCGTGTACCGAAA‑ACAGCAACGGTAAAA‑TGAA‑CAGGAAGGTAAATTCCAGCGGTTCGGTAATCCCCCCCAG‑CGTGGGGTTT‑AGGG‑TGGGGGGAATCAGC‑AGGCCC‑GCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGA‑TTTTTGAGTTA‑CCGTGCAGGGCAAAA‑CCGCCTTCCGGGAACAACGATTTCAGCg < 1:441009/247‑1 (MQ=255) cGACATTGAGGCCGCCAGCACGGCGTGTACCGCAA‑ACAGCAACGGTGAAA‑TGAA‑CAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAG‑CATGGCGGTT‑TAAG‑TTGCGGGGATCAGC‑AAAACC‑CCCCC‑CTTTACCCCATTTTTCGGGGCGGCAGGGAAGTACATCGCTAAAGAAA‑AACCCCCGGCACCCAAAAA‑TTTTT‑AGTTA‑CCGGGCCGGGGAAAA‑CCGCCCTCCCGGaaaaaacatttcaacg > 1:584311/1‑237 (MQ=255) gacatggggcccgccAGCCCGGCTTGTACCGCAA‑ACAG‑AAGCGTGAAA‑AGAAACGGGAGGGAAAAGCCCAGCGGTTGGATAATTCC‑ACCAGGCATGGGTGTT‑AGGG‑TGGGGGAAACCAGA‑AACCCC‑GCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGA‑TTTTTGAGTTA‑CCGTGCAGGGCAAAA‑CCGCCTTCCGGGAACAACGATTTCAGCgg < 1:718609/241‑1 (MQ=255) atgggggCCCCCAGCACGGCGTGTACCGAAA‑ACAGCAACGGGAAAA‑TAAA‑CAGGAAGGAAACCTCCGGGGGTTCGGTAATCCCCACCAG‑CATGGGGGTA‑AGGG‑TGGCGGGAACCAGC‑AAGCCC‑GCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGA‑TTTTTGAGTTA‑CCGTGCAGGGCAAAA‑CCGCCTTCCGGGAACAACGATTTCAGCGGCTc < 1:792948/246‑1 (MQ=255) ccgccAGCACGGCGTGTACCGCAA‑ACAGCAACGGTGAAA‑TGAA‑CAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAG‑CATGGCGGTT‑TAGG‑GTGCGGGGATCAGCCAAACCC‑CCCCC‑TTTAACCCACTTTTCCGG‑GCGGCAGGGAAGGAAATCGCTAAAAAAA‑AGCCCCCGGCACCCAAAAA‑TTTTT‑AGTTAACCG‑GCAGGGGAAAA‑ACACCCTCCGGGGACAACGATTTTTGCGGCTccgggccc > 1:475638/1‑244 (MQ=255) cGCAA‑ACAGCAATGGAGAAA‑TAAA‑CGGGAAGGTAAATTCCAGCGGTACGTTAATCCCCCCCAG‑CAGGGCG‑TT‑AAGGGTTGCGGGAATCAGC‑AAGCCC‑GCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGA‑TTTTTGAGTTA‑CCGTGCAGGGCAAAA‑CCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAt < 1:115803/251‑1 (MQ=255) aaCGGTGAAA‑TGAA‑CAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAG‑CATGGCGGCT‑AAAG‑GTGTGGGGAACAGG‑AAAGCC‑GGCAC‑TTTTTACCGATTTTTCGGGGCGGCAGGGGAGGACATCGCTAAAGAAAATGCC‑ACGGCAAC‑AAAGAATTTTTG‑GTTA‑TCGTGCAGGGAAAAAACCGC‑TTCCGGGAACCACGAATTCAGCGGACCGGGGCTCGAAATGAACTCTTGAAGAtaaggcgcccc > 1:778615/1‑241 (MQ=255) aCGGTGAAA‑TGAA‑CAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAG‑CATGGCGGTT‑AAAG‑GTGCGGGGATCAGC‑AAACCC‑GCCAC‑TTTTTCCCCATTTTTCGGGGGGGCCGGGAAATACCTCGCTAAAGAAA‑AGCCCACCGCACC‑AAAAA‑TTTTTGAGTTT‑ACGTGCAGGGCAAAA‑ACGCCCTCCGGGAAAAAAGATTTCAACGGGCCGGCGCTCCAACTGAACTCTTGCAGAtagctcgcccac > 1:325157/1‑240 (MQ=255) agggtaaaaggAA‑CAGAAAGGTAAATCCCAGCGGTTCGGTAATCCCCCCCAG‑CAGGGGGGTA‑AAGG‑TTGGGGGAATCAGA‑AAGCCC‑GCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGA‑TTTTTGAGTTA‑CCGTGCAGGGCAAAA‑CCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAg < 1:67058/241‑1 (MQ=255) tGAAA‑TGAA‑CAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAG‑CATGGCGGAT‑AAGG‑ATGCGGGGAACAGC‑AAGCCC‑GCCAC‑TTTTTCCCGATTTTCCGGGGCGGCAGGGAAGGACAACGCTCAAGAAA‑TACCCACGGCCCC‑CAAGA‑TTTTTGAGTTA‑ACGTGGAGGGGAAAA‑CCGCCCTCCGGGAACAACGATTTCAGCGGCGCGGCGCTCCAAATGAACTCTCTCACATGGTGCGcccccgacc > 1:59860/1‑245 (MQ=255) gtaaaagaaaagggAAGGTAAACTCCAGGGGTTCGTTAATCCCCAACAG‑CAGGGGGGTT‑AGGG‑GGGGGGAAAGCAG‑‑AAGCCCCGCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGA‑TTTTTGAGTTA‑CCGTGCAGGGCAAAA‑CCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTAc < 2:667204/239‑1 (MQ=255) aGGTAAATTCCAGCGGTTCGGTAATTCCCACCAG‑CATGGCGGTT‑AAAG‑GTGCGGGGATCAGC‑AAGCCC‑GCCCC‑TTTTACCCGATTTTTCCGTGCGGCAGTGGAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGA‑TTTTTTAGTTA‑CCGTGCAAGGCAAAA‑CCCCCCTCCGGGGACAACGATTTCAACGGCTCGGCGCTCAAACTGAACTCTTTCTGATACTGCGCCCCGTACATCtataagacgcct > 1:434041/1‑240 (MQ=255) aaatcccggCGGTTCGTTAATACCCCCCAG‑CTGGGCGTTT‑AGGG‑TGGGGGGAACCAGC‑AAGCCC‑GCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGA‑TTTTTGAGTTA‑CCGTGCAGGGCAAAA‑CCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCaa < 1:362015/243‑1 (MQ=255) aaTTCCAGCGGTTCGGTAATTCCCACCAG‑CATGGCGGTT‑AAGG‑GTGCGGGGATCAGC‑AAGACC‑GCCAC‑TTTTACCCGCTTTTTCGGTGCGGCAGGGAAGGACATCGCTAAAGAAA‑TACCCACCGCACC‑AAAGA‑TTTTTTAGTTT‑CCGTGCAGGGCAAAA‑ACGCCTTCCGGGGACAACGATTTCAACGGCTCGGGGCGCAAACTCAACCCTTGCAGGTGCTGCGCCCCGTACATCTTAATTCCGCCtttcac > 1:814026/1‑247 (MQ=255) ggcgttacgttAATCCCCACCAG‑CAGGGCGGTT‑AAGG‑TTGCGGGAACCAGC‑AAGCCC‑GCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGA‑TTTTTGAGTTA‑CCGTGCAGGGCAAAA‑CCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGc < 1:382568/241‑1 (MQ=255) ggtgcggaattcccccccaG‑CAGGGCGGTT‑AGGG‑TGGTGGGAACAAGC‑AAGCCC‑GCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGA‑TTTTTGAGTTA‑CCGTGCAGGGCAAAA‑CCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCgg < 1:430565/235‑1 (MQ=255) ggTTCGGTAATTCCCACCAG‑CATGGCGGTT‑AAGG‑TTGCGGGAATCAGC‑AAGCCC‑GCCCC‑TTTTACCCGATTTTTCGGGGCGGCAGTGGAGTACATCGCTTAAGAAA‑ATCCCCCGGCACC‑AAAGA‑TTTTTTAGTTA‑CCGGGCAGGGCAAAA‑CCGCCTTCCGGGAACAACGAtt > 1:37196/1‑172 (MQ=255) ggTT‑AAGG‑ATGCGGGAATCAGC‑AAGCCC‑GCCAC‑TTTTACCCGATTTTCCGGTGCGGCCGGGAAGTACATCCCTAAAGACA‑TACCCACGGCCAC‑AAAGA‑TTTTTGAGTTA‑CCGTGAAGGGCAAAA‑CCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCCGATGCTGCGCACAGCACATCTGAATGCCCCCTTCAACACATGCCGGACCAAAAATAAACTGTCCGTAGATgg > 1:112203/1‑250 (MQ=255) cGGGGATCAGC‑AAGCAA‑GCCAC‑TTTTACACGATTTTACGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCAACGGCACC‑AAAGA‑TTTTTGAGTTA‑CCGTGCAGGGCAAAA‑CGGCCTTCCGGGAACCACGATTTCAGCGGCTCGGCGCGCACACTGACCTCTTGCCCATGCTGCGCCCACTACCGCTGCATACCGCCTTCAACAGCGGCCCGCCCAAAGATCCACTGTCCGTAGATGCAgtgtggtaa > 2:264940/1‑244 (MQ=255) | AATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCAT‑GTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAA‑ACAGCAACGGTGAAA‑TGAA‑CAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAG‑CATGGCGGTT‑AAGG‑TTGCGGGAATCAGC‑AAGCCC‑GCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAA‑TGCCCACGGCACC‑AAAGA‑TTTTTGAGTTA‑CCGTGCAGGGCAAAA‑CCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAA > NC_000913/3862428‑3862913 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |