Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 2,375,739 | A→T | 20.0% | L53L (CTT→CTA) | trpD ← | anthranilate phosphoribosyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 2,375,739 | 0 | A | T | 20.0% | 37.8 / 6.1 | 20 | L53L (CTT→CTA) | trpD | anthranilate phosphoribosyltransferase |
Reads supporting (aligned to +/- strand): ref base A (5/11); new base T (2/2); total (7/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.87e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01 |
GCAACAATTGCCGAGGCCGTTGAGATATTAAAAGTGGAAATACCGTCTCCCCCTGTTCCGCATGTATCAACAATATCAGGAAGTCCATCGACTGTAAGAGCGTGTGCCCGCATTGCCTTCACAAAACCCGCAAGCTCTTCTGGCGTCTCCCCCCGATGAGCAAGAATTGACAATATCCCCCCCATTTCAGAAGGAGTCATTTCCGCTGCCATCATCATATTCATCAGCGTTTCAGCCTCACCGGCAGT > BS168_reference_genome/2375608‑2375855 | gCAACAATTGCTGATGCCGTTGAGATATTAAAAGTAGAAATACCGTCTCCTCCGGTTCCGCATGTATCAACGACATCAGGAAGTCCATCGACTGTAAGAGCGCGTGCCCGCATCGCCTATACAAAGCCCGTTAGCtctt > 1:94939/1‑139 (MQ=255) gCCGAGGCCGTTGAGATATTAAAAGTGGAAATACCGTCTCCCCCTGTGCCGCATGTATCAACAATAACCGGAATTCCCTCGACTGTATGAGCGGGGGCCCGCATTGCCTTCACAAAACCCGCAAGCTCTTCTGGCGtc < 8:157311/138‑1 (MQ=255) ccGAGGCCGTTGAGATATTAAAAGTGGAAATACCGTCTCCCCCTGTTCCGCAGGTATCAACAATATCAGGAAGTCCATCGACTGTAAGAGCGTGTGCCCGCATTGCCTTCACAAAACCCGCAAGCTCTTCTGGCGtctc < 4:297771/139‑1 (MQ=255) gCCGTTGAGATATTAAAAGTGGAAATACCGTCTCCCCCTGTTCCGCATGTATCAACAATATCAGGAAGTCCATCGCCTGTAAGAGCGTGTGCCCGCATTGCGTTCACAAAACCCGCAAGCTCTTCTGGCGtctgggcgg > 8:32214/1‑133 (MQ=255) agaTATTAAAAGTGGAAATACCGTCTCCCCCTGTTCCGCATGTATCAACAATATCAGGACGTCCATCGACTGTAAGAGCGTGTGCCCGCATTGCCTTCACAAAACCCGCAAGCTCTTCTGGCGTCTCCCCCCGATGAGc < 5:35474/139‑1 (MQ=255) ttAAAAGTGGAAATACCGTCTCCCCCTGTTCCGCATGTATCAACAATATCAGGAAGTCCATCGACTGTAAGAGCGTGTGCCCGCATTGCCTTCACAAAACCCGCAAGCTCTTCTGGCGTCTCCCCCCg < 6:288056/128‑1 (MQ=255) ttAAAAGTGGAAATACCGTCTCCCCCTGTTCCGCATGTATCAACAATATCAGGAAGTCCATCGACTGTAAGAGCGTGTGCCCGCATTGCCTTCACAAAACCCGCAAGCTCATCTGGCGTCTCCCCCCg > 5:288056/1‑128 (MQ=255) aaatGTAGAAATCCCGTCTCCTCCGGTTCCGCAGGTAGCAACGACATCAGGAAGGCCATCGACTGTAAGAGCGCGTGCCCGCATCGCCTTTACAAAGCCCGTTAGCTCTTCTGGCGTCTCTCCCCGATGAGCAAGAAt < 2:94939/138‑1 (MQ=255) ccGTCTCCCCCTGTTCCGCATGTATCAACAATATCAGGAAGTCCATCGACTGTAAGAGCGTGTGCCCGCATTGCCTTCACAAAACCCGCAAGCTCTTCTGGCGTCTCCCCCCGATGAGCAAGAATTGACAATATcccc < 7:32214/138‑1 (MQ=255) gTTCCGCATGTATCAACGACATCAGGAAGTCCATCGACTGTAAGAGCGCGTGCCCGCATCGCCTTTACAAAGCCCGTTAGCTCTTCTGGCGTCTCTCCCCGATGAGCAAGAAt > 8:205998/1‑113 (MQ=255) gTTCCGCATGTATCAACGACATCAGGAAGTCCATCGACTGTAAGAGCGCGTGCCCGCATCGCCTTTACAAAGCCCGTTAGCTCTTCTGGCGTCTCTCCCCGATGAGCAAGAAt < 7:205998/113‑1 (MQ=255) aaaTATCAGCAAGTCCATCGACTGTAAGAGCGTGTGCCCGCATTGCCTTCACAAAACCCGCAAGCTCTTCTGGCGTCTCCCCCCGATGAGCAAGAATTGACAATATCCCCCCCATTTCAGAAGGAGTCATTTCCGCTGc < 3:95365/138‑1 (MQ=255) tatCAGGAAGTCCATCGACTGTAAGAGCGTGTGCCCGCATTGCCTTCACAAAACCCGCAAg > 3:135157/1‑61 (MQ=255) tatCAGGAAGTCCATCGACTGTAAGAGCGTGTGCCCGCATTGCCTTCACAAAACCCGCAAg < 4:135157/61‑1 (MQ=255) tCCATCGACTGTAAGAGCGTGTGCCCGCATTGCCTTCACAAAACCCGCAAGCTCTTCTGGCGtc > 4:72191/1‑64 (MQ=255) tCCATCGACTGTAAGAGCGTGTGCCCGCATTGCCTTCACAAAACCCGCAAGCTCTTCTGGCGtc < 3:72191/64‑1 (MQ=255) tCGACTGTAAGAGCGTGTGCCCGCATTGCCTTCACAAAACCCGCAAGCTCTTCTGGCGTCTCCCCCCGATGAGCAAGAATTGACAATATCCCCCCCATTTCAGAAGGAGTCATTTCCGCTGCCATCATCATATtcatca < 6:281663/139‑1 (MQ=255) gTAAGAGCGTGTGCCCGCATTGCCTTCACAAAACCCGCa > 2:120444/1‑39 (MQ=255) gTAAGAGCGTGTGCCCGCATTGCCTTCACAAAACCCGCa < 1:120444/39‑1 (MQ=255) gagCGTGTGCCCGCATTGCCTTCACAAAACCCGCAAGCTCTTCTGGCGTCTCCCCCCGATGAGCAAGAATTGACAATATCCCCCCCATTTCAGAAGGAGTCATTTCCGCTGCCATCATCATATTCATCAGCGTTTCAGc < 5:105430/139‑1 (MQ=255) tgtgCCCGCATTGCCTTCACAAAACCCGCAAGCTCTTCTGGCGTCTCCCGGCGAGGAGGAGGATTGGACATTATCGCCCGCGTTGCAGAAGGGGTCATTTCCGCTGGCATGAGCTGAGTATTCAGCGGTTCAGGCTCAc > 6:239144/1‑139 (MQ=255) gCATTGCCTTCACAAAACCCGCAAGCTCTTCTGGCGTCTCCCCCCGATGAGCAAGAATTGACAATATCCCCCCCATTTCAGAAGGAGTCATTTCCGCTGCCATCATCATATTCATCAGCGTTTCAGCCTCACCGGCAGt < 3:119169/139‑1 (MQ=255) | GCAACAATTGCCGAGGCCGTTGAGATATTAAAAGTGGAAATACCGTCTCCCCCTGTTCCGCATGTATCAACAATATCAGGAAGTCCATCGACTGTAAGAGCGTGTGCCCGCATTGCCTTCACAAAACCCGCAAGCTCTTCTGGCGTCTCCCCCCGATGAGCAAGAATTGACAATATCCCCCCCATTTCAGAAGGAGTCATTTCCGCTGCCATCATCATATTCATCAGCGTTTCAGCCTCACCGGCAGT > BS168_reference_genome/2375608‑2375855 |
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: ‑ |