| 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 | 92.3% | 23.1 / ‑4.0 | 14 | intergenic (+74/‑90) | trpS/zapE | tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein |
| Reads supporting (aligned to +/- strand): ref base . (0/1); major base C (3/9); minor base A (1/0); total (4/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.08e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01 | |||||||||||
| * | NC_002947 | 1,499,498 | 0 | T | G | 100.0% | 37.4 / NA | 14 | 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 (4/10); total (4/10) | |||||||||||
GGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGT‑T‑TCGGGCACG‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAA > NC_002947/1499378‑1499621 || ggCCATTCAGTGCGTGAACGAATCGTTGGAACCACAGCAGGACTGAGGCTGCCAACCGGTGTTAGCGGGCGTTCCCGCGAAAAAGGCGCCGCGGTGCATGGCACGGGCTTTGCGGGTGT‑TCGGGGGCACGCCCGCTCCCACAGACTGa < 6:176355/148‑1 (MQ=255) agagCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGCAACAGGCGCCGCGGTGGATGGCCCCCGCGTTGGCCGTGT‑TAGCGGGCCCGCCCGCGCCCCCAGACGGCTTGccct > 3:111276/1‑145 (MQ=21) gCACAACAGCAGGACTGCGGCTGCCAAAGGGTGGTAGCGGGCGGGCCCGCTAAACAGGCGCCGCGGTGCATGGCACGGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACa < 2:125656/147‑1 (MQ=21) cTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTGC‑GCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttt < 6:91148/148‑1 (MQ=21) gggtGCCAAACGGGGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGGGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAATCGGGCCTATCGCTACAGTGACGGCCCGtttttt < 4:36906/143‑1 (MQ=21) ggCTGCCTAACGGTGGGTGGGGGCGTGCCCGCGAAACAGGCGCCGCGGTGAATGGCACCGGCTTTGCGGGTGT‑TCGCGGGCAGGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGTCGGCCCGtttttt < 4:34456/146‑1 (MQ=255) gAGCGGGCGTGCCCGCTAAACAGGCGCCGCGTTGAATGGCACGGGCTTTGCCGGTGT‑TCGCGGGCAGGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAAt < 2:101386/148‑1 (MQ=25) ggCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa < 1:102017/148‑1 (MQ=37) tGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCGGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGTCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTcc < 4:111276/148‑1 (MQ=25) ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGGTTTTTGTTGCCTCGCTAACGAATCATGACTccc > 5:112332/1‑146 (MQ=25) aGGCGCCGCGGTGCATGGCACCGGCTTTGCGGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTa < 6:27750/148‑1 (MQ=37) tGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATc < 4:54331/148‑1 (MQ=38) gCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGTGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCGAACGAATCATGACTCCCTTAGAACGATATCAAGCAGATc > 6:23183/1‑147 (MQ=37) tggcaccGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCACTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa > 4:26856/7‑148 (MQ=37) || GGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGT‑T‑TCGGGCACG‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAA > NC_002947/1499378‑1499621 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |