| 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 | 89.3% | 71.3 / 4.7 | 28 | intergenic (+83/‑81) | trpS/zapE | tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein |
| Reads supporting (aligned to +/- strand): ref base . (3/0); new base C (11/14); total (14/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.22e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01 | |||||||||||
| 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: ‑ |