Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,890 | 0 | C | A | 54.2% | 5.8 / 33.1 | 24 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base C (8/3); new base A (0/13); total (8/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.24e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.12e-03 |
CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAA‑GAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCT > NC_000913/357740‑358032 | cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTCTTCCGCCGcaac > 1:145987/1‑151 (MQ=255) aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc < 2:407876/151‑1 (MQ=255) aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc < 1:59416/151‑1 (MQ=255) tCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctctc > 1:137985/1‑151 (MQ=255) cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc > 1:86572/1‑151 (MQ=255) aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc > 1:398592/1‑150 (MQ=255) tCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTgg > 2:233793/1‑151 (MQ=255) cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg > 1:76401/1‑150 (MQ=255) tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg > 2:538552/1‑140 (MQ=255) gCGGCTCCGGTGTTTCTTTCCCCCCAACAC‑CCCCGCTTTTACCCACCCCATTGTG‑AAAAAATCCTTCCGCGGCGGTTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAA‑GAGAAATAGCtt < 1:371240/150‑1 (MQ=255) gcctccggtgttccATTCCCG‑CAACAACCCGCGCTTTTACCCAAC‑CCATGTGAAAAAAAGGCTTTCGGGGGGGCTGTCTGCCAGATATTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAA‑GAGAAATAGc < 2:811163/133‑1 (MQ=255) gCTACCGTTATTCATTCCCCCCAACAC‑CCGCGCTTTTACCCAACCCCATGTG‑AAAAAATGCTTTCGCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAA‑GAGAAATa < 1:314871/143‑1 (MQ=255) ttttttccccccAACAC‑CCGCGCTTTTATCCACCGCCATGTA‑AAAAAATGCTTTCGCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAA‑GAGAAATAGCTTTAAGCCCGTCAt < 1:318057/139‑1 (MQ=255) aacac‑cCGCCCTTTTACCCACCCCCATGTA‑AAAAAATCCCTTCGGGGCGGCTTTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAA‑GAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAAt < 1:660397/150‑1 (MQ=255) ccGCGCTTTTACCCACCGCCTTGTA‑AAAAAATCCTTCCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAA‑GAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct < 1:825190/149‑1 (MQ=255) tttATCCAACGCCATGTG‑CAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGA‑GAGAAAAAAATTTTAAACCCCCCTGTTTTGTGGGAATCGcttctg > 1:536549/1‑146 (MQ=25) ccAACGCCATGTG‑CAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGGAAAAAAAATTTAAAAACCCCCTTGTT‑TGTGGGAATCGCTGCTGGGAACaa > 1:576555/1‑148 (MQ=25) gtaaaaaaaTCCCTCCGGGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAA‑GAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa < 1:352799/146‑1 (MQ=255) aaaaaTGCCTCCGCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAA‑GAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTg < 1:292437/150‑1 (MQ=255) aaaTGCCTCCGCGGGGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAA‑GAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg < 1:53797/151‑1 (MQ=255) ccTCCGGGGGGGCTGTCGGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAA‑GAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt < 2:783667/151‑1 (MQ=255) ggcggtgttCTGCCGGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAA‑GAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTc < 1:554944/142‑1 (MQ=255) ccgccgCACAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAA‑GAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGc < 1:43915/150‑1 (MQ=255) cgccacAAAAAACCCTCCGTTCTCTCCAGTAGTGGCGGGGCAA‑GAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCt < 2:482782/151‑1 (MQ=255) ccgcAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAA‑GAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc < 2:250548/145‑1 (MQ=255) | CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAA‑GAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCT > NC_000913/357740‑358032 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |