Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,890 | 0 | C | A | 45.7% | 22.5 / 19.4 | 35 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base C (16/3); new base A (0/16); total (16/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.14e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01 |
CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGA > NC_000913/357740‑358038 | cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac > 1:674367/1‑151 (MQ=255) gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca < 2:35841/150‑1 (MQ=255) aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa > 2:1057237/1‑150 (MQ=255) aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa > 2:325787/1‑150 (MQ=255) gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑ccg < 1:203072/151‑1 (MQ=255) cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACACCCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGGCGCAACAAAA‑ccg > 1:22498/1‑150 (MQ=255) gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑ccgcc > 1:387783/1‑151 (MQ=255) gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑ccgcc > 2:1180561/1‑151 (MQ=255) cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTtc > 2:267239/1‑151 (MQ=255) aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCagt > 2:601283/1‑149 (MQ=255) tCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTgg > 1:976714/1‑151 (MQ=255) ccccagccttgtgCCGGCGCCTCCCGTGTTTCTTTCCCGCAACCACCCGCCCTTTTACCCACCCCCTTGTCAAAAAATCCTTCC‑GCGGCGGCTTTCCGCCAGTTTTTTTCCCCCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑gggg < 1:480884/142‑1 (MQ=255) gccttgtgTCAGCGTCTCCCGTGTTTCATTCCCCCAAACACCCGCGCATTTACCCACCGCCATGTGAAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCa < 2:820958/144‑1 (MQ=255) gACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑ggggg > 2:1171124/1‑146 (MQ=255) gACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑ggg > 1:423933/1‑145 (MQ=255) gCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑gggggtagggg > 1:46469/1‑135 (MQ=255) gcccccctgtaTTCTTTCCCGCCACCACCCGCGCTTTTACCCACCGCCAGGTACAAAAATGCCTCC‑GCGGCGGTTGTCCGCCAGCTTTTTTTCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCttt < 2:948479/139‑1 (MQ=255) tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGTCGTTCTCTCCAGTAGTGGC‑GGGTAGAGG‑GAAATAgttgtattac > 1:846512/1‑142 (MQ=21) aCCGTGTTCCTTTCCCCCCAACACCCGCGCTTTTACCCAACGCCATGTAAAAAAAGGCCTCC‑GGGGGGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAg < 2:548822/148‑1 (MQ=255) tgtttCTTTCCCCCCAACACCCGCGCTTTTACCCACCCCCTTGTGAAAAAGTCCCTCG‑GGGGCGTCTGTCGCCCAGCTTTTTTCCCCCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCt < 2:1184059/136‑1 (MQ=255) ttCATTCCCCCCAACACCCGCTCATTTACCCACCGCATTGTGAAAAAATGCCTCC‑GCGGCGGCTGTCTGCCAGTTTTTTTTCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCAt < 2:438928/149‑1 (MQ=255) cATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAAAA‑GAAAAAGTTTTataaccaagattt > 2:788738/1‑138 (MQ=25) cGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGGAAAAAACTTTAAAcccccc > 2:217537/1‑137 (MQ=25) aacaacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGA‑GAAAAAGAATTAAACCCCTCCTGTTTTGTgg > 2:826271/1‑146 (MQ=37) ccccccgcattgtcaaaaaaTGCTGCC‑GCGGCGGCTTTCTCCCAGTTTTTTTCCCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg < 1:853757/136‑1 (MQ=255) ccccaatttgaaaaatTCCCTGC‑GGGGCGCTTTTCTGCCGGTTTTTTCCCCGCCGCAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg < 1:418149/136‑1 (MQ=255) accgctattcgaaaaaTTCCTGCGGCGGCGGCTGCCCCCCAGTTTTTTTCCCACCGAAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGt < 1:30783/130‑1 (MQ=255) gCCATGTGAAAAAATCCCTCC‑CAGCCGTGTGTCCGCCTGATATTTTCCCCCCACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACaa < 2:271472/143‑1 (MQ=255) aTGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGG‑GAAAAAAAATTAAACCCCCCCTTTTTTTGGGCGATCGCGGCTGGGAACAAagtg > 1:1133250/1‑144 (MQ=25) aaaTGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGA‑GAAAAAACTTTTAAACCCCCCTGTTTTGGGGCGATCGCTGCTGGGGACAATATGGGAAGGGGAgggggg > 2:81653/1‑146 (MQ=255) gcggcgGTTGTCTCCCGGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATc < 1:1009975/147‑1 (MQ=255) ggcggcTTTCTGCCTGTTTTTTGCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGc < 1:17732/148‑1 (MQ=255) gcctgcctgcCCTCTTTTTCTCCGCCTCAAAAATG‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg < 1:408156/148‑1 (MQ=255) ctgacattttttttCCCCCCAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcg < 1:267232/140‑1 (MQ=255) gCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGA‑GAAAAAAATTTTAACCCCGCATGTTTTGTGGGGAACGCTGCTGGGGAAAAAGTGGGAAGGGGAGGGGGGGGGAACAACGCCAGaaaaccgc > 2:299021/1‑144 (MQ=255) ccAGTTTTTTTCCCCCCGCAAAAAAT‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc < 2:781044/147‑1 (MQ=255) tttcCCGCGGAAACAAAG‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt < 1:811268/150‑1 (MQ=255) cccgcccCAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc < 1:376815/148‑1 (MQ=255) ccgccgCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCg < 2:178300/149‑1 (MQ=255) caaAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAAA‑GAAAAAACTTTTAACCCGCCCTGGTTTGTGGGAACCGCTGCTGGGGAAAAAATGGGAAGGGGAGGGGGGCGGAACAACCCCCGAAAAGCGCCTAATGAAGTGGCCTCAAtta > 2:41536/1‑147 (MQ=255) | CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGA > NC_000913/357740‑358038 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |