Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,875 | 0 | A | T | 50.0% | 42.1 / 26.2 | 43 | N272K (AAT→AAA) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base A (15/6); new base T (0/21); total (15/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.10e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTA > NC_000913/357728‑358019 | gAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTAttt > 2:708987/1‑151 (MQ=255) gAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTAttt < 2:306798/151‑1 (MQ=255) cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccg < 2:815440/151‑1 (MQ=255) gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccgc < 1:1081563/151‑1 (MQ=255) ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCcgccgc > 2:1027438/1‑151 (MQ=255) tAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa < 1:832153/149‑1 (MQ=255) gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑ccgcc < 2:262286/151‑1 (MQ=255) tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑ccgccg > 1:1045077/1‑151 (MQ=255) aGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGAGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGt > 1:985238/1‑151 (MQ=255) aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTtctct > 1:1102146/1‑151 (MQ=255) aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTcc > 1:71640/1‑150 (MQ=255) ggCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTcc > 1:244591/1‑149 (MQ=255) cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCagt > 1:200823/1‑150 (MQ=255) ccAGCCGCAGACTTGTGTCAGCGGCTCCCGTATTTCATTCCCGCAAACACCCGCGCATTTACCCACCGCAGTGTGCAAAAATGCTTC‑GGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCagt < 2:529796/145‑1 (MQ=255) ccGTAGACCTGGGTCGGGGCCTCCCGTGTTTCTTTCCCCCCAACCCCCGCGTTTTTCCCCACCCCCATGTAAAAAATTCCCTC‑CGGGGCGGCTGTCTCCCATTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCggg < 1:714356/150‑1 (MQ=255) ccGCGGCCCTGTGTCAGCGGCTCCCGTATTTCATCCCCCCCACCACCCGCGCATTTACCCACCGCCTTGGCCAAAATTCCCTC‑CGGGGGGGTTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCgg < 1:830316/149‑1 (MQ=255) ccGCAGCCCTGGTCCGGCGGCTACCGTATTTCTTTCCCGCCAACACCCGCGCTTTTACCCAACCCATTGTGAAAAAATCCCTT‑CGGGGTGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCggg < 2:1062208/150‑1 (MQ=255) ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCggtg > 1:384474/1‑149 (MQ=255) tgtgTCAGCGCTTCCCGTGTTTCTTTCCCCCCAACACCCGCGCTTTTACCCACCCCAATGTGAAAAAAGGCCTC‑CGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgag‑a < 1:761859/149‑1 (MQ=255) tgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgggggccaaaaaaa > 2:1063213/1‑140 (MQ=255) tgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑aaa > 1:401778/1‑149 (MQ=255) cGGCTCCGGTTATTCATTCCCCCCAACACCCGCGCTTTTACCCACCGCCATGTCAAAAAATCCCTT‑CGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag < 1:576557/140‑1 (MQ=255) gcctcccGTGTTTCATTCCCGCCAACACCCGCGCTTTTACCCACCGCCATGTCAAAAAATCCTTC‑CGCGGGGGCTGTCTGCCAGCTTTTTTCCCCCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑aaa < 2:1097357/137‑1 (MQ=255) ccgttttccATTCCCCCCACCACCCCCGCATTTACCAACCCCCATGTAAAAAAGTGCCTC‑CGCGGCGGCTGTCTCCCACTTTTTTTCCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGc < 1:809948/140‑1 (MQ=255) tgtttaatccccGCCAACACCCGGCCTTTTACCCAACGCCTTGTGCAAAATTGCTTC‑CGCGGGGGCTGTCGGCCAGCTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGt < 2:659308/139‑1 (MQ=255) ttCTTTCCCGCCAACACCCGCGCTTTTATCCAACGCCATGTGAAAAAATGCTTT‑CGCGGCGGCTGTCTGCCAGTTTTTTTTCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtg < 2:881029/150‑1 (MQ=255) ttCATTCCCCCCAACACCCGCGCATTTACCCAACGCCATGTGAAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtg < 2:184412/150‑1 (MQ=255) ttccccccaccaCCCGCCCTTTTACCCACCGCAATGTAAAAAATTGCTTC‑CGGGGCGGTTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtgt < 1:281628/141‑1 (MQ=255) ccccccaccaCCCGCCCATTTCCCCAACGCAATGTAC‑AAAATCCCTTCCGCGGCGGCTGTCTCCCACTTCTTTTCCCGCCGCAACAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtg < 1:397683/140‑1 (MQ=255) acccGCCCTTTTACCCACCGCCTTGTCAAAAAAGGCCTC‑CGCGGCGGTTGTCTGCCAGTTTTTTTCCCCCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCg < 2:968768/147‑1 (MQ=255) tttccccacccccATGTGAAAAAATGCCTC‑CGCGGCGGCTGTCTGCCAGCTTTTTTCCCCCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTa < 2:181612/138‑1 (MQ=255) gtgAAAAAAGCCCTC‑CGCGGGGGCTTTCTGCCAGCTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAGGGa < 2:1045220/150‑1 (MQ=255) gtcaaaaaaTGCCTC‑GGGGGCGGTTTTCTGCCAGTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAggg < 1:434130/145‑1 (MQ=255) gCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑AAAAAAATTTAAAACCCCCATGTTGTGTGGGCATAGCTGCT‑GGGAaaaaaatgggaa > 2:773690/1‑136 (MQ=255) acAAAAATCCCTG‑CGCGGCGTTTCTCCCCGAGTTTTTTTCCCCCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAGGGaa < 1:708247/148‑1 (MQ=255) aaaaaTGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGG‑AAAAGGCTTTTAACCCCCCCCTTTTTTTGGCAATTACTGTG‑GGGAAAAATGTGGggaagggaagg > 2:57288/1‑149 (MQ=255) aaaaaTGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAAATTTAAACCC‑CCCTCTTTTGTGGGGATCGCTGCTGGGTAAAAAGTG‑GGAAAGGGagagg > 2:734337/1‑147 (MQ=255) aaaaaTGCCTC‑GGCGGCGCCTTTCTGCCAGTTTTTTTCCCCCCGAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAGGGAAGt < 2:33307/149‑1 (MQ=255) ggcggcTTTCT‑CCAATTTTTTTTCCCGCACAAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGc < 1:1032896/148‑1 (MQ=255) ggcggcTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCt < 2:1004887/149‑1 (MQ=255) tGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGTGCGAGAG‑AAAAAACTTTTAACCCCTCATTTTTTGTGGGGATCGCTGCT‑GGGGACAAAGT‑GGGAAGGGGAGGGGGGGGGAACAGCCCCCGaa > 1:536154/1‑149 (MQ=255) cTGCCATTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCa < 1:204621/143‑1 (MQ=255) ccAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa < 2:1098103/151‑1 (MQ=255) | GAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTA > NC_000913/357728‑358019 |
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: ‑ |