Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,890 | 0 | C | A | 54.5% | 9.5 / 44.0 | 33 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base C (10/5); new base A (0/18); total (10/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.24e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.57e-05 |
TAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCG > NC_000913/357744‑358029 | tAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑aaa < 2:639455/150‑1 (MQ=255) aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAAc > 1:352699/1‑151 (MQ=255) cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAAccgc < 2:28351/151‑1 (MQ=255) cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGCGATTCATTTCCGCAAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAAccgc > 2:201639/1‑151 (MQ=255) gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTtc > 2:73675/1‑150 (MQ=255) tCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGTCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGAtctctc > 2:409593/1‑151 (MQ=255) gCAATTCAAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCagt > 2:407856/1‑151 (MQ=255) cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCagt > 2:18290/1‑150 (MQ=255) ggcacggtctcccGTGATTCATCCCCCCCAACACCCGCCTATTTTCCCACCGCCAGGTGCAAAACGGCTTC‑CGGGGCGCCTGTCGGCAAGTTTTTTTCCCGCCCAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAg < 2:585305/140‑1 (MQ=255) cAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAAG‑gaaaaaat > 1:626449/1‑146 (MQ=255) ggCGGCAACCGTAATTCATTCCCCCCACCACCCCCTCTTTCATCCCCCGCCGTGTTAAAAAATGCGTC‑CGCGGCGCCTGTCTCCCA‑TTTTTTTTCCCCCCAAACAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCtt < 2:458505/150‑1 (MQ=255) gCGGCTCCGGTGTTTCTTTCCCGCCAACACCCGCGCATTTACCCACCGCCATGTGAAAAAAGGCCTT‑CGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCt < 1:621502/149‑1 (MQ=255) ggCTCCCGTGTTTCTTTCCCCCCAACACCCGCGCATTTACCCACCCCCATGTAAAAAAATCCCTT‑CGGGGCGGCTGTCGGCCAGTTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑aa < 1:313683/141‑1 (MQ=255) ggCTCCCGTGATTCTTTCCCGCCAACACCCGCGCATTTACCCACCGCCATGTGAAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCttt < 1:611007/149‑1 (MQ=255) aTTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAAG‑aaaaaagttttaac > 1:799094/1‑132 (MQ=25) ttccccccaccaCCCGCGCTTTTACCAACCCCCTTGTAAAAAAAGGCCTC‑CGCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtgtt < 1:784427/142‑1 (MQ=255) ttcccGCCAACACCCGCGATTTTCCCCAACGCCATGGAAAAAAATGCCTG‑CGGGGGGGTTTCCGCCCAGTTTTTTTCCCGCCAAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTc < 2:526978/140‑1 (MQ=255) aacaccCGCTCTTTTCCCCAACGCCATGTAAAAAAATGCCTC‑GGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCAAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc < 2:440972/151‑1 (MQ=255) caacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑GAAAAAATTTTAACCCGTCATGTTTTGTggggat > 1:541427/1‑144 (MQ=25) ccccgcTTTTACCCACCGCTTTGTGAAAAAATGCCTG‑CGCGGCGGCTGTCTGCAAGTTTTTTTCCCGCCGAAAC‑AAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg < 2:771707/147‑1 (MQ=255) ccccgcTTTTACCCACCCCCAGGT‑CAAAAAAGCCTTCCGGGGGGGCTGCCTGCCGGTTTTTTTCCCGCCACAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct < 2:736871/146‑1 (MQ=255) ccGCGCTTTTCCCCACCCCCTTGTGAAAAAATCCTTC‑CGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct < 1:350740/149‑1 (MQ=255) cctggtaaaaaaaGGCCTC‑CGGGGGGGCTGCCGGCAATTTTTTTTCCCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAg < 2:795735/143‑1 (MQ=255) attggcaaaaaaTGCCTT‑GGGGGGGGTTGTCTGCAGGTTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa < 1:402524/140‑1 (MQ=255) tgtgCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAACTTTAA‑CCCGCCCAGTTTTGTGGGCATTGCCTCTGGTGAAAAAATGGAAAggg > 1:614400/1‑150 (MQ=255) ttgaaaaatTGCCTG‑CGGGGCGGCTGTCTGCCAGCTTTTTTTCCCCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa < 1:233969/146‑1 (MQ=255) tgcgggggggcTTTCTCCCAGCTATTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt < 1:622457/141‑1 (MQ=255) cgtgtttttttcCCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt < 2:532343/143‑1 (MQ=255) ccattttttttcCCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa < 1:389400/144‑1 (MQ=255) tttttttcCCCCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc < 2:640081/138‑1 (MQ=255) ttcccgccgAAAA‑AAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt < 2:612093/147‑1 (MQ=255) ccgccaCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCg < 2:238745/149‑1 (MQ=255) cccgaAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc < 2:611257/145‑1 (MQ=255) | TAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCG > NC_000913/357744‑358029 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |