Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_002947 | 1,499,508 | +C | intergenic (+85/‑79) | 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,506 | 1 | . | C | 91.7% | 58.8 / 1.3 | 24 | intergenic (+83/‑81) | trpS/zapE | tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base . (2/0); new base C (12/10); total (14/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.93e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01 |
AAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCC‑GGTGTT‑TCGGGCA‑CG‑C‑CGC‑TCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAC‑GGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACG > NC_002947/1499365‑1499644 | aaGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCA‑cg‑c‑ccg < 5:163613/147‑3 (MQ=25) ggCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTAGCGGGCA‑CGCC‑CGC‑TCCCACAGACTGAt > 2:44738/1‑149 (MQ=25) gagCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCA‑CGCC‑CGC‑TCCCACAGACTGATTGCCAtt > 5:189673/1‑147 (MQ=21) tGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCTAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCA‑CGCC‑CGC‑TCCCACAGACTGATTGc < 1:323676/138‑1 (MQ=21) tGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCC‑CGCC‑CGC‑GCCCACAGACAGATTGc > 2:323676/1‑138 (MQ=18) aaCGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGGCCGCGAAACCGGCGCCGCGGTGCATGGCACCCGGGTTGCC‑GGTGTTCGCGGGGC‑CGCC‑CGC‑GGCCCCAGACGGATTGGCATTTAGCgggga > 1:315193/1‑147 (MQ=255) gACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCCCAACCGGCGCCCCCGGGCATGGCACCCGCGTTGGC‑GGTGTTAGCGGGCACCC‑C‑CGC‑GCCCCCACACAGAATTCCAATAAGCCGGCCCATCGCGACACTGGCGGGCCCGtt > 8:283524/1‑149 (MQ=255) gACTGAGGCTGCAAAACTGTGGGAGCGGGCGTGTCCGCGAAACATGCGCCGCGGTGCATGGCACAGGCTTTGCC‑GGTGTTCGCGGGCA‑CG‑CACGC‑TCCCAAAGACAGATTGCCATTAAGCGGGCCTATCGCTACAGTGAC‑GGACCGgtt > 8:247449/1‑149 (MQ=18) ggCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCA‑CGCC‑CGC‑TCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAC‑GGCCCGTTTTttgtt < 1:44738/149‑1 (MQ=21) gCCAACCGGGGGGAGCGGGCGTCCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCA‑CGCC‑CGC‑TCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAC‑GGCCCGTTTTTTGTTGCCt < 7:283524/149‑1 (MQ=21) cGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCA‑CGCC‑CGC‑TCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAC‑GGCCCGTTTTTTGTTGCCTCGCTaa < 5:187522/149‑1 (MQ=21) atcgggAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCA‑CGCC‑CGC‑TCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAC‑GGCCCGTTTTTTGTTGCCTCGCTAAc < 7:152715/146‑1 (MQ=21) gCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCA‑CGCC‑CGC‑TCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAC‑GGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt < 3:37650/149‑1 (MQ=25) gggCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCA‑CGCC‑CGC‑TCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAC‑GGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa > 7:71432/1‑149 (MQ=25) gCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCA‑CGCC‑CGC‑TCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAC‑GGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACt > 5:51365/1‑149 (MQ=25) gTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCA‑CGCC‑CGC‑TCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAC‑GGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTcc < 6:189673/149‑1 (MQ=25) cccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCA‑CGCC‑CGC‑TCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAC‑GGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt > 7:201223/1‑148 (MQ=37) aGGCGCCGCGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCA‑CGCC‑CGCGTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAC‑GGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTa > 1:303119/1‑149 (MQ=37) ggtgccgcGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCA‑CGCC‑CGC‑TCCCACAGACTGATTGCCATTAAGCGGGc > 1:74105/4‑80 (MQ=14) ggtgccgcGGTGCATGGCACCGGCTTTGCC‑GGTGTTCGCGGGCA‑CGCC‑CGC‑TCCCACAGACTGATTGCCATTAAGCGGGc < 2:74105/77‑1 (MQ=14) cgggCTTTGCCGGGTTTC‑TCGGTCA‑CGCC‑CGC‑TCCCACAGACTGATTGCCTTTAAGCGGGACTATCTCTACAGTGAC‑TGCCAGTTTTTTTTTGCCTCGCTAACGAATCATTACTCCCTTAGAACGCTATCAATCAGATCTGAAACGTcc < 7:247449/147‑1 (MQ=255) cc‑GGTGTTCGCGGGCA‑CGCC‑CGC‑TCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAC‑GGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaaa < 4:368876/135‑1 (MQ=255) cc‑GGTGTTCGCGGGCA‑CGCC‑CGC‑TCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAC‑GGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaaa > 3:368876/1‑135 (MQ=255) ttcgcGGGCA‑CGCC‑CGC‑TCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAC‑GGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGa > 7:299358/5‑149 (MQ=255) cgcGGGCA‑CGCC‑CGC‑TCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAC‑GGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGAcg > 8:247636/3‑149 (MQ=255) | AAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCC‑GGTGTT‑TCGGGCA‑CG‑C‑CGC‑TCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAC‑GGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACG > NC_002947/1499365‑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: ‑ |