Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,890 | 0 | C | A | 41.2% | 17.5 / 17.2 | 34 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base C (14/6); new base A (0/14); total (14/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.05e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCT > NC_000913/357742‑358032 | gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca < 1:635849/150‑1 (MQ=255) aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa < 1:24643/151‑1 (MQ=255) aggCGTAGCGCATCAGGCAATTCCGGCCGCGGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CACCCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac < 2:362842/144‑1 (MQ=255) gTAGCGCATCAGGCAATTCCAGCGGCAGACCTGTGTCAGCGGCGACCGTAATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAGGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑ccgcc < 1:857634/151‑1 (MQ=255) aGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑ccgccg > 2:843443/1‑150 (MQ=255) cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTcc > 2:1146354/1‑151 (MQ=255) cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCagt > 2:808767/1‑150 (MQ=255) tCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTgg > 2:755213/1‑151 (MQ=255) ccgcccccggtcttgtcTCAGCGCCTCCCGTAATTCTTTCCCCCCAA‑CACCCGGTCTTTTACCCAACCCCATGTGAAAAAATGCCTGCGCGGCGGTTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTg < 1:298750/136‑1 (MQ=255) aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgg > 1:252279/1‑151 (MQ=255) ctgatcTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAA‑CAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTATTTTTTAGGCGCAACAAAA‑CCGCCGCTCCCTCCactagcgacatcg > 2:298741/6‑140 (MQ=11) gtCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑gaa > 2:1158957/1‑148 (MQ=255) aGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAagcggaaga > 2:604256/1‑140 (MQ=255) gCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑aaaaaaactt > 1:59811/1‑141 (MQ=255) cggtttttCTTTCCCCCCAA‑CACCCGCGTTTTTACCCACCCCAATGTAAAAAAAGGCCGCGGGGGGGGCGGTCCGCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCg < 2:82945/145‑1 (MQ=255) tGATTC‑TTTCCCCCAACCACCCCCGCTTTTCCCCACCGCCAGGTAAAAAAATCCTTTCGCGGCGGCTGCCGGCAACTTTTTTTCCCCCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAg < 1:223726/144‑1 (MQ=255) cGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGACAAG‑AAAAAACATTAAgtaccccctgtt > 2:801461/1‑133 (MQ=25) caacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCGCGAG‑AAAAAAAATTAAAAACCTCATCTTGTGTggggatggg > 2:78091/1‑144 (MQ=25) ccGCGCTTTTCCCCAACCCATTGTAAAAAAATCCCTTGGCGGCGGCTGTCGGCCGGTTTTTTTCCCCCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct < 2:927817/149‑1 (MQ=255) tttACCCACCGCTATGTCAAAAAAGCCCTGCGCGGGGGCTGCTTGCGACTTTTTTTCCCGCCGAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTaa < 2:356415/150‑1 (MQ=255) aaCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGCACAG‑AAATAAACTTAAACCCCCCCCGTTTTTTGGGCATTGCTGCCGGGAACAAta > 2:473687/1‑147 (MQ=21) aaacggctttCGCGGCGGTTTGCGGCCATGTTTTTTCCCCCCGCAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcg < 1:439165/143‑1 (MQ=255) aaaaTGCCTCGGGGGCGTCTGTCTGCAAGTTTTTTTCCCCCCGCAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcg < 2:43814/151‑1 (MQ=255) cTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAACCTTTAAACCC‑CCCTTTTTTGTGGGAATAGCCGCTGGTGACAAAGTGGGAAGGGGAgggggggg > 1:156939/1‑140 (MQ=255) ggggctttccgcccAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAg < 2:181523/135‑1 (MQ=255) ggcggcTTTCTGCAATTTTTTTTCCCGCCCAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCt < 1:160803/149‑1 (MQ=255) ggggttGTCTGCCCGTTTTTTTCCCGCCGAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa < 1:329074/140‑1 (MQ=255) gcggtTGTTTGCCATTTATTTTGCCGCCCCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTcc < 1:500030/150‑1 (MQ=255) atgtttttttcCCCCCACAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgcg < 2:982734/141‑1 (MQ=255) gtttttttcCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt < 2:1121203/143‑1 (MQ=255) tATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAATTTTAAACCCG‑CCTGTTTTGTGGGGATAGCCGTTGGGGACAAAGTGGGAAGGGGAGAGGGGCGGATGAGCTCCAGaaaccccccaatt > 1:808056/1‑138 (MQ=255) tttttCCCCCCAAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt < 1:1055638/147‑1 (MQ=255) ttttCCCCCCCAAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg < 1:717254/150‑1 (MQ=255) ccgccaCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc < 1:976489/148‑1 (MQ=255) cccccgaAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGcc < 1:422268/148‑1 (MQ=255) ccgcAAAAACG‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCt < 2:252271/149‑1 (MQ=255) | GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCT > NC_000913/357742‑358032 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |