| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 357,890 | 0 | C | A | 42.9% | 16.6 / 14.2 | 28 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
| Reads supporting (aligned to +/- strand): ref base C (10/6); new base A (0/12); total (10/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 8.83e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
ATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGC‑CTTCGCGGCGGCTGTCTGCCAG‑CTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATG > NC_000913/357743‑358037 | aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGC‑CTTCGCGGCGGCTGTCTGCCAG‑CTATTTTT‑CCGCCGCAACaaa < 1:662916/151‑1 (MQ=255) aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGC‑CTTCGCGGCGGCTGTCTGCCAG‑CTATTTTT‑CCGCCGCAACAAAACCGCCGTtctct > 1:496169/1‑151 (MQ=255) cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGC‑CTTCGCGGCGGCTGTCTGCCAG‑CTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTcc > 1:638287/1‑151 (MQ=255) gCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGC‑CTTCGCGGCGGCTGTCTGCCAG‑CTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCagt > 2:550809/1‑151 (MQ=255) ccGCAGACCTGTATCGGGGCTTCCCGTTATTCTTTCCCCCAAACAACCGCCCATTTACCCAACCCC‑AGGTGAAAAAGTCC‑TTTCGCGGGGCTTTTCTGCCAC‑CTTTTTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgggg < 1:693709/151‑1 (MQ=255) tgtCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCC‑ATGTGCAAAAATGC‑CTTCGCAGCGGCTGTCTGCCAG‑CTATTTTT‑TATGCCCAACACAACCGCCGCTCCCTCCAGGAGcggcggagggcacggaa > 2:9797/1‑138 (MQ=1) gtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGC‑CTTCGCGGCGGCTGTCTGCCAG‑CTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggcgaggaaaa > 1:362500/1‑139 (MQ=255) gccggCGTCTACGGTGATTCATTCCCCCAAACACCCGCCATTTTACCCACCCCC‑TTTTAAAAAAACCC‑CTCCGGGGCGCCTGTCTCCAAG‑TTTTTTTC‑CCCCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaaa < 1:759402/144‑1 (MQ=255) ggCTCCCGTGTTTCATTCCCGCCAACACCCGCGATTTTAACAAACGCC‑TTGTCAAAAAAGCC‑CTTCGCGGCGGCTGCCTCCCAG‑CTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTaa < 1:558131/151‑1 (MQ=255) ggCTACCGTGATTCTTTCCCGCCAACACCCGCGCTTTTATCCAACGCC‑TTGTGAAAAAATGC‑CTCCGGGGCGGCTTTCTGCCAG‑TTTTTTTT‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑tt < 1:599320/148‑1 (MQ=255) attcattcCCGCCAACACCCGCGCATTACCCCACCCCC‑ATGTACAAAAATGC‑C‑TCGCCGCCGTCGGCCGCCCGTTTATTTTC‑CCGCCAAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCa < 2:623354/149‑1 (MQ=255) aTTCCCCCCAACACCCGCGCTTTTATCCACCGCC‑ATGTGCAAAAAGGC‑CTCCGCGGCGGCTTTCTGCCAG‑TTTTTTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCAtgt < 1:840357/148‑1 (MQ=255) ccccccAACACCCGCGCTTTTCCCCACCCCC‑ATGTGAAAAAATCC‑TTCCGGGGGGGCTGTCTGCCAG‑TTTTTTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCa < 1:9823/138‑1 (MQ=255) caaacaccCGCG‑ATTTTACCAACCCCAATGTAAAAAAAGCC‑CTCCGGGGGGGCTTTCGGCAAT‑TTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCa < 1:605823/148‑1 (MQ=255) cccGCCCTTTTATCCACCCCC‑ATTTCAAAAAATGC‑CTCCGGGGGGGCTTTCGGCCAT‑TTTTTTTC‑CCCCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCgctg < 1:535652/150‑1 (MQ=255) ccGCGCTTTTACCCACCCCT‑ATGTAAAAAAAGCC‑CTTCGCGGGGGTTTTCTGCAAT‑TTTTTTTC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCgct < 1:211554/149‑1 (MQ=255) ccGCGCATTTATCCAACGCC‑ATGTGCAAAAATGC‑CTTCGCGGCGGCTGTCTGCCAG‑CTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAAAAC‑TTTTAAACCCTCCTGTTTTGTgggga > 1:60655/1‑141 (MQ=37) gcgcATTTATCCAACGCC‑ATGTGCAAAAATGC‑CTTCGCGGCGGCTGTCTGCCAG‑CTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCCAGGGAAAAAACATTTAAACCCG‑CATTTTTTTTGGGGATCgctgct > 2:59478/1‑150 (MQ=25) cgcATTTATCCAACGCC‑ATGTGCAAAAATGC‑CTTCGCGGCGGCTGTCTGCCAG‑CTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAA‑TTTAAAACCCCCCTGTTGTGTGGGTAACGCCGCTgg > 1:656200/1‑151 (MQ=25) cc‑ATGTGCAAAAATGC‑CTTCGCGGCGGCTGTCTGCCAG‑CTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGGAAAAAAG‑TTTTATACCCTCCCGTTGTGTGGCGATCGCTGCTggggaaaaaaggg > 1:452211/1‑136 (MQ=25) aaaaaTCC‑CTGCGCGGCGGCTGTCTGCCAG‑TTTTTTTC‑CCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt < 1:159661/149‑1 (MQ=255) aaaaTGCTCTTCGGGGGGGGTGTC‑GGCAG‑TTTTTTTTCCCGCGAAAACAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa < 2:635028/146‑1 (MQ=255) c‑cTCCGCGGCGGTTGTCCCCCAG‑TTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGa < 2:452157/149‑1 (MQ=255) ccAG‑TTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgcg < 2:621550/148‑1 (MQ=255) ccAG‑TTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa < 2:258805/140‑1 (MQ=255) tttttccccgccgAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAtt < 2:732586/142‑1 (MQ=255) ttt‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAAAAA‑TTTAAAACCCTCATGTTGTGGGGGGAACGCCGGTGGTGACAAAGTGGGAAGGGAAAAGGGGGGGATCAGCCCCACCACCGCGCTCATTgtg > 1:648958/1‑146 (MQ=255) cgccgcAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCt < 2:738565/151‑1 (MQ=255) aacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATg < 2:454680/150‑1 (MQ=255) | ATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGC‑CTTCGCGGCGGCTGTCTGCCAG‑CTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATG > NC_000913/357743‑358037 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |