Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,880 | 0 | T | C | 55.6% | 19.6 / 47.3 | 36 | K271E (AAA→GAA) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base T (12/4); new base C (0/20); total (12/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.45e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.23e-06 |
GTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGT‑GTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCG‑CGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAC > NC_000913/357734‑358012 | gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGT‑GTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCG‑CGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccgc > 2:648504/1‑151 (MQ=255) gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGT‑GTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCG‑CGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccg < 2:689593/150‑1 (MQ=255) ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGT‑GTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCG‑CGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCcgccgc > 1:882256/1‑151 (MQ=255) ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGT‑GTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCG‑CGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaa < 1:165058/151‑1 (MQ=255) tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGT‑GTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCG‑CGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac > 1:17086/1‑150 (MQ=255) cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGT‑GTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCG‑CGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑ccgc < 1:532441/151‑1 (MQ=255) gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGT‑GTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCG‑CGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑ccgcc > 1:902901/1‑151 (MQ=255) cATCAGGCAATTCCAGCCGCAGACCTGT‑GTCAGCGGCTACCGTGATTCATTTCCGCAAACACCCG‑CGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑ccgccg < 1:39546/146‑1 (MQ=255) tCAGGCAATTCCAGCCGCAGACCTGT‑GTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCG‑CGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTtctctc > 1:242575/1‑151 (MQ=255) ccGCAGACCTGT‑GTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCG‑CGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgggg > 1:467253/1‑151 (MQ=255) ccGCAG‑CCTGTGGTCGGCGCCTACCGTGATTCATTCCCCCCAACACCCG‑CCCATTTATCCACCGCCATGTGAAAAAATCCCTG‑CGCGGGGGCTGTCTGCCAGTTATTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCg < 2:207121/148‑1 (MQ=255) gACCTGT‑GTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCG‑CGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCggggggc > 2:210835/1‑146 (MQ=255) ctttt‑gtCAGCGGCTCCCGTGTTTCTTTCCCCCCAACAACCG‑CGATTTTACCCACCCCAATGTGAAAAAATGCCTT‑CGCGGCGGCTGTCGGCCAGCTATTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCa < 1:212348/142‑1 (MQ=255) cgtgt‑gTCAGCGCCACCCGTTTTTCATTCCCCCAAACACCCG‑CCCTTTTCCCCACCCCCTTGTAAAAAAGGGCCTC‑CGGGGGGCCTGCCGGCGAGCTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag < 2:23954/148‑1 (MQ=255) ccTGT‑GTCGGCGCCTCCCGTGTTTCTTTCCCCCAAACACCCG‑CGCTTTTA‑CCCACCCCCTTTTCAAAAATTCCTTCCGCGGCGGCTGTCTGCCAGTTTTTTTCCCCCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCa < 1:205042/146‑1 (MQ=255) t‑gtCGGCGGCTCCCGTGATTCATTCCCGCCAACACCCG‑CGCTTTTACCAACCGCCTTGTGAAAAAAGGCTTC‑CGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgaga < 2:697275/147‑1 (MQ=255) t‑gtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCG‑CGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgggg > 2:468308/1‑140 (MQ=255) gCGGTTACCGTGATTCTTTCCCGCCAACACCCG‑CGCTTTTACCCAACCCCATGTGCAAAAATGCCTT‑CGCGGGGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg < 1:827479/147‑1 (MQ=255) gtctcccGTGTTTCTTTCCCCCCAACAACCG‑CGCTTTCACCCACCGCCATGTGAAAAAATGCCTT‑CGGGGGGGCTTTCTGCCAGTTTTTTTCCCGCGACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCttt < 2:687850/144‑1 (MQ=255) ggCTACCGTGTTTCATTCCCGCCAACACCCG‑CGCTTTTACCCACCGCCATGTGCAAAAATGCCTT‑CGGGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATa < 1:243737/144‑1 (MQ=25) cccGTGATTCTTTCCCCCAAACACCCG‑CCCTTTCACCCAACCCCATGTGAAAAAAGCCCTC‑CGGGGCGGCTGTCGGCAAGTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGccc < 2:161838/150‑1 (MQ=255) ccttaaTTCTTTCCCGCAAACACCCG‑CACTTTTTTCCACCGCCATGTGCAAAAATGCTTT‑CGCGGCGGCTGTCGGCCAGTTTTTTTCCCCCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTaa < 1:13715/141‑1 (MQ=255) cGTGATTTCTTTCCCCCCACCACCGCCCCTTTTACCCACCCCTTTGTGAAAAAATGCCTC‑CGCGGGGGTTGTCTGCCACTTTTTTTCCCGCCGCAAAAAAT‑CCGCCGTTCTCTGCAGTAGTGGCGGGGCAAGAGAAATAGCttt < 2:724157/144‑1 (MQ=255) gTGATTCATTTCCGCCAACAACCG‑CGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTaaatccgcc > 1:936966/1‑144 (MQ=25) ttccccccAACACCCG‑GGCTTTTACCCACCCCCTGTTGAAAAAATGCCTC‑CGCGGGGGCTTTCTCCAAGTTTTTTTCCCCCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttg < 1:676979/143‑1 (MQ=255) aacCG‑CGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGGAAAAAATTTTAACACCTCCTGTtgtgt > 1:769066/1‑141 (MQ=25) cgcATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGCGAGAAAAAAATTTTAAACCCCCACGTTGTGTGGGGAACgctgcttg > 1:491795/1‑149 (MQ=25) cgcATTTACCCACCGCCATGTGAAAAAATCCCTT‑CGGGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgct < 1:120145/149‑1 (MQ=255) cccATGAGCAAAAATCCCTG‑CGGGGCGGCTGTCTGCCGTTTATTTTCCCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaaa < 2:25789/150‑1 (MQ=255) aaaaaTGCTTC‑CGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgc < 1:106243/151‑1 (MQ=255) aaaaaTCCCTG‑CGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg < 2:149487/148‑1 (MQ=255) aaaTGCCTC‑GGCGGTGTCTTTCTCCCATCTTTTTTCCCCCCCAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg < 2:444938/151‑1 (MQ=255) ggctggctGCCAGTTCTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa < 2:223438/143‑1 (MQ=255) gTCTGCCAGTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa < 2:481386/148‑1 (MQ=255) gTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAAAGAAAAAACTTTTAATCCCCCCTGTTTTTTGGGGATCGCTGCTGGGGACAAAGTGGGAAGGGGAGGGGGGGGGATCAACTCCAgca > 1:59054/1‑146 (MQ=255) tttacctctttttttCCCGCCCAACAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg < 2:309362/135‑1 (MQ=255) tGCCAGTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc < 2:893362/146‑1 (MQ=255) | GTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGT‑GTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCG‑CGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAC > NC_000913/357734‑358012 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |