| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 357,873 | 0 | C | T | 45.7% | 27.8 / 19.3 | 35 | S273N (AGC→AAC) | cynR | transcriptional activator of cyn operon; autorepressor |
| Reads supporting (aligned to +/- strand): ref base C (12/7); new base T (0/16); total (12/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.26e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCG > NC_000913/357728‑358016 | gAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAt < 1:974972/149‑1 (MQ=255) aaCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAtttt < 2:709239/151‑1 (MQ=255) aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt > 2:52385/1‑151 (MQ=255) aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt > 1:347186/1‑151 (MQ=255) tCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt > 2:519106/1‑149 (MQ=255) tAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑ccgc < 2:109404/150‑1 (MQ=255) aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑ccgcc < 2:139544/150‑1 (MQ=255) ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑Ccgccgc > 2:116498/1‑151 (MQ=255) gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAccg > 1:362792/1‑151 (MQ=255) gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAccgc < 1:433418/150‑1 (MQ=255) cccATCAGCCAATCCCAGCCCCAGACCTGTTTCAGCGGCTACCGTGATTCATTCCCGCC‑AA‑CACCCGCGCATTTATCCAACGCCATGTGCAAAATTGCCTCGGCGGCGGCTGTCTGCCAGTTATTTTT‑CCGCCGCAACAAAAccgccg < 2:8652/146‑1 (MQ=18) cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTcc > 1:689215/1‑151 (MQ=255) aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg > 2:455066/1‑151 (MQ=255) gCCCCACACTTGTGTCTGCGCCTCCCGTAATT‑ATTTCCCCCAAA‑CAGCCGCCCTTTTTCCCACCTCAATGTAAAAAAATCCTTGCGGGGGGGTTGTCGCCCGGTTTTTTCC‑CCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGc < 1:131827/148‑1 (MQ=255) cGCAGCCTTGTGTCAGGGCCTACCGTGTTTCATTCCCCCC‑AA‑CACCCGCGCATTTACCCACCGCCTTGTCAAAAAATGCTTGCGGGGGGGCTGTCTGCCGGTTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGc < 1:73522/151‑1 (MQ=255) cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg > 2:923837/1‑148 (MQ=255) cgtgtgTCAGCGGCTCCGGTGATTCTTTCCCCCC‑AA‑CACCCGCGCTTTTACCCACCGCCATGTGAAAAAAGGCTTCCGGGGCGGCTTTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag < 1:558397/146‑1 (MQ=255) gtCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg > 2:468941/1‑139 (MQ=255) gtCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAaggg > 1:270815/1‑143 (MQ=255) cctttattctttCCCCC‑AACCACCCGCCCATTTACCCACCCCCTGGGGAAAAAATCCCTCGGCGGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCaa < 2:837379/122‑1 (MQ=255) ttCATTCCCCCC‑AA‑CACCCGCGCTTTTACCCACCGCATTGTGAAAAAATGCTTCCGGGGCGCCTGTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑a < 1:922392/148‑1 (MQ=255) ctttcccccc‑AA‑CACCCGCGTTTTTACCCACCCCCATGGGAAAAATTCCCTTCGCGGCGGTTGTCT‑CCAGTTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑At < 1:267735/139‑1 (MQ=255) cccccc‑AA‑CACCCGCCATTTTACCCACCGCCTTGTCAAAAAATGCCTCGGGGGCGGCTGTCGGCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTtgtgtg < 1:1278190/147‑1 (MQ=255) aa‑caccCGCCCATTTACCCACCCCCATGTAAAAAAAGGCTTCCGGGGCGGCTGTCGGCCAGCTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATc < 1:896938/151‑1 (MQ=255) gcgcTTTTACCAAACCCTATGTAAAAAAATCCTTGCGCGGGGGCTGTCTGCCAGTTTTTTTT‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATCgctg < 1:1196114/148‑1 (MQ=255) caaacccgtgggcaaaaaaTGCCTCCGGGGCGGCTGTCTGCAAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATCGCTGCTGGTAACAAt < 2:1323813/135‑1 (MQ=255) accgcttgggacAAAAATGCCTCGGCGGCGGCTGTCTGCCAGTTATTTTC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATCGCTGCTGGTAACAATGTgg < 2:187722/140‑1 (MQ=255) gtgAAAAAATCCCTTCGGGGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg < 2:979972/149‑1 (MQ=255) ggcaaaaaaTGCCTCCGCGGGGGTTGCCGGCGAGTTTTTTCC‑CCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa < 2:310983/147‑1 (MQ=255) aaaatTGCCTCCGCGCGGGCTTTCTGCCAGTTATTTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa < 1:79548/147‑1 (MQ=255) aaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAGTTTAAAACCCTCAATGTT‑TTTGGGGATAGCTGCTGGTTACAACGTgg > 2:783339/1‑138 (MQ=255) aattGCCTCCGGGGCGCCTTTCCGCCGTTTTTTTTC‑CCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg < 2:24848/148‑1 (MQ=255) ggcggtTGTCTCCCGCTTTTTTTT‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt < 2:1329885/144‑1 (MQ=255) cTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAATTTTAAACCCCC‑ATGTTTTGTGGGGAACGCTGCCGGGGACAAAATGGGAAGGGGAGGGCGGGGGAACAACTCCCGaa > 1:1253634/1‑150 (MQ=255) ccAGCTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgcg < 1:506326/148‑1 (MQ=255) | GAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCG > NC_000913/357728‑358016 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |