| 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 | 93.8% | 84.5 / 0.7 | 32 | 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 (14/16); total (16/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.84e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01 | |||||||||||
| Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
CTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTT‑TCGGGCA‑CG‑CCGCTCCCAC‑AGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCAT > NC_002947/1499391‑1499640 | cTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCA‑CGCCCGCTCCCAC‑AGACTGATTGCCATTAAGCgg < 8:50957/149‑1 (MQ=21)cTGAACGAAGCCTTGGCACCACAGCAGGAATGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCCCGAAAAACGCGCCCCGGTTCATTGCACCGGATTTGCCGGTGTTAGAGGGCA‑CGCCCGCGCCCAC‑ATCCTGAATGCCATTTATCgg > 8:157735/1‑149 (MQ=255) tGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCGTTGCCGGTGTTAGCGGGCA‑CGCCCGCTCCCAC‑AGACTGATTGCCATTAAGCggg > 2:241822/1‑149 (MQ=25) tGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCCTGCCCCCGAAACAGGCGCCGCCGTGCCTGGCCCCCGCGTTGGCGGTGTTAGCGGGCACCC‑CCCCGCCCAC‑AGCCGGATTGCCATTAACCggg > 7:276185/1‑149 (MQ=255) agcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTTCCCGCGAAACAAGCGCCGCCGTGCAAGGCACCCGCTTTGGCCGTGTTAGCGGGCACCC‑CCCCGCCCCC‑CGCCAGAATTCCAATTAGCGGGCCGATCGCGAAAGTggccggc > 6:283306/1‑143 (MQ=255) tGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCA‑CGCCCGCTCCCAC‑AGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt > 1:95204/1‑149 (MQ=21) ggAGGCTGCCTATCGGTGGGTGCGGGGGGGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCA‑CGCCCGCTCCCAC‑AGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGttttt < 7:157735/147‑1 (MQ=21) aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCA‑CGCCCGCTCCCACTTGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt < 1:237935/148‑1 (MQ=21) ggCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCA‑CGCCCGCTCCCAC‑AGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgtt < 2:252945/149‑1 (MQ=21) ggCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCA‑CGCCCGCTCCCAC‑AGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgtt < 5:283306/149‑1 (MQ=21) aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCCCGGTGCATGGCACCGGCGTTGCCCGTGTTAGCGGGCC‑CGCCCCCTCCCAC‑AGACTGATTTCCATTAAGCCGGCCGATCGCTACAGTGACGGCCCGGTTTTTGTTTCCTccc > 6:246906/1‑147 (MQ=18) tGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCA‑CGCCCGCTCCCAC‑AGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGa < 5:246906/149‑1 (MQ=21) gggAGCGGGCGTGCCCGGTAAACAGGCGCGGCTGTGCTTGGCACCGGCTTTGCCGTTGTTCGCGGGCA‑CGCCCGCTCCCAC‑AGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGaa < 5:191783/149‑1 (MQ=25) ggAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCA‑CGCCCGCTCCCAC‑AGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAAt < 3:284854/149‑1 (MQ=25) taGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCA‑CGCCCGCTCCCAC‑AGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATc < 8:19517/148‑1 (MQ=25) aGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTAGCGGGCA‑CGCCCGCTCCCAC‑AGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa > 3:271784/1‑149 (MQ=25) gcgggcggGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTTTTCGCGGGCA‑CGCCCGCTCCCAC‑AGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt < 7:121726/149‑1 (MQ=25) gCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCA‑CGCCCGCTCCCAC‑AGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt < 7:71391/149‑1 (MQ=25) cGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCA‑CGCCCGCTCCCAC‑AGACTGATTGCCATTAAGCGGGc > 6:242733/1‑99 (MQ=14) cGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCA‑CGCCCGCTCCCAC‑AGACTGATTGCCATTAAGCGGGc > 6:238461/1‑99 (MQ=14) gTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCA‑CGCCCGCTCCCAC‑AGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTcc < 5:135106/149‑1 (MQ=25) gCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCA‑CGCCCGCTCCCAC‑AGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt > 8:160113/1‑149 (MQ=37) ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCA‑CGCCCGCTCCCAC‑AGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt > 7:13180/1‑147 (MQ=37) ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCA‑CGCCCGCTCCCAC‑AGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa < 2:95204/149‑1 (MQ=37) gcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCA‑CGCCCGCTCCCAC‑AGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAg > 4:170815/1‑148 (MQ=37) gAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCA‑CGCCCGCTCCCAC‑AGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAAc > 5:308662/1‑149 (MQ=37) ccgcGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCA‑CGCCCGCTCCCAC‑AGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGc < 3:170815/149‑1 (MQ=37) gCATGGCACCGGCTTTGCCGGTGTTCGCGGGCA‑CGCCCGCTCCCAC‑AGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGTCCGGTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCt > 1:134568/1‑148 (MQ=37) caccGGCTTTGCCGGTGTTCGCGGGCA‑CGCCCGCTCCCAC‑AGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa < 6:60238/111‑1 (MQ=37) caccGGCTTTGCCGGTGTTCGCGGGCA‑CGCCCGCTCCCAC‑AGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa > 5:60238/4‑114 (MQ=37) ggCTTTGCCGGTGTTCGCGGGCA‑CGCCCGCTCCCAC‑AGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTccc > 6:44884/1‑148 (MQ=39) tgttcgcGGGCA‑CGCCCGCTCCCAC‑AGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCAt < 7:160113/143‑1 (MQ=255) | CTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTT‑TCGGGCA‑CG‑CCGCTCCCAC‑AGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCAT > NC_002947/1499391‑1499640 |
| 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: ‑ |