Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 803,869 | 0 | G | T | 35.3% | 37.4 / 14.0 | 34 | V123L (GTG→TTG) | ybhJ | putative hydratase YbhJ |
Reads supporting (aligned to +/- strand): ref base G (2/20); new base T (12/0); total (14/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.66e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGC‑ATATTGCGGTCATCC‑ATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCC‑G‑TTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGATGTGGC > NC_000913/803628‑804113 | ctctTAGCATAATACGTCTTGAAATATGGATCGACTTAAAGTTAAGTTTGATTCATTAGCCTCTCACGCTATTACCTTTATTGGTCTTGTACAGACCGCTAAAGCGTCCGGTATGGATCGTTTCCCGCTGCCGCATGTGCTGACCAACTGCCCTAACTCACTCTGCGCCGTCGGCGGAACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGAGGCCCAGCGTTATGGCGG‑CGTTTTtgttgttctcc > 2:469112/1‑243 (MQ=255) gcGCAAAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTGATGGCGG‑GATTTTGGTGCCTCCGC‑at < 2:214090/248‑1 (MQ=255) cGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTAACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGC‑ATAtt < 2:15709/251‑1 (MQ=255) ataCGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGC‑ATATTGCGGTc < 1:442747/251‑1 (MQ=255) cGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGC‑ATATTGCGGTcatc < 1:182123/251‑1 (MQ=255) ccGGAAATATTGATAAACGTAAAATTAAGATTGATTCATTAGCCTTTCCCGATATTCCCTTTGTTGGTATGGGACAGACCGCTAAAGCGCCCGGTATGTAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGGTGAGCCGCGGCCCAGCGTTATGGCGG‑TATTTCTGTTCCTCCGC‑ATATTGCGGTcatcc‑at < 2:459473/251‑1 (MQ=255) cGGAAAGATGGAAAAACTGAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGCCCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCCCTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGC‑ATATTGCGGTcatcc‑atc < 2:232618/251‑1 (MQ=255) atGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑GATTTTTTTGCCTCCCC‑ATATTGCGGTCATCC‑CTCAatatat > 1:452934/1‑251 (MQ=255) ggATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCTGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGC‑ATATTGCGGTTATCC‑ATCAATATATGc < 1:92816/251‑1 (MQ=255) aTTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑GATTTTTTTGCCTCCCC‑ATATTTCGGGCATCC‑ATCCATATATTCGTggggtggtggg > 1:339497/1‑241 (MQ=255) gTTTGATTCATTAGCCTCTCCCGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TGTTTTTTTGCCTCCGC‑ACATTTCGGGTAACC‑ATCAataaattcggggggtgatggaagggg > 1:115796/1‑228 (MQ=255) tGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGCTATGGCGG‑GATTTTTTTGCCTCCCC‑CTATTTCGGCCATCC‑ATCAatatattggtgtggtgattggaggagggg > 1:464424/1‑228 (MQ=255) tatTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTTTGCCTCCCC‑ATATTTCGGGCCTCC‑CTCAAAATATGCGGGGGGTGATGGCGGGCGGCGGGAAAATTATCCTccgggc > 1:373614/1‑244 (MQ=255) gttgGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGCCCACCTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGC‑ATATTGCGGTCATCC‑ATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCacac < 2:276375/251‑1 (MQ=255) ttGTACAGACCGCTAAAGGGCCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCGGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGC‑ATATTGCGGTCATCC‑ATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCC‑G‑TTa < 2:263868/250‑1 (MQ=255) cAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑GATTTTTTTTCCTCCCC‑ATATTTCGGTCATCC‑ATCAATATATTCGTGGGGTGATTGGAGGCGGCGGCAAAATTATCCCCGGGGCAGACAACCACcccccctacggggc > 1:174667/1‑241 (MQ=255) gACCGCAAAGGCGTCCGTTATGGACCGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGC‑ATATTGCGGTCATCC‑ATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCC‑G‑TTACGGTGCAt < 2:371521/251‑1 (MQ=255) cGGTATGGACCGTTTCCCGCTGCCTTATGTGCTGACCAGCTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGC‑ATATTGCGGTCATCC‑ATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCC‑G‑TTACGGTTCTTTAGGGACAATGGCAg < 2:214966/251‑1 (MQ=255) cGGTATGGAACTTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGC‑ATATTGCGGTCATCC‑ATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCC‑G‑TTACGGTGCATTAGGGACAATGGCAg < 2:452934/251‑1 (MQ=255) tGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGAGGCCATATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGTTTTTTTTTGCCCCCCC‑A‑ATTTCGGGCATCCCATCA‑TATATGCCTGAGATGAAGGAAGGGGGCGGCAAAAAGATCCTCGGGGCAGGCAAACAAACCC‑C‑TTACCGGGCATTAGGGAGAATgggcggcggggaggggggcgggggg > 1:305272/1‑228 (MQ=255) ccGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGC‑ATATTGCGGTCATCC‑ATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCC‑G‑TGACGGTGCATTAGGGACAATTGCAGTCGGTGAGGGCGGCGGTGAGtt < 1:390567/251‑1 (MQ=255) cGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGC‑ATATTGCGGTCATCC‑ATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCC‑G‑TTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTg < 1:161314/251‑1 (MQ=255) gTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTTTGCCTCCGC‑ATATTGCGGGCATCC‑ATCAATATATTCGGGGGATGTTGGGAGGGGGCGGCAAAATGATCCTCCGGGCAGAAAACCCCAACC‑G‑GTACGGGGGCTTTGGGGCAATGGCAGTCGGGGGGGGgggcgggggggtgg > 1:138337/1‑243 (MQ=255) gtgCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTTTGGCTCCGC‑CTATTTCGGGCATCC‑CTCAATATATGCGTGAGAAGGTGGCAGGGGGCGGCCAAAAGATTCCCCGGCCAGGCAGCCACACCC‑GCTAACGG‑GCATTTGGGGCAAAGGCAGGCGGGGAggggggggggggggttggga > 1:50546/1‑233 (MQ=255) ccATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGC‑ATATTGCGGTCATCC‑ATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCC‑G‑TTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATg < 1:150161/251‑1 (MQ=255) aaCTCCCTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGC‑ATATTGCGGTCATCC‑ATCAATATATGCTGGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCC‑G‑TTACGGTGCATGAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAGCAGCTGCTTAAGGacac < 2:147016/251‑1 (MQ=255) aaCTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGC‑ATATTGCGGTCATCC‑ATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCC‑G‑TTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGaca < 2:339497/250‑1 (MQ=255) cactcTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGC‑ATATTGCGGTCATCC‑ATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCC‑G‑TTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCt < 1:346782/249‑1 (MQ=255) ggcggcACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGGCTCCGC‑ATATTTCGGTCATCC‑ATCAATATATTGGTGGGATGATGGGAGGCGGGGGCAAAAAGATCCTCGGGGCAGAAAGCCACACCC‑G‑TTACGGGGCATTTGGGACAAAGGGAGTCGGGGGGGGCGGGGGGGGGTTGGGAAAAAAGCTGGTTTATGAAACCTGGGGTATCggatatcc > 1:92378/1‑244 (MQ=255) gTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTTTGCCCCCCC‑CTATTTCGGGCCTCC‑CTCAATATATGCGGGAGATGATGGGAGGCGGCGGCAAAATGATCCCCGGGGCAAAAAGCCCCCCCCCG‑TT‑CGGGGCCTTTGGGACAATGGCCGGCGGGGGGGGGGGCGGGGGGTTGGTAAAACAGCTGCTTTATGACACCTGGGGGATCGACTATCCGGGCGTGGGTGCgggg > 1:404571/1‑249 (MQ=255) ttttGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTTTGCCTCCGCCAT‑TTGCGGTCCTCC‑CTCAATATATGCGTGAGATGCTGGGCGGCGGCGGGAAAATTATCCCCGGGGCAGACAACCCCCCCC‑G‑CTACGGGGGATTTGGGGCACAGGCAGGCGGGGAGGGGGGCGGTGAGTTGGGAAAAAAGCTGCTTTATGACACCTGGGGGATCGACTTTCCCGGCGTGGTTgggggggctcttaccgga > 1:287899/1‑233 (MQ=255) tATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTTTGCCTCCCC‑ACATTGCGGGCATCC‑ATCAATAAATGCGTGGgatgat > 1:182280/1‑79 (MQ=255) tATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGC‑ATATTGCGGTCATCC‑ATCAATATATGCGTGAgatgat < 2:182280/79‑1 (MQ=255) ggcCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGC‑ATATTGCGGGCATCC‑ATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCCCACCC‑G‑TTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGGGAGGTGGGAAAACAGCTGCTTAATGACACCTGGGGTATCGACTATCCGGGCGTGGTTGCGGGGCATCTGACCGGGAAACCTGCGCCGTAtt > 1:177308/1‑250 (MQ=255) aTTTTTGTGCCTCCGC‑ATATTGCGGTCATCC‑ATCAATATATGCGTGAGAGGATGGCAGGCGGCGGCAAAATGATCCTCGGGGCAGACAGCCACACCC‑G‑TTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTTAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTCAGCGGAAAACCTGCGCCGTATGTGGGGCCACAGGATGTGgc < 1:298216/251‑1 (MQ=255) | CTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGC‑ATATTGCGGTCATCC‑ATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCC‑G‑TTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGATGTGGC > NC_000913/803628‑804113 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |