| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 606,927 | C→T | G153S (GGC→AGC) | ybdK ← | weak gamma‑glutamyl:cysteine ligase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 606,927 | 0 | C | T | 96.2% | 79.0 / ‑3.8 | 26 | G153S (GGC→AGC) | ybdK | weak gamma‑glutamyl:cysteine ligase |
| Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (11/14); total (11/15) | |||||||||||
| 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.73e-01 | |||||||||||
CGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTT > NC_000913/606796‑607043 | cGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCTGTgcagc > 2:515007/1‑139 (MQ=255)cGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCTGTgcag > 1:123050/1‑138 (MQ=255)cGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCTGTgcag < 2:123050/138‑1 (MQ=255) cATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCTGTgc > 1:124803/1‑125 (MQ=255) cATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCTGTgc < 2:124803/125‑1 (MQ=255) cATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCTGTGCAGCAAATAAATGGc < 1:391910/139‑1 (MQ=255) tCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCTGTGCAGCAAATAAATGGCGTCAt < 1:515007/139‑1 (MQ=255) gCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCTGTGCAGCAAATAAATGGCGTCATCGCCACTgg > 1:179422/1‑139 (MQ=255) cGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCCTGACAAGCTGTGCAGCAAATAAATGGCGTCATCGCCa < 2:54051/134‑1 (MQ=255) cGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCCTGACAAGCTGTGCAGCAAATAAATGGCGTCATCGCCa > 1:54051/1‑134 (MQ=255) aaaTATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAAGCGTGACAAGCTGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGcc < 1:458681/139‑1 (MQ=255) tATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGc < 2:352015/139‑1 (MQ=255) tCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCTGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAg > 2:549596/1‑100 (MQ=255) tCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCTGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAg < 1:549596/100‑1 (MQ=255) tCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCTGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGa < 2:179422/139‑1 (MQ=255) gCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCTGTGCAGCAAATAAATGGCGTCATCGCCa < 1:78468/77‑1 (MQ=255) gCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCTGTGCAGCAAATAAATGGCGTCATCGCCa > 2:78468/1‑77 (MQ=255) gACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCTGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGt > 1:477469/1‑139 (MQ=255) cggcggAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCTGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGc > 1:246799/1‑81 (MQ=255) cggcggAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCTGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGc < 2:246799/81‑1 (MQ=255) cggcggAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCTGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGt > 2:256755/1‑109 (MQ=255) cggcggAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCTGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGt < 1:256755/109‑1 (MQ=255) gATAAAGTGCGGCACAAATCGTGACAAGCTGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGa > 2:537169/1‑139 (MQ=255) gATAAAGTGCGGCACAAATCGTGACAAGCTGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGa < 2:570931/139‑1 (MQ=255) tAAAGTGCGGCACAAATCGTGACAAGCTGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATa < 1:537169/139‑1 (MQ=255) tGCGGCACAAATCGTGACAAGCTGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGcgtt < 2:477469/139‑1 (MQ=255) | CGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTT > NC_000913/606796‑607043 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |