Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,890 | 0 | C | A | 41.9% | 41.9 / 21.1 | 44 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base C (16/9); new base A (0/18); total (16/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.12e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC > NC_000913/357742‑358028 | gATAAGCGTAGCGCCTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑Gcaacaa < 1:803661/151‑1 (MQ=255) gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑Gcaaca < 1:223995/150‑1 (MQ=255) aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACaaa < 2:41241/151‑1 (MQ=255) aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACaaaa > 1:669994/1‑150 (MQ=255) aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑c < 2:135736/151‑1 (MQ=255) cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑ccgc < 2:970472/151‑1 (MQ=255) tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑ccgcc < 2:833948/150‑1 (MQ=255) cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCagt > 2:1033389/1‑150 (MQ=255) ccAGCCGCAGCCCTGTGTCAGCGGCTACCGTGTTTCATTCCCGCCAACACCCGCGCATTTACCCACCCCCATGTGAAAAAATGCTTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑ccgccg < 1:165693/134‑1 (MQ=18) cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCg > 2:113506/1‑151 (MQ=255) gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCggg > 2:771791/1‑148 (MQ=255) aCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCggggggaga > 2:889914/1‑145 (MQ=255) aCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgggg > 2:244458/1‑145 (MQ=255) ccTGTGTCGGGGGCTCCCGTGTTTCTTTTCCCCCAACACCCGCGTTTTTACCCAACCCCATGGGAAAAAATCCCTCC‑GGGGCGGCTGTCTGCCAGCTATTTTC‑CCGCC‑GCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAag < 2:888673/148‑1 (MQ=255) tgtCAGCGGCTCCGGTGTTTCATTCCCCCAAACACCCGCGCTTTTCCCCACCGCCATGTAAAAAAATGCCTTC‑GGGGCGTCTGTCGCCCAGTTTTTTTC‑CCGCC‑GCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga < 2:880820/145‑1 (MQ=255) gtCGGCGGCTCCGGTAATTCTTTCCCCCCAACACCCGCGCTTTTACCCACCCCAATGTGAAAAAATCCCTCC‑GCGGGGGCTTTCTGCCATTTATTTTC‑CCGCC‑GAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAag < 2:447368/143‑1 (MQ=255) gCGGCTCCCGTTTTTCTTTCCCCCCAACACCCGCGCTTTTACCCACCCCCTTGTAAAAAAATGCCTCC‑GGGGCGGTTTTCTGCCATTTTTTTTC‑CCGCC‑AAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑tt < 2:708350/150‑1 (MQ=255) aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑gaaaaaaatttatacc > 1:1319091/1‑138 (MQ=25) aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAG‑aaaaaaattttaaacc > 2:63188/1‑138 (MQ=25) tatttCATTCCCCCCAACACCCGCGATTTTACCCACCGCCATGTCAAAAAAGGCCTTC‑GCGGCGGCTGTCTGCCAGGTTTTTTC‑CCGCC‑GCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGcc < 2:151544/143‑1 (MQ=255) tCATTCCCCCCAACACCCACGAATTTCTCCACCCCCATGTACAAAAAGGCGTTC‑GCGGGGGTTGTCTGCACGCTTTTTTC‑CCGCC‑ACAACAAAG‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCa < 2:583011/147‑1 (MQ=255) ttcccccccAACACCCGCGCTTTTACCCACCGCCTTGTCAAAAAATCCCTTC‑GGGCGGGTTTTCCGCAACCTATTTTC‑CCGCC‑GCAAAAAAG‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCAtgttg < 2:15508/143‑1 (MQ=255) tcccccccAACACCCGCCCTTTTA‑CCAACCCCATTGTAAAAAAACCCTTCGGGGGGGGTTTTCTCCCAGTTATTTTC‑CCCCC‑AAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCAtgttgt < 1:445121/144‑1 (MQ=255) ccccaaacccccGCGCATTTACCCACCGCCAGGAGAAAAAATCCCTCG‑GGGGCGGCTGTCTGCCGTCTTTTTCC‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTtgtgtg < 2:1275957/140‑1 (MQ=255) ccccaaacaccCGCTCTTTTACCCACCCCATTGTAAAAAAAGGCCTTC‑GGGGCGGTTTTCTCCCAGTTTTTTTC‑CCGCC‑GAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTtgtg < 2:126057/143‑1 (MQ=255) ccccaccaCCCGCGCTTTCATCCCCCGCTGTGTAAAAAAATTCCCGC‑GGGGCGGCGTCCACCCAGCT‑TTTTT‑CCCCC‑GGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTtgtgtg < 2:1307605/147‑1 (MQ=255) aacaccCCCGCATTTACCCACCGCCATGTGAAAAATTGCCTTC‑GCGGCGTCTGTCTCCCATTTATTTTC‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATc < 2:185310/151‑1 (MQ=255) gcATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑AAAAAAA‑TTTTAAACCCTCCTGTTTTGTGGGAAACgc > 2:610625/1‑144 (MQ=25) tttccccACCGCCATGTGAAAAAATCCCTGC‑GGGGGGGCTTTCTCCCAGTTTTTTTC‑CCGCC‑GAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTaa < 2:1225706/145‑1 (MQ=255) tttACCCACCGCCATGTAAAAAAATGCCTCC‑GCGGGGCCTGCCGCCCAGTTTTTTTC‑CCCCC‑GAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa < 2:908709/149‑1 (MQ=255) cAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGAGGCAGAG‑AAAAAAA‑TTTTAGAACCCCACGTTGTGTGGGGATCGGCGCTGGtgaaaaagt > 2:1018711/1‑142 (MQ=21) aaCGCCATGTAAAAAAATCCCTGC‑GGGCCGGTTGTCGCCCATTTTTTTTC‑CCGCC‑GAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgtg < 1:1165073/150‑1 (MQ=255) ccATGTGAAAAAATCCCTCC‑GCGGCGGCTGTCTCCCAGCT‑TTTTC‑CCCCC‑GCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAgg < 2:1279255/151‑1 (MQ=255) tgCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAA‑AAAAAAA‑TTTAAACCCGCCACGGTTGGTGGGGAACCCTGCTGGGGACAAAGAGGGAAggg > 2:486236/1‑148 (MQ=255) aaaaaGCCCTTC‑GGGGCGGCTTTCTGCCTGTTTTTTTC‑CCCCC‑GCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa < 2:106927/147‑1 (MQ=255) aaaaaCGCCTCC‑GTGCGGGCCTTCTGCCA‑TTATTTTT‑CCCCCAACAATAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg < 2:356179/148‑1 (MQ=255) gCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGG‑AAAAAAA‑TTTTAACCCCTCCTGTTTTGTGGGGAACGCTGCTTGGGACAAAGTGGGAAGGGGAgggggggg > 2:205143/1‑142 (MQ=255) ttC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAAA‑TTTAAACCCGCCA‑GTTTTGTGGGGAAAGCTTCTGGGGACAAAGTGGGAAGGGGAggggggggg > 2:420308/1‑139 (MQ=255) gcggtTGTCTCCCAGTTTTTTTC‑CCGCC‑CAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa < 2:998501/151‑1 (MQ=255) gcggcTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAAAG‑AAAAAAA‑ATTTAACCCCCCCTGTTGTGTGGGGATCCCTGCTGGTGAAAAAGTGGGAAGGGGAGGGGGGGGGATCAGATcc > 2:1186128/1‑150 (MQ=255) ggcTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAA‑AAAAAAC‑ATTAAACCCGCCCTGTTTTGTGGGGATCACCGCCTGGAACAAATTGGGAAGGGGAGGgcgggggaagagctgc > 2:224169/1‑136 (MQ=255) gTCTGCAAGTTTTTTTC‑CCCCC‑CAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc < 2:549864/150‑1 (MQ=255) cTGCAAGTTTTTTTC‑CCGCC‑GAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg < 2:1233299/144‑1 (MQ=255) tccAGCTTTTTTTCCCGCC‑GCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc < 1:200379/145‑1 (MQ=255) cAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAAAAACATTTAAACCC‑CCCCGTTTTGTGGGAATCGCTGCCGGTGAAAAAGTGGAAAAGGGAGGGCGGGGGAACAACTcc > 2:394392/1‑138 (MQ=255) tcccgcc‑gAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc < 2:401488/148‑1 (MQ=255) | GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCC‑GCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC > NC_000913/357742‑358028 |
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: ‑ |