Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,571,368 | 0 | A | T | 42.9% | 80.8 / 27.6 | 56 | R69R (CGT→CGA) | yjiR | putative DNA‑binding transcriptional regulator/putative aminotransferase |
Reads supporting (aligned to +/- strand): ref base A (9/23); new base T (24/0); total (33/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.18e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAAT‑AACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTG > NC_000913/4571219‑4571511 | ggCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTAc < 2:913108/151‑1 (MQ=255) gcTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCCGGGGGTTACGtt > 1:1033980/1‑151 (MQ=255) ttCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTTGCGGCCGGGGTTTTCGtt > 1:819785/1‑149 (MQ=255) cATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGAc < 1:209690/151‑1 (MQ=255) cATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGAACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGAc < 1:472317/151‑1 (MQ=255) aTCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCCTTTCGGTGTGCggc > 1:361411/1‑148 (MQ=255) tCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGAT‑CCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCCGGGGTTTTAGTTTTGCGGCaa > 1:501242/1‑150 (MQ=255) cGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACaaaa < 2:819744/151‑1 (MQ=255) tgtgCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCATCCAGCATATCCAGCATCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCCCCGGACGCGTATTCGGCGGTACTGGTGGAGGGGCTTTACGTTGTGCGTCaaa < 2:448332/147‑1 (MQ=255) gtgCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAAT‑AAc < 2:794649/151‑1 (MQ=255) gtgCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCCGGGGTTTTTGTTGTGCGACAaaaa‑aa > 1:1019423/1‑150 (MQ=255) cTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGGTGGGGTTTTTGCTGTGCGGCAAAAA‑AAacc > 2:533458/1‑150 (MQ=255) cTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCCGGGGGTTTCGTTTTGCGACAAAAA‑AAacc > 2:340833/1‑150 (MQ=255) tttGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAAT‑AACCCg < 1:89585/149‑1 (MQ=255) tttGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAAT‑AACCCGa < 1:473448/150‑1 (MQ=255) ttGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCGGGGGGTTTCGGTGTGCGGCAAAAA‑AAacc > 2:385052/1‑147 (MQ=255) gCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAAT‑AACCCGAAcgc < 1:113665/151‑1 (MQ=255) gCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAAT‑AACCCGAAcg < 1:999600/150‑1 (MQ=255) tAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGGTTTTCGTTTTGCGGCAAAAA‑AAACCCaaaccgg > 2:908138/1‑145 (MQ=255) cGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGGCTTACGTTGTGGGGCAAAAA‑AAACCGAACGCgg > 1:354504/1‑147 (MQ=255) aaTGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTTCGTTGTGGGGCAAAAA‑AAACCGAACGCGGCCGCgggg > 1:583444/1‑147 (MQ=255) aTGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTTGCGGCTGGGCGTTTCGTTTTGCGGCAAAAA‑AACCCCCACGCCGCCGCgggggg > 2:1026450/1‑146 (MQ=255) aTGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCCGGGCTTTTCGTTTTGCGGCAAAAA‑AAACCCAACGCGGCTGCGGGGtg > 1:295234/1‑150 (MQ=255) gAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCCGGGGTTTTCGGTGTGCGGCAAAAA‑AAACCCAACCCGGCCGCgggggg > 1:304691/1‑143 (MQ=255) aaCTGTCGCTATGCGCCTCCAGCATATCCAGCGCCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAAT‑AACCCGAACGCGGCTGCGGAGTGATGAg < 2:196866/151‑1 (MQ=255) ccAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTTCGTTTTGCGGCCAAAA‑AAACCCAACGCGGCCGCGGGGGGGTGAGCTTTATCGTCTcc > 2:381659/1‑146 (MQ=255) atatCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAAT‑AACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGAt < 1:434709/151‑1 (MQ=255) tatCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCCGGGGGGTTCGGTGTGGGACAAAAA‑AAACCCAAAGCGGCTGCGGGGGGATGGGCTTCCACCTCTcccaccgctg > 2:61333/1‑140 (MQ=255) atCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGGTTTTCGGTGTGGGGCAAAAA‑AAACCCCAAGCGGCCGCGGGGGGGTGGGCTTCATCGTCTCCCGCGTCCGAt > 2:1036028/1‑149 (MQ=255) ccAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGGGGCCTGGCGTTTCGCTTTGGGGGCAAAA‑AAACCGAAAGCGGCGGCGGGGGGGTGAGCTTCATCGTCCCCCACGTCCGGTAct > 2:480910/1‑149 (MQ=255) aGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGGGTTTCGTTGTGGGACAAAAA‑AAACCCAAAGCGGGTGGGGGGGGGTGAGGTTTATTGTTTCCAGCGGCTGATaaggc > 2:68609/1‑146 (MQ=255) ccTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAAT‑AACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCctgct < 1:263411/150‑1 (MQ=255) cccACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAAT‑AACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGt < 2:473288/151‑1 (MQ=255) tAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGGTTTTCCTTGTGCGGCAAAAA‑AACCCCAAACCGGCCGCGGGGGGATGAGGCTCATCGTCTCCCGCGGCTGATACGCCCGCTTCCCGGGGGTGAt > 1:983740/1‑148 (MQ=255) aaTTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAAT‑AACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACt < 1:487555/150‑1 (MQ=255) aTTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAAT‑AACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGa < 1:542770/151‑1 (MQ=255) cACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGGGTTTCGTTGTGCGGCAAAAA‑AAACCCAACGCGGCGGGGGAGGGATGGGGTTCATCGTCTCCCGCGTCTGATACGCCTGCCGCCCGGGGGTGGTAATGACGcc > 1:901000/1‑150 (MQ=255) aCCGGGCGCTGCACCGGACGCGTCATCGGCGGTTCTGGCGGGCGGGCTTTTCGGTGTGCGGCAAAAA‑AAACCCAAACCGGCCGCGGGGGGGTGAGCTTTCTCGTCTCCCGCGGCTGATACGCCTGCTGCAAGGGGCTGATACTGACGcc > 1:189359/1‑149 (MQ=255) gggCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAAT‑AACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCt < 2:442299/151‑1 (MQ=255) ggCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGGGTTTCGTTTTGCGGCAAAAA‑AAACCCAAACCGGCCGCGGGGGGGTGAGGTTCATCCTCTCCCGCGGCTGATACGCCCCCTGCCCGGGGGCGGTACTGAAGcc > 2:555235/1‑145 (MQ=255) cgcTGCACCGGACGCGTCATCGGCGGTACTGGCGGGTGGGGTTTTCGGTGTGGGACAAAAA‑AACCCCGAAGCGGGTGCGGGGGGGTGAGGCTCCTCGTCTCCCCCGTCTGATACCCCCGCTGCACGGGGCTGGTACTGACGCCCCGCTcc > 1:887285/1‑150 (MQ=255) gCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAAT‑AACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACt < 2:899613/151‑1 (MQ=255) cACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGGTTTTCGTTGTGCGACAAAAA‑AAACCCAACGCGGCTGCGGGGGGGTTAGCTTCATCGTCTCCAGCGTCTGATACGCCCGCTGCACGGGGCTGATACTGACGCCCTGCTCCTGACtt > 1:475291/1‑151 (MQ=255) aCCGGACGCGTCATCGGCGGTACTTGCGGCTGGGGTTTTCGGTGTGGGGCCAAATAAACCCGAA‑GCGGGTTCGGGGGGATGAGGTTTCTTGTCTCCCGCCTCTTATAACCCCGCTGCCCGGGGGGGGTACTGACGCCCTGCTCCTTACtt > 1:105785/1‑150 (MQ=255) aCGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAAT‑AACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTg < 1:183867/151‑1 (MQ=255) cgcgTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAAT‑AACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTgc < 2:492182/151‑1 (MQ=255) cgcgTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAAT‑AACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTgc < 1:266455/151‑1 (MQ=255) ggcggTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAAT‑AACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACg > 1:914363/1‑151 (MQ=255) gcggTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAAT‑AACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACgg < 2:304675/151‑1 (MQ=255) cggTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAAT‑AACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACgg > 1:957843/1‑150 (MQ=255) tACTGGCGGCTGGGCTTTACGTTGTGCGACAAAAT‑AACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAAt > 1:562834/1‑151 (MQ=255) cTGGCGGCTGGGCTTTACGTTGTGCGACAAAAT‑AACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAAttt > 1:254272/1‑151 (MQ=255) cggcTGGGCTTTACGTTGTGCGACAAAAT‑AACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTcc < 1:1022185/151‑1 (MQ=255) ggcTGGGCTTTACGTTGTGCGACAAAAT‑AACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTccc > 2:954998/1‑151 (MQ=255) cTGGGCTTTACGTTGTGCGACAAAAT‑AACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCg > 1:1011432/1‑151 (MQ=255) gggCTTTACGTTGTGCGACAAAAT‑AACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTg > 1:136487/1‑151 (MQ=255) | GGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAAT‑AACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTG > NC_000913/4571219‑4571511 |
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: ‑ |