| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 357,890 | 0 | C | A | 66.7% | 3.8 / 19.4 | 24 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
| Reads supporting (aligned to +/- strand): ref base C (4/4); new base A (0/16); total (4/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.59e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.55e-02 | |||||||||||
AGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCT‑T‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAA > NC_000913/357746‑358035 | 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) agtcttttgtcccgcccAAACAACTCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt < 1:176614/138‑1 (MQ=255) ttttc‑CCCCGGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTaa < 1:402145/145‑1 (MQ=255) tttcccccagacaacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg < 2:81619/139‑1 (MQ=255) tttc‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt < 1:619908/150‑1 (MQ=255) aacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCaa < 1:534726/148‑1 (MQ=255) | AGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCT‑T‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAA > NC_000913/357746‑358035 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |