Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,880 | 0 | T | C | 48.6% | 36.4 / 12.8 | 35 | K271E (AAA→GAA) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base T (12/6); new base C (0/17); total (12/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.97e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGT‑C‑TGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATT > NC_000913/357731‑358022 | cTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGT‑C‑TGCCAGCTATTTTT‑c > 1:260472/1‑151 (MQ=255) aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGT‑C‑TGCCAGCTATTTTT‑ccgccg < 2:915765/151‑1 (MQ=255) aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGT‑C‑TGCCAGCTATTTTT‑ccgccg > 1:500000/1‑151 (MQ=255) gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGT‑C‑TGCCAGCTATTTTT‑Ccgccgc > 1:414511/1‑150 (MQ=255) ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGT‑C‑TGCCAGCTATTTTT‑CCGCCGca > 2:865499/1‑150 (MQ=255) tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGAGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGT‑C‑TGCCAGCTATTTTT‑CCGCCGcaa > 1:886484/1‑149 (MQ=255) gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGT‑C‑TGCCAGCTATTTTT‑CCGCCGcaaca < 2:315062/150‑1 (MQ=255) aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGT‑C‑TGCCAGCTATTTTT‑CCGCCGCAACaaa < 2:741167/151‑1 (MQ=255) aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGT‑C‑TGCCAGCTATTTTT‑CCGCCGCAACaaaa < 1:866002/150‑1 (MQ=255) aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCCCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGT‑C‑TGCCAGCTATTTTT‑CCGCCGCAACAAAAc < 1:707111/151‑1 (MQ=255) aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGT‑C‑TGCCAGCTATTTTT‑CCGCCGCAACAAAAcc > 1:528184/1‑151 (MQ=255) aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGT‑C‑TGCCAGCTATTTTT‑CCGCCGCAACAAAAcc > 2:880994/1‑151 (MQ=255) tcccACCCGCACACTTGTGTCGGGGGCTCCCGTTTTTCTTCCCCCCC‑ACCACCCGCGATTTTACCCAACGCCATGTCAAAAAATCCCTCGGGGGGGGTTGT‑C‑TGCCACTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCagt < 1:784328/145‑1 (MQ=255) ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGT‑C‑TGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg > 1:132932/1‑151 (MQ=255) ccGAAGCCCTGTTTCAGCGGCTCCCGTGATTCATTCCCCCC‑AACACCCGCGCATTTATCCACCCCCATGTCAAAAAATGCCTTCGGGGCGGTTGT‑C‑TGCCAGCTATTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCggg < 2:362566/150‑1 (MQ=255) aGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGT‑C‑TGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg > 1:254731/1‑147 (MQ=255) ggccggcgcCTACCGTGATTCATTCCCGCC‑ACAACCCGGTCATTTACCCACCGCCATGTGAAAAAATCCCTTCGGGGC‑GCTGTCC‑TGCC‑GCTATTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt < 2:6418/145‑1 (MQ=255) tgccaggggctccGTTATTTCTTTCCCCCAAACACCCGCTCATTTACCCAACCCCAGAGAAAAAAAGGCTTTCGCGACGGCTGT‑C‑TCCAAGCTATTTTC‑CCCCGGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaa < 2:802039/137‑1 (MQ=255) gCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGT‑C‑TGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggagagacagaaaaa > 1:835867/1‑135 (MQ=255) ggttccgtatttCCTTTCCCCC‑CACAACCCCCCTTTTTACCCACCCCCAGGGAAAAAAAGCCTTCCGGGCGGGTTT‑CTTGCGGGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgaga < 1:931508/130‑1 (MQ=255) ccGTGTTTCATTCCCCCC‑AAAACCCGCGAATTTCTCCAACGCCATGTCAAAAAATGCTTTCGCGGCGGTTGT‑C‑TCCCAGTTTTTTTC‑CCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTaa < 2:165768/146‑1 (MQ=255) ccGTGATTCATTCCCCCC‑AACACCCGCCCATTTACCCAACGCCGTGTGAAAAAATGCCTCCGCGGCGGCTGT‑C‑TGCCAGCTATTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCg < 2:632854/151‑1 (MQ=255) cGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGT‑C‑TGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGaaaaaacatttaaccc > 2:1017863/1‑136 (MQ=25) tttCCCCC‑AACCCCCGCGTTTTTACCCAACGCCATGTAAAAAAAGCCCTCCGCGGCGGCTGT‑C‑TCCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttgt < 2:840784/150‑1 (MQ=255) ccgcttaTTTATCCACCCCCTTGCGAAAAATTGCCGTCGCGGGGGTTGT‑C‑CGCCGTTTTTTTTC‑CCGCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct < 1:698152/143‑1 (MQ=255) tttACCCACCCCCTTGTAAAAAAATCCTTTCGGGGGGGCTGT‑C‑TGCCAGCTTTTTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTgg < 2:604339/147‑1 (MQ=255) acccccATGTGAAAAAATGCCTTCGGGGGGGCTGT‑C‑TGCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACaa < 2:273425/142‑1 (MQ=255) cccgggtgAAAAAAGCCCTCCGGGGGGGCTTT‑C‑TGCAAGTTTTTTTC‑CCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaaa < 1:783157/146‑1 (MQ=255) cATGTGCAAAAATGCCTTCGCGGCGGCTGT‑C‑TGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAACTTAAACCCCGCATGTTTGGTGGGGGTCGCCGCCGGTGACAATGTGGAGAgg > 2:125111/1‑151 (MQ=25) ggaaaaaaacgtttccgcgGGGGTTTT‑C‑TGCCGGTTTTTTTC‑CCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa < 1:960677/135‑1 (MQ=255) gaaaaaaaaTCCCTTCGCGGCGCTTTT‑C‑GGCAAGTTTTTTTC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa < 1:571137/147‑1 (MQ=255) aaaaaTCCTTCCGGGGCGGCTGT‑C‑TCCCAGCTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg < 2:381267/148‑1 (MQ=255) ctccgcgGCGGCTGC‑C‑TGCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt < 2:358245/147‑1 (MQ=255) ggttGT‑C‑TGCCAGTTATTTTC‑CCGCCACAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa < 1:652324/146‑1 (MQ=255) tGT‑C‑GGCCGGATTTTTTC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGa < 1:299525/148‑1 (MQ=255) ttttttc‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAtt < 1:926752/147‑1 (MQ=255) ttttc‑CCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt < 2:495390/146‑1 (MQ=255) tt‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAAGAAAAAACTTTAAAACCCCCCTTCTGTGTGGGAAACGCTGCTTGGAACAAAATGGGAAGGGGAGTGGGGGGGATCACCCCCACaaaccccccca > 2:582101/1‑133 (MQ=255) | CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGT‑C‑TGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATT > NC_000913/357731‑358022 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |