Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,890 | 0 | C | A | 61.2% | 4.0 / 48.1 | 50 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base C (14/5); new base A (0/30); total (14/36) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.72e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.73e-03 |
GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCG‑GCTG‑TC‑TGCCAGCTATTTTT‑CCG‑CC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTT > NC_000913/357742‑358027 | gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCG‑GCTG‑TC‑TGCCAGCTATTTTT‑CCG‑CC‑Gcaacaa > 1:1145070/1‑151 (MQ=255) aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCG‑GCTG‑TC‑TGCCAGCTATTTTT‑CCG‑CC‑GCAACAAAA‑c > 2:948607/1‑150 (MQ=255) gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCG‑GCTG‑TC‑TGCCAGCTATTTTT‑CCG‑CC‑GCAACAAAA‑ccg > 1:545504/1‑151 (MQ=255) gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCG‑GCTG‑TC‑TGCCAGCTATTTTT‑CCG‑CC‑GCAACAAAA‑ccg > 1:1613410/1‑151 (MQ=255) gCGTAGCGAATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCT‑GCTG‑TC‑TGCCAGCTATTTTT‑TCG‑CC‑GCAACAAAA‑ccg > 1:445459/1‑151 (MQ=255) cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCG‑GCTG‑TC‑TGCCAGCTATTTTT‑CCG‑CC‑GCAACAAAA‑ccg < 1:89342/150‑1 (MQ=255) ggCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCG‑GCTG‑TC‑TGCCAGCTATTTTT‑CCG‑CC‑GCAACAAAA‑CCGCCGTTCTCTCCag > 1:409512/1‑151 (MQ=255) gCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCG‑GCTG‑TC‑TGCCAGCTATTTTT‑CCG‑CC‑GCAACAAAA‑CCGCCGTTCTCTCCagt > 2:546493/1‑151 (MQ=255) ccACCCGCAGCCCGGTGTCGGCGGCTCCCGTGTTTCTTTCCCCCCACCACCCGCGCTTTTACCAACCCCCATGTAAAAAAATCCCTCCGG‑GGCG‑GCTG‑TC‑T‑CCAGGTTTTTTTCCCG‑CC‑GCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGc < 1:220148/151‑1 (MQ=255) ccGCTGATCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCG‑GCTG‑TC‑TGCCAGCTATTTTT‑CCG‑CC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑gggg > 1:109193/1‑151 (MQ=255) ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCG‑GCTG‑TC‑TGCCAGCTATTTTT‑CCG‑CC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑gggg > 1:1406907/1‑151 (MQ=255) gACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCG‑GCTG‑TC‑TGCCAGCTATTTTT‑CCG‑CC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑ggggggagg > 1:915892/1‑146 (MQ=255) tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCG‑GCTG‑TC‑TGCCAGCTATTTTT‑CCG‑CC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCCAGAGaaa > 2:670461/1‑151 (MQ=255) gtgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCG‑GCTG‑TC‑TGCCAGCTATTTTT‑CCG‑CC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGGAGaaa > 1:1620230/1‑151 (MQ=255) gtgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCG‑GCTG‑TC‑TGCCAGCTATTTTT‑CCG‑CC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑ggggggacgggaaa > 1:1730209/1‑141 (MQ=255) tccggCGGCTCCTGTGTTTCTTTCCCCCCACCACCCGCGCTTTTACCAACCGCCATGTGAAAAAATCCCTTCGG‑GGCG‑GCTG‑TC‑T‑CCAGGTTTTTTTCCCG‑CC‑ACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAg < 1:1325786/147‑1 (MQ=255) gtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCG‑GCTG‑TC‑TGCCAGCTATTTTT‑CCG‑CC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGggaaaggaaaaaa > 2:879291/1‑148 (MQ=255) tcggcgcCTACCTTGTTTCTTTCCCCCCAACACCCCCCCTTTTACCCACCCCATTTTGAAAAAATCCTTTCGG‑GGCG‑GCTG‑TC‑GGCCAGTTTTTTTC‑CCG‑CC‑ACAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAg < 2:1593542/147‑1 (MQ=255) tCAGGGGCTACCGTGATTATTTCCCCCCAACACCCGCGCTTTTACCCACCGCCATGTGAAAAAATGCCTCCGG‑GGGG‑GTTG‑TC‑GGCCAGCTTTTTTC‑CCG‑CC‑GCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAg < 2:1075527/150‑1 (MQ=255) cgggggCAACGGTGATTCTTTCCCCCAAACACCCGCGCTTTTACCCACCGCAAGGTGAAAAAATCCTTCCGG‑GGCG‑GCTG‑TC‑TGCCGGTTTTTTTC‑CCG‑CC‑GCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGc < 2:1256241/146‑1 (MQ=255) tcccGTTATTCATTCCCCCCAACACCCGCGTTTTTCCCCACCCCCA‑GTGAAAAAAATCCTTCGG‑GGGGTCCTT‑TC‑TGCAGTCTTTTTCC‑CCG‑CC‑GAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGcc < 1:987096/149‑1 (MQ=255) cacggtgtTTCATCCCCCCCACCACCCGCGATTTTATCCAACGCCATGTGAAAAAGTCCCTCCGC‑GGCG‑GCTG‑TC‑TGCCAGCTTTTTTC‑CCG‑CG‑GAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGcc < 1:1090650/147‑1 (MQ=255) aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCG‑GCTG‑TC‑TGCCAGCTATTTTT‑CCG‑CC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGCAGAGgaaaaaaatttatac > 2:1136716/1‑138 (MQ=25) ctttcccGACAACACCCGCGCTTTACTCCACCCCCATGACAAAAAAGACCTCCGC‑GGCG‑GCTG‑CC‑GCCCGGTTTTTTTC‑CCG‑CC‑GCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttg < 1:1075255/146‑1 (MQ=255) catccccccaaacccccGCACATTTCTCCAACGCCAGGTGCAAAAATCCCTTCGC‑GGGG‑GCTG‑TC‑TGCAAGTTATTTTC‑CCG‑CC‑AAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttg < 1:403137/136‑1 (MQ=255) ccccaaacaccCCCGCATTAACCCACCGCCGGGTACAAAAATCCCTTGCG‑GGGG‑GCTG‑TC‑TGCAAGTTTTTTTC‑CCC‑CC‑CCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTgg < 1:1061975/146‑1 (MQ=255) ccccAACACCCCCGCATTCCTCCAACCCCGGGCAAAAAAATCCCTCCGG‑GGCG‑GCTG‑TC‑GGCCAGTTATTTTC‑CCG‑CC‑AAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGc < 1:1764011/149‑1 (MQ=255) ccaaacaccCGCGATTTTACCCACCCCCATGTGCAAAAAGCCCTCCGG‑GGGG‑GCTG‑TC‑GCCCAGCTTTTTTC‑CCG‑CC‑GCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCa < 2:1275665/148‑1 (MQ=255) caaacaccCGCGATTTTACCCACCGCCATGTCAAAAATTGCCTTCGG‑GGGG‑GCTG‑TC‑TGCCAGTTTTTTTC‑CCC‑CC‑GCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCaa < 2:826089/149‑1 (MQ=255) cacccGCGCTTTTACCCACCGCAATGGGAAAAAATCCCTCCGG‑GGGG‑GCTG‑TC‑TGCCAGCTTTTTTC‑CCG‑CC‑GCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc < 2:1542018/148‑1 (MQ=255) cacccGCGCATTTCCCCACCCCCATGTGAAAAAATGCTTCCGG‑GGGG‑GCTT‑TC‑GCCCAGTTTTTTTC‑CCG‑CC‑GCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc < 2:72773/148‑1 (MQ=255) cattaccccAACGCCATGTGCAAAAATGCCTCCGG‑GGCG‑GCTG‑TC‑TGCCAGCTTTTTTC‑CCG‑CC‑GCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa < 1:152231/144‑1 (MQ=38) ttcccccgaccccATGGGAAAAAATGCCTCCGC‑GGC‑‑GCTT‑TCGGGCGAGATTTTTTC‑CCG‑CA‑GAACAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTAGTGTGGCAATCGCTGCTGGTAAc < 1:916950/140‑1 (MQ=255) acccaccccaaTGTGAAAAAATCCTTCCGC‑GGGG‑GCTG‑TC‑GGCAACTTTTTTTC‑CCG‑CC‑GAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt < 2:1715860/141‑1 (MQ=255) cccAACCCCTTGTACAAAAGTCCCTTCGC‑GGCG‑GCTT‑TC‑TCCCGGTTTTTTTC‑CCG‑CC‑GCAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg < 1:1012489/150‑1 (MQ=255) ccaaccccctggcaaaaaaaGCCCTGCGC‑GGCG‑CTTG‑TC‑TGCCGGTTTTTTTC‑CCG‑CC‑GAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg < 1:292997/136‑1 (MQ=255) ccacccctatgtcaaaaatTCCCTGCGC‑GGCG‑GCTG‑TC‑TGCCAGTTTTTTTA‑CCG‑CC‑GCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt < 1:976903/137‑1 (MQ=255) cccgtttcAAAAAAGCCTTTGG‑GGGG‑GTTGTTC‑TGCCATTTTTTTTT‑CCG‑CC‑GAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaaa < 2:889788/144‑1 (MQ=255) atgtaaaaaaaTCCTTGCGG‑GGGG‑GCTG‑TC‑GGCCAGTTTTTTTC‑CCG‑CC‑GAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg < 2:1093109/145‑1 (MQ=255) tgaaaaaaaTGCCCTGGGGGGCG‑GCTG‑TC‑GGCCAGTTTTTTTC‑CCG‑CC‑GAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa < 1:768995/147‑1 (MQ=255) gtcaaaaaaTGCCTCCGC‑GGCT‑GTTG‑TC‑T‑CCAACTTTTTTT‑CCC‑CCGGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa < 2:1144492/146‑1 (MQ=255) ctttttcccactttttttc‑CCG‑CC‑ACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc < 2:672838/137‑1 (MQ=255) cgcatggCCCGCTATTTGC‑CCG‑CC‑AAAAAAAAT‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGGGGGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc < 1:652104/145‑1 (MQ=255) ccggCCAGGTATTTTC‑CCG‑CC‑ACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcg < 2:369335/147‑1 (MQ=255) tttttttc‑CCG‑CC‑GAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg < 2:1423202/148‑1 (MQ=255) cttttttc‑CCC‑CC‑GAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg < 2:1494916/148‑1 (MQ=255) ttttttccccg‑cc‑aAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg < 1:1397620/143‑1 (MQ=255) tttt‑CCC‑CC‑CAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt < 2:24324/151‑1 (MQ=255) tttc‑CCG‑CC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt < 1:813270/151‑1 (MQ=255) ttt‑CCGCCC‑GAAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt < 2:638152/151‑1 (MQ=255) ttcccg‑cc‑cCAAGAAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt < 2:1507096/148‑1 (MQ=255) | GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCG‑GCTG‑TC‑TGCCAGCTATTTTT‑CCG‑CC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTT > NC_000913/357742‑358027 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |