Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_002947 | 1,499,498 | 1 bp→CG | 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 | 91.4% | 85.5 / ‑3.7 | 37 | intergenic (+74/‑90) | trpS/zapE | tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base . (0/1); major base C (12/20); minor base A (3/0); total (15/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.95e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
* | NC_002947 | 1,499,498 | 0 | T | G | 100.0% | 102.9 / NA | 37 | 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 (15/22); total (15/22) |
CCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGT‑T‑TCGGGCA‑CG‑‑‑CC‑GCTCCCACAG‑ACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGA > NC_002947/1499369‑1499642 || ccGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGCAACAGGCGCCGCGGTGCATGGCACCGGCGTTGCCGGTGT‑TCGCGGGCACCC‑‑‑CC‑GCTcccc > 6:427297/1‑148 (MQ=25) aTGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCA‑CGC‑‑CC‑GCTCCCACAG‑Ac < 4:113173/149‑1 (MQ=21) gCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCA‑CGC‑‑CC‑GCTCCCACAG‑ACTGAtt < 6:307836/149‑1 (MQ=21) ttCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCA‑CGC‑‑CC‑GCTCCCACAG‑ACTGATTGCCa < 2:47151/149‑1 (MQ=21) ccTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAACGCGCCGCGGTGCAAGGCACCGGCCTTGCCCGTGT‑TCGCGGGCA‑CGC‑‑CC‑CCGCCCACCG‑ACGGATTGCCATTAAGCg > 3:322412/1‑149 (MQ=21) ccTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCCCGAAACACGCGCCCCCGTGCATGGCACCGGCGTTGCCGGTGT‑TAGCGGGCC‑CGC‑‑CC‑GCGCCCACAA‑ACTGATTGCCAATAACCg > 8:382161/1‑149 (MQ=255) aaCGAAGCCTTGGCACCACAGCAGGACTGCGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTGC‑GCGGGCA‑CGC‑‑CC‑GCTCCCACAG‑ACTGATTGCCATTAAGCGGGcc < 7:382161/149‑1 (MQ=21) gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAGGCGCCCCCGTGCAAGGCACCCGCATTGCCGGTGT‑TAGCGGGCACCC‑‑‑CC‑CCACCCACAC‑ACAGAATTCCAATAATCggggcaa > 6:376745/1‑144 (MQ=16) aaGCCTTGGCAACACAGCAGGACTGAGGCTGCCAAACGGTGGGAGAGGGCCTAACCCCGAAACACGCGCCCCCGTGAAATGCAACAGCATTGCCGGTGT‑TCGCGGGCA‑CGC‑‑CC‑GCTACCAAAA‑ACTGGTTGCCATTAATCGAGCCAATc > 6:112311/1‑149 (MQ=255) aGCCTGGGCACCACAGCAGGCCTGAGGCTGCCAACCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCA‑CGC‑‑CC‑GCTCCCACAG‑ACTGATTGCCATTAAGCGGGCCTATCg < 3:44334/149‑1 (MQ=18) gCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGGGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCA‑CGC‑‑CC‑GCTCCCACAG‑ACTGATTGCCATTAAGCGGGCCTATCGc < 7:9466/149‑1 (MQ=21) cagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCA‑CGC‑‑CC‑GCTCCCACAG‑ACTGATTGCCATTAAGCGGGc < 2:66569/127‑1 (MQ=18) cagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TAGCGGGCA‑CGC‑‑CC‑GCTCCCACAG‑ACTGATTGCCATTAAGCGGGc > 1:66569/1‑127 (MQ=18) aGGCTGCCAAACGGTGGGAGCGGGCGTGGCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCA‑CGC‑‑CC‑GCTCCCACAG‑ACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttg < 5:376745/148‑1 (MQ=21) cTGCCAAACGGTGGTTGCGGGCGTTCCCGCTAAACAGGCGCCGCTGTGCATGGCACGGGCTTTGCGGGTTT‑TCGCTGGAA‑CGC‑‑CC‑TCTCCCACAG‑ACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGc < 1:371091/149‑1 (MQ=255) gCCAAACGGTGGGTGGGGGCGTGCCCGCTAAACAGGCGCGGCGGTGCTTGGCACCGGCTTTGCCGGTGT‑TCGGGGGCA‑CGC‑‑CC‑GCTCCCACAGCACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt < 1:97892/143‑1 (MQ=21) ggggggCGCGGGCGTGCCCGCGAAACAGGGGCCGCGGTGCATGGGACCGGGTTTGCCGGTGT‑TGGCGGGCA‑CGCCCCC‑GCTCCCACAG‑ACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttt < 1:323251/132‑1 (MQ=21) gggCGTGCCCGCTAAACTGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCA‑CGC‑‑CC‑GCTCCCACAG‑ACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa < 3:349226/149‑1 (MQ=37) gggCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCA‑CGC‑‑CC‑GCTCCCACAG‑ACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa < 4:322412/149‑1 (MQ=37) cccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCA‑CGC‑‑CC‑GCTCCCACAG‑ACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATc < 2:2673/138‑1 (MQ=25) cccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCA‑CGC‑‑CC‑GCTCCCACAG‑ACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATc > 1:2673/1‑138 (MQ=25) cccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCA‑CGC‑‑CC‑GCTCCCACAG‑ACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt < 8:359184/140‑1 (MQ=25) cccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCA‑CGC‑‑CC‑GCTCCCACAG‑ACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt > 7:359184/1‑140 (MQ=25) ccGCTAAACTGGCGCCGCGGTGCATGGCCCCGGCTTTGCGGGTGT‑TCGCGGGCA‑CGC‑‑CC‑GCTCCCACAG‑ACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa < 2:144332/149‑1 (MQ=37) aCAGGCGCCGCGGTGCATGGCACCGGCTGTGCCGGTGT‑TCGCGGGCA‑CGC‑‑CC‑GCTCCCACAG‑ACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCt < 5:427297/149‑1 (MQ=37) cAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCA‑CGC‑‑CC‑GCTCCCACAG‑ACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTa > 5:427143/1‑149 (MQ=37) cAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCA‑CGC‑‑CC‑GCTCCCACAG‑ACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTa > 1:347758/1‑149 (MQ=37) gcgcCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCA‑CG‑‑‑CCAGCTCCCACAG‑ACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc > 6:392407/1‑148 (MQ=38) cgcgGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCA‑CGC‑‑CC‑GCTCCCACAG‑ACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCa < 6:427143/149‑1 (MQ=37) ggcaccGGCTTTGCCGGTGT‑TCGCGGGCA‑CGC‑‑CC‑GCTCCCACAG‑ACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAAt > 2:212861/6‑111 (MQ=25) ggcaccGGCTTTGCCGGTGT‑TCGCGGGCA‑CGC‑‑CC‑GCTCCCACAG‑ACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAAt < 1:212861/106‑1 (MQ=25) caccGGCTTTGCCGGTGT‑TCGCGGGCA‑CGC‑‑CC‑GCTCCCACAG‑ACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGt > 1:102602/4‑149 (MQ=38) ccGGCTTTGCCGGTGT‑TCGCGGGCA‑CGC‑‑CC‑GCTCCCACAG‑ACTGATTGCCAGTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTcc > 7:67810/2‑149 (MQ=38) cGGCTTTGCCGGTGT‑TCGCGGGCA‑CGC‑‑CC‑GCTCCCACAG‑ACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAt > 1:191376/1‑137 (MQ=38) cGGCTTTGCCGGTGT‑TCGCGGGCA‑CGC‑‑CC‑GCTCCCACAG‑ACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAt < 2:191376/137‑1 (MQ=38) ccGGTGT‑TCGCGGTCA‑CGC‑‑CC‑GCTCCCACAG‑ACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtctt < 3:251264/149‑1 (MQ=255) ccGGTGT‑TCGCGGGCA‑CGC‑‑CC‑GCTCCCACAG‑ACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtctt < 2:347758/149‑1 (MQ=255) ttcgcGGGCA‑CGC‑‑CC‑GCTCCCACAG‑ACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGa < 8:67810/145‑1 (MQ=255) || CCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGT‑T‑TCGGGCA‑CG‑‑‑CC‑GCTCCCACAG‑ACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGA > NC_002947/1499369‑1499642 |
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: ‑ |