| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | AE015451 | 1,499,506:1 | +C | 100% | intergenic (+83/‑81) | trpS → / → zapE | Tryptophan‑‑tRNA ligase/putative enzyme with nucleoside triphosphate hydrolase domain |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | AE015451 | 1,499,506 | 1 | . | C | 94.1% | 49.4 / ‑2.6 | 17 | intergenic (+83/‑81) | trpS/zapE | Tryptophan‑‑tRNA ligase/putative enzyme with nucleoside triphosphate hydrolase domain |
| Reads supporting (aligned to +/- strand): ref base . (0/1); new base C (3/13); total (3/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01 | |||||||||||
| Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
CCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGT‑GTT‑T‑CGGGCAC‑G‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGG > AE015451/1499288‑1499695 | ccGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GTTCG‑CGGGCAc‑g‑ccc < 1:489518/223‑2 (MQ=255) tCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GTTCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc < 1:357684/235‑1 (MQ=39) cGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GTTCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc < 1:194281/234‑1 (MQ=39) aaTGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GTTCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc < 2:224791/221‑1 (MQ=38) aTGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GTTCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAc < 1:421873/235‑1 (MQ=38) ggCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GTT‑CGCGGGCAC‑GCCCGCTCCCACAGAc < 2:207563/145‑1 (MQ=21) tGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GTTCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa < 2:352601/153‑1 (MQ=25) cGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GTTCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTcc > 1:517682/1‑194 (MQ=38) cccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTGGCCGGT‑GTTCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGc < 1:74756/217‑1 (MQ=39) gcgcCGCGGTGCATGGCACCGGCTTTGCCGGT‑GTTCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTCTCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa < 2:168247/138‑1 (MQ=37) gcgcCGCGGTGCATGGCACCGGCTTTGCCGGT‑GTTCG‑CGCGCAC‑GCCCGCTCGCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGt < 1:449548/165‑1 (MQ=38) cgccgcGGTGCATGGCACCGGCTTTGCCGGT‑GTTCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTgg < 1:429131/234‑1 (MQ=255) gccgcGGTGCATGGCACCGGCTTTGCCGGT‑GTTCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt < 1:115073/126‑1 (MQ=25) ggcaccGGCTTTGGCGGT‑GTTCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTTCCTTAGAACGCTATCAAGCAGATCTGATAc > 1:49683/6‑149 (MQ=37) ggcaccGGCTTTGCCGGT‑GTTCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCa > 1:466658/6‑139 (MQ=38) gCTTTGCCGGTGGTC‑G‑CGGGCACTC‑CCGCTCCCACAGTCTGTTTGCCAGTAAGCGGGCCTATCGCTACCGTGTCGGCCCGTTTTTTGTTGCCTCGCTAACGAAt < 2:14853/104‑1 (MQ=255) ttGCCGGT‑GTTCG‑CGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTg < 2:460432/175‑1 (MQ=255) gggCAC‑GCCCGCTCCCATAGACTGATTGCCATTAAGCGGGCCTATCGCTACACTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa < 2:253992/89‑1 (MQ=255) | CCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGT‑GTT‑T‑CGGGCAC‑G‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGG > AE015451/1499288‑1499695 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |