| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 357,890 | 0 | C | A | 48.5% | 16.2 / 19.4 | 33 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
| Reads supporting (aligned to +/- strand): ref base C (8/9); new base A (0/16); total (8/25) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.68e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01 | |||||||||||
GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGT‑C‑TGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC > NC_000913/357742‑358028 | 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 > 2:880994/1‑151 (MQ=255) aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGT‑C‑TGCCAGCTATTTTT‑CCGCCGCAACAAAAcc > 1:528184/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) tcccgccgCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc < 1:1019566/148‑1 (MQ=255) ccgcccAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc < 2:1040438/148‑1 (MQ=255) | GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGT‑C‑TGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC > NC_000913/357742‑358028 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |