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 | 76.5% | 27.9 / 9.8 | 17 | intergenic (+83/‑81) | trpS/zapE | tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base . (3/1); new base C (4/9); total (7/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.50e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-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: ‑ |