| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,161,142 | 0 | C | A | 26.5% | 39.4 / 10.4 | 34 | Q26K (CAA→AAA) | fabR | DNA‑binding transcriptional repressor FabR |
| Reads supporting (aligned to +/- strand): ref base C (21/4); new base A (1/8); total (22/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.12e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTG > NC_000913/4160903‑4161382 | ttgTCTGCTTCTTACTTTTGACCTGCCGTTTTATTGGCTTTTTATCGACGATACTCCCCGTTTAAGCGGCAGGTTTTCGATGTCCGTAAAAGAACCGCCCCAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCCAGGATTCAGACATCGTTATGGGCGTAAGAGCGCAACAAACAGaaaa > 2:425779/1‑250 (MQ=255) gATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCa > 1:436627/1‑250 (MQ=255) gattctctCCATTTAAGCGTCAGGTTTCCGCTGTTAGTAAAAGAACCTACAAAAGAATTGTAGTAAATTTGTTTTATTGCGTTACTGTTATTTCACATTTCTGGAGCAATCCAGTATTTTCTTTCCGTGGTATAGTGCCAGCAGTACTGTTGGCAAGGATTCAGACTTCGTGAGGGGCGTAAGCGCCCAAAAAAAAGAAAAAC‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAAc < 2:247032/247‑1 (MQ=255) gaTACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGGGGAAGCCGCCTTTAGCCAATTAAGTGCTgcac > 2:112756/1‑248 (MQ=255) gaTACTCTCCGATTAAGAGGCAGGTTTCCGCTGTCCGTAAACGAACCTGCCAAAGAATTGCAGTACATATGTTTTATTGCGCTACAGTTCATTCACAATCCTGGAGCAATCCAGAATGTTAATTCTCTGGTATAGTGACAGCCGTACTTTTGGCACGGCTTCAGACATCGTGATGGGCGTAAGAGCGCAAAAAAAAGAACAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAAc > 2:179630/1‑251 (MQ=255) tACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCACAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTCAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGACGCCGCATTTCTCCAATTAAGTGCTGAACGc > 1:157590/1‑251 (MQ=255) ctcCGCTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAAACGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCt > 2:435888/1‑251 (MQ=255) gTACGTAAAAGAACCGCCAAAAGAATTGCAGTAAATTTGTTTTTTTGCTTTCCCGTTCATTCAAAATACGGGAGAAATCCAGTATGTTCATTCTCGGTTAAAGTCCCAGCATTCCTTTTGGAAAGGATTCAGACATCGGGATGGGCGTAAGAGCGCAAAAAAAAAAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTg < 2:358562/251‑1 (MQ=255) tcccgttcggAACCGGCCAAAGAATTGCATTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCATTATTTCCTTTCTCTGGTATAGGGCCAGCAGTTCTTTGGGCAAGTATTCAGACATCGGGATGGGCGTAGGACCGCAACAAAAAAAAAAAA‑TCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTg < 1:191526/241‑1 (MQ=255) gTAAAAGAACCGGCCCAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGgcgc > 1:370920/1‑251 (MQ=255) gAACCGACCAACGAATTGCAGTAAATATGTTTTATTGCGTTCCCGTTCATTCACACTACTGGAGCAATCCAGTATGTTCATTCTCTGGTCTAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTg > 2:191526/1‑241 (MQ=255) cGGCCAAATAATGGCAGTAAATATGGTTGATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATTTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACTTCTTGATGGGCGAAAGAGCGCAGCAAAA‑GAAAAAAGCCCGCCGTTCGCTGTTGGAAGCCGCATTTAGCCAATTAAGGGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCggg < 1:65604/251‑1 (MQ=255) ggCCAAAGAATTGCAGTAACTATGTTTTATTGCGTTACCGTTCATTCACAATGCTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGcaac > 1:155770/1‑144 (MQ=255) ggCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCATACATCGTGATGGGCGTAAGAGCGcaac < 2:155770/144‑1 (MQ=255) aatgtgcataaaaTATGTTTTATTGCTTTCCCTTTCATTCACAATCCGGGACAAATCCAGAAGGTTAATTCTTTGGTTTAGGGCCGGCAGTTCTTTGGGAAAGGATCCAGACACCGTGAGGGGTTAAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTc < 2:378969/241‑1 (MQ=255) gTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACctct > 1:136572/1‑251 (MQ=255) gTTACCGTTCATTCACAATACTGGAGCAATCAAGAAGTTTCATTCTCTGGTATATTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGGGATGGGCGTAAGAGCGCAACAAAAAGAAAAAC‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCcgc < 1:302567/251‑1 (MQ=255) tcaCAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGa > 1:153033/1‑250 (MQ=255) agaaaTCCATTATTTTCATTCTCGGGTAGAGGGCCAGCAGTTCTTTTGGAAAGGATCCAGACATCGGGATGGGCGAAAGAGCGCAAAAAAAAAAAAAAC‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGAcc < 2:370920/248‑1 (MQ=255) tATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAg > 1:285597/1‑167 (MQ=255) gTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTtgatg > 2:161633/1‑248 (MQ=255) ttCATTCTCTGGTATAGGGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGACGTGGCGCGTGAAGCGGGCATTGCTCCCACCTATTTTTATCGGCAATTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATgag > 2:187751/1‑249 (MQ=255) ggTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATg > 1:402900/1‑251 (MQ=255) aggggCCAGCAGTACTTTTGGGAAGGATGCAGACCGCGTGATGGGCGTAAAGGCGCAAAAAAAAAAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTAc < 1:28036/247‑1 (MQ=255) aGTGCCAGCCGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACg > 2:156844/1‑251 (MQ=255) aGTGCCAGCAGTACTTTTGGCAAGGATTAAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCCGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACg > 2:322140/1‑251 (MQ=255) gCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGcc > 1:383147/1‑250 (MQ=255) ccAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCa > 2:283337/1‑250 (MQ=255) gcagTACTTTTGGCAAGGATTCAGACATCGTGGTGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCt > 1:264372/1‑211 (MQ=255) gcagTACTTTGGGCAAGGATTCAGACATGGGGGGGGGCGTAAGAGCACAAAAAAAAGAAAAAC‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCt < 2:264372/211‑1 (MQ=255) gacAAGGATTCAGACATCGTGAGGGGCGTAAGAGCGCAAAAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGgc < 1:342750/249‑1 (MQ=255) aGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGt > 2:226550/1‑250 (MQ=255) aGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGt > 1:150081/1‑250 (MQ=255) aTCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGc > 1:220695/1‑251 (MQ=255) aaggagCGCAAAAAAATGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTGCGCCAGCCTGATTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGGCCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAgt < 1:226291/248‑1 (MQ=255) aaGAGCGCAAAAAAAAAAAAAAC‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTGTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGGATCGCCAAAGGCGGGAgtg < 2:402900/251‑1 (MQ=255) | TTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTG > NC_000913/4160903‑4161382 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |