Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,890 | 0 | C | A | 42.3% | 34.2 / 11.1 | 26 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base C (12/3); new base A (0/11); total (12/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.87e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCT‑GCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAA > NC_000913/357742‑358035 | gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca < 2:1029326/150‑1 (MQ=255) aTAAGCGTAGCGCATCAGGCAATTCCTGCCGCAGTCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa < 1:967669/150‑1 (MQ=255) aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa > 1:908025/1‑150 (MQ=255) aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc > 2:431597/1‑151 (MQ=255) aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc < 2:1029039/151‑1 (MQ=255) gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgcc > 1:956022/1‑151 (MQ=255) ctAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc < 1:793542/146‑1 (MQ=255) tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgccg > 1:785286/1‑151 (MQ=255) aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc > 2:790741/1‑150 (MQ=255) cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg > 1:134257/1‑151 (MQ=255) gtccTGTGTCCGCGGCTCCCGTGATTCTTTCCCCCAAACCACCGCTATTTTACCCAA‑CGCGATGGGCAAAAAAGCGTT‑GGGGGCGGTTGTCTGCAACTTATTTTCCGGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑Aag < 2:117488/148‑1 (MQ=255) ccTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg > 2:733225/1‑144 (MQ=255) tgtgTCAGCGGCTCCCGTGTTTCTTTCCCGCCAACACCCGCGCTTTTACCCAC‑CGCCTTGTGAAAAAATGCCTT‑CGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAgag‑a < 2:261310/149‑1 (MQ=255) tgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGGAGTGGCggggggag > 1:68381/1‑140 (MQ=255) tgacggCGGCTCCGGTGTTTCTTTCCCGCCAACACCCGCATTTTTCCCCAC‑CCCCATGTGAAAAATTGCTTC‑CGGGGCGGCTGTCTGCCGGTTTTTTTCCCCCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCgg < 1:9844/133‑1 (MQ=255) ggCTCCCGTATTTCATTCCCGCAAACACCCGCGCATTTACCCAC‑CGCCTTGTGAAAAAATGCCTC‑CGCGGCGGTTGTCTGCCGGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGAG‑AAATAGCtt < 2:999392/148‑1 (MQ=255) aTTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAAGAG‑AAAAAATTTAAAGccgtgactcttt > 1:788053/1‑141 (MQ=25) tttACCCAA‑CGCAATGAAAAAAAACCCCTC‑CGGGGCGTCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCT‑GCT‑GGt < 1:255967/148‑1 (MQ=255) ttttCCAACCGCTATGTAAAAAAACGCTTT‑CGGGAGGGCTTTTTGCAAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGct‑gct‑g < 2:839742/143‑1 (MQ=255) ccAC‑CGCCATGGAAAAAATTGCCTC‑CGGGGGGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCT‑GCT‑GGTAACa < 2:372900/147‑1 (MQ=255) gtaaaaaaaTCCCTC‑GGCGGCGCCTTTTTGCAATCTTTTTGCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCT‑GCT‑GGTAACAATGTGGAAAg < 1:751062/143‑1 (MQ=255) taaaaaaTGGCTTCCGGGGCGGCTGTCTGCAAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGCGTGGCAATCGCT‑GCT‑GGTAACAATGTGGAAAGGGaa < 2:664413/148‑1 (MQ=255) gCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGG‑AAGAG‑AAAAAAATTTTAGCCCCCCATCTTTTGTGGGGAACCCT‑GCTGGGGAACAATGGGGGAAggaagggggggg > 1:255091/1‑139 (MQ=255) ccTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGG‑AAGAG‑AAAAAGATTTAAAACCCTCATTCTTTGTGGCTATCGCTCGCT‑GTAAAGAATTGGGGAAGGGagggcgggcg > 1:719707/1‑140 (MQ=255) gggggcttccTGCCATTTTTTTTCCCCCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCT‑GCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCt < 2:24668/140‑1 (MQ=255) ggcTGTCTGCCATCTATTTTCCCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCT‑GCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa < 2:956036/149‑1 (MQ=255) aGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑GAGAGAAAAAAAATTTAAACCC‑TCCTGTTTTTTGGGAATAGCT‑GCT‑GGGGACAAAATGGAGAGGGGAGGGCGGCGAAAAAGCCCCCGAACCCCGCCaa > 2:26525/1‑150 (MQ=255) gtggtttgcccgcagAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGGTTGGCAATCGCT‑GCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt < 1:27151/141‑1 (MQ=255) tttCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑GAGGG‑AAAAAACCTTATGCCCCTCATGTTGTGTTGGAAACGCT‑GCT‑GGTAACAATGTGGGAAGGGGAGGGGGGCGGAACAGCTCCCCCACCGCGCCaaatttg > 2:86141/1‑143 (MQ=255) agaacaaAACCGCCGTCCTCTCCAGTAGTGGCGGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAAGCGCT‑GCT‑GGTAACAATGTGGAAAGTGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCGTCaa < 1:192266/148‑1 (MQ=255) aacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCT‑GCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCa < 2:27866/147‑1 (MQ=255) | GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCT‑GCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAA > NC_000913/357742‑358035 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |