Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,880 | 0 | T | C | 43.8% | 36.1 / 11.1 | 32 | K271E (AAA→GAA) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base T (12/6); new base C (0/14); total (12/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.12e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCAATCGCTGCTGGT‑AACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTA > NC_000913/357731‑358019 | cTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTAttttt < 1:675467/150‑1 (MQ=255) cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTccg < 2:597493/151‑1 (MQ=255) ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGca > 2:733787/1‑150 (MQ=255) gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCACACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACATCCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGcaaca < 1:146131/150‑1 (MQ=255) aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGcaacaa < 1:747260/148‑1 (MQ=255) aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc > 1:314963/1‑151 (MQ=255) aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc < 1:555643/151‑1 (MQ=255) aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc < 1:385957/151‑1 (MQ=255) cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCTGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgc > 1:811103/1‑151 (MQ=255) tCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctctc > 1:759031/1‑151 (MQ=255) aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt > 2:485783/1‑149 (MQ=255) aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtagt > 1:594350/1‑151 (MQ=255) aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtagt > 2:329554/1‑151 (MQ=255) cccccccccaccccTGTGTCGGCGCTCACCTTAATTCTTTCCCCCA‑AACACCCGCGCATTTACCAACCCCCATGTGAAAAAATCCCTCCGGGGCGG‑CTGTCTGCAAGTTTTTTTCCCGCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTg < 1:735079/138‑1 (MQ=255) ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg > 1:214886/1‑150 (MQ=255) cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg > 2:700539/1‑150 (MQ=255) agcccTGTGTCGGCGCCTACCGTTTTTAATTCCCGCA‑AACACCCGCGCTTTTACCCAACCCCAGGTGAAAAAATGCCTCCGCGGGGC‑CTTTCTCCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCaa < 2:759126/147‑1 (MQ=255) tCAGCGGCTCCGGTG‑TTTATTTCCCCCAAACACCCGCGCTTTTACCAAACGCAGTGTCAAAAAATCCTTTCGCGGGGG‑CTGTCTGCAATGTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaaa < 2:382761/147‑1 (MQ=255) gccccccGTGTTTCATCCCCCCC‑AACACCCGCGCTTTTCACAACCGCATTGTGAAAAAATGCCTGCGCGGCGG‑CTGTCGGCAAGTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgaga < 1:774401/135‑1 (MQ=255) ccGTGATTCATTCCCCCC‑ACCAACCGCGCTTTTACCCACCGCCATGTCAAAAAATCCCTGCGGGGCGG‑CTGTCGCCCAGCTTTTTTCCCGCCCCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAg < 2:83495/147‑1 (MQ=255) ttCATTCCCCCC‑AACACCCGCCATTTTACCAAACGCCATGTACAAAAAGCCCTCCGCGGCGG‑CTTCCTGCCGGTTTTTTTCCCCCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑t < 1:225395/149‑1 (MQ=255) cc‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCTGAGAGAAAAAATTTAGTCCGTCA‑TGTtgtgt > 1:203953/1‑146 (MQ=25) accacccgcatttttACCCAGCGCCAGGTGAAAAACTGCTTGCGGGGGGT‑CTGTCGCCCAGTTTTTTTCCCCCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCa < 1:65712/136‑1 (MQ=255) ccGCGCTTTTCCCCACCCCCATGCGAAAAACTCCTTCCGGGGCGG‑CTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCAATc < 1:228436/146‑1 (MQ=255) tttACCCAACGCCAGGACAAAAATTTCCTCGGGGGCGG‑TTTTCTGCCGGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCAATCgctgct < 1:785538/145‑1 (MQ=255) ccccATGTAAAAAAATCCCTCCGCGGCGG‑CTGTCTCCCAGTTTTTTTCCCGCCACAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCAATCGCTGCTGGT‑AACAATGTGGaa < 1:853991/149‑1 (MQ=255) aTGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAACTTTTAACCCCTCA‑TTTTTTGTTGGGATAGCTGCCGGTGAACaaatgggaaggg > 2:794437/1‑141 (MQ=255) aaaaaaTCCTTCCGCGGGGG‑CTGTCGGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCAATCGCTGCTGGT‑AACAATGTGGAAAGGGAAGTg < 2:116116/150‑1 (MQ=255) aaaaaTGCTTCGGCGGCGG‑TTGTCTGCCAGCTTTTTTTCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCAATCGCTGCTGGT‑AACAATGTGGAAAGGGaa < 2:159876/147‑1 (MQ=255) aaaactgctgcgcgGGGGCCTGTCTCCCCGTTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCAATCGCTGCTGGT‑AACAATGTGGAAAGGGa < 1:729017/136‑1 (MQ=255) ttCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTTAACCCCTCACTGTT‑TGTTGGCAGCACTGCTGGT‑TACAATGGGGGAAGGGGAgggaggggg > 2:410447/1‑139 (MQ=255) tcggaaggtttttttcCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCAATCGCTGCTGGT‑AACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcg < 1:385731/139‑1 (MQ=255) gcTGCCGGTTTTTTTCCCGCCAAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCAATCGCTGCTGGT‑AACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGa < 1:444663/145‑1 (MQ=255) ccAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCAATCGCTGCTGGT‑AACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa < 1:285797/151‑1 (MQ=255) | CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCAATCGCTGCTGGT‑AACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTA > NC_000913/357731‑358019 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |