Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 4,581,138 | C→T | 6.5% | G109S (GGT→AGT) | hsdS ← | type I restriction enzyme EcoKI specificity protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,581,138 | 0 | C | T | 6.5% | 178.8 / 2.5 | 62 | G109S (GGT→AGT) | hsdS | type I restriction enzyme EcoKI specificity protein |
Reads supporting (aligned to +/- strand): ref base C (32/25); new base T (2/2); total (35/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01 |
CAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACA > NC_000913/4581001‑4581280 | caaatcaaaGCTTGCCGGCTTAATATTATTTATATTTGCACCGGCAGAAAGTGATGAAATTGTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCaaa > 1:945583/1‑143 (MQ=255) caaatcaaaGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCaaa > 1:1006075/1‑143 (MQ=255) aaatcaaaGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAACCACAGAGGCTATCAGTTTTTCAGGACGTAATACACCGCaaaa > 2:899755/1‑143 (MQ=255) aatcaaaGCTTGCCGGCTTAATATTATTAATATTTGCAGCAGCAGAAAGTGATGAAATTTTGTTTCGAGAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAAt < 2:1040310/143‑1 (MQ=255) aatcaaaGCCTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAAt < 1:379017/143‑1 (MQ=255) atcaaaGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATg < 1:975082/143‑1 (MQ=255) tcaaaGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATgc > 1:76282/1‑143 (MQ=255) gCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGGAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACt > 2:864993/1‑143 (MQ=255) cGGCTTAATATTATGAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAACCCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTAc < 2:1207294/143‑1 (MQ=255) cTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAGTATAAGTGTTTCAGGACGTAATACACCGCAAAATGCGCCGACACTACAtt > 1:723827/1‑143 (MQ=255) aaTATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCaa > 1:1012323/1‑143 (MQ=255) tattatTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAAt < 1:954223/143‑1 (MQ=255) ttaATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAg > 2:1146868/1‑143 (MQ=255) ttaATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAg < 2:732322/143‑1 (MQ=255) taATATTTGCACCAGCAGAAAGTGATGAAATTCTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGa < 1:495502/143‑1 (MQ=255) atatTGGCACCAGTAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATg < 2:744626/143‑1 (MQ=255) cagcagAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACTGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGAtgtgtg > 2:1097681/1‑141 (MQ=255) cagcagAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACTGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGAtgtgtg > 2:1096742/1‑141 (MQ=255) cagcagAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCg < 2:788596/143‑1 (MQ=255) gcagAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGa < 2:293246/143‑1 (MQ=255) gcagAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGa < 2:338956/143‑1 (MQ=255) aaaGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTa < 1:125403/143‑1 (MQ=255) aaaGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCCGATGTGCGGAttga > 2:1259868/1‑141 (MQ=255) gTGATGACATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCa > 1:908017/1‑143 (MQ=255) ttttGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGAt < 1:1008568/143‑1 (MQ=255) tttGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGAtt < 1:1072870/143‑1 (MQ=255) gTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGc < 1:516446/143‑1 (MQ=255) ttCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCCGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTc > 2:634698/1‑143 (MQ=255) gcGATAAAGAGAAGATTTTGTGAAATGAGAAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTcc < 2:1213887/142‑1 (MQ=255) gATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCt > 1:643578/1‑143 (MQ=255) gATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCt < 2:1034748/143‑1 (MQ=255) tctaaaagaagaTTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGAATTTACCAACTACGGATTTGCTCCCtg > 2:141348/7‑143 (MQ=255) tAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCtga > 2:14020/1‑143 (MQ=255) aaaGAGAAGAGTTCGTGAAATGAGCCATAGAACCAGAAAATATAAGTGTCTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGAGTTACCAACTACGGATTTGCTCCCtgat < 2:950212/143‑1 (MQ=255) agagAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCtgatga > 1:1019324/1‑143 (MQ=255) agaagaTTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACa > 1:582558/1‑143 (MQ=255) aagaTTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACAtt > 1:250368/1‑143 (MQ=255) agaTTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTg > 1:546665/1‑143 (MQ=255) aTTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTCGAAGCTGATGTGCGGATTTACCAACTACGGCTTTGCTCCCCAATGACAtcaca > 2:616588/1‑139 (MQ=255) ttttGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCaa > 2:399799/1‑143 (MQ=255) tggggAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAAGACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATa < 2:316913/139‑1 (MQ=255) ggtgAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATaa < 1:490571/142‑1 (MQ=255) aaTGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATCGCAATACCAata > 1:229878/1‑143 (MQ=255) cAATAAATCCATAAAATATAAGTTGTGCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCa < 2:303189/143‑1 (MQ=255) aaTAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAg < 2:1012323/143‑1 (MQ=255) aTAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAgg > 1:1137690/1‑143 (MQ=255) tAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGga < 1:919821/143‑1 (MQ=255) tAAAACCAGAAAATATAAGTTTTTCAGGAAACAAAAGACAGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGga > 1:878929/1‑143 (MQ=255) cAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATAtt > 1:298192/1‑143 (MQ=255) aaaaTATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGAGGGGTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATATTCAGGGGATAttttt > 1:654545/1‑143 (MQ=255) atatAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGCGACATTTTTTGa > 2:1083412/1‑143 (MQ=255) attttttCAGGACGTAATACACCGCATATTGCGCCGAAAATACAGTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGActttc < 2:140687/141‑1 (MQ=255) tttttCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGActttctt > 1:185962/1‑143 (MQ=255) tttttCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGActttctt < 1:629411/143‑1 (MQ=255) tttCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATGTTCAGGAGATCTCTTTTGACTTTCTTTa > 2:554675/1‑143 (MQ=255) aGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACaa > 1:1113251/1‑143 (MQ=255) ggACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAg > 1:1148981/1‑143 (MQ=255) cGTAATACACCGCAAAATGCGACGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGAATTGCACCCTGATGACATTGCAATAACAATAACGTCAGGAGATATTTTTAGACTTTCGTTaaggagaat > 2:1057318/1‑136 (MQ=255) aTACACTGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGTGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTa < 1:1097681/143‑1 (MQ=255) aTACACTGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGTGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTa < 1:1096742/143‑1 (MQ=255) acacCCCAAAATGCGCCGAAACTACCTTCAAATGGTAGCTGCTGATGTGCGGAATTACCAACTACAGATTTCCTCCCTGATGACATTGCAATAACAATATCCTCAAGAGATACTTTTTGACTTTCTTTTACAAAATTTTTAgg > 2:1213470/1‑143 (MQ=255) cacCGCAAAATGCGCCGAAACTACATTCAAATGGTAGCTGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAGTATCTTCAGGAGATATTTTTTGACTTTCTTTCACAACATTTATAgta > 2:593249/1‑141 (MQ=255) cGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGaa < 1:468019/141‑1 (MQ=255) cGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGaa > 2:468019/1‑141 (MQ=255) cGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACa > 1:890100/1‑143 (MQ=255) cGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGCGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTCCAATAACAATATCCTCAGCAGACATTTTTTCACTTTCTTTAACAAGATCCTTAGGaccc > 2:325675/1‑140 (MQ=255) | CAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACA > NC_000913/4581001‑4581280 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |