| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | BS168_reference_genome | 1,747,662 | A→G | 6.8% | V180V (GTA→GTG) | dapG → | aspartokinase I (alpha and beta subunits) |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | BS168_reference_genome | 1,747,662 | 0 | A | G | 6.8% | 182.1 / 3.8 | 59 | V180V (GTA→GTG) | dapG | aspartokinase I (alpha and beta subunits) |
| Reads supporting (aligned to +/- strand): ref base A (28/27); new base G (2/2); total (30/29) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01 | |||||||||||
GGTAGTCGTCGCAGGCTTCCAGGGCGCAACTGAAAAGGGTGATACAACAACAATCGGCAGAGGGGGCAGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTCCTGTTGTGACGTATACGGAAATCTGCAACCTGGCGTATCAGGGGGCCAAGGTCATTTCTCCGCGTGCTGTGGAAATCGCAATGCA > BS168_reference_genome/1747527‑1747799 | ggTAGTCGTCGCAGGCTTCCAGGGCGCAACTGAAAAGGGTGATACAACAACAATCGGCAGAGGGGGCAGCGATACGTCTGCAGCCGCACTCGGCGCCGCTGTTGATGCGGAATACCTCGATATTTTCACCGATGTAg > 8:290306/1‑137 (MQ=255) tAGTCGTCGCAGGCTTCCAGGGCGCAACTGAAAAGGGTGATACAACAACAATCGGCAGAGGGGGCAGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAgg < 7:21752/139‑1 (MQ=255) tcggcgCAGGCTTCCAGGGCGCAACTGAAAAGGGTGATACAACAACAATCGGCAGAGGGGGCAGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGt < 7:290306/139‑1 (MQ=255) gtcgCAGGCTTCCAGGGCGCAACTGAAAAGGGTGATACAACAACAATCGGCAGAGGGGGCAGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGtga > 3:199191/1‑139 (MQ=255) gcgtCCAGGGCGCAACTGAAAAGGGTGATACAACAACAATCGGCAGAGGGGGCAGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGActgc < 2:66208/136‑1 (MQ=255) gCTTCCAGGGCGCAACTGAAAAGGGTGATACAACAACAATCGGCAGAGGGGGCAGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCTCCGATGTAGAAGGCGTGATGActgc < 8:78070/139‑1 (MQ=255) gCTTCCAGGGCGCAACTGAAAAGGGTGATACAACAACAATCGGCAGAGGGGGCAGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGActgc < 4:199191/139‑1 (MQ=255) gCTTCCAGGGCGCAACTGAAAAGGGTGATACAACAACAATCGGCAGAGGGGGCAGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGActgc < 4:284365/139‑1 (MQ=255) cTTCCAGGGCGCAACTTAAAAGGGTGATACAACAACAATCGGCAGAGGGGGCAGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGActgct < 2:36690/139‑1 (MQ=255) cgcAACTGAAAAGGGTGATACAACAACAATCGGCAGAGGGGGCAGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGACTGCTGACCCaa > 1:36920/1‑137 (MQ=255) cgcAAATGAAAAGGGTGATACAAGAACAATCGGGAGATGGGGCAGCGATCCGTCAGCAGCCGCACTCGGAGCAGCTGTTGATGCGGATAACATCGTTATTTTCACCGATGAAGAAGGCGTGATGCCTGCTGACCCaa > 4:9371/1‑137 (MQ=255) cTGAAAAGGGTGATACAACAACAATCGGCAGAGGGGGCAGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTgg > 3:21758/1‑138 (MQ=255) aaGGGTGATACAACAACAATCGGCAGAGGGGGCAGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGaa > 1:181513/1‑138 (MQ=255) aaGCGTGATACAACAACAATCGGCAGAGGGGGCAGCGATACGTCTGCAGACGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTAACCGATGTAGAAGGAGTGATGACAGCTGACCCAAGAGTGGTACaa > 2:964/1‑138 (MQ=255) ggTGATACAACAACAATCGGCAGAGGGGGCAGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGActgc > 8:101293/1‑116 (MQ=255) ggTGATACAACAACAATCGGCAGAGGGGGCAGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGActgc < 7:101293/116‑1 (MQ=255) gATACAACAACAATCGGCAGAGGGGGCAGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCa > 6:256510/1‑138 (MQ=255) tACAACAACAATCGGCAGAGGGGGCAGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTGGAAGGCGTGATGActgct > 5:236798/1‑112 (MQ=255) tACAACAACAATCGGCAGAGGGGGCAGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTGGAAGGCGTGATGActgct > 5:291027/1‑112 (MQ=255) tACAACAACAATCGGCAGAGGGGGCAGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTGGAAGGCGTGATGActgct < 6:236798/112‑1 (MQ=255) tACAACAACAATCGGCAGAGGGGGCAGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTGGAAGGCGTGATGActgct < 6:291027/112‑1 (MQ=255) caacaaTCGGCAGAGGGGGCAGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCt > 6:172668/1‑139 (MQ=255) caacaaTCGGCAGAGGGGGCAGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCt > 3:214163/1‑139 (MQ=255) caacaaTCGGCAGAGGGGGCAGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGAGGTAGAAAATGCAAAACCGCt > 7:25510/1‑139 (MQ=255) acaaTCGGCAGAGGGGGCAGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTc < 4:119710/139‑1 (MQ=255) acaaTCGGCAGAGGGGGCAGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTc < 8:41819/139‑1 (MQ=255) acaaTCGGCAGAGGGGGCAGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTc > 3:265738/1‑139 (MQ=255) cAGAGGGGGCAGCGATTCGGCTGCAGCCGCACTCGGATCAGCTGTTGATGCGGATTACATCGATATTTTCACCGATGTAGAAGGCGTGATGACTGCTGACCCAATAGTGGTAGAAAATGCAAAACCGCTTCCTGTtgtg < 8:25510/139‑1 (MQ=255) cAGAGGGGGCAGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTCCTGTtgtg < 6:105468/139‑1 (MQ=255) gCAGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGACTGCTGACCAAAGAGTTGTAGATAATGCAAACCCGCTTCCTGTTGTGACGTATAc > 6:27303/1‑139 (MQ=255) aGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTCCTGTTGTGACGTATACgg < 4:214163/139‑1 (MQ=255) aGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTCCTGTTGTGACGTATACg > 1:145732/1‑138 (MQ=255) aGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTCCTGTTGTGACGTATACg > 3:256889/1‑138 (MQ=255) aGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTCCTGTTGTGACGTATACg > 2:244599/1‑138 (MQ=255) tACGTCTGCAGCCGCACTAGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGTCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTCCTGTTGTGACGTATACGGAAATc > 6:154166/1‑139 (MQ=255) gCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAAcc < 6:88021/104‑1 (MQ=255) gCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAAcc > 5:88021/1‑104 (MQ=255) gCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTCCTGTTGTGACGTATACGGAAATCTGCAACCTgg < 4:21758/139‑1 (MQ=255) cGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGGGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTCCTGTTGTGACGTATACGGAAATCTGCAACCTGGCg < 8:32336/139‑1 (MQ=255) cACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTCCTGTTGTGACGTATACGGAAATCTGCAACCTGGCGTa > 4:266285/1‑139 (MQ=255) cGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGACTGCTGACCCTAGAGTGGTAGAAAATGCAAAACCGCTTCCTGTTGTGACGTATACGGAAATCTGCAACCTGGCGTATCa < 1:964/138‑1 (MQ=255) ggAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTCCTGTTGTGACGTATACGGAAATCTGCAACCTGGCGTATCAg < 3:101012/138‑1 (MQ=255) gATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTCCTGTTGTGACGTATACGGAAATCTGCAACCTGGCGTATCAGGGGGCCAAGGTCa > 5:81143/1‑139 (MQ=255) gATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTCCTGTTGTGACGTATACGGAAATCTGCAACCTGGCGTATCAGGGGGCCAAGGTCa > 8:237871/1‑139 (MQ=255) ggAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTCCTGTTGTGACGTATACGGAAATCTGCAACCTGGCGTATCAGGGGGCCAAGGTCATTtc < 5:6394/138‑1 (MQ=255) ggAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTCCTGTTGTGACGTATACGGAAATCTGCAACCTGGCGTATCAGGGGGCCAAGGTCATTtc > 6:6394/1‑138 (MQ=255) aTACATCGATATTTTCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTCCTGTTGTGACGTATACGGAAATCTGCAACCTGGCGTATCAGGGGGCCAAGGTCATTTCTcc < 7:237871/138‑1 (MQ=255) tattatCACCGATGTAGAAGGCGTGATGACTGCTGACACAAGAGTGGTAGAAAATGCAAGACCGCTTCCTGTTGTGACGTATACGGAAATCTGCAACCTCGCGTATCAGGGGGCCAAGGTCATTTCTCCGCGTGCtgtg > 4:187835/1‑139 (MQ=255) ttttCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTCCTGTTGTGACGTATACGGAAATCTGCAACCTGGCGTATCAGGGGGCCAAGGTCATTTCTCCGCGTGCTGTgg > 8:291250/1‑138 (MQ=255) ttttCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTCCTGTTGTGACGTATACGGAAATCTGCAACCTGGCGTATCAGGGGGCCAAGGTCATTTCTCCGCGTGCTGTGGa < 4:265738/139‑1 (MQ=255) ttttCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTCCTGTTGTGACGTATACGGAAATCTGCAACCTGGCGTATCAGGGGGCCAAGGTCATTTCTCCGCGTGCTGTGGa < 3:295210/139‑1 (MQ=255) tCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTCCTGTTGTGACGTATACGGAAATCTGCAACCTGGCGTATCa > 2:292275/1‑103 (MQ=255) tCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTCCTGTTGTGACGTATACGGAAATCTGCAACCTGGCGTATCa < 1:292275/103‑1 (MQ=255) gATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTCCTGTTGTGACGTATACGGAAATCTGCAACCTGGCGTATCAGGGGGCCAAGGTCATTTCTCCGCGTGCTGTGGAAATCGCaa > 4:112830/1‑139 (MQ=255) tGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTCCTGTTGTGACGTATACGGAAATCTGCAACCTGGCGTATCAGGGGGcc < 1:163986/103‑1 (MQ=255) tGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTCCTGTTGTGACGTATACGGAAATCTGCAACCTGGCGTATCAGGGGGcc > 2:163986/1‑103 (MQ=255) tGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTCCTGTTGTGACGTATACGGAAATCTGCAACCTGGCGTATCAGGGGGCCAAGGTCATTTCTCCGCGTGCTGTGGAAATCGCAATg < 2:36920/139‑1 (MQ=255) tGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTCCTGGTGTGACGTATACGGAAATCTGCAACCTGGCGTATCAGGGGGCCAAGGTCATTTCTCCGCGGGCTGTGGAAATCGCAATg < 2:181513/139‑1 (MQ=255) tAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTCCTGTTGTGACGTATACGGAAATCTGCAACCTGGCGTATCAGGGGGCCAAGGTCATTTCTCCGCGTGCTGTGGAAATCGCAATGCa < 2:145732/139‑1 (MQ=255) | GGTAGTCGTCGCAGGCTTCCAGGGCGCAACTGAAAAGGGTGATACAACAACAATCGGCAGAGGGGGCAGCGATACGTCTGCAGCCGCACTCGGAGCAGCTGTTGATGCGGAATACATCGATATTTTCACCGATGTAGAAGGCGTGATGACTGCTGACCCAAGAGTGGTAGAAAATGCAAAACCGCTTCCTGTTGTGACGTATACGGAAATCTGCAACCTGGCGTATCAGGGGGCCAAGGTCATTTCTCCGCGTGCTGTGGAAATCGCAATGCA > BS168_reference_genome/1747527‑1747799 |
| 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: ‑ |