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 | 93.8% | 89.3 / 0.6 | 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 (15/15); total (17/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.86e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01 |
AAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGT‑‑GTT‑TCGGGCA‑CG‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCG > NC_002947/1499365‑1499649 | aaGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCA‑cg‑cccgc < 3:250800/148‑4 (MQ=25) cccGTGATGCGGCCATTCAGAGCCTGCACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAATCGGTGTGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCA‑CGCCCGCTccc < 5:175979/149‑1 (MQ=21) cccGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCA‑CGCCCGCTccc < 2:300409/149‑1 (MQ=25) gTGATGCGGCCATCCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCA‑CGCCCGCTCCcaca > 6:275120/1‑149 (MQ=21) agagCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCCGT‑‑GTTCGCGGGCA‑CGCCCGCGCCCACAGACTGATTGCCaat > 6:183341/1‑146 (MQ=21) ccTCAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCA‑CGCCCGCTccc < 3:350572/127‑1 (MQ=18) aGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCGTTGCCGGT‑‑GTTCGCGGGCA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg > 6:285171/1‑149 (MQ=21) ggCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACa < 5:275120/148‑1 (MQ=21) gCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGt < 1:52367/149‑1 (MQ=21) acaGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCCCGGTGCAAGGCACCGGCATTTCCGGT‑‑GTTATCGGGCACCC‑CCCCTCCCACAGCCAGATTGCCAATAAGCGCGCCAATACCTACAGAGAccg > 6:113594/1‑147 (MQ=16) aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCAGGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgt < 5:113594/149‑1 (MQ=21) gCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGGGCATGGCACCGGCGTTGCCGGT‑‑GTTCGCGGGCACCC‑CCGCTACCCCAGACTGATGTCCATTAAGCGGGCCTAGCGCTACAGTGACGGCCCGTTTGTTGTTGCCt > 1:70685/1‑149 (MQ=18) aaaGGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGc < 7:326145/149‑1 (MQ=21) aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCCGTGCATTGCACCCGCTTTGCCGGTGAGTA‑GCGGGCA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGGTACAGTGACGGCCCGGTTTTTGTTGGCTCg > 1:307937/1‑149 (MQ=21) cgggcgggCCCGCTAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATg < 2:307937/149‑1 (MQ=25) cccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt > 1:147046/1‑148 (MQ=37) ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt > 6:282803/1‑147 (MQ=37) gcgAAACAGGCGCCGCGGTGCATTGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa > 2:345167/1‑147 (MQ=37) gcgAAACAGGCGCCGCGGTGCATTGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa < 1:345167/147‑1 (MQ=37) gcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAg > 7:158436/1‑148 (MQ=37) gAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAAc < 8:158436/149‑1 (MQ=37) gAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAAc < 5:119545/149‑1 (MQ=37) aGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTAt < 5:282803/149‑1 (MQ=37) ggCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc < 2:70685/149‑1 (MQ=37) gcgGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg > 5:234819/1‑146 (MQ=37) gCCGGT‑‑GTTCGCGGGCA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtct > 4:179837/1‑149 (MQ=255) ccGGT‑‑GTTCGCGGGCA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtctt > 4:89159/1‑149 (MQ=255) gtgttcgcGGGCA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCa < 5:285171/142‑1 (MQ=255) gtgttcgcGGGCA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCa < 4:117411/142‑1 (MQ=255) tcgcGGGCA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGACTCGCTAACGAATCATGACTCACTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGAc > 4:91102/4‑149 (MQ=255) cGGGCA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGAcgcg > 6:192033/1‑149 (MQ=255) cGGGCA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGAcgcg > 4:2554/1‑149 (MQ=255) gCA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACgcggcg > 3:94048/1‑149 (MQ=255) | AAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGT‑‑GTT‑TCGGGCA‑CG‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCG > NC_002947/1499365‑1499649 |
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: ‑ |