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 | 90.9% | 51.8 / ‑3.5 | 22 | intergenic (+57/‑107) | trpS/zapE | tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base C (2/0); new base A (8/12); total (10/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.95e-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% | 64.3 / NA | 22 | 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 (10/12); total (10/12) |
ACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑CACGGCTTTGCCGGTGT‑T‑TCGGGCACG‑C‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGC > NC_002947/1499341‑1499612 || aCGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCttt < 1:136591/148‑1 (MQ=255) gCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCg > 8:31988/1‑148 (MQ=255) gCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCg < 6:72774/148‑1 (MQ=255) cGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGGCCGCGAAACACGCCCCGCGGTGCATGGC‑CCCCGCGTTGCCGgt > 5:41813/1‑148 (MQ=255) gCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCGgtg < 8:49097/148‑1 (MQ=255) cAGTTCGACAACGCTGAAGCCCGTGAGGCGGCCATTCAGAGCCTGAACGACGCGTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCGgtgt < 1:180028/148‑1 (MQ=38) 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 > 5:89203/1‑145 (MQ=37) cgcgAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCGGTGT‑TCGCGGGCACGCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTccc > 7:31576/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) || ACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑CACGGCTTTGCCGGTGT‑T‑TCGGGCACG‑C‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGC > NC_002947/1499341‑1499612 |
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: ‑ |