| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | BS168_reference_genome | 2,375,679 | A→G | 15.9% | I73I (ATT→ATC) ‡ | trpD ← | anthranilate phosphoribosyltransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | BS168_reference_genome | 2,375,679 | 0 | A | G | 15.9% | 57.9 / 2.1 | 25 | I73I (ATT→ATC) ‡ | trpD | anthranilate phosphoribosyltransferase |
| Reads supporting (aligned to +/- strand): ref base A (8/13); new base G (2/2); total (10/15) | |||||||||||
| 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 | |||||||||||
AGCGCTTCCGCTTTTAGAAGAGACAGAGCGATTGCCATGCTTAGCGATTTTCGCACCAGCTGCCGAGGCAACAATTGCCGAGGCCGTTGAGATATTAAAAGTGGAAATACCGTCTCCCCCTGTTCCGCATGTATCAACAATATCAGGAAGTCCATCGACTGTAAGAGCGTGTGCCCGCATTGCCTTCACAAAACCCGCAAGCTCTTCTGGCGTCTCCCCCCGATGAGCAAGAATTGACAATATCCCCCCCATTTCAGAAGGAGTCATTTCCGCTGC > BS168_reference_genome/2375541‑2375816 | aGCGCTTCCGCTTTTAGAAGAGACAGAGCGATTGCCATGCTTAGCGATTTTCGCACCAGCTGCCGAGGCAACAATTGCCGAGGCCGTTGAGATATTAAAAGTGGAAATACCGTCTCCCCCTGTTCCGCATGTATcaaca < 1:29528/139‑1 (MQ=255) ccGCTTTTAGAAGAGACAGAGCGATTGCCATGCTTAGCGATTTTCGCACCAGCTGCCGAGGCAACAATTGCCGAGGCCGTTGAGATATTAAAAGTGGAAATACCGTCTCCCCCTGTTCCGCATGTATCAACAATATCAg > 7:103205/1‑139 (MQ=255) ccGCTTTTAGAAGAGACAGAGCGATTGCCATGCTTAGCGATTTTCGCACCAGCTGCCGAGGCAACAATTGCCGAGGCCGTTGAGATATTAAAAGTGGAAATACCGTCTCCCCCTGTTCCGCATGTATCAACAATATCAg > 5:281663/1‑139 (MQ=255) aagaGACAGAGCGATTGCCATGCTTAGCGATTTTCGCACCAGCTGCCGAGGCAACAATTGCCGAGGCCGTTGAGATATTAAAAGTGGAAATACCGTCTCCCCCTGTTCCGCATGTATCAACAATATCAGGAAGTCCATc < 8:103205/139‑1 (MQ=255) gagaCAGAGCGATTGCCATGCTTAGCGATTTTCGCACCAGCTGCCGAGGCAACAATTGCCGAGGCCGTTGAGATATTAAAAGTGGAAATACCGTCTCCCCCTGTTCCGCATGTATCAACAATATCAGGAAGTCCATCGa < 2:207257/139‑1 (MQ=255) gagaCAGAGCGATTGCCATGATTAGCGATTTTCGCACCAGCTGCCGAGGCAACAATTGCCGAGGCCGTGGAGATATTAAATGTGGAAATACCGTCTCCCCCTGGTCGGCATGTATCAACAATATCAGGAAGTCCATCGa < 5:203205/139‑1 (MQ=255) agCGATTGCCATGCTTAGCGATTTTCGCACCAGCTGCCGAGGCAACAATTGCCGAGGCCGTTGAGATATTAAAAGTGGAAATACCGTCTCCCCCTGTTCCGCATGTATCAACAATATCAGGAAGTCCATCGACTGTAag > 1:287268/1‑139 (MQ=255) aGCGATTTTCGCACCAGCTGCCGAGGCAACAATTGCCGAGGCCGTTGAGATATTAAAAGTGGAAATACCGTCTCCCCCTGTTCCGCATGTATCAACAATATCAGTAAGTCCATc < 2:34318/114‑1 (MQ=255) aGCGATTTTCGCACCAGCTGCCGAGGCAACAATTGCCGAGGCCGTTGAGATATTAAAAGTGGAAATACCGTCTCCCCCTGTTCCGCATGTATCAACAATATCAGGAAGTCCATc > 1:34318/1‑114 (MQ=255) 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 < 7:255786/91‑1 (MQ=255) gTGGAAATACCGTCTCCCCCTGTTCCGCATGTATCAACAATATCAGGAAGTCCATCGACTGTAAGAGCGTGTGCCCGCATTGCCTTCACaa > 8:255786/1‑91 (MQ=255) ccGTCTCCCCCTGTTCCGCATGTATCAACAATATCAGGAAGTCCATCGACTGTAAGAGCGTGTGCCCGCATTGCCTTCACAAAACCCGCAAGCTCTTCTGGCGTCTCCCCCCGATGAGCAAGAATTGACAATATcccc < 7:32214/138‑1 (MQ=255) gTTCCGCATGTATCAACGACATCAGGAAGTCCATCGACTGTAAGAGCGCGTGCCCGCATCGCCTTTACAAAGCCCGTTAGCTCTTCTGGCGTCTCTCCCCGATGAGCAAGAAt < 7:205998/113‑1 (MQ=255) gTTCCGCATGTATCAACGACATCAGGAAGTCCATCGACTGTAAGAGCGCGTGCCCGCATCGCCTTTACAAAGCCCGTTAGCTCTTCTGGCGTCTCTCCCCGATGAGCAAGAAt > 8:205998/1‑113 (MQ=255) aaaTATCAGCAAGTCCATCGACTGTAAGAGCGTGTGCCCGCATTGCCTTCACAAAACCCGCAAGCTCTTCTGGCGTCTCCCCCCGATGAGCAAGAATTGACAATATCCCCCCCATTTCAGAAGGAGTCATTTCCGCTGc < 3:95365/138‑1 (MQ=255) | AGCGCTTCCGCTTTTAGAAGAGACAGAGCGATTGCCATGCTTAGCGATTTTCGCACCAGCTGCCGAGGCAACAATTGCCGAGGCCGTTGAGATATTAAAAGTGGAAATACCGTCTCCCCCTGTTCCGCATGTATCAACAATATCAGGAAGTCCATCGACTGTAAGAGCGTGTGCCCGCATTGCCTTCACAAAACCCGCAAGCTCTTCTGGCGTCTCCCCCCGATGAGCAAGAATTGACAATATCCCCCCCATTTCAGAAGGAGTCATTTCCGCTGC > BS168_reference_genome/2375541‑2375816 |
| 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: ‑ |