Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,880 | 0 | T | C | 57.9% | ‑0.8 / 19.6 | 20 | K271E (AAA→GAA) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base T (4/4); new base C (0/11); total (4/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.81e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.79e-02 |
GTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCG‑CCAA‑CAA‑CCGCGCATTTATCCAAC‑GCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑C‑CGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTG > NC_000913/357734‑358023 | gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCG‑CCAA‑CAA‑CCGCGCATTTATCCAAC‑GCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑c‑cgc < 2:591778/151‑1 (MQ=255) gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCG‑CCAA‑CAA‑CCGCGCATTTATCCAAC‑GCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑C‑CGCCGCAACAAAAccgcc < 2:29102/151‑1 (MQ=255) gCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCG‑CCAA‑CAA‑CCGCGCATTTATCCAAC‑GCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑C‑CGCCGCAACAAAACCGCCGTTCTCTCCagt > 2:547480/1‑151 (MQ=255) tcccACCCGCAGCCCTGTGTCGGCGCCTCCCGTGTTTC‑TTT‑CCC‑CCAA‑CAACCCGCGCTTTTACCCAAC‑GCCATGTCAAAAAATCCTTT‑CGCGGCGGCTTTCTCCCAGCTTTTTTC‑C‑CGCCGCAAAAAAACCGCCGTTCTCTCCagta < 1:486966/146‑1 (MQ=255) ccAGCCGCAGACTTGTGTAGGCGGCTCCCGTGATTCTTTC‑CCC‑CAAA‑CAC‑CCGCGATTTTACCCACC‑GCCATGTGAAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGTTTTTTTC‑C‑CGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTg < 1:601154/149‑1 (MQ=255) ccGCAGACTTGTTCCAGCGCCTACCGTGTTTCTTTC‑CCG‑CAAC‑CAC‑CCCCGAATTTACCCACC‑CCCTTGGGAAAAAATGCCTT‑CGCGGGGGCTGTCTGCCAGTTTTTTTC‑C‑CGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgg < 1:472741/149‑1 (MQ=255) cGCAGCCCTGTGTCAGCGGCTCCCGTGATTCTTTC‑CCC‑CCAA‑CAC‑CCGCCCTTTTACCCACC‑GCCATGTGAAAAAATCCCTC‑CGGGGCGGCTGTCGGCCGCCTATTTTC‑C‑CGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGc < 1:453450/151‑1 (MQ=255) gggtCGGCGGCACCCGT‑ATTTATTTCCCC‑CCAA‑CAC‑CCGCGCTTTAACCCAAC‑GCAATGT‑CAAAAATGGCTTCCGGGGCGGCTGTCTGCCAGTTTTTTTT‑C‑CGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaaa < 1:547525/148‑1 (MQ=255) tCAGCGCCTCCGGTGTTTC‑TTT‑CCC‑CCAACCAC‑CCGCCCATTTCCCCCAC‑CCTTTTTGAAAAACTCCCTG‑CGGGGGGGTTCTCCGCCAGTTTTTTTC‑C‑CGCCAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCaa < 2:315094/141‑1 (MQ=255) gcgcCTCCCGTAATTCTTTC‑CCC‑CCAA‑CAC‑CCCCGCTTTTACCCACC‑GCCTTGTAAAAAAATGCCTC‑CGCGGCGGCTGTCGCCCAGTTTTTTTC‑C‑CGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaa < 2:816482/143‑1 (MQ=255) gCTACCGTGATTCATTT‑CCG‑CCAA‑CAA‑CCGCGCATTTATCCAAC‑GCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑C‑CGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGgaaaaaatttta > 1:803257/1‑141 (MQ=25) cTACCGTGATTCATTT‑CCG‑CCAA‑CAA‑CCGCGCATTTATCCAAC‑GCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑C‑CGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGagggcgagaa > 2:222198/1‑135 (MQ=255) caCCGTTATTTGTTT‑CCCCCCAA‑CAC‑CCGGGCATTTCT‑CAAC‑GCCAAGGACAAAAATGCCTT‑GGGGGGGGGCGCTAGCCCGGTTGTTGACC‑CGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGc < 2:546411/150‑1 (MQ=255) cccccAA‑CAC‑CCGCGC‑TTTCTCCCACGGCAATGTGAAAAATTCCCTC‑CGGGGCGTCTTTC‑CCCTGATTTTTTT‑CACGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgtg < 2:530433/147‑1 (MQ=255) caaccgctgtgtgAAAAAACGCTTT‑CGCGGCGTCTGTCTGCCATTTTTTTTC‑C‑CCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg < 1:376743/141‑1 (MQ=255) aaTGTGAAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGTTTTTTTC‑C‑CGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAgg < 1:152752/150‑1 (MQ=255) aaaaaggcttgggcggcgGTTGTCGCCCAGCTTTTTTC‑C‑CGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTg < 1:50201/138‑1 (MQ=255) aaaaTTCCTTCCGGGGCGGTTTTCTCCCAGTTTTTTTC‑C‑CCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt < 1:302722/149‑1 (MQ=255) tGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑C‑CGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGGAAAAAACATTTAACCCCCCATTTTTTTTGGGGAACGCCGCTGGGGAAAAAGTGGGAAGGGGAgggggggg > 2:54284/1‑143 (MQ=255) tCTGCAAGTTATTTTC‑C‑CGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc < 2:803190/149‑1 (MQ=255) ttttttc‑C‑CCCCCCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg < 2:583093/148‑1 (MQ=255) | GTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCG‑CCAA‑CAA‑CCGCGCATTTATCCAAC‑GCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑C‑CGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTG > NC_000913/357734‑358023 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |