| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 372,640 | 0 | T | G | 28.8% | 112.9 / 11.7 | 52 | L176V (TTA→GTA) | mhpD | 2‑hydroxypentadienoate hydratase |
| Reads supporting (aligned to +/- strand): ref base T (8/29); new base G (15/0); total (23/29) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.09e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TTGAATGGGTACTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCGCT > NC_000913/372503‑372779 | ttGAATGGGTACTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTa > 2:413643/1‑140 (MQ=255)ttGAATGGGTACTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTa > 1:980431/1‑140 (MQ=255) aaTGGGTACTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGAc < 1:1066787/140‑1 (MQ=255) aaTGGGTACTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGAc < 1:431181/140‑1 (MQ=255) cTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGaaaaa < 2:152930/140‑1 (MQ=255) ccGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTg < 1:413643/140‑1 (MQ=255) ccGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTTGACCTGAAAAAATg > 1:281003/1‑140 (MQ=255) ccGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTg > 1:1235978/1‑140 (MQ=255) cGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTgc < 2:857441/140‑1 (MQ=255) gcgcTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTgcgc > 1:1183785/1‑140 (MQ=255) aGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCC‑ATGAAGa < 1:137561/140‑1 (MQ=255) ggtggGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCC‑ATGAAGATGa < 2:415788/140‑1 (MQ=255) gtggGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCC‑ATGAAGATGAc > 2:478657/1‑140 (MQ=255) tggGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCC‑ATGAAGATGAcg < 2:1235978/140‑1 (MQ=255) ggggAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCC‑ATGAAGATGAcgc > 2:638180/1‑140 (MQ=255) gggAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCC‑ATGAAGATGAcgcg < 1:638180/140‑1 (MQ=255) gggAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCC‑ATGAAGATGAcgcg < 1:478657/140‑1 (MQ=255) ggAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGGAACGGCCGGCGGGGTTAGACCTGAAAAAAAGCTCC‑ATGAAGATGACGCGt > 2:760751/1‑140 (MQ=255) ggAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGt > 1:567279/1‑140 (MQ=255) aGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGtaa < 1:733084/140‑1 (MQ=255) aGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGtaa > 1:507846/1‑140 (MQ=255) cGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGtaataa < 2:567279/140‑1 (MQ=255) gCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAAc > 2:1048342/1‑140 (MQ=255) cATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAACg < 2:1183785/140‑1 (MQ=255) cATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAAATGCGCC‑ATGAAGATGACGCGGAATAACg > 2:506623/1‑140 (MQ=255) cATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGGCCGGCGGAACGTCCGGCGGGGGTAGAACTGAAAAAATGCGCC‑ATGAAGATGAAGCGgaataaag > 2:45172/1‑138 (MQ=255) ttCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAACGaa < 1:1048342/140‑1 (MQ=255) tCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAACGAag < 1:760751/140‑1 (MQ=255) gTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAACGAAGAGGTTTCTa < 1:328997/140‑1 (MQ=255) cAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGgcg < 2:980431/140‑1 (MQ=255) tttGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAAAGAAGAGGGTTCTAGCGGGGGCgg > 2:189604/1‑140 (MQ=255) tttGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCCGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCC‑ATGAAGAAGACGCGTAATAAAGAAGAGGGTTTTAACGGGGGCgg > 2:328997/1‑140 (MQ=255) cGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCC‑CTGAAGATGACGCGTAATAACGAAGAGGGTTCTAGCGGGCGCGGCcgcg > 2:458876/1‑140 (MQ=255) gATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGa < 1:1099869/140‑1 (MQ=255) gATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGa < 2:608549/140‑1 (MQ=255) gATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGa < 2:624787/140‑1 (MQ=255) gTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCt > 2:727759/1‑140 (MQ=255) ggCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTACCGCGCAACGTTCGGCGGGGGTTCAGCTGAAAAAATGGGCCAATGAA‑AAGAAGCGGCATAAAAGAGAGGGTTTTAGCGGGCGCGGGCGCGAAtacactg > 2:1248316/1‑134 (MQ=255) gATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATg < 2:1064103/131‑1 (MQ=255) gATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAACGAAGAGGGTTCTAGCGGGCGCGGCAGCGAATg > 1:1064103/1‑131 (MQ=255) gCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCt < 1:1248316/140‑1 (MQ=255) tgtATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTc > 2:1109682/1‑140 (MQ=255) gTCATCGGCGGTCCGGCGCCACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAAAGAAGAGGGTTCTAGCGGGCGCGGCAGCGAAAGCCTGGGACATCCCCTTAATGCGGCCGGctggcc > 2:700781/1‑139 (MQ=255) gTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAAATGCGCC‑ATGAAGATGGCGGGGAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATTGCTGGGGCATCCGGTTAATGCGGCCGTctgggc > 2:327680/1‑138 (MQ=255) cATCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAAATGCGCC‑CTGAAAAAGACGCGTAAAAAAGAAGAGGGTTCTAACGGGCGCGGCCGCGAAAGCCTGGGGCAACCGCTTAATGCGGCCGCctggctgg > 2:1252729/1‑140 (MQ=255) aTCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCtggctggc > 2:985464/1‑140 (MQ=255) tCGGCGGTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCC‑CTGAAGATGACGCGTAATAACGAAGAGGGTTCTAGCGGGGGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTctggctggga > 2:478462/1‑138 (MQ=255) cggcggTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCa < 2:1075095/139‑1 (MQ=255) cggcggTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCa > 1:1075095/1‑139 (MQ=255) cggcggTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAACGAAGAGGGTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCAc > 2:1157109/1‑140 (MQ=255) gcggTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTctggctg < 1:742771/134‑1 (MQ=255) gcggTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGc > 1:868106/1‑140 (MQ=255) gcggTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTctggctg > 2:742771/1‑134 (MQ=255) gcggTCCGGCGCAACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAACGAAGAGGGTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGc > 2:961984/1‑140 (MQ=255) cGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATg < 1:506623/140‑1 (MQ=255) ggCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATgg < 1:985464/140‑1 (MQ=255) gcgcAACGTCCGGCGGGGTTAGACCTGAAAAACTGTGCC‑ATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGc < 2:868106/140‑1 (MQ=255) gcgcAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGc < 1:700781/140‑1 (MQ=255) cAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAg < 1:189604/140‑1 (MQ=255) ggTTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCg > 1:743150/1‑140 (MQ=255) ttAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCGCt < 1:109705/140‑1 (MQ=255) | TTGAATGGGTACTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCC‑ATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCGCT > NC_000913/372503‑372779 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |