Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,890 | 0 | C | A | 50.0% | 15.8 / 24.3 | 34 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base C (14/3); new base A (0/17); total (14/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.77e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01 |
GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCC > NC_000913/357742‑358031 | gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca < 2:161219/150‑1 (MQ=255) aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaaa < 2:122449/150‑1 (MQ=255) gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtct > 1:617630/1‑151 (MQ=255) cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑g > 1:788951/1‑151 (MQ=255) ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑ggg > 2:389249/1‑150 (MQ=255) cgcgacccTGTGTCGCCGGCTACCTTGTTTCTTTCCCGCAACCACCCGCGCATTAACCAACCCCCAGGTCAAAAAATCCCTT‑CGGGGCGGCTGCCGCCCAGTTTTTTTCCCCCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGc < 2:407285/145‑1 (MQ=255) gACCTGTGTCAGCGGCTCCGGTGTTTATTTCCCGCCAACACCCGCCCTTTTCCCCACCCCAATGTGAAAAAATCCCTC‑CGCGGCGGCTGTCTCCCGGTTTTTTTCCCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gggg < 1:196613/146‑1 (MQ=255) cttttgtCAGGGGCTACCGTGATTCTTTTCCCCCAACACCCGCGCTTTTACCCAACGCCTTGTGAAAAAATGCCTT‑CGCGGGGGCTTTCTGCCAGTTATTTTTCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAaga‑g < 2:758470/146‑1 (MQ=255) gtgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGCAGA‑G‑aaa > 2:407010/1‑150 (MQ=255) gtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gggg > 1:845445/1‑139 (MQ=255) cAGCGGCCACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAACGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGGAACAAAACCGCCGTTCTCTCCAGTACTGGC‑GGTCGAAga‑g‑acgaaga > 1:392705/1‑144 (MQ=21) cGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAaga‑gggaaaagct > 1:364790/1‑133 (MQ=255) cGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGC‑G‑aaaaaaa > 2:612424/1‑136 (MQ=255) ttCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGGG‑G‑AAAAAACTTTAAACCCgcc > 2:259955/1‑141 (MQ=25) aacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCCAGA‑G‑GAAAAAATTTTATCCAGTCAtgttgttttggggac > 1:619834/1‑139 (MQ=25) ccGCGCTTTCCTCCACCGCCATGTAAAAAAATCCCTC‑CGCGGCGGCTGTCTGCCAGTTTTTTTTCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct < 1:99916/149‑1 (MQ=255) tttccccACCGCCATGTGAAAAATTGCTTT‑CGGGGCGGTTGTCTGCGAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa < 2:766787/144‑1 (MQ=255) tttACCCACCCCCATTTAAAAAAAGCCCTC‑CGGGGCGGCTTTCGCCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGt < 2:145730/148‑1 (MQ=255) ccAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGA‑GAAAAAAGCTTTAAACCC‑TCCTGTTTTGTGGGGATCggctgtgggtaaa > 2:420469/1‑135 (MQ=25) aaaaaTCCTTC‑CGCGGCGGTTGTCTCCCACCTTTTTTTCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg < 2:439233/148‑1 (MQ=255) aaaaaTCCCTC‑CGCGGGGCCTGTCTCCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt < 2:3894/149‑1 (MQ=255) aaaaGCCTTCCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcg < 2:243066/151‑1 (MQ=255) aaaTGCCTC‑CGCGGCGCCTTTCGCCCACTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg < 2:294702/151‑1 (MQ=255) aaTGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGACGA‑G‑AAAAAGATTTTAACCCCCCATGTTGTGTGGCCATTGCCGCTGGTGAAAAAGTGggaagggaagggggg > 2:406343/1‑145 (MQ=255) aaTGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggc < 2:618075/151‑1 (MQ=255) cgcgGCGGCTGTCTCCCAGTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGc < 1:215249/151‑1 (MQ=255) gcggcgGCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGc < 2:306914/150‑1 (MQ=255) cggcggCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGA‑G‑AAAAAAATTTTAACCCCTCCTTTTGTGTGGGGATCCCTGCTGGTGAAAAAGTGGGAAGGGGAgggggggggagc > 2:215273/1‑135 (MQ=255) gggggcttcctcccAGTTTTTTTCCCGCCCAAAAAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGc < 1:673827/136‑1 (MQ=255) ggcggcTTTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTc < 2:458246/150‑1 (MQ=255) ggcTGTCTCCCAGTTTTTTCCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg < 1:660753/150‑1 (MQ=255) cccgttaaTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa < 2:362261/137‑1 (MQ=255) ccAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc < 2:756897/144‑1 (MQ=255) tATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGG‑AAAAAATTTTAAGCCC‑TCCTGTTGTTTGGGGAACGCTGCTGGGGAAAAAGGGGGAAGGGGAGGGGGGGGGAGCAGCCCCCCAACCCCGCCAAtt > 2:348568/1‑149 (MQ=255) tttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt < 1:844925/150‑1 (MQ=255) cccccacaaaaaAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGcc < 2:701197/145‑1 (MQ=255) | GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCC > NC_000913/357742‑358031 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |