| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NZ_CP009273 | 36,297 | Δ1 bp | coding (1528/1554 nt) | caiC ← | crotonobetaine/carnitine‑CoA ligase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 36,294 | 0 | T | . | 91.7% | 82.2 / 3.5 | 24 | coding (1531/1554 nt) | caiC | crotonobetaine/carnitine‑CoA ligase |
| Reads supporting (aligned to +/- strand): ref base T (1/1); new base . (9/13); total (10/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01 | |||||||||||
CATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGC > NZ_CP009273/36173‑36421 | cATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTA‑TTTCCCCGAGCAATTAc < 1:549245/138‑1 (MQ=255)cATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTA‑TTTCCCCGAGCAATTACg > 2:505094/1‑139 (MQ=255) tttGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTA‑TTTCCCCGAGCAATTACGTgg < 1:290008/138‑1 (MQ=255) gggACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGCTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTA‑TTTCCCCGAGCAATTTCGTGGCAGATCTTTTCTg < 2:273146/139‑1 (MQ=255) caacaTCCCTGCAAAAAATGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTA‑TTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCa < 2:195818/139‑1 (MQ=255) aTCCCTGCAAAAAATGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTA‑TTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGAt > 1:324641/1‑138 (MQ=255) tGCAAAAAATGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTA‑TTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGAt > 1:459236/1‑125 (MQ=255) tGCAAAAAATGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTA‑TTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGAt < 2:459236/125‑1 (MQ=255) atatTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTA‑TTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAAtt < 1:505094/139‑1 (MQ=255) tgATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTA‑TTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTc > 2:577695/1‑139 (MQ=255) aGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTA‑TTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGtt > 1:405829/1‑130 (MQ=255) aGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTA‑TTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGtt < 2:405829/130‑1 (MQ=255) cTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTA‑TTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCg > 1:190422/1‑138 (MQ=255) gCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTA‑TTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGaa < 2:317097/138‑1 (MQ=255) ggTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTcccc > 2:91741/1‑51 (MQ=255) ggTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTcccc < 1:91741/51‑1 (MQ=255) ggTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTA‑TTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAAtt > 2:132124/1‑139 (MQ=255) ccTGCTGTTTCATTTATTTCAGATTCTTTCTAATTA‑TTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCg < 2:407979/139‑1 (MQ=255) ctgTTTCATTTATTTCAGATTCTTTCTAATTA‑TTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTc < 2:204022/138‑1 (MQ=255) tttCATTTATTTCAGATTCTTTCTAATTA‑TTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGAtct > 2:350160/1‑67 (MQ=255) tttCATTTATTTCAGATTCTTTCTAATTA‑TTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGAtct < 1:350160/67‑1 (MQ=255) tttatttCAGATTCTTTCTAATTA‑TTTCCCCGAGc > 2:511536/1‑35 (MQ=255) tttatttCAGATTCTTTCTAATTA‑TTTCCCCGAGc < 1:511536/35‑1 (MQ=255) ctttctAATTA‑TTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGc < 2:298798/138‑1 (MQ=255) | CATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGC > NZ_CP009273/36173‑36421 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |