| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 357,890 | 0 | C | A | 48.3% | 17.1 / 18.2 | 29 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
| Reads supporting (aligned to +/- strand): ref base C (12/3); new base A (0/14); total (12/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.05e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTA > NC_000913/357747‑358019 | gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAACAAAAccg > 2:8227/1‑151 (MQ=255)gCGTAGCGCATCAGCCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCCCGCATTTATCCAACGCCATGCGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCCCCGCAACAAAAcc > 1:233872/1‑150 (MQ=255) gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAACAAAAccgcc > 1:780471/1‑151 (MQ=255) ggCATTTCCAGCCGCAGCCCTGTTTCGGGGGCTCCCGTGTTTCTTTCCCCCCAACACCCGCCCTTTTACCCACCCCCTTGTGAAAAAATG‑CCTT‑GGCGGCGGTTTTCTCCC‑ATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCag < 2:648146/151‑1 (MQ=255) ccGCAGCCTTGTGTCAGCGGCTCCGGTGTTTCTTTCCCCCCAACACCCGCGCATTTACCCACCGCCATGTGAAAAAATG‑CCTT‑CGGGGGGGCTGTCTCCC‑AGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCagtagt < 2:145602/144‑1 (MQ=255) ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg > 2:404092/1‑151 (MQ=255) tctcccgatctcggCGGCTCCGGTGATTATTTCCCGCAAACACCCCCTCATTTATCCACCGCCAGGTGAAAAAATG‑CCTT‑CGCGGCGGCTGTCGGCC‑AGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTg < 2:291020/129‑1 (MQ=255) cTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCagg > 2:722709/1‑145 (MQ=255) tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg > 2:406207/1‑142 (MQ=255) tgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAGGAGaaa > 2:582323/1‑149 (MQ=255) cAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTg > 1:291000/1‑131 (MQ=255) aGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCC‑AGCTA‑TTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAACAGaaaaaaactt > 2:920843/1‑144 (MQ=255) ggCTCCGGTGTTTCTTTCCCCCCAACACCCGCGCATTTACCCACCGCAATGTGAAAAAATC‑CCTT‑CGCGGCGGCTGTCGGCC‑AGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGc < 2:383537/146‑1 (MQ=255) ctttgtttCTTTCCCGCCAACACCCGCGCTTTTACCCACCCCCAGGTGAAAAAATC‑CCTC‑GGGGGGGGCTGTCTCCC‑ATCTTTTTTCCCGCCACAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCtt < 1:82146/137‑1 (MQ=255) cGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGGGAAAAAGGTTTAAAccc > 2:958366/1‑149 (MQ=255) gTGATTCTTTCCCCCCAAACCCCGCGCTTTTACCCACCCCAATGTGAAAAAATC‑CCTC‑CGGGGCGCCTTTCTCCC‑AGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGt < 2:872343/150‑1 (MQ=255) ttCTTTCCCGCCAACACCCGCGCTTTTACCCACCGCCTTGT‑CAAAAATT‑CCTTCCGCGGCGGTTGTCTGCC‑AGTTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGt < 2:116817/146‑1 (MQ=255) ttCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGaaaaaaactttaaacc > 2:606923/1‑131 (MQ=25) accacccGCGCATTTACCCACCCCCATGTGAAAAAATG‑CCTC‑CGCGGCGGCTGTCTGCC‑AGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAAt < 2:809971/145‑1 (MQ=255) acCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAAAAATAAATTTTAAACCGCACTTTTGTGTGGCGATCgc > 1:95770/1‑149 (MQ=25) cccAACGCCATGAAAAAAAATG‑CCTC‑CGGGGCGTCTGTCGCCC‑AGTTTTCTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt < 1:832176/149‑1 (MQ=255) cccaccccaatgtcAAAAATGCCCTC‑CGGGGCGGCT‑TCCGCCAAGTTTTTTTCCCCCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACaa < 1:375821/136‑1 (MQ=255) ccacccccttttcaaaaaaTC‑CTTC‑CGCGGCGTCTTTCTCCC‑CGTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc < 1:410443/132‑1 (MQ=255) cacccccttgtaaaaaaaTC‑CCTC‑CGGGGCGGTTGTGTCCC‑GGTTTTTTTCCCGCAGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg < 1:603782/136‑1 (MQ=255) cccctggaaaaaaaaTC‑CTTC‑CGGGGCGGCTGTCCGCC‑AGTTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgtg < 1:717219/138‑1 (MQ=255) cctggtgAAAAATTC‑CCTC‑CGCGGCGGTTGTCCCCC‑ACTTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaa < 1:158096/145‑1 (MQ=255) tGTCTCCC‑AGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATc < 2:79190/139‑1 (MQ=255) cc‑AGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa < 1:688169/151‑1 (MQ=255) ttttttgCCGCCGCAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa < 1:790066/144‑1 (MQ=255) | GCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTA > NC_000913/357747‑358019 |
| 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: ‑ |