Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,880 | 0 | T | C | 44.2% | 58.9 / 20.4 | 52 | K271E (AAA→GAA) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base T (23/6); new base C (0/23); total (23/29) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.64e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGA > NC_000913/357730‑358024 | aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTAttttt > 2:927869/1‑151 (MQ=255) aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTAttttt > 2:656118/1‑151 (MQ=255) cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTcc < 1:421444/150‑1 (MQ=255) ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCcgccgc > 1:652724/1‑151 (MQ=255) gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCGca < 2:313630/151‑1 (MQ=255) ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCGca > 1:1075262/1‑150 (MQ=255) ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCGca > 2:884057/1‑150 (MQ=255) ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCGca < 2:1091368/150‑1 (MQ=255) cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac > 1:935568/1‑151 (MQ=255) gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGCCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca < 2:667789/150‑1 (MQ=255) gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca > 2:555429/1‑150 (MQ=255) aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa > 2:291608/1‑150 (MQ=255) tAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa > 1:395012/1‑150 (MQ=255) aaGCGGAGATCATCAGGCAATTCCAGCCCCAGGCCTTTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc < 2:511056/151‑1 (MQ=18) gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtct > 2:816823/1‑151 (MQ=255) cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCag > 1:1030471/1‑149 (MQ=255) attcccccccGCAGCCTTGTCTCGGCGCCTACGGTTTTTCATTCCCCCC‑AACACCCCCCCTTTTACCCACCGCC‑TTTTAAAAAAATCCTTCCGC‑GGCGCCTTTCTCCCATTTTTTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCag < 2:688152/140‑1 (MQ=255) ccAGCCGCGGCCTTTTGTCGGCGCCCCCGGTGATCCTTTTCCCCC‑AACACCCGCTATTTTACCCACCGCG‑ATGGAAAAAAATGCTTCCGC‑GGCGGCTGTCCGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTg < 2:636503/149‑1 (MQ=255) ccAGCCGCAAACTTGTGCCAGGGCCTCCGGTGTTTCTTTCCCGCA‑AAAACCCGCGCATTTA‑CCAACCCCGGTGTGAAAAAATCCCTCCGC‑GGCGGCTGTCTGCAAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCagtagt < 1:548180/148‑1 (MQ=255) aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg > 1:811073/1‑151 (MQ=255) aGACCTGTGTCCGCGGCTACCGTGTTTCATTCCCGCC‑AACACCCGCGCTTTTACCCAACCCC‑TTTTAAAAAATTGCCTCCGC‑GGCGGTTGTCTGCCGGCTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgg‑gg < 1:829976/147‑1 (MQ=255) aGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg‑gg > 1:547509/1‑147 (MQ=255) gccttgttTCGGCGGCTCCCGTGTTTCTTTTCCGCA‑ACCACCCGCGATTTTATCCACCTCA‑TTGTTAAAAAGTCCCTCCGG‑GGCGCCGGTCTGCCATTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAag < 1:308815/146‑1 (MQ=255) aCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg‑gg > 2:521682/1‑145 (MQ=255) ttttccggcgcCTCCCGTATTTCTTTCCCCCC‑AACACCCCCCCTTTTCCCCACCCCC‑TTGTAAAAAACTCCCTCCGC‑GGCGGTTTTCTCCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAgaga < 2:521458/142‑1 (MQ=255) tttcgggggCTACCGTTATTCTTTCCCCCA‑AACACCCGCGCTTTTCCCCACCGCA‑ATGTCAAAAAACCCCTCGGG‑GGCGGTGGTCGGCCACTTTTTTTCCCGCCACAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGaa < 1:775326/141‑1 (MQ=255) tgtCGGCGGC‑ACCGGTATTTAATTCCCCCAAACACCCGCGTATTTACCCAACCCA‑ATGTACAAAAATGCTTTGGC‑GGCGGCTGTCTGCCATTTTTTTTCCCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAAt < 1:272983/150‑1 (MQ=255) gccggCGGCTACGGTTTTTCTTCCCCCCC‑AACACCCGCCCTTTTACCCACCCCC‑ATGTGAAAAATTCCCTCCGG‑GGGGGCTGTCGGCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGc < 2:451251/136‑1 (MQ=255) cAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGGGAGGGaaa > 2:355818/1‑146 (MQ=255) gCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGGGAGGGaaaaaaattt > 2:1102774/1‑141 (MQ=25) cTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTTCGC‑AGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGGGAGAGaaaaaaactt > 1:632429/1‑140 (MQ=25) tACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGAGGCAaggaacaaaaacttac > 2:240183/1‑135 (MQ=25) tACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGGAGGAGaaaaaaactt > 2:1052696/1‑139 (MQ=255) gATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGGCAGAGaaaaaaattttaaaccc > 1:591254/1‑133 (MQ=25) ttcttccccccc‑ATCACCCGCGCATTAATCCAACACC‑ATCTGAAAAAATG‑CTGCGCGGGCGTTTATCACCCATTTTTGTTTCCGCCCCAACAAAGCCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAG‑CTTTAAGCCCGTCAt < 1:704211/139‑1 (MQ=255) ttcccGCA‑AACACCCGCGCTTTTACCCACCGCC‑ATGTAAAAAATTGCCTTGGG‑GGCGGCTGTCTGCCGGTTTTTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAG‑CTTTAAGCCCGTCAtgttg < 1:1086462/146‑1 (MQ=255) ccGCC‑AACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGGGAGAGAAAAAA‑AATTAAAACCGTCATGTTGTGTgg > 2:1031655/1‑151 (MQ=25) caaacaccCGCCCTTTAACCCACCCCC‑ATTAGCAAAATTCCCTTCGC‑GGCGGCTTTCTCCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGTAATAG‑CTTTAAGCCCGTCATGTTGTGTGGc < 2:388299/147‑1 (MQ=255) c‑aacaaCCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGGAAAAGAATTTAGGCCCgcctggt > 2:475887/1‑135 (MQ=25) acccGCGATTTTACCCACCGCC‑GTTTAAAAAAGTCCCTCCGC‑GCCGGTTTTCTCCCTGTTATTTTCCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTtgtg < 1:25973/138‑1 (MQ=255) ccgcctttttACCCACCGCC‑ATGGAAAAAAATGCCTGCGG‑GGCGGCTGTCTCCCACTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCgctg < 2:97961/143‑1 (MQ=255) ccGCGCATTTACCCACCCCA‑ATGTAAAAAAATCCCTCCGC‑GGGGGCTGTCTCCCAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCgc < 1:1093902/148‑1 (MQ=255) cgcTTTTACCCACCCCC‑ATGTAAAAAAATGCTTCCGC‑GGCGTTTCCCGCCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCgctgc < 1:888313/148‑1 (MQ=255) aTTTATCCAACGCC‑ATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGGGAGAGAAAAAA‑CATTAAAACCCTCCCGTTGTGTGGGGATCgctgcttg > 1:16556/1‑146 (MQ=25) ccACCGCC‑TTGGGAAAAATTCCCTCCGC‑GGCGGCTGTTTGCCATTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt < 1:763301/149‑1 (MQ=255) accccc‑ATGTGAAAAAATGCCTTCGG‑GGGGGCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTgg < 1:291620/147‑1 (MQ=255) cc‑ATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGGAAGAGAAAAAA‑CTTTTAACCCCCCCTGTTGTGTGGGGATCGCCGCTgggtaaaaaatgggaa > 1:281252/1‑136 (MQ=25) ccatggaaaaaaaTGCCTCCGC‑GGCGGTTGTCGCCCAGTTTTTTTCCCGCGACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAg < 2:729331/143‑1 (MQ=255) cgatgtaaaaaaaGCCCTCGGGGGCGGTTGTCTCCCATTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAg < 2:1101658/144‑1 (MQ=255) tgtgAAAAAATGCTTTCGG‑GGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg < 1:281888/150‑1 (MQ=255) aaaaaTGCCTCCGG‑GGCGGCTGTCTGCAAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa < 1:1047391/147‑1 (MQ=255) ccatttttttttcCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc < 1:240193/138‑1 (MQ=255) cAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGa < 2:287126/141‑1 (MQ=255) gCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCCAGAGAAAAAA‑CTTTTAACCCCGCCCGTTGTTTGGGGATTGCCGCTGGGGACAATGTGGGAAGGGGAGGGGGGGGGATAAGCTCCCAAAACgcgcccaat > 1:1090809/1‑146 (MQ=255) tATTTTCCTCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAtt < 2:456992/149‑1 (MQ=255) tCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGGCAAAGAAAAAA‑GTTTTAGCCCCTCCTGTTGTGTGGCAATCGCTGCTGGTAAAAATGAGGGAAAGGAAGGGGGGCGGATCAGCACCCGAACCCCGCCAATTGa > 1:1065843/1‑145 (MQ=255) | ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGA > NC_000913/357730‑358024 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |