| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,383,722 | 0 | G | A | 29.6% | 42.1 / 10.8 | 27 | F80F (TTC→TTT) | ycjW | LacI family putative transcriptional repressor |
| Reads supporting (aligned to +/- strand): ref base G (17/2); new base A (0/8); total (17/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.03e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCAATATCATAAATAGTAG > NC_000913/1383509‑1383954 | gCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGGGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGgtt > 2:385879/1‑239 (MQ=255) atatCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGaaa > 1:289458/1‑238 (MQ=255) tatTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCAt > 1:504289/1‑239 (MQ=255) aTCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTATCAGCCGTAGTTTTTGTCGAACCGGTTGAAAttccactaatgt > 2:451245/1‑228 (MQ=255) tCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGc > 1:70988/1‑239 (MQ=255) cGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGt > 1:12621/1‑211 (MQ=255) aCTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTc > 1:526013/1‑239 (MQ=255) ggcaatcTGAGCAATATATTCCTTTTAATCATTGTCGTCTTATTTTTGTTTTCGATTTTTGGCATATCTTCGGCCGGTTGGTGCAACGTTTAAAATATCACGTCGAACCTTTTTTCTTCAGCTTGGGCGGTGATGGCATGCAAAACTTCAGAAAAAAACGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGc < 2:504289/233‑1 (MQ=255) tggactatattaATGCTTTAAATCATTTTCTGTTTATTTTTGTTTTCCATTTTTTGAAGTTCTTCTGCCGAATTGGGCGACGTTTGTAATACCACGTCAAACCCTTCTTCTGCAGCTTGGGCGGTGAGGGAATGCAAACCTTCAGAAAAAAACGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCtt < 2:526013/233‑1 (MQ=255) tGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCGGCTTGTTTGGCTGAt > 1:655500/1‑237 (MQ=255) ttATTTTTGTGTTCGTATTTTTGCAGATTGTCTGCGGGTTGGTGCGGCGTTTCAAATATCACGTCAAAACCTTCTTCTTCAGCTTTGGCGTTGATGGCATGAAAAACTTAAAAAAAAAACGTATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCGGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCg < 2:655500/239‑1 (MQ=255) tctGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGGGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGat > 1:434032/1‑239 (MQ=255) cGGATGGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTGCAGTTTTGGCGGTGATGGCATGCAAAACTTAAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAg < 1:385879/239‑1 (MQ=255) tttGAAATCACACGGAGAACCTTTCTTCCTGCGTTTGGGCGGGAATGGAAAGAAAAACTTAAAAAAAAACCGGATCACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTcgcg < 2:558772/205‑1 (MQ=255) tttGAAATATCGCGTCAAAACCTTCTTCTTCAGTTTTGGCGGTGATGGCATGCAAAACTTAAGAAAAAAGCGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGc < 1:359834/239‑1 (MQ=255) cACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCggg > 1:558772/1‑205 (MQ=255) tCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATg > 1:320534/1‑171 (MQ=255) aaaCCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGt > 1:651055/1‑209 (MQ=255) ttGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCaa > 2:261453/1‑238 (MQ=255) gTGATGACAGGAAAACCTTAAGAAAAAAACGGATGCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAAACACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGc < 2:681405/123‑1 (MQ=255) gATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCaaa < 1:24848/84‑1 (MQ=255) gATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCaaa > 2:24848/1‑84 (MQ=255) cATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATAtt > 1:87118/1‑175 (MQ=255) gCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTa > 1:681405/1‑123 (MQ=255) aaacCTTAAGAAAAAACCGGTTTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCAATATCATaaa < 1:451245/238‑1 (MQ=255) tgaagagAACAACTGATTACCCGCCGTAGTTGTGGTCGAACGGTGAGAAATCCCCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCAATATCATAAAtagtag < 2:104408/236‑1 (MQ=255) cAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCgggg > 2:361107/1‑156 (MQ=255) cAGAAAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCgggg < 1:361107/156‑1 (MQ=255) | GCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGATTGTGCGACGTTTGTAATATCACGTCGAAACCTTCTTCTTCAGCTTTGGCGGTGATGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCAATATCATAAATAGTAG > NC_000913/1383509‑1383954 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |