| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 357,890 | 0 | C | A | 45.7% | 31.7 / 16.3 | 35 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
| Reads supporting (aligned to +/- strand): ref base C (16/3); new base A (0/16); total (16/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.14e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01 | |||||||||||
GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGG‑CAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC > NC_000913/357742‑358028 | gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac < 1:423317/149‑1 (MQ=255) tAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa > 1:421422/1‑150 (MQ=255) tAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa > 1:589605/1‑150 (MQ=255) aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaaa > 1:679869/1‑150 (MQ=255) cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccg < 1:535179/150‑1 (MQ=255) gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgc < 2:172550/150‑1 (MQ=255) aGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGt > 1:799171/1‑151 (MQ=255) cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc > 1:164129/1‑151 (MQ=255) cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt > 2:508389/1‑150 (MQ=255) cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg > 1:286505/1‑151 (MQ=255) cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg > 2:349044/1‑151 (MQ=255) cccGCAGACTTGTGTCAGCGCCTCCCTTGTTTCTTCCCCCCCCACACCCGCTTTTTTACCCACCGCAATGTAAAAAATTCCCTC‑CGGGGCGTCTTTCGGCCGTTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCg < 1:411626/148‑1 (MQ=255) ccccAGACTTGTGCCGGCGGCTCCCTTGATTCTTTCCCGCAAACACCCGCGCATTTATCCACCCCCTGGTGAAAAATTGCCTT‑CGGGGGGGCTGCCTGCCACCTATTTTCCCGCGGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg < 2:1190860/146‑1 (MQ=255) cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg > 1:643905/1‑148 (MQ=255) gcctggtgtCGGCGCCTCCGGTGATTCTTTCCCGCCAACACCCGCTCTTTTACCCACCCCCATGCAAAAAACTCCCTC‑CGGGGCGGGTGTCGCCCCGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCggg < 1:848170/140‑1 (MQ=255) ccTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg > 1:373149/1‑144 (MQ=255) tgtcggggcctcccGTGTTTCTTTCCCGCAAACACCCGCATATTTCTCCACCCCCTTGTAAAAAAATCCTTT‑CGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGG‑CAAGA‑Gaa < 2:722801/136‑1 (MQ=255) aGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggaaaggaaaaaaa > 1:779391/1‑136 (MQ=255) cGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGG‑GGAGA‑Gaaaaaaatttt > 2:153221/1‑143 (MQ=255) gcctcccGTAATTCATTCCCCCAAACACCCGCGCATTTACCCAACGCCATGT‑CAAAAACTCCTTCCGCGGGGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGG‑CAAGA‑GAAATAg < 2:1127902/140‑1 (MQ=255) cccccAACACCCGCCCTTTTCCCCACCGCCTTGTAAAAAAATGCTTC‑CGGGGCGGCGTTCTGCCAGCTTTTTTCCCGCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGG‑CAAGA‑GAAATAGCTTTAAGCCCGTCATGTtgtgtg < 2:509746/147‑1 (MQ=255) aacacccccccttttACCCAACGCCTGGAAAAAAAATCCCTC‑GGCGGCGGCTTTCTGCAGGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGG‑CAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc < 2:936221/139‑1 (MQ=255) aacaacCGCGCTTTTATCCACCGCTATGTGAAAAAATGCCTG‑CGCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGG‑CAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc < 2:896495/151‑1 (MQ=255) aacaacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGTCAGGG‑GAAAAAATTTTAATCCCGTCGttttgttggg > 1:619248/1‑143 (MQ=25) cgaaTTTACCCACCGCCATGTGAAAAAATGCCTT‑CGCGGGGCCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGG‑CAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTgg < 2:236284/148‑1 (MQ=255) cATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGG‑GGAGG‑GAAAAAAATTTATAACCCTCCTTTTTTGTGGGGATCgctgggggg > 1:1216039/1‑145 (MQ=21) aacgccgtgacaaaaaaTGCTGG‑CGCGGGGGCAGGCGCCCGGCTTTTTTCCCCCCAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGG‑CAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg < 1:1229565/136‑1 (MQ=255) aaCGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGG‑GGAGA‑GAAACAAATCTAAACCCCGCATTTTTTGTGGGCGTCGCTGCTggggaaaaaatg > 1:598594/1‑140 (MQ=25) cccatggaaaaaaaaTCCCTC‑CGCGGGGGTTGTCGGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGG‑CAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACa < 1:634640/133‑1 (MQ=255) cggtgtgAAAAAATGCCTG‑CGCGGGGTGTTTCGCCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGG‑CAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAg < 2:1246767/148‑1 (MQ=255) aTGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGG‑GGAGAGGAAAAAATTTTAAACCC‑CCCTGTTTTGTGGGAAACGCCGCTGGTTAAAAAATGGGAAgg > 2:506544/1‑150 (MQ=255) ggttttcgcccAGCTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGG‑CAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAg < 1:341403/135‑1 (MQ=255) tgcttgcaattttttttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGG‑CAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa < 1:1130418/127‑1 (MQ=255) ccattttttttcCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGG‑CAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCgcgc < 1:1127492/142‑1 (MQ=255) ttttttcCCCCCGCAAAAAAACCGCCGTTGTCTCCAGTAGTGGCGGGG‑CAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg < 2:695639/148‑1 (MQ=255) ccgccaAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGG‑CAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc < 2:1117558/148‑1 (MQ=255) | GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGG‑CAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC > 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 < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |