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 | 81.2% | 31.4 / 5.2 | 17 | intergenic (+74/‑90) | trpS/zapE | tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base . (1/2); new base C (7/6); total (9/8) | |||||||||||
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% | 45.4 / NA | 17 | 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 (9/8); total (9/8) |
CCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑CACGGCTTTGCCGGTGT‑T‑TCGGGCACG‑C‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACG > NC_002947/1499369‑1499644 || ccGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCCGGCCCCCGCAACACGCGCCGCCGGGCAAGGCACCCGCGTTTGCCGGGTT‑G‑GCCGGCAcg‑c‑ccgcgac > 8:12122/1‑141 (MQ=255) aTTCAGAGCCTGAACGAAGCCTTGGCACCACCGCAGGCCTGCGGCGGCCAAACGGGGGGCGCGGGCGGGCCCGCGAAACAGGCGCCGCGGGGCAGGGC‑ACCGGCTTTGCCGGTGTGC‑GCGGGCACGCC‑CGCTCCCACAGACTGATTGc < 3:77382/148‑1 (MQ=21) tGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTACCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGGGCATGGC‑ACCGGCTTTGCCGGTGT‑TCGCGGGCACGCC‑CGCTCCCACAGACTGATTGCCATTAAGCgg < 4:141722/148‑1 (MQ=21) accacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCGGTGT‑TCGCGGGCACGCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTg < 6:28800/148‑1 (MQ=21) aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCGGTGT‑TCGCGGGCACGCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTttg > 3:131519/1‑148 (MQ=21) ggCTGCCAAACGGGGGGCGCGGGCGGGCCCGCTAAACTGGCGCCGCGGGGCAGGGC‑ACCGGGTTTGCCGGTGTGC‑GCGGGCACGCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTttgt < 5:121023/148‑1 (MQ=25) tGCCAAACGGTGGGAGCGGGCGTGGCCCCGCAACCCGCGCCCCGGTGGATGGC‑ACCCGCGTTGCCCGTGT‑TAGCGGGCCCGCC‑CCCGCCCCCACACGGATTGGCAATAAGCCGGCCGATAGCGACAGGGCCCGCCCCGGTTTttgtgg > 2:20509/1‑146 (MQ=255) gcgggcggGCCCGCGAAACAGGCGCCGCGGTGCTTGGC‑ACCGGCTTTGCCGGTGT‑TCGCGGGAACGCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa < 2:118341/148‑1 (MQ=25) ggCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCGGTGT‑TCGCGGGCACGCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa < 2:11781/148‑1 (MQ=37) ccGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCGGTGT‑TCGCGGGCACGCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTccc > 2:155036/1‑146 (MQ=37) ccGCGAAACAGGCGCCGCGGTGCATGGC‑ACAGGCTTTGCGGGTGT‑TCGCGGGCACGCC‑CGATCCCACAGACTGATTGCCATTAAGCTGGCCTATCGCTACAGTT‑ACGGCCCGTTTTTTGTTGCCTCTCTAACGAATCATGACTCCCt < 4:135919/147‑1 (MQ=25) cgcgAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCGGTGT‑TCGCGGGCACGCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTccc > 7:31576/1‑145 (MQ=37) cgcgAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCGGTGT‑TCGCGGGCACGCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTccc > 5:89203/1‑145 (MQ=37) gccgcGGTGCATGGC‑ACCGGCTTTGCCGGTGT‑TCGCGGGCACGCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCaa > 1:37324/1‑148 (MQ=37) cgcgGTGCATGGC‑ACCGGCTTTGCCGGTGT‑TCGCGGGCACGCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGc < 4:131519/148‑1 (MQ=37) gcgGTGCATGGC‑ACCGGCTTTGCCGGTGT‑TCGCGGGCACG‑CACGCTACCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTTTGTTGCCTAGCTAACGAATCATGACTCCCTTAGAACGCTATCaa > 2:104701/1‑145 (MQ=37) ccGGTGT‑TCGCGGGCACGCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtct > 4:83942/1‑148 (MQ=255) gcGGGCACGCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGAcg > 3:172901/2‑148 (MQ=255) || CCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑CACGGCTTTGCCGGTGT‑T‑TCGGGCACG‑C‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACG > NC_002947/1499369‑1499644 |
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: ‑ |