| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 357,890 | 0 | C | A | 58.6% | ‑2.4 / 26.5 | 29 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
| Reads supporting (aligned to +/- strand): ref base C (7/5); new base A (0/17); total (7/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.07e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.21e-03 | |||||||||||
ATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTG‑TCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACC > NC_000913/357743‑358040 | aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTG‑TCTGCCAGCTATTTTT‑CCGCCGcaacaa > 2:667458/1‑150 (MQ=255) aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTG‑TCTGCCAGCTATTTTT‑CCGCCGCAACaaaa < 1:344411/150‑1 (MQ=255) aGCGTAGCGCATCAGTCCATTCCACCCGCAGACCTGTGCCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTG‑TCTGCCAGCTATTTTT‑CCGCCGCAACAAAAcc < 2:565990/151‑1 (MQ=255) gTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTT‑CGCAGCGGCTG‑TCTGCCAGCTGTTGTT‑TCGACGGAACAAGAcc > 2:702085/1‑148 (MQ=14) gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTG‑TCTGCCAGCTATTTTT‑CCGCCGCAACAAAAccgcc > 2:649064/1‑151 (MQ=255) cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTG‑TCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTtc > 1:646977/1‑151 (MQ=255) gCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTG‑TCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCagt > 1:401337/1‑151 (MQ=255) gCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTG‑TCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCagt > 2:677882/1‑151 (MQ=255) aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTG‑TCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gg > 1:810992/1‑151 (MQ=255) ccGCAGCTCTGTGTCAGCGGCTCCCGTGTTTCTTTCCCGCCACCACCCGCCCTTTTACCCACCCCCTTGTACAAAAATCCTTC‑CGGGGCGGCTG‑TCGCCCAGTTTTTTTT‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑ggg < 1:246805/150‑1 (MQ=255) gtCGGCGGCTCCCGTACTTCTTTCCCCCCAACACCCGCGCTTTTACCCACCGCCTTTTGAAAAAATGCCTG‑CGCGGCGGCTG‑CCTGCCAGCTATTTCG‑CCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCa < 1:232998/141‑1 (MQ=255) ggCCACGGTGATTCATTCCCCCCAACACCCGCGCTTTTACCAACCGCCATGTGCAAAAATGCCTC‑CGCGGGGGCTG‑CCTGCCAGTTTTTTTC‑CCGCAGAAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTa < 2:370155/150‑1 (MQ=255) cggTGATTCTTTCCCCCAAACACCCGCCCTTTTCCCCACCCCTTTGTGAAAAAATCCCTC‑CGCGGCGGCTT‑CCTGCCGCCTATTTTC‑CCGCCAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGcc < 1:460831/147‑1 (MQ=255) ccgtatttCTTTCCCGCCAACACCCGCGAATTTACCCACCCCCATGTGAAAAAATCCCTC‑GGCGGCGGCTG‑TCTGCCATCTTTTTTC‑CCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTaa < 1:660354/140‑1 (MQ=255) ccGTGTTTCTTTCCCCCCAACACCCGCGCATTTACCCAACCCCATGTGAAAAAATGCCTC‑CGCGGCGGCTG‑TCTGCCAGCTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCg < 1:70850/151‑1 (MQ=255) ttccccccAACACCCGCGCTTTTCCCAACCCCCATGT‑AAAAAAAGCCTTCCGGGGCGGCTT‑TCTGCCAGCTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgtt < 1:702308/142‑1 (MQ=255) ccccgcATTTACCCACCCCCATGTGAAAAAAGGCTTG‑CGCGGCGGCTT‑TCTGCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct < 1:599236/146‑1 (MQ=255) cccacccccattacAAAAATCCCTT‑CGGGGGGGGTGTTCGCCCTGTTTTTTTC‑CCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg < 2:701030/138‑1 (MQ=255) cgcccATGTGCAAAAACGCCTC‑CGCGGCGTTTT‑TC‑GCCCGTTTTTTTTCCCGCCGCAACAAATCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa < 2:16107/148‑1 (MQ=255) cGCTATGTGAAAAAATGCTTT‑GGGGGGGGTTG‑CCTGCCAGTTTTTTTC‑CCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt < 2:805084/145‑1 (MQ=255) cGCCTTGTAAAAAAATCCCTT‑CGCGGCGGCTG‑TCGCCCATTTTTTTTC‑CCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa < 1:204171/150‑1 (MQ=255) aaaTGCTTT‑CGGGGCGGCTG‑TCTGCCAGTTTTTTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg < 1:460032/151‑1 (MQ=255) aaaCGCCGT‑CGGGGCGGCTT‑TGGGCAGGCTTTTTGC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg < 2:198804/151‑1 (MQ=255) aaTGCCTT‑CGCGGCGGCTG‑TCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAAAGAAAAACCTTTATTTCC‑TCCTGTTGTGTTGCGATCGCCGCTTGGATACATGTGGGAAGGGGAgggggg > 1:422480/1‑145 (MQ=255) cctcggcggcgCCTG‑TCTGCCATTTATTTGC‑CCCCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggcg < 2:296418/143‑1 (MQ=255) ccattttttttc‑CCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc < 1:563226/140‑1 (MQ=255) cggttgttttc‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa < 2:651676/133‑1 (MQ=255) ccatttttttggccgccaAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt < 2:30864/120‑1 (MQ=255) tttc‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt < 1:786282/150‑1 (MQ=255) ccgccgAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGcc < 2:670863/151‑1 (MQ=255) ccacaaaaaAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTc < 2:359457/147‑1 (MQ=255) caacaaAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAAAAACTTTAAAACCCTCCTGTTTTGTGGGGAACGCCGCTGGTAACAATATGGGAAGGGGAGGGGGGGGGATCAGCCCCCGAACCGCGCCAATTGAGTTCGCCTCCAtt > 1:222444/1‑150 (MQ=255) acaaACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAAt < 2:630808/148‑1 (MQ=255) caaAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGAGAAAAAAATTTTAACCCCGCCTGTTTTGTGGGGAACGCTGCTGGGGAAAAAATGGGAAGGGGAGGGGGGGGGAACAGCCCCACAAACGCGCCAATTTAGTTTGCCTccaatacc > 1:602814/1‑144 (MQ=255) | ATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTG‑TCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACC > NC_000913/357743‑358040 |
| 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: ‑ |