Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,111,429 | 0 | T | G | 39.7% | 92.8 / 10.4 | 58 | G496G (GGT→GGG) | metL | bifunctional aspartate kinase II/homoserine dehydrogenase II |
Reads supporting (aligned to +/- strand): ref base T (15/20); new base G (23/0); total (38/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.13e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑TGTGGTGGACAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGACGCCAGCCGCGCGTTAGCCTTCTTCAACGAT > REL606/4111354‑4111504 | ccGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGg‑t < 2:5760375/76‑1 (MQ=255) ccGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGg‑t < 2:2180417/76‑1 (MQ=255) ggtGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGg‑tgt < 1:1121619/74‑1 (MQ=255) gTTGTCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGg‑tgt < 1:733981/76‑1 (MQ=255) tGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑Ggtgg > 1:595912/1‑76 (MQ=255) ggCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑Tgtggt < 2:5290267/76‑1 (MQ=255) ggCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGGGCTGGCAGG‑Tgtggt < 2:3886536/76‑1 (MQ=255) gCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑Tgtggtg < 2:5753515/76‑1 (MQ=255) cTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGCGTTTGTGCTGGCAGG‑TGtggggg > 2:828041/1‑76 (MQ=255) cTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑TGtggtgg < 1:251822/76‑1 (MQ=255) cTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑GGtggggg > 1:5404941/1‑76 (MQ=255) cTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑GGgggtgg > 1:5240296/1‑76 (MQ=255) tGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑TGTGGTGGa < 1:2222401/76‑1 (MQ=255) tGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑GGTGGTGGa > 2:3620483/1‑76 (MQ=255) ggAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑TGTGGTGGAc < 1:2473575/76‑1 (MQ=255) ggAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑TGTGGTGGAc < 1:6207538/76‑1 (MQ=255) ggAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑TGTGGTGGAc < 2:420226/76‑1 (MQ=255) gAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑TGTGGTGGACa < 2:1579419/76‑1 (MQ=255) gAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑TGTGGTGGACa < 2:4018490/76‑1 (MQ=255) gAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑GGTGGGGGGCa > 2:3311648/1‑76 (MQ=255) gAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAAG‑TGTGGTGGACa < 1:5907226/76‑1 (MQ=255) aaCTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑TGTGGTGGACAg < 1:2301105/76‑1 (MQ=255) aaCTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑TGTGGTGGACAg < 1:4081134/76‑1 (MQ=255) aaCTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑GGGGGGGGGCAg > 2:477592/1‑76 (MQ=255) aCTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑GGGGGTGGGCAgc > 2:2660977/1‑76 (MQ=255) cTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑TGTGGTGGACAgcc < 1:2043569/76‑1 (MQ=255) cTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑GGTGGTGGGCAgcc > 2:4271823/1‑76 (MQ=255) gTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑TGGGGGGGACAaccac > 1:2230247/1‑74 (MQ=255) ttCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑TGTGGTGGACAgccgcc < 2:1958951/76‑1 (MQ=255) ttCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑GGTGGTGGGGAgccgcc > 2:3519786/1‑76 (MQ=255) tCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑TGTGGTGGACAGccgccg < 1:1102628/76‑1 (MQ=255) tCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑GGTGGGGGGCAGccgccg > 1:2230621/1‑76 (MQ=255) tCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑GGGGGGGGGCAGccgccg > 1:4086202/1‑76 (MQ=255) tCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑GGGGGGGGACAGgcgccg > 2:5147035/1‑76 (MQ=255) cGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑TGTGGTGGACAGCcgccgc > 2:4702768/1‑76 (MQ=255) cccGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑TGTGGTGGGGAGCCGCCGCAg > 2:1178516/1‑76 (MQ=255) cccGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑TGTGGTGGGCAGCCCCCGCAg > 2:2558404/1‑76 (MQ=255) ccGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑GGTGGTGGGCCGCCGCCGCAGc > 1:4107969/1‑76 (MQ=255) ccGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑GGTGGTGGACAGCcgccgccgc > 1:2569247/1‑76 (MQ=255) cGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑GGTTGTGGACCGCCGCCGCAGcc > 2:5131595/1‑76 (MQ=255) cGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑GGTGGTGGGCAgccgccgccgcc > 2:356423/1‑76 (MQ=255) ggcagagcaCGCTTTCGGCACGTCCCGGCTTTGAGTTTGTGCTGGCAGG‑TGTGGTGGACAGCCGCCGCAGCCTGtt < 1:5813005/75‑1 (MQ=255) cagagcaCGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑TGGTGAGGACAGCCGCCTCAGTCTCTTgt > 1:1338452/1‑75 (MQ=255) cagagcaCGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑GGTGGTGGGCAGCCGCCGCAGCCCGTTGa > 2:4482130/1‑76 (MQ=255) agagcaCGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑TGTGGTGGACAACCGCCCCCGCCTGTTTAg > 2:2160682/1‑76 (MQ=255) accaCGCTTTCGGCACGTACCGGCTTTGATTTTGTGCTGGCAGG‑TGTGGTGGACAGCCGCCGCAGCCTGTTGAGCt < 1:137887/74‑1 (MQ=255) gcaCGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑TGTGGTGGACAGTCGCCGCAGCCTGTTGAGCTa < 2:4021436/76‑1 (MQ=255) aCGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑TGTGGTGGGCAACCGCCGCCACCTGTTGAGCTTTg > 2:4551580/1‑76 (MQ=255) cGCTTTCGGCACGTCCGGGCTTTGAGTTTGTGCTGGCAGG‑TGTGGGGGACAGCCGCCGCAGCCTGTTGAGCTATGa < 1:881607/76‑1 (MQ=255) cGCTTTCGGCACGTCCCGGCTTTGAGTTTGTGCTGGCAGG‑TGTGGGGGACAGCCGCCGCAGCCTGTTGAGCTATGa < 1:5939160/76‑1 (MQ=255) cGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑GGTGGTGGACAGCCGCCGCAGCCCGTTGAGCTATGa > 2:1683585/1‑76 (MQ=255) gCTTTCGGCCCGTCCCGGCTTTGAGTTTGTGCTGGCAGG‑TGTGGTGGACAGCCGCCGCAGCCTGTTGAGCTATGAc < 1:5819785/76‑1 (MQ=255) gCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑GGTGGTGGACAGCCGCCCCATGCCGTTTAGCTAttta > 2:4071620/1‑73 (MQ=255) gCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑GGGGGGGGGCCGCCGCCGCCGCCTGTTGCGCTATggc > 1:4822922/1‑74 (MQ=255) cTTTGGGCGCGACCCGGTTTTGGGTTTGTGCTGGCGGG‑GGGGGTGGACAGCCGCCGCAGCCTGTTGAGCTATGACg < 2:1366152/76‑1 (MQ=255) cTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑TGTGGGGGGGCATCGCCGCTGCCTGTTGTGCCATGAcc > 2:5933664/1‑75 (MQ=255) cTTGCGGCACGTGCCGGTTTTGATTTTGTGCTGGCGGG‑TGGGGTGGACAGCCGCCGCAGCCTGTTGAGCTATGACg < 1:4781308/76‑1 (MQ=255) tttCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑TGTGGGGGACCGGCACCGCAGCCCGTTGAGCTAAGAcag > 2:1110651/1‑74 (MQ=255) tttCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑GGGGGTGGACAGCCGCCGCAGCCTGTTGAGCTATGACgg > 2:6058853/1‑76 (MQ=255) tttCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑GGGGGGGGGAAGGCCCCGCAGCCCGTTTGGCCTTGACgg > 2:2464534/1‑76 (MQ=255) ttCGGCACGTACCGGCTTTGATTTGGTGCTGGAAGG‑TGTGGGGGACAGCCGCCGCAGCCTGTTGAGCTATGACggg < 1:586105/76‑1 (MQ=255) ttCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑TGTGGTGGGCAGCCGCCGCAGCCTGTTGGGCTATGAAggg > 1:1988256/1‑76 (MQ=255) ttCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑GGTGGGGGACAGCCGCCGCAGCCCGGTGGGCTATGACggg > 1:275460/1‑76 (MQ=255) ggCCCGTCCCGGTTTTGAGTTTGTGCGGGCGGG‑TGTGGGGGACAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTg < 2:4627904/76‑1 (MQ=255) cACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑GGTGGTGGACAGCCGCCGCAGCCTGTTGAGCTATGAAGGGCTggg > 1:2615361/1‑75 (MQ=255) aCGTCCCGGTTTTGATTTTGTGCTGGCAGG‑TGTGGTGGCCAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGAc < 1:1984874/76‑1 (MQ=255) aCGTCCCGGCTTTGAGTTTGTGCTGGCAGG‑TGGGGTGGACAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGAc < 1:2556026/76‑1 (MQ=255) cGTACCGGCTTTGAGTTTGTGCTGGCAGG‑GGTGGTGGACAGCCGCCGCAGCCTGTTTAGCTATGTCGGGGTGGGCg > 1:4977716/1‑76 (MQ=255) cGTACCGGCTTTGAGTTTGTGCTGGCAGG‑GGTGGTGGACAGCCGCCGCAGCCCGTTTTGCTATGACgggcggggcg > 2:5372864/1‑76 (MQ=255) cGTACCGGCTTTGAGTTTGTGCTGGCAGG‑GGTGGGGGACAGGCGCCGCAGCCTGTTTAGGTATGACAGGCgggggg > 2:609355/1‑73 (MQ=255) gTACCGGCTTTGAGTTTGTGCTGGCAGG‑GGGGGTGGGCAGCCGCCGCCGCCTGTTGAGCTATGACGGGGTGGagcc > 2:1390983/1‑73 (MQ=255) tcccGGTTTTGAGTTGGGGC‑GGCGGGTTGGGGGGGACAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGACGcc < 1:4906865/74‑1 (MQ=255) tACCGGCTTTGAGTTTGTGCTGGCAGG‑GGTGGGGGACAGCCGCCGCAGACTGTTGAGCTATGACCGGCTGGGCGcc > 1:5960045/1‑76 (MQ=255) cccGTCTTTGGGTTTGGGCGGGAAGG‑TGTGGGGGACAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGACGCCa < 1:1782541/75‑1 (MQ=255) aCCGGTTTTGAGTTTGTGCGGGGGGG‑TGTGGTGAACAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGACGCCa < 1:5642225/76‑1 (MQ=255) aCCGGCTTTGAGTTTGTGCTGGCAGG‑GGTGGTGGGCAGCCGCCGCCGGCTGTTGAGCTATGACGGGGTGGGCGccc > 2:5088371/1‑75 (MQ=255) ccGGTTTTGATTTTGTGCGGGCAGG‑TGGGGGGGACAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGACGCCAg < 1:309172/76‑1 (MQ=255) ccGGTTTTGAGTTTGGGCGGGCAGG‑TGTGGGGGCCAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGACGCCAg < 1:138742/76‑1 (MQ=255) ccGGCTTGGAGTTTGTGCTGGCGGG‑TGTGGTGGACAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGACGCCAg < 1:4597682/76‑1 (MQ=255) gCTTTGAGTTTGTGCTGGCAGG‑TGGGGGGGGCAGCCGGCGCAGCCTGCTGAGCTATGACGGGCTGGGagccagccc > 1:3277340/1‑75 (MQ=255) cTTTGAGTTTGTGCTGGCAGG‑GGTGGTGGGCAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGGCGCCcgccgc > 2:4903942/1‑76 (MQ=255) tttGAGTTTGTGCTGGCAGG‑GGTGGTGGACAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGGCGCCAGCcgcg > 1:4974803/1‑76 (MQ=255) tttGAGTTTGTGCTGGCAGG‑GGTGGGGGACAACCGCCGCAGCCCGTTTAGCTATTGCGGGCGGGGCGCCCGCcgcg > 1:6157616/1‑76 (MQ=255) tGAGTTTGTGCTGGCAGG‑GGTGGGGGACAGCCGCCGCAGCCTGTTGAGCTATGAAGGGCTGGGCGGCCACcgcgcg > 2:2036581/1‑76 (MQ=255) aGTTTGTGCTGGCAGG‑GGGGGTGGACAGCCCCCGCAGCCTGTTGAGCTATGACGGGGCGGAAGCCACCcccgccgt > 2:164865/1‑73 (MQ=255) tgtgCTGGCAGG‑GGTGGTGGACAGCCGCCGCAGCCTGTTTAGCTATGACGGGCTGGACGCCAGCCGCGCGTTAGcc > 2:6090598/1‑76 (MQ=255) tgCTGGCAGG‑TGTGGGGGACAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGGCGCCCGCCGCGCGTTAGCctt > 1:674323/1‑76 (MQ=255) gCTGGCAGG‑GGTGGGGGACAGCCGCCGCCGCCCGTTTAGCTCTGACGGGGCGGACGCCCGCCGCGCGGTAGCcctc > 2:2244950/1‑76 (MQ=255) cTGGCAGG‑TGTGGTGGACAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGACGCCAGCCGCGCGTTAGCcttct > 1:4609740/1‑76 (MQ=255) cTGGCAGG‑TGTGGTGGACAGCCGCCGCAGCCTGTTGAGCTATGACGGGCCGGACGCCAGCCGCGCGTTAGCcctct > 1:696162/1‑76 (MQ=255) tGGCAGG‑TGTGGGGGACAGCCGCCGCAGCCTGTTGTGCTATGACGGGCCGTGCCCCAGCCGCGCGTTATGcttctc > 2:2248522/1‑75 (MQ=255) ggCAGG‑TGTGGTGGACAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGACGCCAGCCGCGCGCTAGCCttcttc > 1:4192386/1‑76 (MQ=255) ggggg‑TGGGGGGGCCAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGACGCCAGCCGCGCGTTAGCCTTCTTCa < 1:2923666/73‑1 (MQ=255) ggtg‑tggGGTGGCCCCCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGACGCCAGCCGCGCGTTAGCCTTCTTCaa < 1:2892968/73‑1 (MQ=255) cAGG‑TGTGGTGGACAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGACGCCAGCCGCGCGTTAGCCTTCTTCaa > 1:4140789/1‑76 (MQ=255) aGG‑TGTGGTGGACAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGGCGCCAGCCGCGCGTTAGCCTTCTTCAAc > 1:1504286/1‑76 (MQ=255) aGG‑TGTGGTGGACAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGACGCCAGCCGCGCGTTAGCCTTCTTCAAc > 2:5463807/1‑76 (MQ=255) aGG‑TGTGGTGGACAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGACGCCAGCCGCGCGTTAGCCTTCTTCAAc > 2:4253236/1‑76 (MQ=255) aGG‑TGTGGTGGACAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGACGCCAGCCGCGCGTTAGCCTTCTTCAAc > 1:3678197/1‑76 (MQ=255) tgtgGTGGACAGCCGCCGCAGCCTGTTGCGCTATGACTGGATGGACGCCAGCCGCGCGTTACCCTTCTTCAACGAt > 2:4323387/1‑76 (MQ=255) | CCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGTTTGTGCTGGCAGG‑TGTGGTGGACAGCCGCCGCAGCCTGTTGAGCTATGACGGGCTGGACGCCAGCCGCGCGTTAGCCTTCTTCAACGAT > REL606/4111354‑4111504 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |