Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,510,395 | 0 | C | A | 40.0% | 89.8 / 32.8 | 56 | S131R (AGC→AGA) | ydcR | putative DNA‑binding transcriptional regulator and putative aminotransferase |
Reads supporting (aligned to +/- strand): ref base C (28/5); new base A (0/22); total (28/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.21e-11 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CGAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTAGAGGCATTA > NC_000913/1510249‑1510542 | cGAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCaaaa < 1:99567/151‑1 (MQ=255) gagaTGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑c > 2:469308/1‑151 (MQ=255) gagaTGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑c < 2:953835/151‑1 (MQ=255) agaTGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑cc > 1:707828/1‑151 (MQ=255) aTGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGc < 1:305398/151‑1 (MQ=255) aaGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCcac > 1:680238/1‑151 (MQ=255) gCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCcac > 2:757277/1‑149 (MQ=255) aGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCg > 1:38617/1‑151 (MQ=255) atCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCg > 1:740570/1‑151 (MQ=255) cAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTg > 2:855929/1‑151 (MQ=255) cAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTC‑CCCGA‑CCCGCGACTTTCCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGAAAAG‑CCGCCACAGCGATGAGCg < 1:137881/149‑1 (MQ=255) acTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTg > 2:130434/1‑151 (MQ=255) ttttttttttGATTTGTTGCGGCCCAGTCGCGACCCGTGGGTC‑TTCCCTTTTCCCCCGCCTTTCCCCGC‑CCCGCTATTTTCCCCCCCCC‑AAAAAAAAACCGGCTCGCGGGCGGGGGTAAAAAAAA‑CCGCCACAGCGATGAGCGt < 1:601257/139‑1 (MQ=255) ttATTTTTTTGATAGGTGGCAGCCCGGCCGCGCCCCGTGTGTCG‑TCCGTTTTCCCCCGCCTTTCCCCGC‑CCCGCGATTTCCCCCCCCCC‑AACAACAAACCCGCTCGCTGGCGAAGGAAAGAAAAA‑CCGCCACAGCGATGAg < 1:886062/141‑1 (MQ=255) ttATATTTTTAATATGTTGCAGGCCAGTCGCGATCCGTCGGTCTTTCCTTTTCCCCCGGCCTTC‑CCCGC‑CCCGCGATTTTCCCCCCCCC‑AACAACAAACCCGCTCGCTGGCGCAGGTAAGAAAAA‑CCGCCACAGCGATGAGCGTGATTg < 2:247929/149‑1 (MQ=255) tttttttGATTTGTTGCGGGCCGGTCGCGACCCGTCGGTCTTTCCTTTTGCCTCGGCTTTC‑CCCGC‑CCCGCGCTTTTCCCCCCCCC‑AACAACAAACCCGCTCGTTGGCGCAGGAAAGAAAAA‑CCGCCACAGCGATGAGCGTGATTg < 1:659192/144‑1 (MQ=255) tttttGATTTGTGGCAGGCCAGGCGCGCCCCGCGGTTCGTTCGGTTTCCCGCGCCTTTC‑CCCGA‑CCCGCG‑CTTTCCCCCCCCCAAACAACAAACCCGCTCGCTGGCGCAGGAAACCAAAT‑CCGCCACAGCGATGAGCGTGATTGaaaa < 2:613704/148‑1 (MQ=255) ttttgGATTTGTGGCGGGCCAGGCGCGTCCCGTCGTTCGTCCCGTTGCCCCCGCCCTTC‑CCCGC‑CCCGCGACTTTCCCCCCTCC‑AACAACAAACCCGCTCGCGGGCGCGGGAAAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGaaaa < 2:382237/148‑1 (MQ=255) ttttGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACtt > 2:943703/1‑150 (MQ=255) ttttGATATGTTGCAGCCCGGCCGCGATCCGGGGGTCTTTCTTTTTCCCCGGCCCTTC‑CCCGC‑CCCGCCTCTTTCCCCCCCCC‑AAAAAAAAACCCCCCCGCTGGCGCAGGTAAGAAAAC‑CCGCCACAGCGATGAGCGTGATTGAAAACtt < 1:316134/150‑1 (MQ=255) atGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGcc > 2:782535/1‑149 (MQ=255) tGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAgg > 2:822317/1‑151 (MQ=255) tGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGAGTGATTGAAAACTTACCGCCAgg > 1:850575/1‑151 (MQ=255) tGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGaa > 1:263547/1‑150 (MQ=255) aGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGc > 2:678032/1‑151 (MQ=255) ggccggccGCGACCCGCGGGTC‑TTCCCTTTTCCCCCGCCTTTCCCCGC‑CCCGCGTTTTTCCCCCCCCA‑AAAAACAAACCCGCCGGCGGGCGAAGGAAAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACg < 2:21365/149‑1 (MQ=255) ccacccgcgccccGTCGTTCGTCCCTTTGCCCTCGGCTTTC‑CCCGC‑CCCGCGTTTTTCCCCCCCCC‑ACCAACAAACCCGCTCGCTGGCGGAGGAAAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGc < 2:241028/137‑1 (MQ=255) aGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAAc > 1:680196/1‑151 (MQ=255) gcgATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCg > 1:715351/1‑151 (MQ=255) gccacccGTCGGTCTTTCGTTTGGCCTCGGCTTTC‑CCCCC‑CCCCCGCTTTTCCCCCCCCC‑AACAACAAACCCCCTCGCGGGCGCGGGAAAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGc < 1:830892/145‑1 (MQ=255) gcgccccttggttcttCCCTTTGCCCCCGCCTTTCCCCGC‑CCCGCGTTTTTCCCCCCCCC‑AACAACAAACCCGCTCGCGGGCGAAGGAAAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACt < 2:409808/134‑1 (MQ=255) gATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTc > 1:219915/1‑151 (MQ=255) gtcccttcgttcttCCCTTTGCCCCCGCCTTT‑CCCCACCCCGCGCTTTTCCCCCCCCC‑AACAACAAACCCCCCCGCGGGCGCAGGAAAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACt < 1:761777/134‑1 (MQ=255) ccGTGGGCCTTTCCGTTT‑CCCCCGCCTTT‑CCCCG‑CCCCCCGCTTTTCCCCCCCCAAACAACAAACCCCCCCGGGGGCGCAGGAAAAAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTg < 1:993116/144‑1 (MQ=255) gtcggtcgTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCt > 1:743935/1‑151 (MQ=255) ggtcgTTCCGTTTGCCTCGGCCTTT‑CCTGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTAtt > 1:761692/1‑151 (MQ=255) cctttcgccgcgCCTTTCCCCGC‑CCCGCGCCTTTCCCCCCCCA‑AAAAACAAACCCGCTCGCGGGCGCAGGAAAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTa < 2:671055/131‑1 (MQ=255) gcccTTC‑CCCGC‑CCCGCG‑CTTTTCCCCCCCCAAAAAACAAACCCGCCCGGGGGCGAGGAAAAAAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCt < 2:97755/147‑1 (MQ=255) ctttccccgc‑cccgcGTTTTTCCCCCCCCC‑AAAAACAAACCCGCTCGCGGGCGGAGGAAAAAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGcc < 1:831445/141‑1 (MQ=255) ctttccccgc‑cccgcGCTTTTCCCCCCCCC‑AACAACAAACCCCCTCGCGGGCGCAGAAAAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCtt < 1:757156/143‑1 (MQ=255) ttccccgc‑cccgcGTTTTTCCCCCCCCA‑AAAAACAAACCCGCTCGCGGGCGGGGGAAAAAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACa < 2:606687/146‑1 (MQ=255) cccGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCAcc > 2:368598/1‑151 (MQ=255) cgACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCAttt > 1:245437/1‑151 (MQ=255) cTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTtctc > 1:326265/1‑151 (MQ=255) tttcCCCCCCCA‑AACAACAAACCCCCCCGCGGGGGCGGGTAAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCAt < 2:126953/145‑1 (MQ=255) ccccc‑AACAACAAACCCCCCCGCTGCCG‑AGGGAAACAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGaa < 2:931363/148‑1 (MQ=255) ccTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCtgatga > 1:175967/1‑150 (MQ=255) c‑aacaaCTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTg > 1:439043/1‑151 (MQ=255) aacccccccgctGGCGAAGAAAAAAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACt < 1:782414/143‑1 (MQ=255) aaCCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTg > 1:973437/1‑151 (MQ=255) ccctcgggggCGAAGGAAAAAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCg < 2:570332/143‑1 (MQ=255) ccctGCCGCAGGTCAATAAAT‑CCGCCACAGCGATGAGCGTGATTGAAATCTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGc < 1:382226/149‑1 (MQ=255) ggCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGtt > 1:303650/1‑150 (MQ=255) ggCGAAGGTAAAAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGtt < 2:1069621/150‑1 (MQ=255) gcAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTag > 2:144598/1‑149 (MQ=255) aaaacAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTAGAGGCAtt < 2:884168/147‑1 (MQ=255) aGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTAGAGGCATTa > 2:377316/1‑150 (MQ=255) | CGAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTAGAGGCATTA > NC_000913/1510249‑1510542 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |