Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,890 | 0 | C | A | 42.9% | 39.2 / 23.4 | 42 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base C (17/7); new base A (0/18); total (17/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.26e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01 |
AAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC > NC_000913/357745‑358028 | atgtgtcgggcATCAGGCAATTCCAGCCGCGGACATGTGTCAGCGGCTACCGTGATTCATTTCCTCCAACAACCGCGTCGTCATCCAACCCCCTGTCCCATAATGCCCT‑CTCGGCGGCT‑GTCTTCCAGCTATTTTT‑CCGCCGCAACaaac < 2:97376/141‑2 (MQ=17) aaGCGTAGCGCATCAGGCAATTCCAGCCGAAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCACCAAAAc > 2:808749/1‑151 (MQ=255) tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAccgcc > 2:8948/1‑150 (MQ=255) aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTcc > 2:336118/1‑150 (MQ=255) ggaaaccccACCCGCAGCCCGGTGTCCGCGGCTCCGGTGATTAATTCCCGCCAGCACCCGCGTATTAACCCAGCGCCA‑GGGCAAAAAACCCCT‑CGGGGGGGGTGGTCTCCCGGCTTTTTTC‑CCCCCGCAACAAAACCGCCGTTCTCTCCa < 2:9119/143‑1 (MQ=255) ggCAATTCCAGCCGCAGCCTTGTGTCGGCGGCTACCGTGTTTCATTCCCGCCAACACCCGCGCTTTTACCCAACGCCATGTGAAAAATTGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCag < 1:927432/151‑1 (MQ=255) cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCagt > 1:422261/1‑150 (MQ=255) ttCCAGCCCAATCCCTTTGCCAGCGGCTCCGGTGTTTCTTTTCCGCCAACACCCCCGCTTTTTTCCAACGCCATGTACAAAATTCCTTG‑CGGGGGGGCT‑GTCTGCCAGTTTTTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCagtagt < 1:848376/150‑1 (MQ=255) tcccACCCGCAGCCTGGTGTCAGCGGCTACGGTGATTCTTTCCCCCCACCAACCGCGCATTTACCCACCGCCATGTGCAAAAATGCCTT‑CGGGGGGGCT‑GTCTGCCAGCTATTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCagtagt < 1:944520/148‑1 (MQ=18) tCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGGGCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTgg > 2:35156/1‑151 (MQ=255) aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg > 1:362075/1‑151 (MQ=255) ccGCAGCCCTGTGTGAGGGGCTCGTGGTTTTCTTTCCCGCAACCACCCGCGCTTTAACCCACCCCCTTGTGC‑AAAATCGCTT‑CGGGGGCGGT‑GGCCGGCAGTTTTTTTTGCCGCCGCAACAAACCCGCCGTTCTCTCCAGTAGTGGc < 2:270161/147‑1 (MQ=255) ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCT‑GTATGGCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCagtag > 1:377940/1‑143 (MQ=255) cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg > 2:670049/1‑148 (MQ=255) tgtttaggCGGCTCCCGTGATCCTTTCCCGCCAACACCCGCGTTTTTACCCACCCCCTTGTGAAAAAATCCTTT‑CGGGGCGGCT‑GTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga < 1:60502/140‑1 (MQ=255) gtaaGCGGCTCCTGTGTTTCATTCCCCCCAACACCCGGTTATTTAACAAACCCCATGTGAAAAATTGCTTC‑CGGGGCGGCT‑CTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAAt < 1:927840/146‑1 (MQ=255) ggCTCCCGTGTTTCTTTCCCCCAAACACCCGCGCTTTTACCCACCCCCGTGTCAAAAAATCCCTT‑CGCGGCGGCT‑GTCTGCCGGTTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑aaa < 1:1047150/142‑1 (MQ=255) ggCTACCGTGATTCTTTCCCCCAAACACCCGCCTTTTTACCAACCGCCATGTG‑AAAAATTCCTTCCGGGGCGGCT‑GTCGGCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTaa < 1:463497/151‑1 (MQ=255) cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑aaaaaaagt > 1:240837/1‑140 (MQ=255) tACCGTGATTCATTTCCGCCAGCAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGTGGACGAG‑AGAAAG‑CGTTTAgac > 2:205751/1‑149 (MQ=21) ccGTGATTCATTCCCCCCAACACCCGCGCATTTATCCACCCCCTTGTAAAAAAATGCCTC‑CGGGGCGGCT‑GTCTGCCAGTTATTTTC‑CCCCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTAAGc < 1:695044/148‑1 (MQ=255) tCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg > 2:808750/1‑113 (MQ=255) ccccaaacaacCGCCCTTTAACCCACCGCCATGT‑CAAAAATGGCTTCCGCGCCGGC‑‑GTCTGCCAGGTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTAAGCCCGTCATGTtgtg < 2:707340/143‑1 (MQ=255) ccGCCAACACCCGCGCTTTTACCCACCCCCATGTGAAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCg < 1:808796/112‑1 (MQ=255) ccGCCAACACCCGCCCTTTTACCAACCGCTATGTGAAAAAGTGCCTC‑CGCGGCGGCT‑GTCTGCCGGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTgg < 2:62156/151‑1 (MQ=255) acaacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAG‑GAAAAA‑CATTAAACCCCTcc > 2:230993/1‑134 (MQ=25) acccAACACCGTGTACAAAATTGCCTC‑CGCCGCGTCT‑GTTTGCCATCTTTTTTC‑CCGCCGCAAAAAAACCGCCTTTCTCTCCAGTAGTGGCGGGGCAAGAT‑AAATAG‑CTTTGAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc < 2:760515/146‑1 (MQ=255) cacccccATGTAAAAAAATCCCTC‑GGGGGGGGCT‑GTCTGCCATTTATTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgtg < 2:72720/146‑1 (MQ=255) aCGCCATTTACAAAAGTGCCTT‑CGCGGGGGGT‑GTCTGCGTGTTTTTTTC‑CCCCCACAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt < 2:1059963/148‑1 (MQ=255) cGCCATGTGCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAA‑CTTTAAACCC‑CCCTGTTTTGTGGGAAACGCtggtggggaaaa > 1:611827/1‑138 (MQ=25) cATGTGCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑AAAAAGAATTTAAACCC‑CCCTTTTGTGTGGGGATCGCTGCTGGGGACAAAGTgg > 1:353179/1‑146 (MQ=25) tgggaaaaatTCCCTC‑CGGGGCGGCT‑GTCGGCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg < 1:1020648/145‑1 (MQ=255) tgCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCCAAAGAAAAAAC‑CTTTAAACCC‑CCATGGTGTGTGGCAATCGCTCCTggggaaaaagtgggaa > 1:525807/1‑131 (MQ=255) gCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAAAG‑AAAAAG‑CATTAAAAAGGGCGTGCTGTTTGGCCATCGCCTCTGGGGACAATGGGGGAAgg > 1:18427/1‑146 (MQ=255) aaaaaTCCCTC‑CGCGGCGGCT‑TTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt < 1:336143/149‑1 (MQ=255) ccAGCTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg < 2:654314/141‑1 (MQ=255) aTTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAAAG‑AAAAAA‑AATTTATACCCCCATTTTTTGTGGGAATAGCTGGTGGGGACAATAGGGAAAGGGGAGTGGGGGGGATCAACCCCCGAAACccgcccatt > 1:929622/1‑143 (MQ=255) ttttc‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg < 2:1013718/150‑1 (MQ=255) ttttc‑CCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAtt < 1:180035/147‑1 (MQ=255) ttttc‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt < 1:912511/146‑1 (MQ=255) ttgccgcggCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt < 2:715880/147‑1 (MQ=255) ccgccgAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc < 1:543420/148‑1 (MQ=255) | AAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCT‑GTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC > 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 < 24 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |