Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,880 | 0 | T | C | 50.0% | 20.3 / 13.5 | 28 | K271E (AAA→GAA) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base T (12/2); new base C (0/14); total (12/16) | |||||||||||
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 |
TAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCA‑TGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCC‑GCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAG > NC_000913/357735‑358025 | tAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCA‑TGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTccgcc > 1:227066/1‑151 (MQ=255) ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCA‑TGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCC‑Gcaa < 2:521094/151‑1 (MQ=255) aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCA‑TGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCC‑GCAACAAAAc > 1:701323/1‑151 (MQ=255) cATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCA‑TGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCC‑GCAACAAAACCGCCGTtctc > 2:579320/1‑151 (MQ=255) tCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCA‑TGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCC‑GCAACAAAACCGCCGTtctctc > 1:72053/1‑151 (MQ=255) cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCA‑TGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCC‑GCAACAAAACCGCCGTTCTCTcc > 1:680592/1‑151 (MQ=255) ccGCTGCCCGGTGTCGGCGTCTCCCGTGTTTCTTGCCCGCCAACACCCGCTCTTTTCCCCAC‑CGCAA‑TGGACAAAAAGGCCTCCGCGGCGCTTGTCGGCGA‑TCTTTTTTCCCCCC‑AAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgggg < 2:5075/151‑1 (MQ=255) cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCA‑TGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTCTTTTTCCGCC‑GCAACAAAACCGCCGTTCTCTCCAGTAGTGGAgggg > 2:431060/1‑148 (MQ=255) gcctttggtCAGCGGCTACGGGTTTTTATTCCCGCAAACACCCGC‑CCTTTCCCCAACCGCAG‑TGTGAAAAAATGCCCCGGGGGGGGCTGTCGGCGA‑GATTTTTGCCCGCC‑GAAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg < 2:317348/137‑1 (MQ=255) aCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCA‑TGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCC‑GCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg > 2:288079/1‑145 (MQ=255) ggCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCA‑TGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCC‑GCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGggaaaggaaaaaaatt > 1:658782/1‑142 (MQ=255) ccGTGTTTCTTTCCCCCAAACACCCGCGCTTTTACCCAC‑CGCCA‑TGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCG‑GTTTTTTTCCCGCC‑GAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGc < 1:680620/148‑1 (MQ=25) tGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCA‑TGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCC‑GCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGaaaaaaattttaaacc > 1:607957/1‑134 (MQ=25) tattcttccccccaaacaccCGCTCTTTTCCCCAA‑CGCCA‑TGTGCAAAAGCCCCTCCGGGGCGGCGGTCTGCCA‑GCTATTGTCCCGCG‑GCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCa < 2:572864/136‑1 (MQ=255) ccGCCAACACCCGCCCATTTACCCAA‑CGCAA‑TGTAAAAAAAGGCCTCCGCGGCGGCTGTCTGCCA‑GCTTTTTTCCCGCC‑GCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTtgtgtg < 1:124819/150‑1 (MQ=255) gCCAACAACCGCGCATTTATCCAA‑CGCCA‑TGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCC‑GCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGGGAGAAAAAA‑GTTTTAGGCCGCCCTGTTGTgggggg > 1:418748/1‑149 (MQ=25) ccAACAACCGCGCATTTATCCAA‑CGCCA‑TGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCC‑GCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGGAAAAA‑ATTTaaagccccccagt > 1:737141/1‑130 (MQ=25) caaccacccGCTCATTTCCCCAC‑CGCCA‑TGAAAAAAATTCCCTCCGCGGCGGCGGTCTACCA‑TTTTTTTCCCCGCA‑AAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCa < 2:227020/143‑1 (MQ=255) tttACCCAA‑CGCCA‑TGTGAAAAAATCCTTCCGCGGCGGCTGTCTGCCA‑GCTTTTTTCCCGCC‑GAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTaa < 1:486689/150‑1 (MQ=255) ttaacccAG‑CGCCA‑TGTGAAAAAAGCCCTGCGGGGGGGTTTTC‑GCCAGGTTATTTTCCCGCC‑ACAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTaa < 2:202815/145‑1 (MQ=255) ttcatcAA‑CCCCATTGTGAAAAAATCCCTCCGGGGCGGCTCTC‑GCCC‑GGTTTTTTTCCCCCGGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATc < 2:696274/134‑1 (MQ=255) aa‑CGCAA‑TGTGAAAAAATGCTTTCGCGGCGGCTTTCTCCCA‑GTTTTTTTCCCGCC‑GCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAAtgtg < 1:85798/150‑1 (MQ=255) cgcggttcacAAAAATCCCTGCGGGGGGTCTGTCGGCGT‑GTTTTTTTCCCGCC‑AAAACAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGa < 2:535933/141‑1 (MQ=255) cattgtgAAAAAAGCCTTTCGGGGGGGCTTTCTCCCA‑GTTTTTTTCCCGCC‑GAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAg < 1:841424/148‑1 (MQ=255) aaaaaTGCCTTCGCGGCGGTTGTCTCCCA‑GTTTTTTTCCCGCC‑GCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAg < 1:285316/148‑1 (MQ=255) gCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCC‑GCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAGCCTTTAATCCCG‑CCCGTTTTGTGGGCATAGCTGCTGGGTGAAAATG‑Gggaagggaaggggggg > 1:393016/1‑142 (MQ=255) ggcggcTTTCCGCCT‑GTTTTTTCCCCGCC‑AAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa < 2:418646/146‑1 (MQ=255) gCCA‑GTTTTTTTCCCGCC‑GCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc < 2:225829/146‑1 (MQ=255) caattttttttcCCGCC‑GAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTCGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc < 2:143612/140‑1 (MQ=255) ttttcCCGCC‑ACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg < 1:579459/150‑1 (MQ=255) | TAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCA‑TGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCC‑GCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAG > NC_000913/357735‑358025 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |