| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 357,890 | 0 | C | A | 62.5% | ‑3.9 / 36.3 | 32 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
| Reads supporting (aligned to +/- strand): ref base C (9/3); new base A (0/20); total (9/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.84e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.99e-05 | |||||||||||
TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGT‑GTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCG‑CGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAC > NC_000913/357741‑358012 | 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) | TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGT‑GTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCG‑CGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAC > NC_000913/357741‑358012 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |