Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_002947 | 1,499,480 | 2 bp→AC | intergenic (+57/‑106) | 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,480 | 0 | C | A | 95.5% | 125.0 / ‑6.2 | 44 | intergenic (+57/‑107) | trpS/zapE | tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base C (1/1); new base A (13/29); total (14/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.40e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
* | NC_002947 | 1,499,481 | 0 | A | C | 100.0% | 131.6 / NA | 44 | intergenic (+58/‑106) | trpS/zapE | tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base A (0/0); new base C (14/30); total (14/30) |
GGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTT‑TCGGG‑‑‑CACG‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAG > NC_002947/1499344‑1499611 || ggCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCg < 6:302991/149‑1 (MQ=255) ggCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCg < 3:144286/149‑1 (MQ=255) gCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGCAACAGGCGCCGCGGTGCATGGCCCCGGCGTTGCCgg > 4:88206/1‑149 (MQ=255) ccGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgt < 3:385992/149‑1 (MQ=255) cGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgtg < 4:302192/149‑1 (MQ=255) cAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGCACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGCGCGGGCGTGCCCGCGAAACAGGCGCCGCGGGGCATGGCACCGGCGTTGCCGGGGtt > 6:105635/1‑149 (MQ=255) cAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgtgt > 5:115743/1‑148 (MQ=255) cGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcg < 4:261195/148‑5 (MQ=255) cGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcg < 4:11108/148‑5 (MQ=255) cGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CAcg‑c < 2:37615/149‑1 (MQ=37) gCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCCTTAACTCAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CAc < 1:137929/149‑1 (MQ=21) cTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CAcg‑ccc < 3:143505/149‑2 (MQ=37) gCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTcc > 7:342386/1‑149 (MQ=25) cGGCCATTCAGAGCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CAcg‑cccg > 5:399210/1‑133 (MQ=25) cGGCCATTCAGAGCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CAcg‑cccg < 6:399210/135‑3 (MQ=25) ggCCATTCAGAGCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CAcg‑cccgctc < 7:376523/137‑6 (MQ=21) gCCATTCAGAGCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGAtt < 5:314660/149‑1 (MQ=21) cATTCATAGCCTGAAGGAAGCCTTGGC‑‑‑ACCACAGCTGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CAGGCCCGCTCCCACAGACTGATTGc < 6:115743/149‑1 (MQ=18) agagCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgt < 7:318743/109‑1 (MQ=255) agagCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgt > 8:318743/1‑109 (MQ=255) agagCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCAGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CAcg‑cccgct < 1:311837/128‑5 (MQ=21) agCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgg > 1:422006/1‑110 (MQ=255) agCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgg < 2:422006/115‑6 (MQ=255) cTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCgg > 7:321540/1‑149 (MQ=21) gAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCg > 4:157098/1‑100 (MQ=255) gAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCg < 3:157098/100‑1 (MQ=255) aaCGAAGCCTTGGC‑‑‑ACCACACCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGTGc < 2:63656/148‑1 (MQ=21) gC‑‑‑ACCACAGCAGGACTGAGGCTGCCAACCGGTGGGAGCGGGCGTGCCCGCTAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGt < 3:133593/149‑1 (MQ=21) ccacAGCAGGACTGAGGCTGCAAAAGGGGGGCTGCGGGCGTGCCCGCTAAACTGGCGCCGCGGTGCAGGGCACGGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAc < 4:414339/149‑1 (MQ=21) acaGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACgg > 5:427971/1‑149 (MQ=21) aGGCGGCCAAACGGGGGGAGCGGGCGGCCCCGCGAAACAGGCGCCGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt < 3:179131/147‑1 (MQ=21) gCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgtt < 2:256361/148‑1 (MQ=21) cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGc < 8:342386/149‑1 (MQ=21) gCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgtgt < 3:274761/69‑1 (MQ=255) gCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgtgt > 4:274761/1‑69 (MQ=255) gggAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGaa < 3:310004/149‑1 (MQ=25) gCCCGCCAAACCGGCGGCGGGGTGCCCGGCCCGGGGTTTGGGGGTGTG‑GCGGGGCCCCCC‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGAcccct < 2:194897/149‑2 (MQ=255) aCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCt < 6:427971/149‑1 (MQ=37) aCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCt < 8:321540/149‑1 (MQ=37) cAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTa > 3:346028/1‑149 (MQ=37) gcgcCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACa > 6:292970/1‑90 (MQ=17) gcgcCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACa < 5:292970/90‑1 (MQ=17) gccgcGGTGCATGGCACGGGCTTTGCCGGTGTTCGGGGG‑‑‑AACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg < 1:399416/149‑1 (MQ=255) ccgcGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg > 5:392424/1‑148 (MQ=37) || GGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTT‑TCGGG‑‑‑CACG‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAG > NC_002947/1499344‑1499611 |
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: ‑ |