| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 357,875 | 0 | A | T | 44.8% | 29.5 / 13.0 | 29 | N272K (AAT→AAA) | cynR | transcriptional activator of cyn operon; autorepressor |
| Reads supporting (aligned to +/- strand): ref base A (10/6); new base T (0/13); total (10/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.14e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
AACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCG‑TGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGC > NC_000913/357729‑358004 | aaCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGTCTACCG‑TGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAtttt < 1:127498/151‑1 (MQ=255)aaCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCG‑TGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttt < 2:234571/150‑1 (MQ=255) gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCG‑TGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTc < 2:802434/148‑1 (MQ=255) ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCG‑TGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCcgccgc > 1:643509/1‑151 (MQ=255) tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCG‑TGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac > 1:86845/1‑150 (MQ=255) aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCG‑TGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa > 2:759071/1‑150 (MQ=255) gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCG‑TGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑ccgcc < 2:482796/151‑1 (MQ=255) cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCG‑TGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTcc > 2:311803/1‑151 (MQ=255) ccACCCGCAGCCCTGTG‑CCGCCGCCACCGTTGTTTTTTTCCCCCCAACACCCGCATTTTTCCCCACCCCCATGTGAAAAAACCCCTCCGGGGCGGCTGCCCCCCAGTTTTTTTCCCCCCGCAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGc < 1:819509/151‑1 (MQ=255) gtgccggCGGCTCCCG‑TGTTTCTTTCCCCCCACCACCCGCTCTTTTCTCCACCGCCTTGGAAAAAAATCCCTGGGGGGCGGCTGTCGGCCAGTTTTTTTCCCGCCAAAACAAAC‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑Gaa < 1:663304/143‑1 (MQ=255) tgtCAGCGGCTACCG‑TGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑gg > 1:12291/1‑138 (MQ=255) aGCGGCTACCG‑TGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGA‑Gaaaaaaaat > 2:830408/1‑145 (MQ=255) gcgcCTCCCG‑TGATTCATTCCCGCCAACACCCGCGCTTTTACCCACCGCCTTGTCAAAAATTGCCTTCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCtt < 2:709592/150‑1 (MQ=255) aacctttaTTCTTTCCCCCCAACACCCGCGCATTTCCCCACCCCCTTGTGAAAAAAGCCCTCCGCGGCGGTTTTCGGCCAGTTTTTTTCCCGCCGCAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAaga‑g < 1:317380/129‑1 (MQ=255) gtatttCTTTCCCCCCACCACCCGCCCATTTACCCACCCCAATGTAAAAAAATCCCTCCGCGGCGGCTGTCGGCCAGTTTTTTTCCCCCCAAAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGc < 1:381450/135‑1 (MQ=255) ttCTTTCCCCCCAACACCCGCGCTTTTCCCCACCCCCGGTTGAAAAAATGCCTTCGGGGGGGCTTTCCCCCATCTTTTTTCCCGCCGCAACAAAC‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCa < 2:170845/148‑1 (MQ=255) ccttcccccccAACACCCGCGCATTTACCCAACGCCATTTGAAAAAAGCCCGTCGCGGCGGCTGTCGGCCATCTTTTTTCCCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCAtgt < 1:569327/141‑1 (MQ=255) ccccccAACACCCGCGTTTTTACCCACCCCCTTGTGAAAAAATGCCTCCGCGGCGTCTTTCTCCAAGTTTTTTTCCCGCCACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTtgtgtg < 2:529352/147‑1 (MQ=255) cGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCGAGA‑GAAAAAGCtcttattcccctatg > 1:858423/1‑128 (MQ=25) ccccgcTTTTATCCACCCCCATGTGAAAAAATCCTTTCGGGGGGGCTGTCTGCCATTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct < 2:869325/146‑1 (MQ=255) cATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGA‑GAAAAAAATTTAATCCCGCCCTTTTGTTTgggggtcactgccgg > 1:672954/1‑136 (MQ=25) tttACCCACCGCCATGTGCAAAAATGCTTTCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCATCAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGt < 2:100332/148‑1 (MQ=255) acgcctgggcaaaaatTGCCTCCGCGGCGGTTTTCTGCCGGTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTgg < 1:574429/139‑1 (MQ=255) tgtgCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGGAAAAAATTTTAAACCC‑TCCTGGTTTGTGGGGAACGCTGCTGGGAACaaaaggggaaagg > 2:436347/1‑139 (MQ=255) aaaacTCCCTCCGCGGGGCCTGTC‑GCCAAGTTTTTTCCCCCCCAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACa < 1:305758/133‑1 (MQ=255) aaaaaTCCCTCCGCGGCGGCTGTCTCCCAGTTTTTTTCCCCCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg < 2:148687/148‑1 (MQ=255) aaaTGCCGTCGGGGCGGCTTTCTGCCGGCTATTTTCCCCCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg < 1:758866/151‑1 (MQ=255) ccTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAAAG‑GAAAAAACTTTAA‑CCCCCCATGTTTTTGGGGGATCGCCGCTGGGAGAAATATGGGAAGGGGAgggggggg > 2:546616/1‑141 (MQ=255) ggcggttttcccccAGTTTTTTTCCCGCCGAAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGc < 1:311711/136‑1 (MQ=255) | AACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCG‑TGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGC > NC_000913/357729‑358004 |
| 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: ‑ |