| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 3,271,600 | (C)5→6 | coding (171/1146 nt) | garK ← | glycerate kinase I |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,271,595 | 1 | . | C | 85.7% | 41.4 / 0.9 | 14 | coding (176/1146 nt) | garK | glycerate kinase I |
| Reads supporting (aligned to +/- strand): ref base . (0/2); new base C (6/6); total (6/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.73e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAG‑CCC‑C‑CTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGC > NC_000913/3271458‑3271711 | cgcgTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCTGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAG‑c < 2:74326/139‑1 (MQ=255) gcgTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAG‑cc < 2:86317/139‑1 (MQ=255) tttttCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAG‑ccc‑c‑c < 2:559278/139‑1 (MQ=255) tCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGACGTGACGTTCAG‑ccc‑c‑cctgg > 2:474714/1‑135 (MQ=255) cAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCATTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAg < 2:392025/113‑1 (MQ=255) cAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAg > 1:392025/1‑113 (MQ=255) cAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCC‑C‑CTGGGTGGCTGCAATCATCGc > 2:213408/1‑139 (MQ=255) cTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCC‑C‑CTGGGTGGCTGCAATCATCGCTTCCAcc > 1:156773/1‑139 (MQ=255) gccgccgccATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCC‑C‑CTGGGTGGCTGCAatcat < 2:103644/126‑1 (MQ=255) gccgccgccATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCC‑C‑CTGGGTGGCTGCAatcat > 1:103644/1‑126 (MQ=255) cgccATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAG‑ccc‑c‑cctggg < 1:398776/108‑6 (MQ=255) cgccATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAG‑ccc‑c‑cctggg > 2:398776/1‑103 (MQ=255) aaCGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCC‑C‑CTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCa < 1:69477/139‑1 (MQ=255) cgcgGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCC‑C‑CTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAAc < 1:327565/139‑1 (MQ=255) ttttGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCC‑C‑CTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGaa < 1:63542/139‑1 (MQ=255) ttctGCCATCGCCCGAGACGCCCCACCTGGCATTCACTCTCTCCCCGAGCGGACCTGTACTCCACACCTGACGTTCAG‑CCCGC‑CTGGGTGGCTGCAATCCTGGCTTCCACAGGGCCTTCGCCACCGTCGGCAACCGGaa < 2:409773/136‑1 (MQ=255) tGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCC‑C‑CTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAg > 2:116450/1‑139 (MQ=255) aCTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAG‑CCC‑CGCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGa < 2:156773/139‑1 (MQ=255) tctcGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCC‑C‑CTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATc > 2:324104/1‑139 (MQ=255) gTAACCCAGGCGTGACGTTCAGCCCC‑C‑CTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTgc < 1:116450/139‑1 (MQ=255) aCGTTCAGCCCC‑C‑CTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAAc > 2:600965/1‑76 (MQ=255) aCGTTCAGCCCC‑C‑CTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAAc < 1:600965/76‑1 (MQ=255) | CGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAG‑CCC‑C‑CTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGC > NC_000913/3271458‑3271711 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |