Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 40,278 | C→T | V37I (GTC→ATC) | caiA ← | crotonobetaine reductase subunit II, FAD‑binding |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 40,278 | 0 | C | T | 98.0% | 159.0 / ‑4.6 | 51 | V37I (GTC→ATC) | caiA | crotonobetaine reductase subunit II, FAD‑binding |
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (23/27); total (23/28) | |||||||||||
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.97e-01 |
CCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGG > NC_000913/40155‑40411 | ccAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCa < 1:97675/139‑1 (MQ=255) tccaCACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCaaa > 1:223520/1‑139 (MQ=255) tccaCACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCaaa > 2:665686/1‑139 (MQ=255) acacGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATa < 1:828907/139‑1 (MQ=255) ggcgAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTcc > 1:78137/1‑139 (MQ=255) cAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCt < 2:223520/139‑1 (MQ=255) cccGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCa < 2:165096/139‑1 (MQ=255) cgcgTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTc < 1:366522/105‑1 (MQ=255) cgcgTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTc > 2:366522/1‑105 (MQ=255) cgcgTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCaaa < 1:414797/111‑1 (MQ=255) cgcgTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCaaa > 2:414797/1‑111 (MQ=255) cgTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAg > 1:234398/1‑139 (MQ=255) cgTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAg > 2:603043/1‑139 (MQ=255) gTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGt > 2:86328/1‑139 (MQ=255) cGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATAcc > 2:559104/1‑139 (MQ=255) cGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATAcc > 2:157848/1‑139 (MQ=255) tCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACaa < 2:234398/139‑1 (MQ=255) aGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCa > 1:480356/1‑78 (MQ=255) aGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCa < 2:480356/78‑1 (MQ=255) aGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAAc < 1:405588/139‑1 (MQ=255) aGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAAc < 1:86328/139‑1 (MQ=255) agcagGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGt < 2:897618/61‑1 (MQ=255) agcagGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGt > 1:897618/1‑61 (MQ=255) agcagGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCa < 2:756627/71‑1 (MQ=255) agcagGCTGTCGATACCCACATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCa > 1:756627/1‑71 (MQ=255) gcagGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTg < 2:1039491/139‑1 (MQ=255) caggcTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTc < 2:305374/68‑1 (MQ=255) cagGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTc > 1:305374/1‑72 (MQ=255) ccATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCg < 2:487276/83‑1 (MQ=255) ccATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCg > 1:487276/1‑83 (MQ=255) tatCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTGGCTGGCCATCAGTTCGCGGATACCGGCGACaaa > 2:1020857/1‑112 (MQ=255) tatCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAg < 1:34837/91‑1 (MQ=255) tatCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAg > 2:34837/1‑91 (MQ=255) tatCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACaaa < 1:1020857/112‑1 (MQ=255) tatCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATttaaattaaa > 2:586672/1‑139 (MQ=255) ccAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATAc < 1:696168/96‑1 (MQ=255) ccAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATAc > 2:696168/1‑96 (MQ=255) cAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCAt > 2:326876/1‑81 (MQ=255) cAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCAt < 1:326876/81‑1 (MQ=255) cAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATc < 2:296543/135‑1 (MQ=255) cAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATc > 1:296543/1‑135 (MQ=255) cTTTGACAAAACGTTCCGGGTAGATGTTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCa < 2:689088/139‑1 (MQ=255) cTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCa < 1:157848/139‑1 (MQ=255) tttGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCaa < 1:603043/139‑1 (MQ=255) tttGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCaa < 2:78137/139‑1 (MQ=255) aaaaCGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCtt < 1:665686/139‑1 (MQ=255) aaaCGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCCCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTg > 1:79192/1‑139 (MQ=255) ccGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATAtttt > 1:182775/1‑139 (MQ=255) gggTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTgg < 1:170925/139‑1 (MQ=255) gTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTgggg < 1:948578/139‑1 (MQ=255) gTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTgggg < 1:831751/139‑1 (MQ=255) | CCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGG > NC_000913/40155‑40411 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |