Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,880 | 0 | T | C | 47.4% | 38.9 / 13.0 | 39 | K271E (AAA→GAA) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base T (18/2); new base C (0/18); total (18/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.69e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGT > NC_000913/357733‑358026 | cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTccg > 2:86869/1‑151 (MQ=255) aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTccgccg > 2:942897/1‑151 (MQ=255) cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac > 1:133534/1‑151 (MQ=255) cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgc > 1:375463/1‑151 (MQ=255) aGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGt > 1:535071/1‑151 (MQ=255) cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtc > 1:349805/1‑151 (MQ=255) gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtct > 2:541577/1‑151 (MQ=255) attcccAGCCGCAGCCTTGTGTCAGCGGCTACCGTGATTCTTTCCCCCAAACACCCGC‑GCTTTTATCCACCGCCTTGTGAAAAAATCCTTTC‑GCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCagta < 2:763287/146‑1 (MQ=18) tcccAGCCGAAGCCCTGTGTCAGCGCCTACCGTAATTCATTCCCCCCAACACCCGC‑GCTTTTATCCAACCCCTTGTGAAAAAATGCCTTC‑GGGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCagtagt < 2:163983/148‑1 (MQ=255) tCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTgg > 1:45126/1‑151 (MQ=255) ccAGCCGCAGACTTGTGTCAGCGGCTCCCGTGATTCATTCCCCCCAACACCCGC‑GCTTTTACCCAACGCCATGTGAAAAAATGCCTCC‑GGGGCGGCTGTCTGCCAGCTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCagtagt < 2:296652/148‑1 (MQ=255) ccGGAGCCCTGTGTCTGGGGCTCCCGTGATTCATTCCCGCCAACACCCGC‑GCATTTACCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgg < 2:268957/149‑1 (MQ=18) ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg > 2:421789/1‑151 (MQ=255) gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg > 1:715467/1‑147 (MQ=255) gccttgtgTTAGGGGCAACCGT‑ATTCCATTCCGCCAACTACCGCGGCTTTTACCCAACCCAATGTGAAAAAATCCTTTC‑GCGGCGGCTGTCGCCCAGCTATTTCACCCCCCAAACAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCa < 2:402837/144‑1 (MQ=255) cgtggccggCGACTCCCGTGATTCTTTCCCGCAAACACCCGC‑CCATTTACCCACCGCATTGTGCAAAAAGCCCTCG‑GCGGGGCTTGTCCGCCGGTTTTTTTCCCCCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga < 1:42925/140‑1 (MQ=255) ccttttttCGTTCCCGCCACCAACCGC‑GCATTTACCCAGCCCTATGTGCAAAAAT‑CCTTCGGCGGGGTCTTCCCGCTATCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAg < 2:188540/141‑1 (MQ=255) ccGTTATTCATTCCCGCCAACACCCGC‑GATTTTATCCAACGCCATGTGAAAAAATGCCTTC‑GGGGCGGCTGTCTCCCAGCTATTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCg < 2:316783/151‑1 (MQ=255) cGTGATTCATTTCCGCCAACAACCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGG‑G‑aaaaaaactttaaaccc > 1:526064/1‑136 (MQ=25) ttCATTTCCGCCAACAACCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGA‑A‑AAAAACCATTaacacc > 2:433699/1‑140 (MQ=25) ttcccGCAACCACCCGC‑CCATTTACCCACCCCCATGGCAAAAAATGCCTCC‑GCGGCGGTTGTCTGCCAGTTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAtgttgt < 1:391667/147‑1 (MQ=255) tCCGCCAACAACCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAAAAAATTTTAACCCGCCAAGTtttgtg > 2:650954/1‑151 (MQ=25) ccAACACCCGC‑GCTTTTACCCACCCCCATGTGAAAAAATCCCTTC‑GGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCaa < 1:65122/151‑1 (MQ=255) accacccGC‑GCTTTTCCCCACCGCTATGGGAAAAAATCCCTTC‑GGGGCGGCTTCCTGCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc < 1:606031/146‑1 (MQ=255) accccc‑gcTTTTACCCAACGCCATTTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGTTATTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct < 2:919631/146‑1 (MQ=255) aacCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCCGA‑G‑GAAAAGCTTTAAAAGCGGCCTGTTGTGTgggga > 2:706203/1‑143 (MQ=25) acCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGG‑AAAAAGCTTAAACCCCgcacatttt > 2:322970/1‑130 (MQ=25) cgc‑gCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAAA‑G‑AAAAAAATTTTAACCCCCCACGTTTTGTGGGAATc > 2:605593/1‑145 (MQ=25) cGCCATTTGCAAAAATCCTTTC‑GCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa < 2:28460/150‑1 (MQ=255) cGCCATGTGAAAAAAGGCTTTC‑GCGGCGGCTGCCTGCCGACTTTTTCCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa < 2:873466/150‑1 (MQ=255) aTGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGA‑GAAAAAAATTTTAAACCCG‑CATGTTTTGGGGGAAACGCTGCCGGGGACaaaatgggaaagg > 2:511680/1‑140 (MQ=255) tgtgAAAAAATGCCTGC‑GGGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAgg < 2:375506/149‑1 (MQ=255) gtgAAAAAATCCCTCG‑GGGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa < 1:942787/150‑1 (MQ=255) tgaaaaatTGCTTTC‑GCGGGGGCTGTCTGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa < 2:749586/146‑1 (MQ=255) tgCAAAAATCCTTTC‑GCGGGGGCTGTCTGCCAGTTATTTTTCCCCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa < 1:779015/150‑1 (MQ=255) aaaaaTGCCTCC‑GCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg < 2:535129/148‑1 (MQ=255) gcggcgGCTGTCTCCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAg < 1:749937/149‑1 (MQ=255) tGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑G‑AAAAAACTTTTAAACCCCCATGTTTTGTGGGGATAGCTGCTGGTGAAAAAGTGGGAAGGGGAGGGGGGCGGATCAGCCCCAGAAACccgc > 2:15777/1‑151 (MQ=255) cccAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt < 1:63488/150‑1 (MQ=255) ccattttttttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc < 1:312226/137‑1 (MQ=255) ttttcCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg < 2:334764/150‑1 (MQ=255) ttCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGTGCCAGG‑G‑AAAAAACATTAAACCCCCCCTGTTGTGTGGCGATAGCCGCTGGTGACAATTTGGGAAGGGGAGTGCGGGGGGTCAGCCCCACAACCCCGCCAAATGaat > 1:502373/1‑146 (MQ=255) | CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGC‑GCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGT > NC_000913/357733‑358026 |
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: ‑ |