Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,499,498 | 1 bp→CG | 100% | 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.7% | 57.2 / 4.9 | 31 | intergenic (+74/‑90) | trpS/zapE | tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base . (0/0); major base C (10/15); minor base A (6/0); total (16/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.77e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
* | NC_002947 | 1,499,498 | 0 | T | G | 100.0% | 85.8 / NA | 31 | 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 (16/15); total (16/15) | |||||||||||
Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
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: ‑ |