Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,890 | 0 | C | A | 42.4% | 35.2 / 12.3 | 33 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base C (13/6); new base A (0/14); total (13/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.04e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCG‑CGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTC > NC_000913/357744‑358033 | tAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCG‑CGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACaaa > 2:114900/1‑150 (MQ=255) aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCG‑CGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAc > 2:254684/1‑151 (MQ=255) gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCG‑CGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAccg < 2:530718/151‑1 (MQ=255) cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCG‑CGCATTCATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCtgt > 1:218229/1‑150 (MQ=255) cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCG‑CGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg > 2:598763/1‑148 (MQ=255) gACCTGTTTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCG‑CGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg > 1:352057/1‑145 (MQ=255) tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCG‑CGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGA‑Gaa > 1:945995/1‑150 (MQ=255) gtgtCAGCGGCTCCCTTGTTTCTTTCCCCCCACCACCCG‑CGCATTTACCCAACCCTATTTG‑AAAAATGCCTTGCGGGGGGGCTGCCT‑CCCGCTTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAga‑ga < 1:27110/148‑1 (MQ=255) gtCAGCGCCTCCCGGAATTCTTTCCCCCCAACACCCG‑CGCATTTACCCACCCCCATGTGAAAAAATCCCTT‑CGCGGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑Gaaa < 1:930958/148‑1 (MQ=255) tCAGCGGCTACCGTGATTCATTTCCGCCAACAACCG‑CGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCCCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggagg > 2:855325/1‑138 (MQ=255) gCTACCGTGATTCATTTCCGCCAACAACCG‑CGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAAA‑GAAAAAACtttt > 1:1014623/1‑148 (MQ=25) ccGTGTTTCATTCCCCCCAACACCCG‑CGCTTTTACCCACCCCCATGTCAAAAATTCCCTT‑CGGGGCGGCTGCCGGCCAGTTTTTTTC‑CCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCg < 2:25903/151‑1 (MQ=255) tatttCATTCCCGCCAACACCCG‑CGAATTTATCCAACCCAATGTCAAAAAATGCCTT‑CGGGGCGGCTTTCTGCCAGTTTTTTCC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGt < 1:628959/146‑1 (MQ=255) tGATTCTTTCCCCCCAACACCCG‑CCCTTTTACCCACCGCCTTGTACAAAAATCCCTC‑CGCGGCGGCTGTCGCCCGGTTTTTTCC‑CCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCtt < 2:911987/140‑1 (MQ=255) ttCCCCCCACCACCGCCGCTTTCACCCACCCCCTTGTGAAAAAATGCCTC‑CGCGGCGGTTTTCTCCCAGTTATTTTC‑CCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCAtgt < 2:415092/147‑1 (MQ=255) caacacCCG‑CGCTTTTACCAACCGCAATGTAAAAAAATGCCTG‑CGCGGCGGCTGTCGGCCGGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGc < 2:528243/148‑1 (MQ=255) acaacCG‑CGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGGAAAAAGATTTATGCCC‑TCCTGTtgtgttggg > 2:199453/1‑143 (MQ=25) ccACCGCCAGGCCCAAAAATCCCTT‑CGCGGCGCCTCTCCGCCAGGTTTTTGC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt < 2:455919/151‑1 (MQ=255) ccAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCTCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg < 2:45837/150‑1 (MQ=255) cGCCTTGTAAAAAAATGCCTC‑GGGGGCGGCTGCCTGCCAGTTTTTTTC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa < 1:218110/150‑1 (MQ=255) aaaaaTGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑GAAAAAGATTTTAACCCGCCCCGTTTGGTGGGAAGCGCTGTTGGGGAAAAAGTGGGAAGGGGAgggg > 1:570680/1‑148 (MQ=255) aaaaaGGCCTC‑GGGGGCGGCTTTCGGCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg < 2:1160507/148‑1 (MQ=255) tGCCTT‑CGCGACGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑GAAAAAAATTTTAACCCCTCCTGTTGTTTGGGGATTGCCGCCGGTTACAAAATGGAAGGGGGAgggggggg > 2:541937/1‑143 (MQ=255) cggcggCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑GAAAAAACATTAAAACCGCCATGTTTTGTGGGAATTCCTGCTGGGAAAAAAGCGGAAAGGGGAggggggggggata > 2:287769/1‑135 (MQ=255) ggcggcTTTCTGCCACTTATTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTcc < 1:305022/151‑1 (MQ=255) gcgggttGTCTCCCGGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGc < 2:341877/142‑1 (MQ=255) tGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑AAAAAAGAGTTTATACCGCGAGGTTGTGTGGGAGTAGCTGATGGTAGAAAAATGATAAAGGGAGAGGGGCGGATAAGCTTCAGCACCccgc > 1:196576/1‑151 (MQ=255) ccAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa < 1:854815/151‑1 (MQ=255) ccAGTTATTTTC‑CCGCCCAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc < 1:371252/145‑1 (MQ=255) cAGTTATTTCC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa < 2:56973/150‑1 (MQ=255) ttttc‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa < 1:1076809/149‑1 (MQ=255) ccgccgCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCg < 1:597946/149‑1 (MQ=255) ccgcAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTc < 1:15456/150‑1 (MQ=255) | TAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCG‑CGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTC > NC_000913/357744‑358033 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |