| 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 | 85.7% | 31.5 / 1.3 | 15 | 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 (5/7); total (6/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.05e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01 | |||||||||||
| * | NC_002947 | 1,499,498 | 0 | T | G | 100.0% | 36.2 / NA | 15 | 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 (6/9); total (6/9) | |||||||||||
AGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGT‑T‑TCGGGCACG‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAC > NC_002947/1499399‑1499631 || aGCCTTGGCACCACAGCAGGACTGAGGCGGCCAAACGGGGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATc < 6:126238/148‑1 (MQ=21) ttGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTTCCGGTTT‑TAGCGGGCACGCCCGCTACCACAGCATGATTTGCATTAAGCGGGCCTAtaggt > 4:104747/1‑143 (MQ=18) ggAGGCTGCCAAACGGGGGGAGCGGGCGGGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTGC‑GCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGttttt < 7:68955/147‑1 (MQ=25) gCTGCCAAACGGTGGGCGCGGGCGGGCCCGCGAAACCGGCGCCGCGGGGCAGGGCACCGGCTTTGCCGGTGTGC‑GCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgtt < 1:21800/148‑1 (MQ=25) cGGGCGTGCCGGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt < 4:69775/148‑1 (MQ=25) ggCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa < 2:2555/148‑1 (MQ=37) cgcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTccc > 7:141634/1‑145 (MQ=37) cATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCt > 7:99683/1‑147 (MQ=38) tggcaccGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa < 3:150193/142‑1 (MQ=38) ggcaccGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaaa > 4:74849/6‑148 (MQ=38) ggcaccGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaaa < 8:99683/143‑1 (MQ=38) accGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGt < 5:41043/146‑1 (MQ=38) cGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTcc > 6:74716/1‑148 (MQ=39) cGGCTTGGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTcc < 8:66783/148‑1 (MQ=39) tttGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAc > 6:133506/1‑148 (MQ=39) || AGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGT‑T‑TCGGGCACG‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAC > NC_002947/1499399‑1499631 |
| 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: ‑ |