Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,890 | 0 | C | A | 48.4% | 18.4 / 23.1 | 31 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base C (11/5); new base A (0/15); total (11/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.77e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCAATCGCTGCTGGT‑AACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCA > NC_000913/357742‑358034 | 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:385957/151‑1 (MQ=255) aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc < 1:555643/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 > 2:329554/1‑151 (MQ=255) aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtagt > 1:594350/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) ccccccgcAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCAATCGCTGCTGGT‑AACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCg < 1:310911/146‑1 (MQ=255) ccgaAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCAATCGCTGCTGGT‑AACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCa < 1:264622/151‑1 (MQ=255) | GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCA‑TGTTGTGTGGCAATCGCTGCTGGT‑AACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCA > NC_000913/357742‑358034 |
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: ‑ |