Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,875 | 0 | A | T | 44.4% | 51.3 / 12.2 | 36 | N272K (AAT→AAA) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base A (16/4); new base T (0/16); total (16/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.07e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACC > NC_000913/357735‑358013 | tAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCCACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑ccgcc > 2:240065/1‑151 (MQ=255) ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑Ccgccgc < 1:706442/151‑1 (MQ=255) ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGca > 1:665604/1‑150 (MQ=255) aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCTTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAAcc < 2:728787/151‑1 (MQ=255) gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAAcc < 2:786506/150‑1 (MQ=255) aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTtctct > 1:17788/1‑151 (MQ=255) aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTtctct > 1:476869/1‑151 (MQ=255) cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTtctct > 2:734458/1‑149 (MQ=255) attcccAGCCGCAGCCCTGTGTCAGCGGCTCCCGTGATTCTTTTCCGCCAACACCCGCGCATTTACCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTTTCTGCCAGTTTTTTT‑C‑CCGCCGCAAAAAAACCGCCGTTCTCTCCagtag < 1:181352/147‑1 (MQ=21) gCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg > 1:430604/1‑151 (MQ=255) ccccAGACCTTTTTCCGCGCCTACCGTGATTCATTTCCCCCAACCCCCGCGCATTTTTCCAACGCCTGGTGCAAAAATCCCTTCCC‑GGGGTCTGTCGTCCAGCTTTTTT‑C‑CCGCCGCAACAAAACCGCCGTtc < 2:815484/130‑1 (MQ=255) cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggtg > 1:678193/1‑146 (MQ=255) cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg > 2:984825/1‑147 (MQ=255) aGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg > 2:832828/1‑147 (MQ=255) cccTGTGCCAGCGGCTCCGTTTTTTCTTTCCCGCAAACACCCGCCGATTTACCCAACCCCATGCACAAAATGCCCTCCGG‑GGCGGCTGTCTGCCAGTTTTTTT‑C‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga < 2:508847/149‑1 (MQ=255) aCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggga > 1:309704/1‑145 (MQ=255) cttgtgTCGGCGGCTCCCGTATTTCATTCCCGCCAACACCCGCGCTTTTATCCACCCCCATGTGAAAAAATCCCTTCGC‑GGCGGCTGTCTGCCAGTTTTTTT‑C‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag < 1:585679/146‑1 (MQ=255) ccTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg > 2:834207/1‑143 (MQ=255) tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggagggaa > 2:959260/1‑142 (MQ=255) tgtgCCGGCGGCTACCGTATTTCTTTCCCCCCAAAACCCGCGCTTTTACCCAACCCATTGTGAAAAAAGGCTTCCGG‑GGCGGCTGTCTGCCAGTTTTTTT‑C‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga < 1:263130/147‑1 (MQ=255) ctggATCCATTCCCGCCAACCCCCGCTCATTTCCCCAACCCCGTGTAAAAAAATGCTTTCGG‑GGCGGCTTTCTGCCAGTTTTTTT‑C‑CCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTa < 2:890070/141‑1 (MQ=255) ttcccGCAAAAACCCGCGCTTTTACCCACCGCAATGTGAAAAAATGCCTCCGC‑GGGGGTTGTCTGCAAGTTTTTTT‑T‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCAtgttg < 1:34897/146‑1 (MQ=255) cccccaaacaccCGCGCTTTTACCCAACCCAATGTGAAAAAATGCTTCCGC‑GGCGGCTGTCGGCCAGTTTTTTT‑C‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTtgtg < 1:841971/143‑1 (MQ=255) ccGCAAACACCCGCGCATTTACCCAGCCCCTTGTAAAAAACTGCCTTCGG‑GGGGGCGTTCCGCCAGCTTTTTT‑C‑CCGCCACAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCAtgttgt < 1:529052/147‑1 (MQ=255) aacaccCGCGCATTTATCCAACCCCAGGTGCAAAAATGCCTTCGG‑GGGGGCTGTCTGCCAGCTTTTTT‑C‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGc < 1:218647/147‑1 (MQ=255) ccGCGCTTTTCTCC‑ACCCCATGGAAAAAAAGTGCTTCGCGGGCGGCTGTC‑CCCAACTTTTTT‑TCCCGCAAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct < 2:560931/149‑1 (MQ=255) cgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑GGAAAAAATTTTAACCCCTCCTGTTGTGTGGGAGTCGCTGCTgg > 1:69191/1‑151 (MQ=25) cAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAA‑GAAAAAAATTTTAACCCCCCCTGTTTTGTGGGGAACGCTGCTGGGGACaa > 1:156832/1‑147 (MQ=25) ggcggtTGCCTGCCATATTTTTT‑G‑CCGCAAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTc < 2:971473/150‑1 (MQ=255) ggcggcTGTCTGCAGGTTTTTTT‑C‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa < 1:834432/146‑1 (MQ=255) ggcggcTGTC‑GCCAATTTTTTTCT‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGc < 1:485496/148‑1 (MQ=255) ggcggcGGTCTCCCAGTTTTTTT‑C‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa < 1:593768/146‑1 (MQ=255) gcggcTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGGAAAAAATTTTAAACCCCCCCAGTT‑TGTGGGAAACGCTGCTGGGGAAAAAATGGGAAGGGGAGTGGGGGGGAACAACTcc > 1:778456/1‑150 (MQ=255) ggcTGTCTGCCAGTTATTTT‑C‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa < 1:835601/149‑1 (MQ=255) tcggCCAGTTTTTTT‑C‑CCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa < 2:464625/141‑1 (MQ=255) cTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑GAAAAAAATTTTAAACCCCCCTGTTGTGGGGGAAACGCTGCTGGGTAAAAAATGGGAAGGGGAGGGGGGGGGAACAGCTCCAGaaac > 1:641567/1‑146 (MQ=255) | TAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACC > NC_000913/357735‑358013 |
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: ‑ |