| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | BS168_reference_genome | 2,375,710 | T→C | 15.9% | H63R (CAC→CGC) | trpD ← | anthranilate phosphoribosyltransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | BS168_reference_genome | 2,375,710 | 0 | T | C | 15.9% | 54.0 / 5.4 | 26 | H63R (CAC→CGC) | trpD | anthranilate phosphoribosyltransferase |
| Reads supporting (aligned to +/- strand): ref base T (7/14); new base C (2/2); total (9/17) | |||||||||||
| 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.07e-01 | |||||||||||
AGCGATTTTCGCACCAGCTGCCGAGGCAACAATTGCCGAGGCCGTTGAGATATTAAAAGTGGAAATACCGTCTCCCCCTGTTCCGCATGTATCAACAATATCAGGAAGTCCATCGACTGTAAGAGCGTGTGCCCGCATTGCCTTCACAAAACCCGCAAGCTCTTCTGGCGTCTCCCCCCGATGAGCAAGAATTGACAATATCCCCCCCATTTCAGAAGGAGTCATTTCCGCTGCCATCATCATATTCATCAGCGTTTCAGCCTCAC > BS168_reference_genome/2375583‑2375848 | aGCGATTTTCGCACCAGCTGCCGAGGCAACAATTGCCGAGGCCGTTGAGATATTAAAAGTGGAAATACCGTCTCCCCCTGTTCCGCATGTATCAACAATATCAGGAAGTCCATCGACTGTAAGAGCGTGTGCCCGCa < 3:93599/137‑1 (MQ=255) gATTTTCGCACCAGCTGCCGAGGCAACAATTGCCGAGGCCGTTGAGATATTAAAAGTGGAAATACCGTCTCCCCCTGTTCCGCATGTATCAACAATATCAGGAAGTCCATCGACTGTAAGAGCGTGTGCCCGCATTGcc < 1:133969/139‑1 (MQ=255) ttCGCACCAGCTGCCGAGGCAACAATTGCCGAGGCCGTTGAGATATTAAAAGTGGAAATACCGTCTCCCCCTGTTCCGCATGTATCAACAATATCAGGAAGTCCATCGACTGTAAGAGCGTGTGCCCGCATTGCCTTca > 4:119169/1‑139 (MQ=255) gCACCAGCTGCCGAGTCATCAATTGCCGAGGCCGTTGAGATATTAAAAGTGGATATACCGTCTCCCCATGTTCCGCATGTATCAACAATATTAGGAAGTCCATCGACTGTTAGAGCGTGTGCCCGCATTGCCTTCACaa < 8:83441/139‑1 (MQ=255) 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) gTGGAAATACCGTCTCCCCCTGTTCCGCATGTATCAACAATATCAGGAAGTCCATCGACTGTAAGAGCGTGTGCCCGCATTGCCTTCACaa > 8:255786/1‑91 (MQ=255) gTGGAAATACCGTCTCCCCCTGTTCCGCATGTATCAACAATATCAGGAAGTCCATCGACTGTAAGAGCGTGTGCCCGCATTGCCTTCACaa < 7:255786/91‑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) | AGCGATTTTCGCACCAGCTGCCGAGGCAACAATTGCCGAGGCCGTTGAGATATTAAAAGTGGAAATACCGTCTCCCCCTGTTCCGCATGTATCAACAATATCAGGAAGTCCATCGACTGTAAGAGCGTGTGCCCGCATTGCCTTCACAAAACCCGCAAGCTCTTCTGGCGTCTCCCCCCGATGAGCAAGAATTGACAATATCCCCCCCATTTCAGAAGGAGTCATTTCCGCTGCCATCATCATATTCATCAGCGTTTCAGCCTCAC > BS168_reference_genome/2375583‑2375848 |
| 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: ‑ |