Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,880 | 0 | T | C | 44.4% | 28.5 / 10.5 | 27 | K271E (AAA→GAA) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base T (13/2); new base C (0/12); total (13/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.98e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAG > NC_000913/357733‑358025 | cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccg < 1:798818/151‑1 (MQ=255) tAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccgcc > 1:119651/1‑151 (MQ=255) gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac < 2:296858/149‑1 (MQ=255) tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgccg > 1:1033444/1‑151 (MQ=255) tAGCGCATCAGGCAATACCAGCCGCAGACCTGGGGCAGCGGCGACCGGGATTCATTTCCGAC‑AA‑CAA‑ACGCGCATTTATCCAACGACATGTGCACAAATGCCTT‑CGCGGCGGCTGTCTGCAAGCTATTTTTCCGCCGCACCAAAAccgccg > 1:849939/1‑151 (MQ=255) gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtc > 1:216541/1‑150 (MQ=255) aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg > 2:43617/1‑151 (MQ=255) gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg > 1:770305/1‑149 (MQ=255) gccttgtgTCGGCGCCTCCGGTTATT‑ATTTCCCCCAAA‑CAC‑CCGCCCTTTTACCCACCCCCATGTCAAAAAATGCTTC‑GGCGGGGGCTGTCTGCCGGTTTTTTTCCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAag < 2:303657/146‑1 (MQ=255) gtgtCGGCGGCGACCGTGTTTCTTCCCCCCC‑AA‑CAA‑CCGCGCATTTCCCCAACGCGATGTCAAAAAATGCCTC‑CGCGGCGGCTGTCTCCCAGCTTTTTTCCCCCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaaa < 1:888834/150‑1 (MQ=255) tCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggcgggagaaaa > 2:699384/1‑138 (MQ=255) aGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggcgaggaaaaaactt > 2:76980/1‑136 (MQ=255) gtctcccGTGTTTC‑TTTCCCCC‑AA‑CAACCCGCCCTTTTACCCAACCCCTTGTCAAAAAATCCTTT‑CGCGGCGGCTTTCTCCCAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCtt < 2:852225/143‑1 (MQ=255) ggCTACCGTGATTCATTTCCGCC‑AA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGgaaaaaaatt > 2:758301/1‑142 (MQ=25) tGATTCATTTCCGCC‑AA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAGCttttaaaccc > 2:327848/1‑141 (MQ=255) tgtttaTTTCCCCC‑AACCAC‑CCGCCCTTTTCCCCAACGCCATGTG‑AAAAAAGCCTTCCGGGGGGGCTGTCTCCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTaa < 2:919841/138‑1 (MQ=255) cccccc‑ac‑caC‑CCGCGCTTTTCCCCAACGCCATGGAAAAAAGTCCCTC‑CGCGGCGCCTGTCTGCGAGCTGTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtg < 1:7904/145‑1 (MQ=255) cccccgattttccccACCGCCAGGTGAAAAAACGCCTC‑CGGGGGGGCTTTCTGCCAGCTTTTTTCCCGCCACAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg < 1:734657/138‑1 (MQ=255) cgcacttttACCAAACCCCATGTGAAAAAATGCTTT‑CGGGGCGTCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgct < 2:840398/145‑1 (MQ=255) gcgcATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGCGAGAGAAAGATGTAACACCGTGTTGTTGTGTggtgagcata > 1:384876/1‑139 (MQ=25) ttcccccacccccATGTAAAAAAGTGCCTC‑CGGGGCGCCTTTCGCCCAGTTTTTTTCCCGCCACAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGctgctg < 1:356888/135‑1 (MQ=255) ttaacccACCGCAATGTAAAAAATTCCCTC‑CGGGGCGCCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGctgctg < 2:1003851/141‑1 (MQ=255) acccAACCCCATGTGAAAAAATGCTTC‑CGCGGCGGCTTTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACaa < 2:419089/148‑1 (MQ=255) cAACCCCAGGTGAAAAAATCCCTG‑CGCGGCGGCTTTCTCCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACa < 2:366743/146‑1 (MQ=255) aaCGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAATTTTAAACCGCCACGTTGTGTTGGAAACgccccttgttaaa > 2:197281/1‑134 (MQ=21) cAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAAGAAAAAACTTTTAACCCCTCCTGTTTTGTGGGGATCGCTGCTGGGGACAAAGTGGGAAAGGGag > 2:588503/1‑148 (MQ=255) aaagtcccTG‑CGGGGCGCCTGTCGGCCGGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGTgcg < 2:159535/145‑1 (MQ=255) ttttcCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg < 2:854897/150‑1 (MQ=255) | CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAG > NC_000913/357733‑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: ‑ |