| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,658,167 | 0 | G | T | 35.3% | 23.2 / 12.5 | 34 | L33F (TTG→TTT) | ynfE | putative selenate reductase YnfE |
| Reads supporting (aligned to +/- strand): ref base G (11/11); new base T (12/0); total (23/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.65e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCG‑AAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAG‑CGATGAGTACGGCAACCATC‑AGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTTGAATTACCCAAT > NC_000913/1657920‑1658415 | acaacaTAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGAGGGCATCAGCAGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGCCGCTGGCTGCAGGCGGTTTTTTTTTG‑CCg > 2:48953/1‑251 (MQ=255) aTCCCACCTTTTTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGGAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTACAAAAATGACCGAATGGTGGGCCTCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGATCTCTGGCGCTGGCTGCAGGCGGTGTTT‑TTTGTCCGttttaacc > 2:293494/1‑246 (MQ=255) aTCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGGTTTTTTTTT‑CCGTTTTCCCt > 2:201582/1‑251 (MQ=255) taCTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCCTGCGCAATgcag < 1:234618/250‑1 (MQ=255) cATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCCTGCGCAATgcagcag < 1:423159/249‑1 (MQ=255) aTACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGCACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTATCCTtcgccctttagctgc > 2:112257/1‑238 (MQ=255) tCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCCTGCGCAATGCAGCAGCAGCGGGACAACAAGCCCGCG‑AAAAAGAGGTc > 1:155590/1‑251 (MQ=255) tAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATATACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTT‑CCGGTTTCCCTCCGCACAGCAGCCGCAAGGGTGCAAACGGCACGCG‑AAACAGTGGTCTgg > 2:138846/1‑251 (MQ=255) aaaTGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTG‑CCGTTTATCCTGCGCCATGCAGCCGCAGCGGGACAACAGGCCCCCG‑AGAAAGGGGGCTGGGGGGccc > 2:191745/1‑250 (MQ=255) ggTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTT‑CCGTTTACCCTGCGCAATGCAGCAGCCGCGGGACACCAGGCCCCCGGAAAAAG‑GGGCTGGGGGGCCTGTTcc > 1:453536/1‑251 (MQ=255) aTTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCG‑AAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGt < 2:17408/249‑1 (MQ=255) cTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTT‑CCGGTTTCCCTGCGCAAAGCAGCAGCAGCGGGACAAAAGGCCCGCG‑GAAAAGGGGGCTGGGGGGCCTGTTTCGTCCACTGTGGTAGGCGCCGTGGACt > 1:155302/1‑251 (MQ=255) taACTATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTT‑CCGTTTAACCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCG‑AAAAAATGGTCTGGGGTGGCTGTTTCGTCAACTGTGGGAGCCGCTGTGCACTTCGTCTTAATGTTTAAAATAATGAAGTg > 1:458404/1‑251 (MQ=255) tAAAGAACATGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCG‑AAAAAGTGGTCGGGGGTGCCTGTTCCGGCAACTGTGGTATCCGCTGTGCACTTc < 1:375289/219‑1 (MQ=255) aaaGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTT‑CCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCG‑AAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTTAAGATAATGAAGTGACCTggg > 1:114011/1‑251 (MQ=255) aaaGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCG‑AAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCt < 1:11169/209‑1 (MQ=255) aaCAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTTCCCTGCGCAATGCAGCAGCAGCGGGACAACAAGCCCGCG‑AAAAAGTGGGCTGGGGGGGCTGTTCCGGCCACTGTTGGAGCCCCTGTGCACTTCGGCCACATGTTAAAAATAAAGaaagggcctggggggg > 2:396296/1‑237 (MQ=255) aCAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGCTTTTTTTTT‑CCGCTTAACCCGCGCAATGCCGCAGCAGCGGTATAACCAGCCCGCG‑AAAAAATGGTCTGGGGGGCCTGTTCCGTCCACTGTGGGAGCCCCTGGGCCCCTTGGCTTCAAGTTTAAGATAATTAAGTGCCCTggggggg > 1:365088/1‑246 (MQ=255) aaGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCG‑AAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGaaa < 2:453536/250‑1 (MQ=255) agtgagtTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCG‑AAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGAc < 1:123464/251‑1 (MQ=255) tATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCCTGCGCCATGCAGCAGCAGCGGGACACCAAGCCCGCG‑GAAAAATGGTCTGGGGGGGCTGTTCCGGCAACAGTGGGTGCCGCTGGCCACTTTCTCTACCTTTTAAAAATAATGAAGGGGCCCGGGGGGGAACTGACAAAAccc > 1:443553/1‑250 (MQ=255) aTGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTT‑CCGTTTACCCTGCGCAATGCAGCAGCAACGGTACAACCAGCCCGCG‑GAAAAGGGGGCTGGGGTGCCCGTTCCGTCCACTGGGGGAGGCGCTGTGCCCCTCGGCTACCTGTTAAAGATAATGAAGTGACCTGGGGGGAAACTGGCAAAACCgg > 1:354570/1‑251 (MQ=255) ccAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTG‑CCGTTTAACCTGCGCAATGCAGCAGGAGCGGTACAACAGGCCCGCG‑GAAAAGAGGTCTGGGGGGCCTGTTCCGTCGACCATGGGAGTCGCTCTGCAATTCGGCTACATGTTTAAGAAAATGAAGAGGCCCGGGTGGGAACTGAAAATACCAGCaac > 1:391474/1‑249 (MQ=255) aaaTGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTT‑CCGGTTACCCTGCGCCAGGCCGCAGCCGCCGTACAACAAGCCCGCG‑AAAAAGGGGGCTGGGGGGCCCGGTCCGTCAACCGTGGTAGGCGCTGTGCACTTTGGGTACCTGTTTAAGATAATGAAATGACCTGGGTGGGAACTGACAAAACCCGCAgggttt > 1:49654/1‑246 (MQ=255) aaTGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACCGGCCCGCG‑GAAAAGCGGTCTGGGGTGCCTTTTCCGCCAACTGTGGTAGCCGCTGTGCACTTCGTCTCCATGTTAAAGATAATGAAATGACCTGGGGGGAAACTGGCAATTACGGCAG‑CGATTAGGACGGc > 1:159471/1‑250 (MQ=255) ggtggGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTT‑CCGTTTACCCTGCGCCAAGCAGCCGCCGCGGTACAATCGGCCCCCGGCAACAG‑GGGCTTGGGGGCCTGTTCCGGCAACTGTGGGCGCCGCTGTGCAACTCGTCTTCATTTTTCAGGAAATGAAGTGACCTGGGGGGAAACTGACAATACCGGCAGGCGA‑GAGTACGGCAAcc > 2:419035/1‑251 (MQ=255) aTCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGATTTTTTTTT‑CCGTTTACCCTTCGCAATGCAGCAGCCGCGGGACAACAAGCCCGCG‑AACAAGTGGTCTGGGGTGGCTGGTCCGTCAACTGTGGGTGCCGCGGTGCCCTTCGGCTACATGTTTAAGATAATGAAGAGACCTGGGGGGGAACTGACAATACCGGGAG‑AGAAGAGTACgggcaacacccggg > 2:12419/1‑239 (MQ=255) cacTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCG‑AAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAG‑CGATGAGTACGGCAACCATC‑AGGTACGCGCCTGTTTg < 1:191745/251‑1 (MQ=255) gcgATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTT‑CCGTTTTCCCTGCGCAAAGGAGCAGCAGCGGGACAACAGGCCAGCG‑GAAAAGGGGGCTGGGGGGGCTGTTTCGTCAACTGTGGTAGCCGCTGTGCACTTCGGCTACATGTTAAAGATAATGGAGGGACCTGGGGGGGAACTGACAAAACCGGCAG‑GGAAGGGGACGGGAACCATC‑AGGGACGCGCCCGTTTGCGCGGGCGCGCCATccg > 1:15352/1‑251 (MQ=255) gTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCG‑AAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAG‑CGATGAGTACGGCAACCATC‑AGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGt < 2:155302/251‑1 (MQ=255) ttCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTG‑CCCTTTAACCCGCGCAATGCAGCCGCAGCGGAACAACAGGCCCCCG‑AAAAAGTGGGCTGGGGGGCCTGTTCCGTCAACTGTGGTTGCCGCTGTGCCATTCGCCTTCATGTTTAAGATAATGAAATGACCCGGGGGGGAAATTACAATACCGGCAG‑CGATGAGTACGGCAACCATCCGGGCAC‑CGCCCGTTTGCGCGGTCGCCCCCTCCGCCCGCGt > 1:77771/1‑250 (MQ=255) gggCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCCGGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCG‑AAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAG‑CGATGAGTACGGCAACCATC‑AGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAat < 2:365088/250‑1 (MQ=255) tgtgggcttgtttattttG‑CCGTTGTCCCTTCGCAAGGCAGCAGCCGCGGTACAAAGGGCCCACG‑AAAAAGTGGTGTGGGGTGTCTGTTCCGACACCTTTGGTGGCCTCTGTTCACTGCTTCTCCAGTTTAAAGCTAATGAAGTGACCTGGGTGGTATCTGACAATACCGGCAG‑CGATGAGTACGGCACCAATC‑AGGTACGCGTCTGTTAGCGCGGTCGCTCCATCCGCCGGCGTATTAATCAGCCCGAt < 2:189094/235‑1 (MQ=255) ttG‑CCGTTTACCCTGCGAAATGCAGCAGCCGCGGTACAACAGGCCCGCG‑AAAAAGTGGTCTGGGGTGCCTGTTCCGTCCACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAG‑CGATGAGTACGGCAACCATC‑AGGTACGCGCCTGTTTGCGCGGTCGCGCCATCCGCCGGCGTGTTAATCATCCCGATCGCTTGAATTACCCAAt > 2:254183/1‑251 (MQ=255) | ACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTG‑CCGTTTACCCTGCGCAATGCAGCAGCAGCGGTACAACAGGCCCGCG‑AAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAG‑CGATGAGTACGGCAACCATC‑AGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTTGAATTACCCAAT > NC_000913/1657920‑1658415 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |