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 | 88.9% | 18.0 / NA | 10 | 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 (5/3); minor base A (1/0); total (6/4) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
* | NC_002947 | 1,499,498 | 0 | T | G | 100.0% | 23.5 / NA | 10 | 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/4); total (6/4) |
CCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCAC‑GGCTTTGCCGGTGTT‑TCGGG‑CACG‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCG > NC_002947/1499380‑1499629 || ccATTCTTAGCCTGAACGAAGCCTTGGCACCACAGCAGTACTGTGGCTGCCAAACGTTGGTAGCGGGCGTGCCCGCTAAACAGGCGCCGCGGTGCATGGCACG‑GGCTTTGCGGGTGTTCGCGGG‑CACGCCCGCTCCCACAGACTGAtt < 4:53782/148‑1 (MQ=21) gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGGCCGCGCAACCCGGGCCGCGGGGCATGGCACC‑CGCGTTGGCCGTGTTCGCGGG‑CCCGCCCGCGCCCACAGCCGGATTGGCCTTAAGCggggcga > 3:28587/1‑144 (MQ=255) agcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCAGGGCACC‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGc < 7:25786/148‑1 (MQ=21) gCCAAACGGTGGCCGCGGGCGGGCCGGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTGGGCCGGGGTTCGCGGG‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGcc < 3:106263/148‑1 (MQ=25) aaaCGGTGGGAGCGGGCGTGCCCGCGAAACACGCGCCCCCGTGCATTGCAAC‑GGCGTTTCCCGTGTTAGCGGG‑CACGCCCCCGCCCACAGACAGATTGCCATTAAACCGGCCTATCGCGAAACTGAACGCCCCTTTTTTTTTTCCTCg > 8:42723/1‑148 (MQ=12) ggggggggCCGGCTAAACTGGCGCGGGGGGGC‑TGGCCCCGGGTTTTGCGGGTGTG‑GGGGGCCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt < 2:143406/141‑1 (MQ=255) aaCAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCGTATCGCTACAGTGACGGCGCGTTTTTTGTTGCCTCGCTAACGAATCCTGCCGCCCTTAGAACg > 3:27585/1‑148 (MQ=25) gTGCATGGCACC‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAt > 2:115448/1‑148 (MQ=38) gTGCATGGCACC‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAt > 2:97939/1‑148 (MQ=38) gCTTTGCCGGTGTTCGCGGG‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCg > 6:20681/1‑148 (MQ=39) || CCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCAC‑GGCTTTGCCGGTGTT‑TCGGG‑CACG‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCG > NC_002947/1499380‑1499629 |
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: ‑ |