| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 357,890 | 0 | C | A | 58.3% | 4.9 / 30.2 | 24 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
| Reads supporting (aligned to +/- strand): ref base C (6/4); new base A (0/14); total (6/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.56e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-05 | |||||||||||
AAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAA > NC_000913/357745‑358035 | aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAAc < 1:498715/151‑1 (MQ=255)aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAAc < 2:510061/151‑1 (MQ=255) gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAACCGCCGTtct > 1:261760/1‑151 (MQ=255) aTCAGCAAATTCCAGCCGCAGACCTGTGTCGGCGGCTACCGTGATTATTTCCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAAccgcc < 1:564272/144‑1 (MQ=18) ccgcagtcttggtccAGCGCCAACGGTGTTTCATTCCCCCCAAACACCGC‑TATTTAACCCACCGCTTTTTAAAAAATTCCTTCCGCGG‑CGGTTGCCT‑CCATTTTTTTTTCCCCCCCC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑gggg < 1:629474/137‑1 (MQ=255) ccGCAGACCTGTGTCGGGGCCTCCCGTGATTCTTTCCCGCCAACACCCGCCCATTTCCCCACCGCCATGTGCAAAAATGCCTT‑CGCGG‑CGGCTGTCTGCCAGTTTTTTTC‑CCGCCGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑gggg < 2:314030/151‑1 (MQ=18) cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGC‑ggtcg > 1:552138/1‑146 (MQ=255) tgtCGGCGGCTACCGTGTTTCTTTCCCGCCACCACCCGCGCATTTACCCACCCCCATTTGAAAAAATGCCTT‑CGCGC‑CGGCTGTCGGCCAGCTTTTTTC‑CCCCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCaa < 2:226565/143‑1 (MQ=255) cggtgtTTCATTCCCCCCAACACCCGCGCTTTTACCCACCGCCATGTGAAAAAATCCCTT‑CGCGG‑CGGTTGTCGGCCAGTTTTTTTC‑CCGCCGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAg < 2:475736/145‑1 (MQ=255) ttccccccAACACCCGCGCTTTTACCCAACGCAATGTGAAAAAATGCTTC‑CGCGG‑CGGCTGTCTGCCGGTTTTTTTC‑CCCCCGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttg < 2:102494/143‑1 (MQ=255) cccccaaacaccCGCGCTTTTACCCACCGCAATGTGAAAAAAGTCTTT‑GGCGG‑CGGCTGCCGGCAACCTTTTTTC‑CCGCCGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttg < 2:678071/141‑1 (MQ=255) cccccacaacccGCGCTTTTACCCACCGCCATTTGCAAAAAGCCCTC‑CGCGG‑CGGTTGTCT‑CCCGCTTTTTTTCCCGCCAA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgtg < 1:569719/140‑1 (MQ=255) acaccCGCGCTTTTAACAAACCCATTGTGAAAAAATGCCTT‑CCGGG‑GGGTTGTCTCCCA‑CTTTTTTT‑CCCCCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc < 2:525280/150‑1 (MQ=255) accggcgcATTCATCCAACACAATGTAAAAAACTGCCTC‑GGGGG‑GGGCTGTCGCCCTGTTTTTTTC‑CCGCCAC‑AAAAAAGCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct < 2:183159/147‑1 (MQ=255) aacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGGGAAAAAAATTTTAACCCCGCATGGTTTGTggggatc > 1:635511/1‑143 (MQ=21) aTCCAACGCCATGTGCAAAAATGCCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGCAGAGGAAAAAATTTTAACCCCCCCTGTTGTGTGGGAATCGCCGCTGGtgaaaa > 2:639140/1‑145 (MQ=25) catgatacaaaaTGTCCT‑CGCGGCCGGCTGTCTGCCAGTTATTTTC‑CCGCCAC‑AAAAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGaa < 1:257826/141‑1 (MQ=255) aaaaaaTCCCTC‑GGGGG‑CGGTTTTCGGCCAGTTTTTTTC‑CCGCCCA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGt < 1:164377/149‑1 (MQ=255) cgcgG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAAGGAAAAAAATTTTAAACCCCCCA‑GTTTTGTGGGGAATGCTGCTTGGGACAAAGTGGGAAAGGGaggggggggggaca > 2:361918/1‑136 (MQ=255) cgcccAGTTTTTTTC‑CCCCCGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc < 1:586694/144‑1 (MQ=255) cTGCCAGTTTTTTTC‑CCGCCAA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcg < 1:386483/149‑1 (MQ=255) ccAGTTATTTTT‑CCGCCGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc < 2:10971/147‑1 (MQ=255) ccAGCTATTTTT‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGAGAAAAAAATTTTAACCCCTCCTGTTTTGTGGGGAACGCCGCTGGGGACAAAGT‑GGGAAGGGGAGGGGGGCGGAACAGCTCCAGAAcccc > 2:94707/1‑145 (MQ=255) ttttttc‑CCGCCGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAtt < 2:265568/147‑1 (MQ=255) aacaaAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCaa < 1:507099/148‑1 (MQ=255) | AAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGG‑CGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAA > NC_000913/357745‑358035 |
| 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: ‑ |