Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,880 | 0 | T | C | 46.7% | 26.2 / 11.4 | 30 | K271E (AAA→GAA) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base T (9/7); new base C (0/14); total (9/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.40e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATT > NC_000913/357730‑358022 | aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt > 1:347186/1‑151 (MQ=255) aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt > 2:52385/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) ttttttc‑CCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAtt < 1:1153583/147‑1 (MQ=255) | ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTC‑ATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATT > NC_000913/357730‑358022 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |