| 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% | 28.0 / 11.7 | 31 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
| Reads supporting (aligned to +/- strand): ref base C (12/6); new base A (0/13); total (12/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.42e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTA > NC_000913/357742‑358019 | gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGcaacaa < 1:303066/151‑1 (MQ=255) aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACaaaa < 2:79035/150‑1 (MQ=255) cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc > 2:87042/1‑151 (MQ=255) aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc > 1:455056/1‑150 (MQ=255) aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt > 1:137844/1‑149 (MQ=255) aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg > 2:487203/1‑151 (MQ=255) gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg > 2:239080/1‑149 (MQ=255) aGACCTGTGTCAGCGGCTCCCGTATTTCATTCCCCCCAACACCCGCGCTTTTACCCACCCCCTTGTAAAAAAA‑TCCCTC‑CGGGGCGG‑CTGTCTGCCAGCTATTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag < 1:155782/151‑1 (MQ=255) gACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg > 1:383414/1‑146 (MQ=255) gtgcCAGCGGCACCCGTGTTTCATTCCCCCAAACACCCGCCCTTTTACCCACCGCC‑TGGTCAAAAA‑AGCCTTCCGCGGGGC‑CTGCCTGCCGGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCggg < 2:224921/140‑1 (MQ=255) gtCAGCGGCTCCCGTGTTTCTTTCCCCCCAACACCCGCGCTTTTACCCACCGCCATGTGAAAAAA‑TGCCTC‑CGGGGCGG‑CTGTCTCCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAg < 1:265079/151‑1 (MQ=255) cggcgcCTCCCGTGATTCTTTCCCCCCAACACCCGCGCTTTTACCCACCGCAATGTAAAAAAA‑TGCTTC‑CGCGGGGG‑CTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATa < 1:87063/146‑1 (MQ=255) tACCGTGATTATTTTCCGGGAACACCCGGCCCTTTATCCACCGCTATCCGCAAAAA‑TTCATC‑CGGGGCGG‑CTGCCAGCCGGTTATCTTCCTGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTaa < 1:521961/148‑1 (MQ=255) ccGTGTTTTATCCCCGCCAACACCCGCGCATTTTTCCACCGCCTTGTGCAAAAA‑TGCTTC‑CGCGGCGG‑CTGTCTGCCAGTTTTTTTTCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGc < 1:198309/148‑1 (MQ=255) ttaTTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGA‑GGAAAAAATTTaaaacccccc > 1:43490/3‑143 (MQ=25) aTTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGA‑Gaaaaaaattttaagcc > 1:219925/1‑132 (MQ=25) ttCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGGAAAAAGATTTAAAccccgc > 1:396416/1‑145 (MQ=25) acccGCGCTTTTA‑CCCACCCCATTTTCAAAAATTGCCTC‑CG‑GGCGGCCTGTCTGCCATTTTTTTTCCCGCCCCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct < 1:596425/149‑1 (MQ=255) acccGCGCATTTACCCACCGCCATGTGAAAAAA‑TGCCTC‑CGCGGCGG‑CTGTCTCCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc < 1:330549/148‑1 (MQ=255) ccGCGCTTTTACCCAACGCTTTGGGAAAAAA‑TGCCTT‑GGCGGCGT‑CTGTCGGCCAGTTATTTTCCCGCGAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc < 2:298847/146‑1 (MQ=255) ccATTTATCAAGCCCAATGTGAAAAAA‑TCCCTT‑CGCGGCGC‑CTGTTGGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGt < 1:400833/150‑1 (MQ=255) ttttcccAACCCCTTTTGAAAAAA‑TGCCTC‑CGCGGCGG‑CTGTCTGCCACTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgct < 1:314820/140‑1 (MQ=255) tttACCCACCGCCGTTTACAAAAA‑TGCCTG‑CGCGGCGG‑CTGTCTCCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgct < 1:352957/145‑1 (MQ=255) acccaccgcattgtgCAAAAC‑TCCCTG‑CGGGGCGG‑CTGTCTCCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACaa < 2:24637/139‑1 (MQ=255) aTCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGA‑GAAAAAACTTTTAAACCCCCATGTTGTTTGGGGATCCCTGCTGGTGACaa > 1:319332/1‑150 (MQ=25) ccacccccttgtaaaaaaa‑TGCCTC‑CGCGGCGG‑CTGTCTCCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt < 1:245710/135‑1 (MQ=255) aaaaa‑TGCCTG‑CGGGGCGT‑CTGTCTGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa < 1:409438/146‑1 (MQ=255) ccTT‑CGCGGCGC‑TTGTCTCCCAGCTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt < 2:592736/151‑1 (MQ=255) ggcgg‑cTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAA‑GAAAAAAATTTAAACCCGCCCTGTTTTGTGGGGATCGCTGCTGGGTACAATGTGGGAAGGGGAGGGGGGGGAAGCAGcccc > 1:551153/1‑151 (MQ=255) gcgg‑cTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑GAAAAAAATTTTAACCCCCCATGTTGTGTGGGGATCGCTGCTGGGGACAATGTGGGAAGGGGAGAGGGGGGGATCaccccc > 1:615198/1‑145 (MQ=255) g‑cTGTCGGCCAGCTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa < 1:193048/151‑1 (MQ=255) ccAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa < 1:419916/151‑1 (MQ=255) | GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTA > NC_000913/357742‑358019 |
| 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: ‑ |