Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,875 | 0 | A | T | 57.1% | ‑2.1 / 24.1 | 21 | N272K (AAT→AAA) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base A (5/4); new base T (0/12); total (5/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.19e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.15e-02 |
CTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCG‑CCAA‑CAA‑CCGCGCATTTATCCAAC‑GCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑C‑CGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACC > NC_000913/357726‑358013 | cTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCG‑CCAA‑CAA‑CCGAGCATTTATCCAAC‑GCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTAt > 1:200050/1‑151 (MQ=255) 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) | CTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCG‑CCAA‑CAA‑CCGCGCATTTATCCAAC‑GCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑C‑CGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACC > NC_000913/357726‑358013 |
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: ‑ |