Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,890 | 0 | C | A | 54.5% | 3.6 / 25.7 | 22 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base C (8/2); new base A (0/12); total (8/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.41e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.58e-02 |
GTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC > NC_000913/357749‑358028 | gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑ccgc > 2:749669/1‑150 (MQ=255) gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGtt > 2:106096/1‑149 (MQ=255) atttcccccggCAGCCCTGTGCCGGCGGCTCCCGTTT‑TTC‑TTTCCCCCAACCAC‑CCGCGATTTTACCCACCCCCTTGTACAAAATTGCTTC‑CGGGGCGGTTGTCTGCCAGTTTTTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCagtag < 2:153865/141‑1 (MQ=255) attcccccccccAGACCTGTGTCGGCGGCTACGGTGA‑TTC‑TTTCCCCCAACCAC‑CCGCGCTTTTACCCACCGCAGTGTGAAAAAATGCTTT‑CGGGGCGGCTTTCTGCCAGTTTTTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCagtag < 2:746321/140‑1 (MQ=255) attcccAGCCGCAGCCTTGTGTCGGGGCCCCCCGTGT‑TTCTTTCCCGCCAA‑CAC‑CCGCACTTTTCCCCAACCCCATGTGAAAAAGTCCCTC‑CGCGGCGGCTGTCTCCCAGCTATTTTCCCCCCGCAAAAAAA‑CCGCCGTTCTCTCCagtag < 2:613064/147‑1 (MQ=255) attccccgccccgACCTGTGGCGGCGGCTTCCGGTA‑TTTATTTCCCCCAA‑AAACCCCCGCTTTTCCCCACCG‑CAAGTTGAAAAATGCCTTCCGGGGCGCCTTTCGCCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCagtag < 1:591658/139‑1 (MQ=255) cccgCCGCAGCCTTGTTCCAGCGCACACCGGGT‑TTCATTCCCCCCAA‑CAC‑CCGCCCTTTACCCCAACGCCATGTGAAAAATGGCTTT‑CGGGGCGGCTGTCTGCCATTTTTTTTCCCGCCGAAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGc < 1:596116/148‑1 (MQ=255) gACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgggggcagg > 2:229730/1‑146 (MQ=255) cccTGTGCCAGGGGCTACCGGGT‑TTCTTTCCCCCAAA‑CAA‑CCGCCCATTTACCCACCCCCGTGGGAAAAAAGGCCTT‑CGGGGGGGCTGTCGGCCAGCTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag < 2:244679/150‑1 (MQ=255) ccTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGccaaaaa > 2:430780/1‑145 (MQ=255) cGGCTACCGTGA‑TTCATTTCCGCCAA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCaacgagcaaaaacttt > 1:675878/1‑137 (MQ=255) ttccGTGATTTCTTTCCCCCCAA‑CAC‑CCGCGCTTTTACCAAACGCTATGTGCAAAATTGCCTC‑CG‑GGCCGTTTTCCGCCGGGTTTTTTTCCCCCGGAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGc < 1:141499/149‑1 (MQ=255) ccGTGA‑TTCATTTCCGCCAA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAaaaaaaaattttaaacc > 2:217178/1‑137 (MQ=25) gA‑TTCATTTCCGCCAA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGGAAAAGATTTTAGccccccc > 2:35545/1‑146 (MQ=25) taTTCCCGCCAA‑CAC‑CGGCGATTTAACCCACCCCCATGTGCAAAATTGCCTT‑CGGGGGGGCGGTCTCCAAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttg < 1:747483/150‑1 (MQ=255) ccacccccttgtcAAAAATTCCTTCCGGGGCGGTTGTCTGCCACTTTTTTTCCCCCCGCAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt < 2:744250/139‑1 (MQ=255) ttGTCTCCCATCT‑TTTTTCCCCCCCAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc < 2:388495/150‑1 (MQ=255) ctttctcccattttttttcCCGCCGAAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc < 2:62492/137‑1 (MQ=255) cTGTCTGCCAGCTTTTTTTCCGCCGCAAAAAAA‑CCGCCGTCCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc < 2:363870/151‑1 (MQ=255) ttgtttccccgcccCAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg < 2:71594/143‑1 (MQ=255) cTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAATTTTAACCCCCCCTTTTTTTTGGGGCTCGCTGCTGGGGACAAAAAGGGGAGGGAAGGGGGGGGGATCACATCCCGAACCccgccacttt > 2:116652/1‑145 (MQ=255) ttcccgccgAAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc < 2:651123/148‑1 (MQ=255) | GTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC > NC_000913/357749‑358028 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |