Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | AE015451 | 1,499,498 | 1 bp→CG | 100% | intergenic (+75/‑89) | 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,497 | 1 | . | C | 88.2% | 45.6 / 1.6 | 17 | intergenic (+74/‑90) | trpS/zapE | Tryptophan‑‑tRNA ligase/putative enzyme with nucleoside triphosphate hydrolase domain |
Reads supporting (aligned to +/- strand): ref base . (0/2); new base C (3/12); 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 = 8.68e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
* | AE015451 | 1,499,498 | 0 | T | G | 94.1% | 53.4 / ‑2.9 | 17 | intergenic (+75/‑89) | trpS/zapE | Tryptophan‑‑tRNA ligase/putative enzyme with nucleoside triphosphate hydrolase domain |
Reads supporting (aligned to +/- strand): ref base T (0/0); major base G (3/13); minor base C (0/1); 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 = 6.97e-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) ||| 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 < 26 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |