Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | AE015451 | 1,499,480 | 2 bp→AC | intergenic (+57/‑106) | 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,480 | 0 | C | A | 93.3% | 34.2 / NA | 15 | intergenic (+57/‑107) | trpS/zapE | Tryptophan‑‑tRNA ligase/putative enzyme with nucleoside triphosphate hydrolase domain |
Reads supporting (aligned to +/- strand): ref base C (1/0); new base A (4/10); total (5/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.33e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01 | |||||||||||
* | AE015451 | 1,499,481 | 0 | A | C | 100.0% | 41.9 / NA | 15 | intergenic (+58/‑106) | trpS/zapE | Tryptophan‑‑tRNA ligase/putative enzyme with nucleoside triphosphate hydrolase domain |
Reads supporting (aligned to +/- strand): ref base A (0/0); new base C (5/10); total (5/10) |
TTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTT‑TCGGGCAC‑G‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGC > AE015451/1499335‑1499648 || ttGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCtt < 1:111602/153‑1 (MQ=255) tCGACAACGCTGAAGCCCGTGGTGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa < 1:159777/237‑1 (MQ=38) gAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa < 1:253742/237‑1 (MQ=255) gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTa < 2:117045/162‑1 (MQ=25) gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGGGCATGGCCCCgg > 1:97551/1‑112 (MQ=255) gagcGAAGCCTTGGCACCACAGCAGGACTGAGGGTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGCGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCgg < 2:175411/144‑1 (MQ=21) gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCttt > 2:122962/1‑96 (MQ=255) 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) gACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTg < 1:348069/72‑1 (MQ=255) gACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTg > 2:348067/1‑72 (MQ=255) ggCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGc < 2:97550/156‑1 (MQ=25) gCCAAACGGTGGGAGCGGGCGTGCCCCCGAAACAGGCGCCGGGGTGCATGGCACCGGCCTTGCCGGTGTTCGCGGGCACGG‑CCGCTccc > 1:98301/1‑89 (MQ=17) || TTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTT‑TCGGGCAC‑G‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGC > AE015451/1499335‑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: ‑ |