| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 357,890 | 0 | C | A | 60.9% | ‑6.5 / 21.4 | 23 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
| Reads supporting (aligned to +/- strand): ref base C (8/1); new base A (1/13); total (9/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.55e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-03 | |||||||||||
GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCG‑CGCATTTATCCAACGCC‑ATGTGC‑AAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAT > NC_000913/357742‑358021 | gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCG‑CGCATTTATCCAACGCC‑ATGTGC‑AAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCG‑caacaa > 2:11820/1‑151 (MQ=255)gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCG‑CGCATTTATCCAACGCC‑ATGTGC‑AAAAATGCCTTCGCGGCGGCTGTCTGCCAGCCATTTTT‑CCGCCG‑caacaa > 2:295898/1‑151 (MQ=255) gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCG‑CGCATTTATCCAACGCC‑ATGTGC‑AAAAATGCCTTCGCGGCGGCTGTCTGCCAGCAATTTTT‑CCCCCA‑CAACAAAAcc > 2:75283/1‑150 (MQ=255) gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCG‑CGCATTTATCCAACGCC‑ATGTGC‑AAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCG‑CAACAAAACCGCCGtt > 1:215870/1‑149 (MQ=255) gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCG‑CGCATTTATCCAACGCC‑ATGTGC‑AAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCG‑CAAAAAAACCGCCGtt > 2:68713/1‑149 (MQ=255) cccccccccAGCCCTGTGCCAGGGGCCCCCGTGTTTC‑TTTCCCCCAACCACCCG‑CCCTTTTACCCACCCCC‑TTGTCA‑AAAAAGGCTTGCGGGGGGGCTGTCTGCCAGTTTTTTTC‑CCGCCG‑CAAAAAACCCGCCGTTCTCTCCAGTAGTgg < 2:524930/143‑1 (MQ=255) ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCG‑CGCATTTATCCAACGCC‑ATGTGC‑AAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg > 1:133044/1‑151 (MQ=255) gACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCG‑CGCATTTATCCAACGCC‑ATGTGC‑AAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCgaga > 2:535976/1‑151 (MQ=255) ccTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCG‑CGCATTTATCCAACGCC‑ATGTGC‑AAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggcaaga > 2:370797/1‑144 (MQ=255) ccTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCG‑CGCATTTATCCAACGCC‑ATGTGC‑AAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGCGGCGGGGccttacg > 2:74812/1‑145 (MQ=255) gtgtCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCG‑CGCATTTATCCAACGCC‑ATGTGC‑AAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCG‑CAACAAAACCGCCGGTCTCTCCAGTAGTGGCGGGGGCAGAGgaaa > 2:569446/1‑150 (MQ=255) gtCAGCGGCTACCGGGATTCATTTCCGCCAA‑CAACCG‑CGCATTTATCCAACGCC‑ATGTGC‑AAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAAATaaaaaa > 2:53832/1‑148 (MQ=255) cggcctcccGTGATTCTTTCCCGCAAA‑CACCCG‑CGTTTTTACCCACCCCC‑ATGTCA‑AAAAGTCCCTTCGGGGGGGTTGTCTGCAATTTTTTTTC‑CCGCCA‑AAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCtt < 1:373432/143‑1 (MQ=255) tatttattccccGCAAA‑CACCCC‑GGCATTTA‑CCAACCCCGATGAGA‑AAAATTCCCTCCGGGGGGGCTGTCGCCCATTTTTTTTC‑CCGCCG‑CAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCa < 1:353923/141‑1 (MQ=255) ttttatacccGCAAA‑CACCCG‑CCCTTTCCCCCACCGCC‑TTGTAA‑AAAAACCCCTCGGCGGCGGTTGTCTGCGATTTTTTTTC‑CCGCCA‑AAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtg < 1:10271/143‑1 (MQ=255) tgctttCCCCCAA‑CCACCGCCCCATTTACCCACCCCC‑ATGAAA‑AAAATTCCCTTCGCGGCGGCGGTCCGCCAGAGTTTTTC‑CCCCCG‑AAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt < 2:375039/148‑1 (MQ=255) tatcccccccAA‑AAACCG‑CTCTTTTCCCCAACCCC‑ATGGAA‑AAAAATGCTTTCGGGGCGGTTGTCGGCCGGTTTTTTTC‑CCGCCC‑CAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttg < 2:572595/143‑1 (MQ=255) cacCC‑CGCCTTTACCCCACCCC‑CTGGGCAAAAAATCCCTCCGCGGCGGTTTTCGGCC‑GGTTTTTTT‑CCCCCGAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc < 1:448186/150‑1 (MQ=255) tttccccAACCCC‑ATGGCC‑AAAAATGCCTCCGGGGGGCCTTTCCGCCAGCTTTTTTC‑CCGGCG‑CAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc < 2:292264/146‑1 (MQ=255) tttaaccaccccctggtgA‑AAAAATCCCTCCGGGGGGGTTGTCTCCCAGCTTTTTTC‑CCGCCA‑CAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc < 2:522571/136‑1 (MQ=255) cccacccacaggagC‑AAAAATCCCTTCGGGGGGGCTGTCTGCCAG‑TATTTTTCCCGCCG‑AAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg < 2:611721/138‑1 (MQ=255) ctccgcgGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCG‑AAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATc < 1:77671/148‑1 (MQ=255) tGTCTGCCGGCTATTTTC‑CCGCCA‑AAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc < 2:648874/150‑1 (MQ=255) gattttttc‑CCGCCC‑AAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt < 2:349581/146‑1 (MQ=255) cTATTTCC‑CCCCCG‑CAAAAAAACCGCCGTTCTCTCCAGGAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTaa < 2:132629/148‑1 (MQ=255) | GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCG‑CGCATTTATCCAACGCC‑ATGTGC‑AAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAT > NC_000913/357742‑358021 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |