Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,499,506:1 | +C | 100% | intergenic (+83/‑81) | 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 | 95.8% | 61.2 / ‑2.9 | 24 | intergenic (+83/‑81) | trpS/zapE | tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base . (1/0); new base C (8/15); total (9/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.75e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
TGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCAC‑GGCTTTGCCGGTGTT‑TCGGGCACG‑CC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATC > NC_002947/1499363‑1499617 | tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACACGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGTTAGCGGGCCcg‑cccg > 5:164114/1‑147 (MQ=25) gCCTTTCAGAGCCTGAACGAAGCCTTGGCACCAGAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGGGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGGTTTGCCGGTGTTCGCGGGCACGCCC‑GCTCCCACAGACTGAtt < 8:356912/149‑1 (MQ=21) gCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGGTTTGCCGGTGTTAGCGGGCCCGCCC‑GCTCCCCCAGACTGAtt > 6:302248/1‑149 (MQ=25) gCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGACACC‑GGCGTTGCCGGTGTTCGCGGGCACGCCC‑GCTCCCACAGACTGAtt < 2:77009/149‑1 (MQ=21) ttCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGTTCGCGGGCACGCCC‑GCTCCCACAGACTGATTGCCa < 3:283052/149‑1 (MQ=21) cAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGTTCGCGGGCACGCCC‑GCTCCCACAGACTGATTGCCAt < 2:31751/148‑1 (MQ=21) agCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGTTCGCGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAg < 1:240597/149‑1 (MQ=21) tGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGTTCGCGGGCACGCCC‑GCTCCCCCAGACTGATTGCCATTAAGCggg > 7:281275/1‑149 (MQ=21) ttGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACCGGCGCCGCGGTGCATGGCACC‑CGCGTTGCCGGTGTTAGCGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCGATCGCTa > 2:329290/1‑148 (MQ=21) ttGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGTTCGCGGGCCCGCCC‑GCTCCCCCAGACTGATTTCCAATAAGCGGGCCTATCGCTa > 3:79953/1‑148 (MQ=21) ccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGTTCGCGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAc < 4:245245/149‑1 (MQ=21) ggCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACATGGTCCGCTGTGCATGGCACC‑GGCTTTGCCGGTGTTCGCGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTTACGGCCCGTTTTttgtt < 6:194891/149‑1 (MQ=21) ccAAACGGTGGGAGCGGGCGTGGCCGCGGAACAGGCGCCGCGGTGGATGGCACC‑GGCGTTGGCGGGGTTAGCGGGGCCGCCC‑GCGCCCCCAGCCGGGTTGCCATTAAGCGGGGCGATCGGGACAGTGACGGGCCGGTGTTTGTTGccgc > 6:108785/1‑147 (MQ=255) aaaCGGGGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGTTCGCGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGc < 7:190922/149‑1 (MQ=21) aaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGTTCGCGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGc > 1:269393/1‑148 (MQ=21) ggTGGGAGGGGGCGGGCCCGCTAAACAGGCGCCGCGGGGCATGGCACG‑GGCTTTGCCGGTGTTCGCGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAAc < 3:381698/149‑1 (MQ=25) ggAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGTTCGCGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAAt < 5:302248/149‑1 (MQ=25) gcgCGGGCGTGCCCGCTAAACCGGCGCGGCGGCGCAGGGCACC‑GGCTTTGCCGGTGTTCGCGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATc < 8:289551/147‑1 (MQ=25) gggCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑TGCTTTGCCGGTGTTCGCGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa < 7:297403/149‑1 (MQ=37) gcgAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGTTCGCGGGCACG‑CCGGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCACCTTTTTTTTTTCCTCGCTAACTAATCATGACTCCCTTAg > 5:328470/1‑148 (MQ=25) aaaCAGGCGCCGCGGTGCA‑GGCCACGGGGTTTGCCGTTGTTCGCGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACg < 5:15275/149‑1 (MQ=255) aCAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGTTCGCGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCt < 6:164114/149‑1 (MQ=37) aGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGTTCGCGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTAt > 3:43523/1‑149 (MQ=37) gccgcGGTGCATGGCACC‑GGCTTTGCCGGTGTTCGCGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg < 7:41217/149‑1 (MQ=37) gTGCATGGCACC‑GGCTTTGCCGGTGTTCGCGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAACCAGATc > 1:350054/1‑149 (MQ=37) | TGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCAC‑GGCTTTGCCGGTGTT‑TCGGGCACG‑CC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATC > NC_002947/1499363‑1499617 |
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: ‑ |