Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 372,640 | 0 | T | G | 34.0% | 70.1 / 20.3 | 47 | L176V (TTA→GTA) | mhpD | 2‑hydroxypentadienoate hydratase |
Reads supporting (aligned to +/- strand): ref base T (8/23); new base G (16/0); total (24/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.52e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CATTGAATGGGTACTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCG > NC_000913/372501‑372777 | cATTGAATGGGTACTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGt < 1:81251/140‑1 (MQ=255) cATTGAATGGGTACTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGt < 1:307591/140‑1 (MQ=255) aTTGAATGGGTACTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGtt < 2:935119/140‑1 (MQ=255) ttGAATGGGTACTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTa > 2:280767/1‑140 (MQ=255) ttGAATGGGTACTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTa > 2:301360/1‑140 (MQ=255) tGAATGGGTACTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAg < 1:665974/140‑1 (MQ=255) gTACTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACC‑TGaa > 1:599811/1‑140 (MQ=255) aCTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGGCC‑TGaaaa > 2:1508056/1‑140 (MQ=255) ctgccGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACC‑TGaaaaa < 1:1402840/137‑1 (MQ=255) tCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACC‑TGAAAAACt > 1:104602/1‑140 (MQ=255) tCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACC‑TGAAAAACt > 2:289601/1‑140 (MQ=255) ggCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACC‑TGAAAAACTgc‑g < 2:168486/140‑1 (MQ=255) gcgcTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGGCC‑TGAAAAACCgc‑gc > 2:58494/1‑140 (MQ=255) cgcTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACC‑TGAAAAACTGC‑Gcc < 1:289601/140‑1 (MQ=255) cgcTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACC‑TGAAAAACTGC‑Gcc > 1:349353/1‑140 (MQ=255) gcTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGTAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACC‑TGAAAAACTGC‑GCCa < 2:547181/140‑1 (MQ=255) gcTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACC‑TGAAAAACTGC‑GCCa < 2:720019/140‑1 (MQ=255) ggAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAAACTTCGGCGGGGTTTAACC‑CGAAAAACTGC‑GCCATGa > 1:1310119/1‑140 (MQ=255) acaggtggGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACC‑TGAAAAACTGC‑GCCATGAAGa < 1:755503/137‑1 (MQ=255) aGTGGTGGGGAGCCGCATTCGCGACTGGTCTATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACC‑TGAAAAACTGC‑GCCATGAAGa < 2:801750/140‑1 (MQ=255) tggGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACC‑TGAAAAACTGC‑GCCATGAAGATGAcg < 1:1065626/140‑1 (MQ=255) gggAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACC‑TGAAAAACTGC‑GCCATGAAGATGAcgcg < 1:768916/140‑1 (MQ=255) ggAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACC‑TGAAAAACTGC‑GCCATGAAGAAGAcgcgg > 2:209388/1‑139 (MQ=255) gAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGta > 1:212965/1‑140 (MQ=255) ccGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACC‑TGAAAAACTGC‑GCCATGAAGATGAcg < 2:478242/133‑1 (MQ=255) ccGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACC‑TGAAAAACTGC‑GCCATGAAGATGAcg > 1:478242/1‑133 (MQ=255) ccGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGtaata > 1:254342/1‑140 (MQ=255) cATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGTAATAACg < 1:280767/140‑1 (MQ=255) cATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGCCCGGCGGGGGTTGACCTTAAAAAACTG‑‑GCCCTGAAGATTACGCGtaaattgt > 2:755503/1‑135 (MQ=255) gTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACC‑TGAAAAACTTC‑GGCCTGAAGATGAAGCGTAATAACGAAGAGGTTTCtc > 1:1230300/1‑139 (MQ=255) cAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCCGCGGGGGTTGACC‑TGAAAAAATGC‑GCCATGAAGATGACGCGTAATAAAGAAGAGGGTTCTTGCGGgcg > 2:1011815/1‑140 (MQ=255) aGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGgcgc < 2:349353/140‑1 (MQ=255) tttGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCgg > 1:1100737/1‑140 (MQ=255) tttGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGACGGTCCGGCGCAACGTCCGGCGGGGTTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCgg > 1:345802/1‑140 (MQ=255) cGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCg < 2:1325947/140‑1 (MQ=255) gATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGa < 2:104602/140‑1 (MQ=255) gATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGa < 2:254342/140‑1 (MQ=255) aTACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGACGCAACGTCCGGCGGGGTTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGaa < 2:440086/140‑1 (MQ=255) tAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATc > 2:1480190/1‑140 (MQ=255) cGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGc > 1:563043/1‑140 (MQ=255) cGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGGAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGc > 1:1422052/1‑140 (MQ=255) ctcctGTGGGGTGTATGTCATCGGCGGTCCGGCGCACCCGCCCGCGGGGGTTGAAC‑TGAAAAACTGCGGCCAT‑AAGATGACGCGCAATTAAGGAGAGGGTTCTCGCCGGCGCCGCACCGAATGCCTGGGACATTCCCttt > 2:762887/1‑139 (MQ=255) ctGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGTAATAACGAAGAGGGTTCTAGCCGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAAtt > 1:85642/1‑139 (MQ=255) tATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTctg > 1:1096501/1‑140 (MQ=255) tATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTctg > 1:365228/1‑140 (MQ=255) tATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTctg > 2:810050/1‑140 (MQ=255) tGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTctggc > 2:965604/1‑140 (MQ=255) gTCATCGGCGGTCCCGCGCAACGTCCGGCGGGGGTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGGAATAAAGAAGAGGGTTTTAGCGGGCGCCGCAGCGAATGCCTGGGGCATCCGCTTAATGGGGCCGGcttggt > 2:708118/1‑136 (MQ=255) cggcggTCCGGCGCAACGTCCGGCGGGGGTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGTAATAACGAAGAGGGTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCAc > 1:994515/1‑140 (MQ=255) cGGCGCAACGTCCGGCGGGGTTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATg < 1:1011815/140‑1 (MQ=255) gcgcAACGTCCGGCGGGGTTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGc < 2:212965/140‑1 (MQ=255) gcgcAACGTCCGGCGGGGTTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGc < 2:1100737/140‑1 (MQ=255) gcgcAACGTCCGGCGGGGTTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGc < 1:209388/140‑1 (MQ=255) gcgcAACGTCCGGCGGGGTTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACACAAAATGGc < 2:1471611/140‑1 (MQ=255) cGTCCGGCGGGGTTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCt > 1:1373009/1‑140 (MQ=255) ggcggGGTTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTg > 2:263733/1‑140 (MQ=255) ggggTTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAAc < 1:1480190/140‑1 (MQ=255) ggTTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCg > 1:1464990/1‑140 (MQ=255) | CATTGAATGGGTACTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACC‑TGAAAAACTGC‑GCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCG > NC_000913/372501‑372777 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |