Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,003,921 | 0 | T | G | 24.2% | 52.8 / 10.7 | 33 | Y145D (TAT→GAT) | xdhC | putative xanthine dehydrogenase iron‑sulfur‑binding subunit XdhC |
Reads supporting (aligned to +/- strand): ref base T (4/21); new base G (8/0); total (12/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.57e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGACAGTTGCTTATACCTTGCCGCCTGGGCTGAAGGAAAAGAGATCCGCACGCTGGAAGGTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCTTTGTCGCTGCACGGGGTAT‑CAGATG‑ATTGTA‑AATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTG‑AATTCAGGCCGG‑ATTCACTGAGGTTATGTGTTTAACAACTCATATTTCTTAATCTTGCGATAGAGCGTAGCAATGCCGATGCCCAGTTCATCAGCAACTTGCTTCTTGCTGTTATGACGTGAAAGCGCCTCGCGGATCATTTGCTTTTCCATCTCCTCCAGCGC > NC_000913/3003671‑3004148 | aGACAGTTGCTTATACCTTGCCGCCTGGGCTGAAGGAAAAGAGATCCGCACGCTGGAAGGTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCTTTGTCGCTTCACGGGGt < 1:232479/251‑1 (MQ=255) aGACAGTTGCTTATACCTTGCCGCCTGGGCTGAAGGAAAAGAGATCCGCACGCTGGAAGGTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGGCGCGGACTGGCGGGAAATCTTTGTCGCTGCACGGGGt < 1:272838/251‑1 (MQ=255) gCTTATACCTTGCCGCCTGGGCTGAAGGAAAAGAGATCCGCACGCTGGAAGGTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGGCGCGGACTGGCGGGAAATCTTTGTCGCTGCACGGGGTAT‑CAgatg < 1:93544/251‑1 (MQ=255) tataCCTTGCCGCCTGGGCTGAAGGAAAAGAGATCCGCACGCTGGAAGGTGACGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCTTTGTTGCTGCACGGGGTAttagctgaat > 1:193789/1‑242 (MQ=255) ataCCTTGCCGCCTGGGCTGAAGGAAAAGAGATCCGCACGCTGGAAGGTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCTTTGTCGCTGCACGGGGGAT‑CAGATG‑ATTg > 1:233319/1‑251 (MQ=255) gggCTGAAGGAAAAGAGATCCGCACGCTGGAAGGTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGGAATATTTTTTGCTGCACGGGGGAT‑AAGATG‑ATTGTA‑AAAACAGTTCTg > 2:93530/1‑251 (MQ=255) cACGCTGGAAGGTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCTTTGTCGCTGCACGGGGTAT‑CATATG‑ATTGTA‑AATACAGTTCTGGATTGCGAGAAAACGAAGTaaa < 1:196026/251‑1 (MQ=255) cGCTGGAAGGTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCTTTGTCGCTGCACGGGGTAT‑CAGATG‑ATTGTA‑AATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAg < 1:97595/251‑1 (MQ=255) cGCTGGAAGGTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATCCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCCTTGACGCTGCCCGGGGGAT‑CAGATG‑ATTGTC‑AATACAGTTTTGGGATGCGCGAAAACGAAGGaaaa > 1:119269/1‑250 (MQ=255) cGCTGGAAGGTGAAGCGAAAGGCGGAAAACTGTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCTTTGTCGCTGCACGGGGTAT‑CAGATG‑ATTGTA‑AATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAg < 2:171459/251‑1 (MQ=255) cTGGAAGGTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCTTTTTCGCTTGACGGGGGGT‑CAAATG‑AT‑GTACAATACAGTTTTGGATTGTGAGAAAAAGAAGAaaaaagg > 1:14937/1‑248 (MQ=255) ggAAGGTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCTTTGTCGCTGCACGGGGTAT‑CAGATG‑ATTGTA‑AATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGata < 2:238188/251‑1 (MQ=255) aGGTGAATCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCTTTGTCGCTGCACGGGGTAT‑CAGATG‑ATTGTA‑AATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATcc < 2:167942/251‑1 (MQ=255) gAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCTTTGTCGCTGCACGGGGTAT‑CAGATG‑ATTGTA‑AATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGc < 2:182045/250‑1 (MQ=255) aGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGGAATCTTTGTCGCTGCACGGGGGGT‑CAGCTG‑ATTGTA‑AATACAGTTTTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGcc > 1:150534/1‑249 (MQ=255) aGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCTTTGTCGCTGCACGGGGTAT‑CAGATG‑ATTGTA‑AATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGcc < 2:150534/249‑1 (MQ=255) gCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCTTTGTCGCTGCACGGGGTAT‑CAGATG‑ATTGTA‑AATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTG‑AATTCAgg < 1:190697/251‑1 (MQ=255) gAAAACTTACTCAGGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTGTGTACGCCGGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCAGTACGGAAATTCGTCGCGGACTGGCGGGAAATCGTTGTCGCTGCACGGGGTAT‑CAGATG‑ATTGTA‑AATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGACATCCGGCCTG‑AATTCAGGCCg < 2:158120/251‑1 (MQ=255) aaaaCTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCTTTGTCGCTGCACGGGGTAT‑CAGATG‑ATTGTA‑AATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTG‑AATTCAGGCCgg < 1:77149/251‑1 (MQ=255) cgttCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTAGTACGCCGGGCCTGAGTATGGCTACCACGGCAATGCGGGCGAGACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCTTTGGCGCTGCACGGGGTAG‑CAGATG‑ATTGTA‑AAGACAGTTCTGGATTGCGAGAAAACGAAGCAAAAGGAGATCCGGCCTG‑AATTCAGGCCGG‑ATTc < 2:247461/249‑1 (MQ=255) cagcagGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCTTTGTCGCTGCACGGGGTAT‑CAGATG‑ATTGTA‑AATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTG‑AATTCAGGCCGG‑ATTCACTGAGGTTAtg < 1:93530/251‑1 (MQ=255) cagcagGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCTTTGTCGCTGCACGGGGTAT‑CAGATG‑ATTGTA‑AATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTG‑AATTCAGGCCGG‑ATTCACTGAGGTTAtg < 2:233319/251‑1 (MQ=255) tATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCTTTGTCGCTGCACGGGGTAT‑CAGATG‑ATTGTA‑AATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTG‑AATTCAGGCCGG‑ATTCACTGAGGTTATGTGTTTaac < 2:86180/250‑1 (MQ=255) cGGCGCAGTGCAGTGCGACTTTTGTACGCCGGGCCTGATTAGGGCTACCCCGGCAATGCTGGCGAACCCACGCGAGAAGCCATTAACCATTACGGAAATTCTGCGCGGACTAGCGTGAAATCTTTGTCGCTGCACGGGGTAT‑CAGATG‑ATTGTA‑AATACAGTGCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTG‑AAGTCAGGCCGT‑ATTCACTGAGGTTATGTGTTTAGCTACTCATATTTc < 2:175768/251‑1 (MQ=255) cgcAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAACTTTTTTGCTGCACGGGGGATACA‑ATGGATT‑TA‑AAAACAGGTCTGGAATGCGAGAAAACGAAGTAAAAAGAGATCCGGCCTT‑AAATCAGGCCGG‑GTTCGCTGAGGGTATgtgtgttacaaacccaactttgtt > 2:6233/1‑231 (MQ=255) cctggcctgATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCTTTTTCGCTGCACGGGGGAT‑CAGATG‑ATTGTA‑AATACAGTTCTGGGTTGCGGGAAAACGGAGGAAAAGGGTATTCCGCCTG‑AATTCAGGCCGG‑GTTTACTGGGGGTATGTGTTTAAACACTCATATTTTTTCATCTTGTGAt > 1:243394/1‑236 (MQ=255) ccgggcctgATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCTTTGTCGCTGCACGGGGTAT‑CAGATG‑ATTGTA‑AATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTG‑AATTCAGGCCGG‑ATTCACTGAGGTTATGTGTTTAACAACTCATATTTCTTAATCTTGCGAt < 2:243394/233‑1 (MQ=255) ctgATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGGAAACTCTTTCGCTGCACGGGGGAT‑TAGATG‑ATTGTA‑AATACAGTTCTGGATTGCGAGAAAACACAGTAAAAGAGTATCCGGCCTG‑GATTTAGGCCGGCGTTC‑CTTAGGGTTTGTGTTTTACCACTCATATTTTTTTATCTTTCGGTTGAGCGGAGACATGCCGATGc > 2:26015/1‑251 (MQ=255) ttATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGGAATCTTTTTTGCTGCCCGGGGTAT‑CAGATG‑GTTGTA‑AATACCGGTCTGGAGTTCGGGAAAAAGAAGTAAAAAGAGATTCGCCCTGTAATTC‑GGGCGG‑GTTCACTGAGGGTTTGTGTTTTACAAAACACATTTTTTTATCTTGTGGTAGAGGGTGGCAATGCCgcaggcccc > 2:290506/1‑243 (MQ=255) aTGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCTTTGTCGCTGCACGGGGTAT‑CAGATG‑ATTGTA‑AATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTG‑AATTCAGGCCGG‑ATTCACTGAGGTTATGTGTTTAACAACTCATATTTCTTAATCTTGCGATAGAGCGTAGCAATGCCGATGCCCAGtt < 1:263838/251‑1 (MQ=255) aCGCGAGAAGCCATTAACCATTACGGAAATTCGGCGCGGACTGGTGGGAAATCGTTGTCGCTGCACGGGGTAT‑CAGAGG‑AGTGTA‑AATACAGTTCGGGATTGCGAGAAAACGAAGTAAAAgg < 1:42635/122‑1 (MQ=255) gagaAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCTTTGTCGCTGCAAGGGGTAT‑CAGATG‑ATTGTA‑AATACAGTTCTGGGTTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTG‑AATTCAGGCCGG‑ATTCACTGAGGTTATGTGTTTAACAACTCATATTTTTTAATCTTTCGATAGAGCGTAGCAATGCCGATGCCCAGTTCATCAGCAACTTGCTTCTTGCTGTTTTGACggg > 1:276003/1‑249 (MQ=255) tACGGAAATTCGTCGCGGACTGGCGGGAAATCTTTGTCGCTGCACGGGGTAT‑CAGATG‑ATTGTA‑AATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTG‑AATTCAGGCCGG‑ATTCACTGAGGTTATGTGTTTAACAACTCATATTTCTTAATCTTGCGATAGAGCGTAGCAATGCCGATGCCCAGTTCATCAGCAACTTGCTTCTTGCTGTTATGACGTGAAAGCGCCTCGCGGAt < 1:201657/250‑1 (MQ=255) aaTTCGTCGCGGACTGGCGGGAAATCTTTGTCGCTGCACGGGGTAT‑CAGATG‑ATTGTA‑AATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTG‑AATTCAGGCCGG‑ATTCACTGAGGTTATGTGTTTAACAACTCATATTTCTTAATCGTGCGATAGAGCGTAGCAATGCCGATGCCCAGTTCATCAGCAACTTGCTTCTTGCTGTTATGACGTGAAAGCGCCTCGCGGATc < 1:302276/245‑1 (MQ=255) gAAATCTTTGTCGCTGCACGGGGTAT‑CAGATG‑ATTGTA‑AATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTG‑AATTCAGGCCGG‑ATTCACTGAGGTTATGTGTTTAACAACTCATATTTCTTAATCTTGCGATAGAGCGTAGCAATGCCGATGCCCAGTTCATCAGCAACTTGCTTCTTGCTGTTATGACGTGAAAGCGCCTCGCGGATCATTTGCTTTTCCATCTCCTCCAgcgc > 2:159822/1‑251 (MQ=255) | AGACAGTTGCTTATACCTTGCCGCCTGGGCTGAAGGAAAAGAGATCCGCACGCTGGAAGGTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCTTTGTCGCTGCACGGGGTAT‑CAGATG‑ATTGTA‑AATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTG‑AATTCAGGCCGG‑ATTCACTGAGGTTATGTGTTTAACAACTCATATTTCTTAATCTTGCGATAGAGCGTAGCAATGCCGATGCCCAGTTCATCAGCAACTTGCTTCTTGCTGTTATGACGTGAAAGCGCCTCGCGGATCATTTGCTTTTCCATCTCCTCCAGCGC > NC_000913/3003671‑3004148 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |