| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,658,167 | 0 | G | T | 39.1% | 9.7 / 13.1 | 23 | L33F (TTG→TTT) | ynfE | putative selenate reductase YnfE |
| Reads supporting (aligned to +/- strand): ref base G (5/9); new base T (9/0); total (14/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.61e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGG‑TACAACAGGCCCGCG‑AAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATC > NC_000913/1657918‑1658399 | aTACAACATAATCCCACCTTATTACTCATACGCTTCTATTGATATGGATTAATAATTCTTAACCCAATCTGGGTCGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGc > 2:214278/1‑251 (MQ=255)aTACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGc < 2:125794/251‑1 (MQ=255) cccACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTg < 1:336961/250‑1 (MQ=255) cATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCATGCGGGTTTTCGTTGCCGTTTACCCGGCGCAATgcagcagca < 2:519501/251‑1 (MQ=255) atatGGATTAATAATTCTTAACCCAAACTGGGTAGACTCCCACTATTGTTAGCGCGCTCAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCAGCAGAACACTCGCTAAATCTACCGCGATCGGTTCTCTGGCGCTGGCTGCAGGCGGCTTTTCTTTTCAGATTCCCCAGCACAATTGAGCAGCAACGG‑Tacaaccggac > 1:323958/1‑246 (MQ=255) atTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGG‑TACAACAGGCCCGCG‑AAAAAGTgg < 1:477258/251‑1 (MQ=255) cAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGG‑TACAACAGGCCCGCG‑AAAAAGTGGTCTGGGGTGc < 2:458187/251‑1 (MQ=255) gggTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCCTTTACCCTGCGCAATGCAGAAGCCGCGG‑GACCACAGGCCCGCG‑AAAAAAGGGTCTGGGGTGGCTGTTc > 1:94085/1‑251 (MQ=255) ggTATACTCCCGCTATTGTTAGCGCGCTAAATATTCAAGATAGAAACTTGTATATAACGATAAAGAACAGGGAGTGAGTAATGTCCAAAAATGAACGAATGGTGGTCATCAGCCGCAGAACACCCGTTAAATCTACCGCGATAGGTTCTCTGGCGTTGGCTGCTGGCGGGTTTTCTTTGCCGTGTACCCTGCGCAATGCAGCATCAGCGG‑TACAACAGGCCCGCG‑AAAAAGTGGTCATGGGTGCCTGTTcc < 2:221855/251‑1 (MQ=255) tctATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTTCCCTGCGCAATGCCGCAGCAGCGG‑TACAACCGGCCCGCG‑GAAAAATGGGCTGGGGTGCCTGTTCCGttcactttgg > 1:27492/1‑242 (MQ=255) ctATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCACAAATGCACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACAGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTGCCGTTTTCCCTGCGGAATACACCCCCA‑CGGCTAAAACCGGCCCGCG‑AAAAAGGGGGATGGGGGGCCTGTTccccctactgtggg > 1:544300/1‑239 (MQ=255) atatTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTTCCCTGCGCAATGCAGCAGCAGCGG‑TACAACAGGCCCGCG‑GAAAAGTGGTCTGGGGTGCCTGTTCCGTCCACTTTGGTAGCCGCTGTGCACTTCGt > 1:151013/1‑251 (MQ=255) cAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGGTTTTCTTTTCCGTTTACCCTGCGCAATGCACAAGCAGAGG‑GACCACAGGACCCCG‑AAAAAGTGGGCTGGGGTGCCTGTTCCGACAACTGTGGGAGCCGCGGTTCAATTCGTCTAcc > 2:464638/1‑250 (MQ=255) gtttATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACTAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGGTTACCCTGCGCAATGCAGCAGCAGCGG‑TACAACAGGCCCGCG‑AAAAAGTGGTCTGGGATGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATaa < 2:495755/250‑1 (MQ=255) tatataACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGCATGGAGGGCATCAGCCGCAGAACACTCGCTACATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCCTTTTCCCTGCGAAATGCAGCCGCAGCGG‑AACAACACGCCCGCG‑CAACAGTGGTCTGGGGGGTCTGGTCCCTCTAATGTGGTAGACGCCGTGCACATCGCCTACATGTTAAAGATAATg > 1:345146/1‑250 (MQ=255) aaGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTTCCGTTTACCCCGCGCAATGCAGCAGCAGCGG‑GACAACAGGCCCGCG‑AAAAAGTGGGCTGGGGGGGCTGTTCCGTCAACTGTGGGAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTgggg > 1:412460/1‑250 (MQ=255) aCAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGG‑TACAACAGGCCCGCG‑AAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGaa < 1:19237/251‑1 (MQ=255) agtgagtTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTTTACCTGCGCAATGCAGCAGCAGCGG‑GACAACAAGCCCGCG‑AAAAAGAGGGCTGGGGGGCCTGTTCCGGCCAATGGGGGAGCCGCTGTGGACTTCGTCTACATGTTAAAGATAATGAAGGGGCCTGGGGGGGAACTGa > 1:57125/1‑250 (MQ=255) agtTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGGTTTTCTTTGCCCTTTAACCTGCGCAATGCAACAGCAGCGG‑GACAACCGGCCCGCG‑AAAAAGGGGGCTGGGGTGGCTGTTCCGTCCACTGTGGTAGCCGCTGTGCACATAGTCTCCCTGTTTAAGATAATGAAGTGACCtcggggtgaacttccaat > 1:375581/1‑233 (MQ=255) ccGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGCTTTTCTTTTCCGTTTTTCCTCCGCAATGGCGCAGCAGGGG‑GGCAAAAGGGCCGCG‑GAAAAGGGGTCTGGGGGGCCTGTTCCGGCAACAGTGGGAGCCGCTGTGGCCCTCGTCTTCATGTTAAAGAAAAAGAAGTGACCTGGGGGGGAACTGCCAATACCGGGAGAGATGGGTAAAGCAAACCTCAGGTAagcg > 2:505909/1‑251 (MQ=255) gAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTTTTTTCCGTTTACCCTGCGCAATGCAGCAGCAGCGG‑TACAACAGGCCCGCGGAAAAGGTGG‑CTGGGGGGCCTGGTTCGTCCACTGTGGTAGCCTCTGTGCACTTCGTCTAAATGTTAAAGAAAAAGAAGGGAACTGGGGGGGAACTGACAAAACCGGGAGCGATTAGTACGGCGACCATctgggacacgcccgt > 1:460576/1‑237 (MQ=255) tgCAGGCGGTTTTTCTTTTCCGTTTTCCCTTCGCAATGCAGCAGCAGCGG‑TACAACAAGCCCGCG‑AAAAAGGGGTCTGGGGGGCCCGTTCCGGCCACTGTGGGAGCCCCTGTGCACTTCGGCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAAATGACCATACCGGCAGCGGTGGGGACGGGAAACATCCGGGACGCGCCCGTTTGCGCCGTCGCTCCCTCCGCCGGccgatcgggagagcggcgg > 2:134079/1‑233 (MQ=255) tgCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGG‑TACAACAGGCCCGCG‑AAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGc < 1:134079/233‑1 (MQ=255) | ATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGTTTTTCTTTGCCGTTTACCCTGCGCAATGCAGCAGCAGCGG‑TACAACAGGCCCGCG‑AAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCACTTCGTCTACATGTTAAAGATAATGAAGTGACCTGGGTGGAAACTGACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATC > NC_000913/1657918‑1658399 |
| 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: ‑ |