| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 357,890 | 0 | C | A | 43.9% | 32.0 / 19.3 | 41 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
| Reads supporting (aligned to +/- strand): ref base C (15/8); new base A (0/18); total (15/26) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 8.85e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
AAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTC > NC_000913/357745‑358033 | 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‑CGGGGCGGGTGTCTGCCAGCTATTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑aaa < 1:233055/150‑1 (MQ=255) gtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggggggggaaaa > 1:14996/1‑139 (MQ=255) cAGCGGCTCCCGTGTTTCTTTCCCCCCAACACCCGCGCTTTTACCCACCCCCGTGT‑AAAAAATGGCTTGCGCGGGGCCTGTCTGCCATCTTTTTTCCCCCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATa < 2:680066/148‑1 (MQ=255) aGCGCCTCCCGTGATTCTTTCCCCCCAACCCCCGCGCTTTAACCCACCGCCTTGTGAAAAAAGGCCTG‑CGCGGCGGCTGTCGGCCAGCTTTTTTCCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAg < 1:256580/148‑1 (MQ=255) cGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGCG‑aaaaaaa > 1:742869/1‑143 (MQ=255) ccTACCGTGATTCTTCCCCGCAAACACCCGCCCTTTTCCCCAACGCCATGTAAAAAAATCCCTC‑CGCGGCGCTTTCCTGCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAg < 1:84641/143‑1 (MQ=255) ccGTGATTCTTTCCCCCCAACACCCGCGCTTTTACCCACCCCCATGTAAAAAAATCCCTG‑CGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCg < 1:160788/151‑1 (MQ=255) gTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGGG‑aaaaaaaaataagcaaac > 1:713579/1‑135 (MQ=25) ttCTTTCCCGCCAACACCCGCGCTTTTACCCACCGCCATGTGAAAAAAGGCCTT‑CGGGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCAtg < 2:151122/150‑1 (MQ=255) ttCTTTCCCGCCAACACCCGCCCTTTTACCCACCCCCATGTGAAAAAATGCCTT‑CGCGGCGGTTGTCCGCCAGTTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCAtg < 1:729441/150‑1 (MQ=255) ttCATTCCCGCAACCACCCGCGCTTTAACCCACCCCCATTTGCAAAAATCCTTC‑CGCGGCGGCTGTCTGCCACCTTTTTTGCCGCGGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCAtg < 2:443619/150‑1 (MQ=255) tCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑AAAAAgaatttaaccc > 1:573985/1‑133 (MQ=25) cccGCGCTTTTACCCACCGCCATGTGCAAAAATGCCTT‑CGCGGCGGTTGTCTGCCAGCTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCgctg < 2:682027/150‑1 (MQ=38) aTTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAA‑AAAAAACATTTAAACCC‑CCATGTTTTGTGGGGATCGCCGCT‑gg > 2:487542/1‑148 (MQ=25) cccAACCCCATGTCAAAAAATCCCTT‑CGCGGCGGCTGCCTGCCAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACaa < 1:562401/148‑1 (MQ=255) aaCCCCATGTACAAAAATCCCTC‑CGCGGGGGCTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAAtg < 1:458175/148‑1 (MQ=255) aCGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAG‑AAAAAAA‑CTTAATCCCGTCATGGTGTGTGGGAATAGCTggt‑ggtgaaaa > 2:741978/1‑140 (MQ=37) gtaaaaaaaTCCCTC‑CGCGGCGGTTGTCTGCCAGCTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGaa < 1:50233/141‑1 (MQ=255) gggCAAAAAGGCTTG‑CGCGGGGGCTGTCGGCCAGCTATTTTCCCGCCACAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGa < 1:409403/148‑1 (MQ=255) tgCAAAAATGCCTC‑CGGGGGGTCTGTCTGCCAGTTATTTTCCCCCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAggg < 1:990987/148‑1 (MQ=255) aaaaaTGCCTC‑GGCGGGGGCTGCCGGCAATTTTTTTTCCCGCCAAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGa < 1:803090/146‑1 (MQ=255) aTGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAG‑AAAAAACATTTAACCCCCCCTGGTTTTGTGGGAA‑CGCTGCT‑GGGAAAAAAGTGGGAAGGGGAgggggg > 1:478774/1‑144 (MQ=255) cctgcccggcgtCTGTCGCCCGGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACATTGTGGAAAGGGAAGTGCGGCGaa < 1:834864/144‑1 (MQ=255) ggcggtTGTCTGCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTcc < 2:356261/151‑1 (MQ=255) ggcgccTGTCGGCCATTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTc < 1:346635/150‑1 (MQ=255) gTCGGCCTGTTATTTTCCCGCCGAAACAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc < 1:403072/150‑1 (MQ=255) tgccaatttttttcCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa < 2:360052/136‑1 (MQ=255) ccgttttttttcCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc < 2:667630/140‑1 (MQ=255) ccAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg < 1:235537/141‑1 (MQ=255) ttttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa < 2:80031/149‑1 (MQ=255) ccgcAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAAA‑TTTAAACCC‑CCCTGTTTTGGGGGAATTGGTGCTGGGTAAAAAGGGGGAAAGGAAAG‑GCGGGGGAACAACCCCACAACCCCGCCCATTTAGTTCGCCTc > 2:990743/1‑150 (MQ=255) | AAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTC > NC_000913/357745‑358033 |
| 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: ‑ |