Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,880 | 0 | T | C | 44.9% | 60.4 / 17.3 | 49 | K271E (AAA→GAA) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base T (24/3); new base C (0/22); total (24/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.18e-11 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGA > NC_000913/357734‑358024 | gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccgc < 1:1162370/151‑1 (MQ=255) aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccgccg > 2:150231/1‑151 (MQ=255) aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccgcc < 2:499723/150‑1 (MQ=255) tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac > 2:1029316/1‑150 (MQ=255) gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca > 2:530222/1‑150 (MQ=255) aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑c > 1:546255/1‑151 (MQ=255) aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑cc > 1:770908/1‑151 (MQ=255) cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑ccgc < 1:23165/151‑1 (MQ=255) tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑ccgccg > 1:227780/1‑151 (MQ=255) cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGtt > 1:1062702/1‑150 (MQ=255) cATCAGGCATTTCCAGCCGCAGACCTGTGTCAGCGGCTCCCGTGATTCATTTCCGCCAACACCCGCGCATTTACCCAACGCAATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAACAAAA‑CCGCCGTtc < 1:424512/149‑1 (MQ=255) cATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTtctc > 1:277923/1‑151 (MQ=255) cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTcc > 2:1049059/1‑151 (MQ=255) aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCagt > 1:245153/1‑149 (MQ=255) aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCagt > 1:1072101/1‑149 (MQ=255) aTTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATATTTTCGCCGGAACAAAA‑CCGCCGctctctaccctag > 1:952123/1‑143 (MQ=11) aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCagtagt > 2:830113/1‑151 (MQ=255) cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGc > 2:356875/1‑150 (MQ=255) aGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑gggggg > 2:678940/1‑147 (MQ=255) cttgttcctgCGGCTCCTGTTTTTCTTTCCCGCAAACACCCGCGAATTCCCCCACCGCCTTGCGAAAAAATC‑CTTG‑CGCGGCGGCTGTCCGCCAGTTTTTTTCCCCCCCCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑gggg < 2:771676/135‑1 (MQ=255) ttgacggcgcccaCCGTG‑TTCCTTTCCCCCAACCACCCTCTCTTTAACCAAGCCCAATTTCAAAAAAG‑CCTTCCGGGGCGGCTGGCTCCCATTTTTTTTCCCTCCCAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAgaga < 1:1026974/136‑1 (MQ=255) gtgtCAGCGGCTACCTTGATTCTTTCCCCCCAACACCCGCCCATTTACCCAACGCCATGTGAAAAAATG‑CCTT‑CGGGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAag < 1:902492/145‑1 (MQ=255) gtCGGCGGCTCCCGTTTTTCTTCCCCCCCAACACCCGCCCTTTTACCCACCGCCATGTAAAAAAATG‑CCTC‑CGGGGCGGCTGTCGCCCAGTTTTTTTCCCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAg < 2:1127381/151‑1 (MQ=255) cGGCTCCCGTGATTCATTCCCGCCACCACCCGCCCTTTCACCCAACCCCATGT‑AAAAAATGCCCTC‑CGCGGCGGCTGTCGGCCGGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCttt < 2:461991/150‑1 (MQ=255) cGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGgaaaaaactt > 2:540759/1‑143 (MQ=25) cggTGATTCTTTCCCCCCAACACCCCCGCTTTTACCCACCGCCTTGTGAAAAAATG‑CCTG‑CGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGcc < 1:265176/147‑1 (MQ=255) ccgtttttCATTCCCCCCAACACCCGCCCTTTTCCCCAACGCCTTGTGAAAAAATC‑CCTT‑CGGGGGGGCTGTCTCCCATTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATa < 1:1181730/133‑1 (MQ=255) gTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGGGAAAAAGAATT‑aaatc > 2:833654/1‑144 (MQ=25) ttcttcccccccAACACCCGCGCTTTTCCCCACCGCAATGAAAAAAAATG‑CCTG‑CGGGGCGGCTGTCCGCCGGTTTTTTTCCCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGc < 2:384/132‑1 (MQ=255) ccccccAACACCCGCGCTTTTCCCCACCCCTTTGTAAAAAAATG‑CCTG‑CGGGGGGGCTGTCTCCCGTTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTtgtgtg < 2:1078359/147‑1 (MQ=255) ccGCCAACACCCGCG‑ATTTTCCCAACCCCTTTTGGAAAAAAGCCCTC‑CGGGGCGGCTGTCCGCCGTCTTTTTTCCCCCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTtgtgtg < 1:168647/150‑1 (MQ=255) acccGCGCATTTACCCAACGCCATGTGAAAAAATC‑CCTT‑CGCGGCGGCTGTCGGCCAGTTTTTTTCCCCCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCg < 1:1057613/147‑1 (MQ=255) ccgccattttACCCAACCCCATGAAAAAAAAGC‑CCTC‑CGGGGCGGTTGTCTCCCGTTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCgctgc < 2:289570/144‑1 (MQ=255) ccGCGCTTTTATCCACCGCCATGTGAAAAAATG‑CCTC‑CGCGGCGGCTGTCGGCCATTTTTTTTCCCCCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCgctg < 1:562994/150‑1 (MQ=255) cgcgCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCGAGAGGAAAAAATTT‑TATCCCCCCCTGTTTTGTGGGGATCggt > 1:109606/1‑146 (MQ=25) cccATTTACCCACCCCCATGTAAAAAAATG‑CCTC‑CGCGGGGGCTGTCTGCCAGTTTTTTTCCCCCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCgctgct < 1:36258/147‑1 (MQ=255) cctttcaccaaccGCTATGTAAAAAAATC‑CCTC‑CGCGGCGGTTGTCGGCCATTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGt < 2:965972/139‑1 (MQ=255) cATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGAGAAAAAAATTT‑TAACCCCTCATGTTGTGTTGGTATCGCCGCTGGtg > 2:598460/1‑150 (MQ=37) ccATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGCAAAAAAAAAAAGTTT‑AAACCCG‑CCTTTTTTGTGGGGATCCGTGCTGGGGACAAAGTgg > 1:281691/1‑147 (MQ=19) gtaaaaaaaG‑CCTTCCGGGGGGGCTGTCTCCCATTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa < 1:1107354/147‑1 (MQ=255) cAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGAAAAAAATTTAAAACCCTCCATGTT‑TGTGGGAATCGCCGCCGGGAAAAAAGTGGGAAggg > 1:2460/1‑146 (MQ=255) aaaaaGG‑CCGT‑CGCGGCGGCTGTCTGCCTGTTATTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg < 2:677445/148‑1 (MQ=255) aaaaAG‑CCTTCCGCGGCGGTTGTCGGCAACCTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgc < 1:336790/151‑1 (MQ=255) aaaTG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGCAGAGAAAAAAATTT‑TAATCCCCCCTTTTGTGTGGGGATTCCTGCTGGTGACAAAGTGGGAAGGGGAgggggg > 2:825293/1‑146 (MQ=255) gggggggcTGTCCCCCATTTATTTTCCCCCCACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggcg < 1:527057/138‑1 (MQ=255) cggcggCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGAAAAAACTTT‑TAAACCCCCCTCTTTTTTGGGGATCGCTGCTGGGTACAAATTGGGAAAGGGAGGgcgggggaaaaaccc > 1:449998/1‑140 (MQ=255) cggcTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGAGAAAAAGATTT‑TAACCCGTCCTGTTGTGTGGGAATAGCTGCTGGTTAAAAAGTGGGAAGGGGAGGggggcggaaaaaccc > 1:1160813/1‑139 (MQ=255) gcttccTGCCGGTTTTTTTCCCCCAGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa < 2:560124/137‑1 (MQ=255) ccgccaagtttttttcCCCCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc < 2:1072092/138‑1 (MQ=255) gcccattttttttCCCGCACAAAAAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgcg < 2:227779/142‑1 (MQ=255) agttATTTTCCAACCGCAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTaa < 1:181568/147‑1 (MQ=255) ttttcCCGCCCAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa < 1:870831/149‑1 (MQ=255) | GTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTT‑AAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGA > NC_000913/357734‑358024 |
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: ‑ |