| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | AE015451 | 1,499,508 | +C | intergenic (+85/‑79) | 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 | 91.7% | 22.3 / ‑2.6 | 12 | intergenic (+83/‑81) | trpS/zapE | Tryptophan‑‑tRNA ligase/putative enzyme with nucleoside triphosphate hydrolase domain |
| Reads supporting (aligned to +/- strand): ref base . (1/0); new base C (1/10); total (2/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.67e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01 | |||||||||||
TCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTT‑TCGGGCAC‑G‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGC > AE015451/1499353‑1499648 | tCGACAACGCTGAAGCCCGTGGTGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa < 1:159777/237‑1 (MQ=38) gAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa < 1:253742/237‑1 (MQ=255) gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTa < 2:117045/162‑1 (MQ=25) gagcGAAGCCTTGGCACCACAGCAGGACTGAGGGTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGCGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCgg < 2:175411/144‑1 (MQ=21) accgcAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATg < 1:297984/236‑1 (MQ=38) agcagGACTGAGGATGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAGCAGGCGCCGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAc‑g‑c < 2:246984/96‑1 (MQ=17) gcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACgcggc < 1:122960/237‑1 (MQ=255) agGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgt > 2:228610/1‑156 (MQ=25) ggCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGc < 2:97550/156‑1 (MQ=25) gCCAAACGGTGGGAGCGGGCGTGCCCCCGAAACAGGCGCCGGGGTGCATGGCACCGGCCTTGCCGGTGTTCGCGGGCACGG‑CCGCTccc > 1:98301/1‑89 (MQ=17) cgccGGTTTTGCCGGTGTTCGCGGGCGC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGcc < 2:98300/94‑1 (MQ=25) tgttcgcGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACgcggc < 1:122963/151‑1 (MQ=255) gttcgcGGGCAC‑GCCCGCGCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc < 2:144693/111‑1 (MQ=255) | TCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTT‑TCGGGCAC‑G‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGC > AE015451/1499353‑1499648 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |