Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,890 | 0 | C | A | 58.6% | ‑3.4 / 25.3 | 29 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base C (9/3); new base A (0/17); total (9/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.20e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.37e-02 |
TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCA‑TTTATCCAACGCCATGTGCAAAAA‑TG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCG > NC_000913/357741‑358029 | tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCA‑TTTATCCAACGCCATGTGCAAAAA‑TG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca > 1:757110/1‑151 (MQ=255) gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCG‑CCAACAACCGCGCA‑TTCCTCCAACGCCATGTGCAAAAA‑TG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa > 2:62887/1‑151 (MQ=255) gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCA‑TTTATCCAACGCCATGTGCAAAAA‑TG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa > 2:614164/1‑151 (MQ=255) gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCA‑TTTATCCAACGCCATGTGCAAAAA‑TG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa > 1:298881/1‑151 (MQ=255) gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCA‑TTTATCCAACGCCATGTGCAAAAA‑TG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa > 2:460670/1‑151 (MQ=255) aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCA‑TTTATCCAACGCCATGTGCAAAAA‑TG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc > 1:675869/1‑150 (MQ=255) gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCA‑TTTATCCAACGCCATGTGCAAAAA‑TG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccg > 1:562583/1‑151 (MQ=255) gCGTAGCGCACCAGGCAATTCCAGCCGCTGACCTGTGTCAGCGGCCACCGTGATTCATTTCCG‑CCATCAACCGCGCA‑TTTATCCAACGCCATGTGCAAAAA‑TG‑CCTTCGCGGCGGCTGTCTGCCATCTATTTTTCCGCCGCCACAAAAccg > 1:34956/1‑151 (MQ=17) aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCA‑TTTATCCAACGCCATGTGCAAAAA‑TG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCa > 1:962407/1‑151 (MQ=255) ggCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCA‑TTTATCCAACGCCATGTGCAAAAA‑TG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcccg > 2:248731/1‑149 (MQ=255) cacccccgttccTGTGTCAGCGTCTCCCGTGTTTCTTACCCG‑CCAACACCCGCACT‑TTCACCCAACCCCCTGTGAAAAAA‑GC‑CCTGCGGGGCGGCTTTCCCCCATTTTTTTTCCCGCCGCAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCg < 1:866704/141‑1 (MQ=255) cacccGCGGCCCGGTGTCCGCGCCTCCGGTGTTTCTTCCCCC‑CCACCACCCGCCTT‑TTTCTCCACCCCCATGTGAAAAAT‑TG‑CCTCCGGGGCGCTTGTCTGCCGGTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCg < 1:468831/148‑1 (MQ=255) gCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCA‑TTTATCCAACGCCATGTGCAAAAA‑TG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg > 1:169523/1‑151 (MQ=255) ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCA‑TTTATCCAACGCCATGTGCAAAAA‑TG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg > 1:563108/1‑151 (MQ=255) gACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCA‑TTTATCCAACGCCATGTGCAAAAA‑TG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGagg > 2:231101/1‑150 (MQ=255) cctggtgtCAGCGCCTCCCGTGATTCTTCCCCC‑CCAAAACCCCCGTT‑TTTACCCACCGCCATGTCAAAAAA‑TC‑CCTCGGGGGGGGCTGTCGCCAAGTTATTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag < 2:544182/147‑1 (MQ=255) cccTGTGACAGCGGTTCCCTTTTTCCATTCCCG‑CCAACCACCGCTCT‑TTTCCCCAACCCCATGTAAAAAAA‑TG‑CCTCGGGGCCGGCGTCCGGCCAGTTTTTTTCCCGCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag < 1:339742/150‑1 (MQ=255) gtgtCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCA‑TTTATCCAACGCCATGTGCAAAAA‑TG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggcaggggaaa > 1:882633/1‑141 (MQ=255) tCAGGGGCTCCAGGGTTTCA‑TTCCC‑CCAACCACCGGGGATTTTATCAACCGCAATGTAAAAAAA‑TG‑CCTCCGGGGCGGTTGTCGGCCAGTTTTTTTCCCGCACAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGc < 2:118236/151‑1 (MQ=255) tACCGTGATTCTTTCCCG‑CCACCACCCGCTCT‑TTTACCCACCG‑CATGGTCAAAAA‑AGCCCTCCGCGGCGGTTGTCGGCAGGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc < 1:923256/151‑1 (MQ=255) aacgtggatCAATTCCGCCCAACACCCGCGC‑‑TTTAACCCACCGCATGGGCAAAAATTG‑CCTGGGGGGCGGTTGTCCGCCATTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc < 2:1039915/143‑1 (MQ=255) ttcttccccccaaacaccCGCGCT‑TTTACCCACCGCCATGTGCAAAAA‑TG‑CCTCCGGGGCGGCTTTCGGCCGGCTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt < 2:124816/139‑1 (MQ=255) ccaccccccGCGCA‑TTTACCAAGCGCTTTGTGAAAAAA‑TC‑CCTCCGCGGGGCTTGTCTGCCATTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCaa < 1:972837/144‑1 (MQ=255) tttACCCAACGCCATGTGCAAAAA‑TG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc < 1:108282/151‑1 (MQ=255) tttACCCAACCCCATTTCAAAAAA‑TC‑CCTCCGGGGCGGTTGTCTGCCAGCTTTTTCCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc < 2:334695/151‑1 (MQ=255) acccAACGCCATGTAAAAAAA‑TG‑CTTTCGCGGCGGCGGTCGGCCGGATATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt < 2:639857/149‑1 (MQ=255) tCCAACCCCTTGGGAAAAAA‑TC‑CCTTCGGGGCGGCTGTCGCCCAGTTTTTTTCCCGCCGAAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg < 2:563347/151‑1 (MQ=255) ccAACGCCATGGCAAAAAA‑TG‑CCTCCGGGGCGGCTTCCTGCCGGTTATTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt < 1:62879/151‑1 (MQ=255) ctctCTGCCAGTTTTTTCCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc < 1:66000/148‑1 (MQ=255) cgccgtcttttttcCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCgcgc < 1:47368/142‑1 (MQ=255) cggtcTATTTTCCGGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCgcgc < 1:232953/145‑1 (MQ=255) tttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt < 1:854892/151‑1 (MQ=255) ttccccgccgAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt < 1:75374/147‑1 (MQ=255) tcccgccaCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCg < 1:79643/149‑1 (MQ=255) | TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCA‑TTTATCCAACGCCATGTGCAAAAA‑TG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCG > NC_000913/357741‑358029 |
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: ‑ |