| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 357,880 | 0 | T | C | 53.8% | 27.4 / 56.2 | 40 | K271E (AAA→GAA) | cynR | transcriptional activator of cyn operon; autorepressor |
| Reads supporting (aligned to +/- strand): ref base T (15/3); new base C (0/21); total (15/25) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.25e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.10e-08 | |||||||||||
AGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGA > NC_000913/357736‑358024 | aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccgcc < 2:648091/150‑1 (MQ=255) ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaa < 2:39649/151‑1 (MQ=255) ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGca > 1:286530/1‑150 (MQ=255) aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc < 1:565513/151‑1 (MQ=255) gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgcc > 1:798586/1‑151 (MQ=255) aGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGt > 2:774691/1‑151 (MQ=255) gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATACCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtc > 1:25661/1‑150 (MQ=255) tCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctctc > 2:798994/1‑151 (MQ=255) aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc > 2:71558/1‑150 (MQ=255) ccAGCCGCAGCCCTGTGTCAGCGGCTCCCGTGTTTCATTCCCCCCAACACCCGCGCATTTATCCACCGCCATGTGAAAAAATGCTTC‑CGCGGCGGCTTTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGc < 2:669285/151‑1 (MQ=255) ccAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTgg > 1:653304/1‑150 (MQ=255) aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg > 1:133839/1‑151 (MQ=255) cccagccctttccacgcgcCTACCGTAGTTCTTTCCCCCCAACACCCGCACATTTCCCCACCGCCTTGGAAAAAAATGCCTC‑CGCGGCGGTTGCCGCCCGGCTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGc < 1:42657/136‑1 (MQ=255) gtgtCAGCGCCTCCCGTGATTCATCCCCCCCAACACCCGCCCTTTTACCCAACGCCATGTCCAAAAATGCCTG‑CGGGGCGGGTGTCTGCCAGCTATTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaaa < 1:233055/150‑1 (MQ=255) gtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggggggggaaaa > 1:14996/1‑139 (MQ=255) cAGCGGCTCCCGTGTTTCTTTCCCCCCAACACCCGCGCTTTTACCCACCCCCGTGT‑AAAAAATGGCTTGCGCGGGGCCTGTCTGCCATCTTTTTTCCCCCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATa < 2:680066/148‑1 (MQ=255) aGCGCCTCCCGTGATTCTTTCCCCCCAACCCCCGCGCTTTAACCCACCGCCTTGTGAAAAAAGGCCTG‑CGCGGCGGCTGTCGGCCAGCTTTTTTCCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg < 1:256580/148‑1 (MQ=255) cGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGCGaaaaaaa > 1:742869/1‑143 (MQ=255) ccTACCGTGATTCTTCCCCGCAAACACCCGCCCTTTTCCCCAACGCCATGTAAAAAAATCCCTC‑CGCGGCGCTTTCCTGCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg < 1:84641/143‑1 (MQ=255) ccGTGATTCTTTCCCCCCAACACCCGCGCTTTTACCCACCCCCATGTAAAAAAATCCCTG‑CGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCg < 1:160788/151‑1 (MQ=255) gTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGGGaaaaaaaaataagcaaac > 1:713579/1‑135 (MQ=25) ttCTTTCCCGCCAACACCCGCGCTTTTACCCACCGCCATGTGAAAAAAGGCCTT‑CGGGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCAtg < 2:151122/150‑1 (MQ=255) ttCTTTCCCGCCAACACCCGCCCTTTTACCCACCCCCATGTGAAAAAATGCCTT‑CGCGGCGGTTGTCCGCCAGTTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCAtg < 1:729441/150‑1 (MQ=255) ttCATTCCCGCAACCACCCGCGCTTTAACCCACCCCCATTTGCAAAAATCCTTC‑CGCGGCGGCTGTCTGCCACCTTTTTTGCCGCGGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCAtg < 2:443619/150‑1 (MQ=255) tCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAgaatttaaccc > 1:573985/1‑133 (MQ=25) cccGCGCTTTTACCCACCGCCATGTGCAAAAATGCCTT‑CGCGGCGGTTGTCTGCCAGCTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCgctg < 2:682027/150‑1 (MQ=38) aTTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAAAAAAAACATTTAAACCC‑CCATGTTTTGTGGGGATCGCCGCTgg > 2:487542/1‑148 (MQ=25) cccAACCCCATGTCAAAAAATCCCTT‑CGCGGCGGCTGCCTGCCAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACaa < 1:562401/148‑1 (MQ=255) aaCCCCATGTACAAAAATCCCTC‑CGCGGGGGCTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg < 1:458175/148‑1 (MQ=255) aCGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAA‑CTTAATCCCGTCATGGTGTGTGGGAATAGCTggtggtgaaaa > 2:741978/1‑140 (MQ=37) gtaaaaaaaTCCCTC‑CGCGGCGGTTGTCTGCCAGCTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaa < 1:50233/141‑1 (MQ=255) gggCAAAAAGGCTTG‑CGCGGGGGCTGTCGGCCAGCTATTTTCCCGCCACAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa < 1:409403/148‑1 (MQ=255) tgCAAAAATGCCTC‑CGGGGGGTCTGTCTGCCAGTTATTTTCCCCCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg < 1:990987/148‑1 (MQ=255) aaaaaTGCCTC‑GGCGGGGGCTGCCGGCAATTTTTTTTCCCGCCAAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa < 1:803090/146‑1 (MQ=255) aTGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAACATTTAACCCCCCCTGGTTTTGTGGGAA‑CGCTGCTGGGAAAAAAGTGGGAAGGGGAgggggg > 1:478774/1‑144 (MQ=255) cctgcccggcgtCTGTCGCCCGGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACATTGTGGAAAGGGAAGTGCGGCGaa < 1:834864/144‑1 (MQ=255) ggcggtTGTCTGCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTcc < 2:356261/151‑1 (MQ=255) ggcgccTGTCGGCCATTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTc < 1:346635/150‑1 (MQ=255) gTCGGCCTGTTATTTTCCCGCCGAAACAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc < 1:403072/150‑1 (MQ=255) tgccaatttttttcCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa < 2:360052/136‑1 (MQ=255) ccgttttttttcCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc < 2:667630/140‑1 (MQ=255) ccAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg < 1:235537/141‑1 (MQ=255) ttttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa < 2:80031/149‑1 (MQ=255) | AGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGA > NC_000913/357736‑358024 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |