Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,890 | 0 | C | A | 52.4% | 17.6 / 35.3 | 21 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base C (10/0); new base A (0/11); total (10/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.84e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.82e-05 |
AAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC > NC_000913/357745‑358028 | aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAc > 1:442555/1‑151 (MQ=255) cTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGagag > 1:665246/1‑149 (MQ=255) tgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGCG‑aaaa > 1:159741/1‑149 (MQ=255) cAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGG‑aaaaa > 1:48590/1‑146 (MQ=255) aGCGGCTACCGTGATTCATTCCCCCCAACACCCGCGCATTTACCCAACGCAATGTGCAAAAATGCCTTCGGGGGGTCTGTCTGCCAGCTTTTTTT‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgag‑a < 1:486819/143‑1 (MQ=21) tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAG‑gaaaaaacttaaaac > 2:161111/1‑139 (MQ=25) ccGTGTTTCATTCCCCCCAACACCCGCGCTTTTCCCCACCCCAAGGTGAAAAAATCCTTTCGCGGCGGCTGTCTGCCAGTTTTTTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGcc < 1:127859/149‑1 (MQ=255) ccGTGATTCATTCCCCCCACCACCCGCGCTTTTACCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTT‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGccc < 1:495704/150‑1 (MQ=37) gtattccttcccccccAACACCCGCGCTTTTACCCAACCCCATTTGCAAAAA‑GCCCTCGGGGGGGGTTGCTGGCCGTTTTTTTTCCCGCCCCAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCg < 2:405949/135‑1 (MQ=255) gTGATTCATTCCCGACACCACCCGCCTATTTACCCACCCCCTTGGGCAAAAAGTCTTGCGGGGCGGCTGTCTGCCATTTTTTTTC‑CCACCAAATAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCttt < 2:170574/142‑1 (MQ=255) ttttCTTTCCCCCCACCACCCGCCCTTTTCCCCAACGCCTTGGGAAAAAATCCCTCCGCGGCGGCTGTCTCCGAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCa < 2:227747/148‑1 (MQ=255) gATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAa‑aaaa > 1:65965/1‑133 (MQ=255) ccacccccATGTGAAAAAATGCCTCGGCGGCGGCTGTCGCCAAGTTTTTTTC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt < 1:469587/145‑1 (MQ=255) tgCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAAATTTAAACCC‑CCCTTTTGTGTGGGAATCGCTGCTGGtggaaatgt > 1:192424/1‑132 (MQ=255) cTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAG‑AAAAAAATTTTAAACCGCCCTGTTGTGTTGGAAACGCTGCTGGGGAAAATAT‑GGGAAGGGGAGGGGGGGGGATCAGcccc > 1:120951/1‑142 (MQ=255) gCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAA‑AAAAAAATTTAAAACCCCCACTTTTTGTGGGAAATGCTGCCGGGGACAAAGTGGGAAAGGGAAGGGGGGGGaacaaaacccacaaac > 1:646819/1‑133 (MQ=255) cAGCTTTTTTC‑CCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc < 1:452029/143‑1 (MQ=255) ttttttttc‑CCGCCCAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTaa < 2:766793/145‑1 (MQ=255) tttttccccgccaCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg < 2:297687/143‑1 (MQ=255) cgccgcAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAG‑AAAAAAATTTTAAACCCCCCTGTTTTGTGGGAATAGCTGCCGGGGAAAAAGT‑GGGAAAGGGAGGGCGGGGGAACAAATCCCAAAACCCCCTAATTGAggt > 1:362344/1‑144 (MQ=255) ccgcAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc < 2:777677/145‑1 (MQ=255) | AAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC > NC_000913/357745‑358028 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |