| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_002947 | 1,499,498 | 1 bp→CG | intergenic (+75/‑89) | trpS → / → zapE | tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 1,499,497 | 1 | . | C | 90.0% | 49.5 / 0.4 | 20 | intergenic (+74/‑90) | trpS/zapE | tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein |
| Reads supporting (aligned to +/- strand): ref base . (0/2); new base C (13/5); total (13/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.11e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| * | NC_002947 | 1,499,498 | 0 | T | G | 100.0% | 51.8 / NA | 20 | intergenic (+75/‑89) | trpS/zapE | tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein |
| Reads supporting (aligned to +/- strand): ref base T (0/0); new base G (13/7); total (13/7) | |||||||||||
AACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGT‑T‑TCGGGCACG‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCC > NC_002947/1499394‑1499638 || aaCGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCCCGCCCGCTCCCACAGACTGATTGCCATTAAGCggg > 5:214545/1‑147 (MQ=21) aCGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGGGTGGCCGTTGT‑TCGCGGGCGCGCCCGGGCCCCCGGACTGTATGCGATTAAGCgggg > 6:57985/1‑146 (MQ=16) gcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCTTGGCACCGGCGTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACg < 5:237849/145‑1 (MQ=21) ggACTGAGGCTGCCAAACGGGGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCg < 3:228914/147‑1 (MQ=21) gCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGCAACAGGCGCCCCGGGGCATGGCACCCGCGTTGCCGGGGT‑TCGCCGGCACGCCCGCACCCACCGACCGGTTGCCCTTTAGCCGGGCTATCGCTTCAAGGTCGGCCCCGGTTttgt > 6:234599/1‑147 (MQ=9) ggTGGGAGCGGGGGTGCCCGGAAAACTGGCGCGGGGGTGCTTGGAACGGGCTTTGCCGGTGTCC‑GCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTa < 5:165867/147‑1 (MQ=25) gTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCTGTGCATGGCACGGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTaa < 1:37185/147‑1 (MQ=21) gCGTGCCCGGGTAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCGGGTGTGC‑GCGGGCACGCCCGCTCGCACAGACTGATTGCGATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGa < 2:74627/147‑1 (MQ=25) cGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGAc > 8:239299/1‑147 (MQ=25) gCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAATCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGACtac > 8:142525/1‑145 (MQ=25) gCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGACTcc > 4:106497/1‑147 (MQ=37) cgcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGACTccc > 5:113856/1‑145 (MQ=37) cgcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGACTccc > 5:18978/1‑145 (MQ=37) cgcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGACTccc > 2:162662/1‑145 (MQ=37) cgccgcGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACCGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGACTCCCTTAGAACGCTAt > 6:13352/1‑146 (MQ=37) caccGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAAc < 1:163909/144‑1 (MQ=38) accGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACTATACATCACTCCCTTAGAACGCTATCAAGCTGATCTGAAACg > 2:116843/3‑147 (MQ=37) ccGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGt > 7:9432/2‑147 (MQ=39) ggCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTcc > 1:252510/1‑147 (MQ=39) tttGCCGGTGT‑TCGCGGGCAAGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCTGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCg < 1:74893/147‑1 (MQ=38) gtgttcgcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACttcttc < 2:252510/140‑1 (MQ=255) tgttcgcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTcc > 8:21250/7‑147 (MQ=255) || AACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGT‑T‑TCGGGCACG‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCC > NC_002947/1499394‑1499638 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |