Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 4,155,390:1 | +A | 100% | intergenic (‑181/+43) | trnY‑Lys ← / ← purA | tRNA‑Lys/adenylosuccinate synthetase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 4,155,390 | 1 | . | A | 94.9% | 127.4 / ‑0.7 | 39 | intergenic (‑181/+43) | trnY‑Lys/purA | tRNA‑Lys/adenylosuccinate synthetase |
Reads supporting (aligned to +/- strand): ref base . (1/1); new base A (15/22); total (16/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch. |
ATCTTGTACTCTATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCAAAC‑AAAA‑ACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAATGCCT > BS168_reference_genome/4155259‑4155514 | atattGTACTCTTTTGAAGCAGTTTCAGTAACAACGAGATTTATAATATCACGAACATATTGAAGTGGCAATAACTTTTTTTGTTGTTTTTCATATTTTTATAGCAGACGAACTTAT‑ATGAGGCTTTCatcc‑aaaa‑aa > 6:540/4‑137 (MQ=255) aTCTTGTACTCTATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCaaac‑aaaa‑aa > 7:194899/1‑137 (MQ=255) ttGTACTCTATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCaaac‑aaaa‑aacaag < 6:166190/139‑6 (MQ=255) gTACTCTATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCaaac‑aaaa‑aacaag < 7:37/137‑6 (MQ=255) gTACTCTATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCaaac‑aaaa‑aacaa < 8:69611/136‑5 (MQ=255) tATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCCTTa > 5:136151/1‑138 (MQ=255) ttGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCCTTaaa > 3:181061/1‑138 (MQ=255) gAAGCAGTTTCCGTAACAACGAGATTTATAATATAACTAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCCTTAAATa < 7:92641/138‑1 (MQ=255) gTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAATCCACTTAAATAGGGGCtt > 4:143981/1‑139 (MQ=255) gTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCCTTAAATAGGGGCtt > 1:205978/1‑139 (MQ=255) gTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCCTTAAATAGGGGCtt > 8:67716/1‑139 (MQ=255) gTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCCTTAAATAGGGGCtt > 6:203844/1‑139 (MQ=255) gTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCCTTAAATAGGGGCtt > 4:109135/1‑139 (MQ=255) gTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCCTTAAATAGGGGCtt > 3:84213/1‑139 (MQ=255) tttCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGGGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCCTTAAATAGGGGCTTg < 4:197202/139‑1 (MQ=255) gTAACAACGAGATTTATAATATCACGAACAACTAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCCTTAAATAGGGGCTTGCAGAc < 6:87253/139‑1 (MQ=255) aacGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCGTTAAATAGGGGCTTGCAGACATAtt > 4:117258/1‑139 (MQ=255) aacGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCCTTAAATAGGGGCTTGCAGAc > 1:216540/1‑134 (MQ=255) aacGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACAAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCCTTAAATAGGGGCTTGCAGAc < 2:216540/134‑1 (MQ=255) ataataTCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTGCAAACAAAAA‑ACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGt < 7:170378/138‑1 (MQ=255) taataTCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTc < 8:194899/139‑1 (MQ=255) atatCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGc < 4:181061/139‑1 (MQ=255) atatCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGc < 1:7101/139‑1 (MQ=255) tatCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCa < 6:43672/139‑1 (MQ=255) tCACGAACAAATAGAAGTGGCAATTACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGGGGCTTTCAAAC‑AAAATACAATCCCCTTAAATAGGGGCTTGCAGACATATTCTTTTTAGTTCGCAc < 4:48695/138‑1 (MQ=255) tCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCCTT‑AATAGGGGCTTGCAGCCATATTCTATTTAGTTCGCACgg > 6:219632/1‑139 (MQ=255) cGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAATCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACTGTa < 2:139681/139‑1 (MQ=255) agcAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTaca < 1:57583/137‑1 (MQ=255) cAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTacaca < 7:11381/139‑1 (MQ=255) gAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCg < 7:67716/138‑1 (MQ=255) gtcAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACAAAATAATAATAAAATTTAAAAAA‑AAAA‑ACAAACACCTTAAATAGGGTATTGCAGACATATTCTATTTATTTCTCACGGTACACACTGCGAAggct > 2:57583/3‑137 (MQ=255) ggCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGc < 4:80740/118‑1 (MQ=255) ggCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGc > 3:80740/1‑118 (MQ=255) aaCTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCCTTAAATAGGGGCTTGCAGAc < 1:186127/94‑1 (MQ=255) aaCTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCCTTAAATAGGGGCTTGCAGAc > 2:186127/1‑94 (MQ=255) ttGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCCTT‑AATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAg < 5:219632/138‑1 (MQ=255) tGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAg < 3:109135/138‑1 (MQ=255) attCTCAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAgaga < 5:33/137‑1 (MQ=255) tcAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGc < 4:193389/81‑1 (MQ=255) tcAGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGc > 3:193389/1‑81 (MQ=255) aGACCAACTTAT‑ATGCGGCTTTCAAACAAAAA‑ACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGa < 4:84213/139‑1 (MQ=255) tGCGGCTTTCAAACAAAAA‑ACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAATGCCt < 6:136151/139‑1 (MQ=255) cTTTCAAACAAAAA‑ACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGc > 4:166061/1‑61 (MQ=255) cTTTCAAACAAAAA‑ACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGc < 3:166061/61‑1 (MQ=255) | ATCTTGTACTCTATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTAT‑ATGCGGCTTTCAAAC‑AAAA‑ACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAATGCCT > BS168_reference_genome/4155259‑4155514 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |