| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 357,880 | 0 | T | C | 63.6% | 3.2 / 21.4 | 22 | K271E (AAA→GAA) | cynR | transcriptional activator of cyn operon; autorepressor |
| Reads supporting (aligned to +/- strand): ref base T (5/3); new base C (0/14); total (5/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.13e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.79e-05 | |||||||||||
GGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCT‑T‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGT > NC_000913/357737‑358026 | ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCT‑T‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑Ccgccgc < 2:112123/151‑1 (MQ=255) ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCT‑T‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGca > 2:634644/1‑150 (MQ=255) aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCT‑T‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAcc < 1:656472/151‑1 (MQ=255) tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCT‑T‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAccgc < 2:677479/149‑1 (MQ=255) cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCT‑T‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTtc > 1:383368/1‑151 (MQ=255) aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCT‑T‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCagt > 2:583657/1‑149 (MQ=255) gCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCT‑T‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg > 1:344387/1‑151 (MQ=255) ccTGTGTCAGGGGCAACCGTGATTCATTCCCGCCAACACCCGCGCTTTTACCCAACGCCATGTGAAAAAATGCCT‑T‑CGCGGCGGCTGTCTGCCAGCTATTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga < 1:478856/149‑1 (MQ=21) gtgtCGGCGGCTCCCGTGATTCATTCCCGCCAACACCCCCGCTTTTACCCACCGCCGTGTG‑AAAAATGGCTGT‑GGCGGCGGCTGTCGGCCAGTTATTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag < 2:16147/147‑1 (MQ=21) ggCTACCGTGATCAATCCCCGCCAACACCCGCCATTTTACCCACCCCCATGTGCAAAATTCCCT‑C‑CGGGGCGGCTTTCTGCAATTTATTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaa < 1:190278/141‑1 (MQ=255) cccGTTATTAATTCCCCCCAACACCCGCGCTTTTATCCACCGCCTTGTCCAAAATCGCCT‑T‑GGGGCCGCCTTCCGCCCAGTTTTTTTC‑CCCCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaaa < 1:231657/137‑1 (MQ=255) gttattcattcCCGCAACCACCCGCTCTTTTATCAACCGCATGGTGAAAAATGGCCT‑C‑CGCGGCGTGTGTCTGCCACTTATTTTC‑CCGCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGt < 1:109059/147‑1 (MQ=255) aTTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCT‑T‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGagaggagaaagtttaat > 1:66920/1‑125 (MQ=25) ccccccAACACCCGCTCTTTTATCCACCGCAAGGTAAAAAAAGGCCT‑C‑CGCGGCGGCTGTCGGCCGGTTTTTTTC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTc < 2:646771/137‑1 (MQ=255) acccGCGCATTTACCCAACGCCTGGTGAAAAAATGCCT‑C‑CGCGGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc < 1:631334/146‑1 (MQ=255) ccGCGCATTTACCCAACCCCTTGTGAAAAAATGCCT‑T‑CGCGGCGGTTGTCTGCCAGCTATTTTC‑CCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg < 2:259316/150‑1 (MQ=255) tttACCCACCCCCTTGTCAAAAAATGCCT‑T‑CGCGGCGGCTGTCGGCCAGCTTTTTTC‑CCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa < 1:246087/149‑1 (MQ=255) ccACCGCCATGTCAAAAAATCCCT‑C‑CGGGGCGGTTGTCTCCAGGATTTTTTC‑CCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg < 2:191904/150‑1 (MQ=255) aaaaaTGCTT‑C‑CGGGGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg < 1:187219/145‑1 (MQ=255) aaaaaGCCCT‑C‑CGCGGCGG‑TTTCCGCCCGTTTTTTTCCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt < 2:248624/149‑1 (MQ=255) aaaTGGCT‑TCCGGGGGGGTTGTCTGCCAGCTTTTTTC‑CCGCCGCAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcg < 1:280335/151‑1 (MQ=255) tgccggcagtcTATTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGa < 2:153173/138‑1 (MQ=255) cgcccattttttttc‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc < 1:634665/140‑1 (MQ=255) ttttc‑CCCCGGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTaa < 1:402145/145‑1 (MQ=255) tttc‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt < 1:619908/150‑1 (MQ=255) | GGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCT‑T‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGT > NC_000913/357737‑358026 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |