Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,880 | 0 | T | C | 65.2% | 2.4 / 43.9 | 46 | K271E (AAA→GAA) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base T (11/5); new base C (0/30); total (11/35) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.27e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.15e-07 |
ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAG > NC_000913/357730‑358025 | aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTAttttt > 2:1218880/1‑151 (MQ=255) cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccg < 1:180967/151‑1 (MQ=255) gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccgc > 1:1254192/1‑151 (MQ=255) aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccgccg > 2:604139/1‑151 (MQ=255) aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccgccg > 2:145744/1‑151 (MQ=255) gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGca < 1:302631/151‑1 (MQ=255) aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑c < 1:527466/151‑1 (MQ=255) tCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTtctctc > 2:1002227/1‑151 (MQ=255) gcaaaTCCCAGCCGAGGCCTTGTGTCAGCGGCTCCCGTGA‑TTCATTCCCGCAAA‑CACCCGCGCTTTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTC‑CGCGGCGGCTGTCTGCAAGTTATTTTTCCGCCGCAAAAAAA‑ccgccg < 1:1185898/138‑1 (MQ=17) tcccACCCGCAGCCTGGTGCCGGGGGCTACGGTGT‑TTATTTCCCGCAAA‑AACCCGCGCTTTTA‑CCCAACCCC‑TTGTAAAAAAATGCTTC‑CGGGGGGGCTGTCGGCCAGTTTTTTTCCCGCCGAAACAAAA‑CCGCCGTTCTCTCCagtagt < 2:1211666/148‑1 (MQ=255) tCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTgg > 1:851415/1‑151 (MQ=255) ccggCCGGAGCCTTGTGCCGGCGGCTACCTTGA‑TTTTTTCCCCCAAA‑CACCCGCACATTTA‑CCCACCCCC‑TTTGGAAAAATTGCCTC‑CGCGGCGGCTGTCTGCAAGTTTTTTTCCCGCCCCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTg < 2:1015942/146‑1 (MQ=255) cAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCg > 2:998059/1‑151 (MQ=255) aGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgg > 2:97630/1‑151 (MQ=255) ccGCAGCCCTGGGCCAGCGCCTACCGTAT‑TTCTTTCCCCCCAA‑CACCCGCGCTTTTA‑CCCACCGCC‑ATGTGAAAAAATCCCTC‑CGGGGCGGTTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCggg < 2:642170/150‑1 (MQ=255) ccGCAGACCTGT‑TCAGCCGGAACCCTTATTTCATTCCCGCAAA‑CACCCGCG‑ATTTACCCCAACGCA‑AGGTAC‑AAAATGCCTTGCGCGGGGGCTGTCTGCAACCTATTTTCCCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgggg < 2:1198551/151‑1 (MQ=255) cgcaacctggtgcCAGCGGCTCCCGTTT‑TTCATTCCCCCCAA‑CACCCGCGCTTTAA‑CCCAACGCC‑TTGTAAAAAAATCCCTC‑CGCGGCGGCTGTCGGCCAGTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgggg < 1:1209561/141‑1 (MQ=255) cgaagcccTGTGTCAGGGGCTCCCGTGT‑TTCTTTCCCGCCAC‑CACCCGCGCTTTTA‑‑CCAACCCCAATGTGAAAAATTGCCTC‑CGGGGCGGCTTTCTGCTAGTTTTTTTCCCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGc < 2:1139967/145‑1 (MQ=255) cggccgtttgccggCGTCTACCGTAA‑TTCTTTCCCCCCAC‑CACCCGCGTTTTTA‑CCCACCCCC‑ATGAGAAAAAAGGCCTC‑CGGGGGGGTTTTCTGCCAGTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGc < 1:922202/136‑1 (MQ=255) gACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgg > 2:187327/1‑144 (MQ=255) cggcgcCTACCGTAA‑TTCTTTCCCCCCAA‑CACCCCCCCTTTTA‑TCCAACCCA‑ATGTGCAAAAATGCCTT‑CGGGGCGGCTGTCGGCCCGTTATTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg < 2:515583/147‑1 (MQ=255) gcctccgttGA‑TTCTTTCCCCCCAA‑CACCCCCGCATTTA‑CCCACCGCC‑ATGTGAAAAAATCCCTC‑CGCGGCGGTTTCCGGCCAGTTTTTTTCCCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATa < 1:115389/136‑1 (MQ=255) gcctcccGTGA‑TTCTTTCCCCCCAA‑CACCCGCGCTTTTA‑CCCACCGCC‑ATGTGAAAAAATGCCTC‑CGGGGGGGCTGTCTGCCAGCTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCt < 2:825477/142‑1 (MQ=255) tACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAGAttt > 1:584735/1‑146 (MQ=255) cccGTGT‑TTCTTTCCCCCCAA‑CACCCGCCCTTTTA‑‑CCAACCCCAAGGTGAAAAAATCCCTC‑CGGGGGGGCTGTCGGCAAGTTATTTTCCCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGccc < 1:139969/150‑1 (MQ=255) cccGTGT‑TTC‑TTTCCCCCAACCACCCGCGCTTTTA‑CCCACCGCC‑TGGTAAAAAAATGCCTC‑CGGGGCGGTTGTCGGCAAGTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGc < 1:1061724/148‑1 (MQ=255) aCCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTaaaacc > 2:556077/1‑147 (MQ=25) ccGTGT‑TTTTTTCCCCCCAA‑CACCCGCGCATTTA‑CCCAACGCC‑ATGTGCAAAAATGCCTC‑CGGGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGccc < 2:159503/150‑1 (MQ=25) ccGTGA‑TTCTTTCCCCCCAA‑CACCCGCGCATTTA‑CCCACCGCC‑TTGTGAAAAAATGCCTT‑CGCGGCGGCTGTCGGCCAGTTTTTTTTCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGc < 2:788332/148‑1 (MQ=255) tGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGgaaaaagttt > 2:592861/1‑134 (MQ=25) gA‑TTCTTTCCCCCCAA‑CACCCGCCTATTTA‑CCCACCCCC‑AGGTCCAAAAAGCCCTC‑CGGGGGGCTTGTCCGCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGt < 1:171552/148‑1 (MQ=255) ttCTTTCCCGCAAA‑CACCCGCGATTTTA‑CCAAACGCC‑ATGTGAAAAATTGCCTC‑CGGGGGGGTTGTCGGCCACTTATTTTCCCCCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtg < 2:452240/150‑1 (MQ=255) ccccccAA‑CACCCCCCCTTTTA‑CCCACCGCC‑ATGTGAAAAAATCCCTC‑CGCGGCGGCTGTCTGCAAGCTTTTTTCCCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgtg < 1:556026/147‑1 (MQ=255) acccGCGCATTTA‑CCCACCGCA‑ATGTGAAAAAATCCCTT‑CGGGGCGGTTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc < 2:845587/146‑1 (MQ=255) tttA‑TCCACCGCC‑ATGAAAAAAAGTCCCTG‑CGCGGGGGTTTCCGCCCAGTTTTTTTCCCCCCACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGt < 1:786247/148‑1 (MQ=255) ccACCGCC‑ATGTGAAAAAATGCCTC‑CGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt < 1:370226/149‑1 (MQ=255) ccAACGCC‑ATGTGAAAAAATGCCTC‑CGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACaa < 2:954857/148‑1 (MQ=255) cGCC‑ATGTAAAAAAATGCCTT‑CGGGGCGGCTTTCTGCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGACAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt < 2:14748/147‑1 (MQ=255) cc‑ATGTGAAAAAAGCCCTC‑CGCGGCGGCTGTCGGCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaaa < 2:449158/150‑1 (MQ=255) aaaaaTGCTTG‑CGTGGCGGATGACTGCAACAAATTTTCCCGCCGAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg < 2:462029/148‑1 (MQ=255) aaaaaGGCCTG‑CGCGGCGGCTTTCTCCCATCTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt < 1:568480/149‑1 (MQ=255) aaaccctcccgcgCCGGTTGTCAGCCGTTTTTTTTCCCCCCCCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg < 1:278901/142‑1 (MQ=255) gggggggcTGTCGGCCGGTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATc < 2:1066980/142‑1 (MQ=255) cggcTTTCCGCCAGCTATTTCCCCGCCACAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa < 1:812295/150‑1 (MQ=255) cTGCCAGTTTTTTTCCCCCCCCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg < 1:1143432/144‑1 (MQ=255) ggccggttATTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc < 1:1218756/139‑1 (MQ=255) ccAGTTTTTTTCCCGCCGAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc < 1:145731/145‑1 (MQ=255) ccactttttttCCCACCCAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa < 1:1002127/145‑1 (MQ=255) tatttttcCCGCCGAAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa < 1:1230355/127‑1 (MQ=255) ttttttcCCGCCAAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTaa < 1:1064137/145‑1 (MQ=255) ttttttcCCCCCGAAAAAAAT‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt < 1:899604/146‑1 (MQ=255) tttttcCCCCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg < 2:551856/148‑1 (MQ=255) ttttcCCGCCGAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa < 2:25172/149‑1 (MQ=255) ttttcCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg < 2:694655/150‑1 (MQ=255) | ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAG > NC_000913/357730‑358025 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |