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 | 94.0% | 61.0 / ‑4.7 | 29 | intergenic (+74/‑90) | trpS/zapE | tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base . (0/2); major base C (13/12); minor base A (2/0); total (15/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.87e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
* | NC_002947 | 1,499,498 | 0 | T | G | 100.0% | 81.4 / NA | 29 | 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/14); total (15/14) | |||||||||||
Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
GCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCA‑CG‑GCTTTGCCGGTGT‑T‑TCGGG‑CA‑C‑G‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCC > NC_002947/1499361‑1499628 || gCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCATGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACCGGCGCCGCGGAGCATGGCAC‑CG‑CCTTTGCCCGAGT‑TCGCGGG‑CC‑c‑g‑cc > 7:169549/1‑149 (MQ=25) gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACACGCGCCGCGGTGCATGGCAC‑CC‑GCATTGCCGGTGT‑TCGCGGG‑CACC‑C‑CCGCACCcacaaa > 1:25577/1‑147 (MQ=21) aTGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTACCCGCGAAACAGGCGCCGCGGTGCATGGCAC‑CG‑GCTTTGCCGGTGT‑TCGCGGG‑CA‑C‑GCCCGCTCCCACAGAc < 1:300308/149‑1 (MQ=21) cATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCAC‑CG‑GCTTTGCCGGTGT‑TCGCGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGc > 2:57712/1‑149 (MQ=21) agagCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCAC‑CG‑GCTTTGCCGGTGT‑TCGCGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCAtt > 2:44506/1‑148 (MQ=21) gagCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGGAACAGGCGCCGCGGGGCATGGCAC‑CG‑GCGTTGCCGGTGT‑TAGCGGG‑CC‑CGG‑CCGCGCCCACAGACTGATTGCCAtt > 6:263940/1‑147 (MQ=255) tGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGGCCGCGAAACACGCGCCCCGGTGGCTGGCAC‑CC‑GCATTGCCGGTGT‑TAGCGGGCCC‑C‑C‑CCGCGCCCACCGACTGAAGGCCAATAAGCcgg > 6:69748/1‑149 (MQ=255) aactatgCCTGGGCACCACTGCTGGCGTGAGGCTGCCAATGGGTGGGAGGGGGCGGGCCCGCTAAACAGGCGCCGGGGTGCATGGCAC‑GG‑GCTTTGCGGGTGT‑TCGGGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGc < 2:25577/142‑1 (MQ=255) aaCGAAGCCTTGGCACCACAACAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCAC‑CG‑GCTTTGCCGGTGT‑TCGCGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGc < 7:112072/148‑1 (MQ=21) ttGGCACCACCGCAGGACTGAGGCTGCCAAACGGTGGGCGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCAC‑CG‑GCTTTGCCGGTGT‑TCGCGGG‑CA‑C‑GCCCGCTCCCACAGa < 1:320691/119‑1 (MQ=17) cTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCAC‑CG‑GCTTTGCCGGTGT‑TCGCGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGttttt < 3:225773/149‑1 (MQ=21) cggcggcCATCCGGGGGGCGGGGGCGGGCCCGCGAAAGAGGCGCCGCGGTGC‑TGGCCACCG‑GCTTGGCCGGTGTGC‑GCGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgt < 5:210750/148‑1 (MQ=21) cggCTGCCAACCGGTGGGAGCGGGCGGCCCCGCGAAACCGGCGCCGCGGTCCAT‑GCCA‑CGTGCTTTGCGGGTGTGC‑GCGGG‑AA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgt < 5:205652/148‑1 (MQ=255) gTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCAC‑CG‑GCTTTGCCGGTGT‑TCGCGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgt < 5:98567/136‑1 (MQ=21) ggAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCAC‑CG‑GCTTTGCCGGTGT‑TCGCGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAAt < 1:57712/149‑1 (MQ=25) gggCGTGCCCGCGAAACAGGCGCGGCTGTGCATGGGAC‑GG‑GCTTTGCCGGTGT‑TCGCGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTGGCTAACGAATCATGa < 8:169549/149‑1 (MQ=25) gggCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCAC‑CG‑GCTTTGCCGGTGT‑TCGCGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa < 1:44506/149‑1 (MQ=37) cccGCGAAACAGGCGCCGCGGTGCATGGCAC‑CG‑GCTTTGCCGGTGT‑TCGCGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt > 6:343045/1‑148 (MQ=37) ccGCGAAACAGGCGCCGCGGTGCATGGCAC‑CG‑GCTTTGCCGGTGT‑TCGCGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa < 7:190422/141‑1 (MQ=25) ccGCGAAACAGGCGCCGCGGTGCATGGCAC‑CG‑GCTTTGCCGGTGT‑TCGCGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa > 8:190422/1‑141 (MQ=25) cgcgAAACAGGCGCCGCGGTGCATGGCAC‑CG‑GCTTTGCCGGTGT‑TCGCGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGaa > 7:118803/1‑134 (MQ=25) cgcgAAACAGGCGCCGCGGTGCATGGCAC‑CG‑GCTTTGCCGGTGT‑TCGCGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGaa < 8:118803/134‑1 (MQ=25) cgcgAAACAGGCGCCGCGGTGCATGGCAC‑CG‑GCTTTGCCGGTGT‑TCGCGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt > 7:98768/1‑146 (MQ=37) gcgAAACAGGCGCCGCGGTGCATGGCAC‑CG‑GCTTTGCCGGTGT‑TCGCGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAg > 7:91456/1‑148 (MQ=37) aGGCGCCGCGGTGCATGGCAC‑CG‑GCTTTGCCGGTGT‑TCGCGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTa > 8:186012/1‑148 (MQ=37) gcgcCGCGGTGCATGGCAC‑CG‑GCTTTGCCGGTGT‑TCGCGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTATAAAGCTAta > 8:258810/1‑147 (MQ=37) cgccgcGGTGCATGGCAC‑CG‑GCTTTGCCGGTGT‑TCGCGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCaa < 8:98768/149‑1 (MQ=37) cgcgGTGCATGGCAC‑CG‑GCTTTGCCGGTGT‑TCGCGGG‑CA‑C‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg > 6:91514/1‑147 (MQ=37) cggcaaagCCGGTGT‑TCGCGGG‑CA‑C‑GCCCGCTCCCACAGACAGCTTGCCATTAAGCGGGCCTCTCGCTACAGTGACGGCCCGGTGTTTGTTGCCTCGCTAACGAATAATGCCTCCCGTAGAACGCTATCAAGCAGATCTGAAACCTccc > 1:307600/8‑149 (MQ=255) || GCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCA‑CG‑GCTTTGCCGGTGT‑T‑TCGGG‑CA‑C‑G‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCC > NC_002947/1499361‑1499628 |
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: ‑ |