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 | 83.2% | 20.2 / ‑3.4 | 19 | intergenic (+74/‑90) | trpS/zapE | tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base . (3/0); major base C (3/10); minor base A (3/0); total (9/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.57e-02 | |||||||||||
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 | 84.2% | 24.1 / 2.7 | 19 | intergenic (+75/‑89) | trpS/zapE | tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base T (0/0); major base G (6/10); minor base C (3/0); total (9/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.67e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01 | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
CCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTT‑T‑CGGGCA‑CG‑CC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCT > NC_002947/1499368‑1499596 ||| cccGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGCAACAGGCGCCGCGGTGCATGGCACCGGCGTTGCCGGTGTTAG‑CGGGCC‑cg‑cccgcgccc > 8:77785/1‑142 (MQ=25) cGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGCGGCTGCCAAACGGTGGGAGCGGGCGTGCCCCCGAAACACGCGCCCCGGTGCAAGGCCCCCGCGTTGCCCGTGTTAG‑CGGGCA‑CG‑cc‑cccccccc > 3:132613/1‑146 (MQ=255) gCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCA‑CGCCC‑GCTCCCACAGACTGATTGCCATTAAGc < 3:49821/149‑1 (MQ=21) aaCGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCA‑CGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGc < 7:77785/148‑1 (MQ=21) aaCGAAGCCTTGCCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCA‑CGCCC‑GCTCCCACAGACTGATTGCCATTAAGCggg < 3:216054/147‑1 (MQ=21) gCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGGGAAACAGGCGCCCCGGGGGATGGCACCGGCTTTGCCGGTGTTAG‑CGGGCC‑CGCCC‑GCGCCCACAGAATGCATGCCATTAAACGGGCCAAtagg > 6:262141/1‑146 (MQ=18) ggCACCACAGCAGGAGTGGGGCTGCCAAAGGGTGGGAGGGGGGGTGCCCGCTAAACTGGCGCGGGGGTGCATGGCACGGGCTTTGCGGGTGTTCG‑CGGGCA‑CGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACa < 4:132613/148‑1 (MQ=21) ggCACCACAGCAGGACTGTGGCGGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCAGGGCACCGGCTTTGCCGGTGTTCG‑CGGGCA‑CGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACa < 3:63644/148‑1 (MQ=21) ggCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCACCC‑CC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACa < 1:242070/148‑1 (MQ=21) gCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCA‑CGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGt < 3:112071/149‑1 (MQ=21) agcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCC‑CGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGcc > 7:235036/1‑149 (MQ=21) aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCA‑CGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt < 1:228925/147‑1 (MQ=21) ccAAACGGTGGGAGCGGGCGTGCACGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCA‑CGCCC‑GCTCCCACAGACTGATTGCCa > 7:310694/1‑103 (MQ=14) ccAAACGGTGGGAGCGGGCGTGCACGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCA‑CGCCC‑GCTCCCACAGACTGATTGCCa < 8:310694/103‑1 (MQ=14) tagggggCGTTCCCGCGAAACAGGCGCCGCGGTGCAGGGCACGGGCTTTGCCGGTGTTCG‑GGGGAA‑CGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATc < 4:147555/145‑1 (MQ=255) GAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCA‑CG‑CCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATc > 8:255317‑M1/1‑149 (MQ=255) GCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCA‑CG‑CCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt > 7:151225‑M1/1‑149 (MQ=255) CGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCA‑CG‑CCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATg > 3:220638‑M1/1‑149 (MQ=255) ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCA‑CGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt > 1:253442/1‑147 (MQ=37) ||| CCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTT‑T‑CGGGCA‑CG‑CC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCT > NC_002947/1499368‑1499596 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |