| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,161,142 | 0 | C | A | 22.2% | 159.0 / 13.8 | 54 | Q26K (CAA→AAA) | fabR | DNA‑binding transcriptional repressor FabR |
| Reads supporting (aligned to +/- strand): ref base C (26/16); new base A (0/12); total (26/28) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.17e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAAT > NC_000913/4161093‑4161190 | cagTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGcaac < 1:6219722/50‑1 (MQ=255)cagTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGcaac < 1:3508665/50‑1 (MQ=255)cagTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGcaac < 1:1575956/50‑1 (MQ=255) agTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGcaac > 1:5852997/1‑49 (MQ=255) gTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGcaaca > 1:5707501/1‑49 (MQ=255) tACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGcaacaa > 1:6506010/1‑49 (MQ=255) tACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGcaacaa > 1:5150588/1‑49 (MQ=255) tACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACaaa < 1:6264287/50‑1 (MQ=255) cTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACaaaaa < 1:5198674/50‑1 (MQ=255) cTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACaaaaa < 1:5929293/50‑1 (MQ=255) cTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACaaaa < 1:3923947/49‑1 (MQ=255) cTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACaaaa < 1:1555377/49‑1 (MQ=255) cTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACaaa > 1:3400528/1‑48 (MQ=255) cTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACaaa > 1:2789527/1‑48 (MQ=255) ttGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAg < 1:753465/48‑1 (MQ=255) ggCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAGAGaaa > 1:5447050/1‑49 (MQ=255) gCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaaaa < 1:5504645/50‑1 (MQ=255) gCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaaa > 1:1486127/1‑49 (MQ=255) gCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaaa > 1:4560605/1‑49 (MQ=255) cAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaaaaa < 1:4859102/50‑1 (MQ=255) cAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaaaaa < 1:5356723/50‑1 (MQ=255) cAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaaaaa > 1:1332859/1‑50 (MQ=255) ggATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑ccc < 1:4273698/50‑1 (MQ=255) ggATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑ccc < 1:5804180/50‑1 (MQ=255) gATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑Cccg < 1:3323685/50‑1 (MQ=255) gATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑Cccg < 1:4142778/50‑1 (MQ=255) tCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑Cccgccg > 1:5572140/1‑50 (MQ=255) cAGACATCGTGATGGGCGTAAGGGCGCAACAAAAAGAAAAAC‑Cccgccg < 1:1329939/49‑1 (MQ=38) aGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGtt > 1:31526/1‑50 (MQ=255) gACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTc > 1:215250/1‑50 (MQ=255) aCAGCGTGAGGGGGGAAAGAGCGCAAAAAAAAAAAAAAA‑CCCGCCGTTCg < 1:1875858/50‑1 (MQ=19) cATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGc > 1:5380301/1‑50 (MQ=255) aTCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCt > 1:1198721/1‑50 (MQ=255) aTCGTGAGGGGCGTAAGAGCGAAAAAAAAAAAAAAAA‑CCCGCCGTTCGCt < 1:6270194/50‑1 (MQ=25) cGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCtgg > 1:6035248/1‑50 (MQ=255) gTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCtggt > 1:5456676/1‑50 (MQ=255) gATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCtggtg > 1:6299285/1‑49 (MQ=255) gATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCtggtg > 1:6138718/1‑49 (MQ=255) gggCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAg > 1:1552183/1‑50 (MQ=255) ggCGTAAGAGCGCAACAAAAAGAAAAAC‑CCCGCCGTTCGCTGGTGGAAg < 1:618522/49‑1 (MQ=255) ggCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGc > 1:5566742/1‑50 (MQ=255) gagTAAGAGAGCAAAAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGcc < 1:3163358/48‑1 (MQ=38) gCGTAAGGGCGAAAAAAAAAAAAAAAC‑CCCGCCGTTCGCTGGTGGAAGcc < 1:2939809/50‑1 (MQ=25) gCGAAAGAGCGAAAAAAAAAAAAAAAA‑CCCGCCGTTCGCTGGTGGAAGcc < 1:4820439/50‑1 (MQ=25) cGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCg > 1:5345959/1‑50 (MQ=255) cGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCg > 1:6632310/1‑50 (MQ=255) cGTAAGAGCGCAAAAAAAAAAAAAAC‑CCCGCCGTTCGCTGGTGGAAGCCg < 1:83412/50‑1 (MQ=38) gTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCg > 1:3746531/1‑49 (MQ=255) tAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCa > 1:6293368/1‑50 (MQ=255) aaGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCAt > 1:1709449/1‑50 (MQ=255) aaGAGCGAAAAAAAAAAAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCAt < 1:1491035/50‑1 (MQ=37) agagCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCAtt > 1:5863751/1‑50 (MQ=255) gagCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCAtt > 1:2698565/1‑49 (MQ=255) gagCGCAAAAAAAAAAAAAAC‑CCCGCCGTTCGCTGGTGGAAGCCGCAttt < 1:2701052/50‑1 (MQ=38) ggcgaAAAAAAAAAAAAAAC‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTa < 1:5614782/49‑1 (MQ=37) agCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCAtt > 1:5575332/1‑48 (MQ=255) agCGCAAAAAAAAAAAAAAC‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTa < 1:4175565/50‑1 (MQ=38) agCGAAAAAAAAAAAAAAAC‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTa < 1:181385/50‑1 (MQ=37) gcgaAAAAAAAAAAAAAAC‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAg < 1:2001935/50‑1 (MQ=37) gcagcaAAAAGAAAAAATGCCGCCGTTCGCTGGTGGAAGCCGCATTTAGc < 1:1195524/50‑1 (MQ=37) gcAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGc > 1:5828966/1‑49 (MQ=255) gcAACAAAAAAAAAAAC‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGcc < 1:5110675/50‑1 (MQ=39) aaacaaAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCa < 1:5804167/49‑1 (MQ=255) aacaaAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCaa > 1:4770435/1‑50 (MQ=255) aacaaAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCaa > 1:6427253/1‑50 (MQ=255) acaaAAAAAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAAt < 1:4718191/50‑1 (MQ=255) | CAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAA‑CCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAAT > NC_000913/4161093‑4161190 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |