| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | BS168_reference_genome | 2,375,733 | A→G | 18.2% | G55G (GGT→GGC) | trpD ← | anthranilate phosphoribosyltransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | BS168_reference_genome | 2,375,733 | 0 | A | G | 18.2% | 43.3 / 5.8 | 22 | G55G (GGT→GGC) | trpD | anthranilate phosphoribosyltransferase |
| Reads supporting (aligned to +/- strand): ref base A (6/12); new base G (2/2); total (8/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.02e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-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: ‑ |