| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | BS168_reference_genome | 3,907,393 | G→C | 30.3% | intergenic (‑79/‑451) | ywcI ← / → vpr | conserved protein of unknown function/extracellular serine protease |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | BS168_reference_genome | 3,907,393 | 0 | G | C | 30.3% | 45.1 / 3.9 | 32 | intergenic (‑79/‑451) | ywcI/vpr | conserved protein of unknown function/extracellular serine protease |
| Reads supporting (aligned to +/- strand): ref base G (14/8); new base C (2/8); total (16/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.38e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTA > BS168_reference_genome/3907260‑3907509 | cccAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCACTCCGCTTCGATTTTGCAGGCACCCTTTTTAGCAGCCCGTTTTCTCAACCGCAGCCCGCAAGAAAACACGGc > 1:188962/1‑139 (MQ=255) atcaAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttca > 1:42155/1‑138 (MQ=255) ccATCCACACCGTTAAATATTCACCCAACCGTTTCATACAAAAACATCCCTCCGCTTCTTTTTGGCAGGCACCCTTTTTCCCAGCCGTTTTTCTCCCCCGCCGCCCGAAACAAACGACGGCCGATTTCTTTTCTCGCaa < 4:190496/139‑1 (MQ=255) ccACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAAc > 6:252260/1‑137 (MQ=255) ccccccGAAAACATACAACCAAACGTTTCAAATAAAAACATCCCTCCGCTGCTTTTTGCCAGGCAGCCTTTTACGCCCCCCGTTTTCTCCGCCGCACCCCGCAACAAAACACGGCCGATTTCTTTTCTCGCAAACTAGc < 6:92756/136‑1 (MQ=255) acCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCgg > 7:42893/1‑139 (MQ=255) cccGTAAACAAACACGCAAACGTTCCAACGAAAAACACCCCTCCGCTGCTTTTTGGCCGGCACCCTTTTAACCAGCCCGTTTACTCAGCCGCAGCCCGCAAGAAACGACGGCCGATTTCTTTTCTCGCAAACTggcggc < 3:340010/139‑1 (MQ=255) cGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTGATTTCTCCCCAACAAGCGCctc > 3:190496/1‑139 (MQ=255) gTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCCGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCCGCCCGCAAGCAAAGACGGCCGATGTCGGTTCTCGCCACCTAGCGGctc > 8:74582/1‑138 (MQ=255) tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCGTTTATAGCAAACTAGCGGctc > 3:56735/1‑137 (MQ=255) aaGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCAGATTTCTTTTCTCGCAAACTAGCGGCTCTAg > 8:191303/1‑138 (MQ=255) acaTACACCCAACCGTTTAATACAAAAACCTCCCTCCGCTTATTTTTGGCCGGCACCCTTTTTCGCAGCCCTTTTTCCCCCCCCCCCCCCCCAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa < 3:7168/137‑1 (MQ=255) aTACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTGTAGCaa > 4:305393/1‑138 (MQ=255) aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaaa > 1:190407/1‑137 (MQ=255) aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAAATAGCGGCTCTAGCaaaa > 4:303278/1‑138 (MQ=255) gCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGAAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt < 7:57598/139‑1 (MQ=255) ccAAACGTTTCATACAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCATTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAACAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt < 7:115766/138‑1 (MQ=255) caaaaGGTTTAATCTAAAAACATCCCCCCGATTCTTTTTGGCCGGCCGCTTTTTTAGCCGCCGTTTTTCACACCCGCACCCCCAAATAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt < 2:190407/137‑1 (MQ=255) caaaaCTTTTAATACAAAAACATCCCTCCGCTTCTTTTTGCCTGCCAGCTTTTTTAGCCGCCGTTTTTCTCCCCCCCACCCCGCAATAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt < 3:303278/137‑1 (MQ=255) aaccttttCATACAAAAACATCCCTCCGATTTTTTTTGGCAGCCAGCCTTTTTAGCAGCCGGTTTTCCCCGCCGCCGCCCGCAATAAAACACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGa < 1:25990/134‑1 (MQ=255) aaaCGTTTCATATAAAAACAGCCCTCCGCTTCTTTTTGGCAGGCTGCCTTTTTAGCCGCCCGTTTTCTCAGCCCCAGCCCGCAATAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGa < 7:788/139‑1 (MQ=255) ccGTTTCATATAAACACATCCCTCCGCTTATTTTTGGCAGCCACCCTTTTTCCCAGCCCGTTTTCTCACCCGCACCCCCAAACAAAACACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct < 7:207929/138‑1 (MQ=255) ccGTTTCATACAAAACCCCCCCGCCGCTGCTTTGTGCCCGGCAGCGTTTTTCGCCCCCCGTTTCCCCCGCCCCCCCCCCAAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct < 1:305433/138‑1 (MQ=255) ccGTTTCATACAAAAACATCCCTCCGCTTCTTTTTGGCCGGCCGCTTTTTTAGCAGCCGGTTTTCTCACCCCCCCCCCGAAATAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct < 5:252260/138‑1 (MQ=255) aCGTTTCATAGAAAAACATCCCTCCGCTTATTTTTGCCAGGCAGCATTTTTCGCCGCCGTTTTTCTCACCCCCCGCCCGCAATAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct < 8:42893/139‑1 (MQ=255) aCGTTTAATACAAAAACATCCCTCCGCTTCTTTTTGGCCGGACGCCTTTTTAGCCCCCCGTTTGCTCCCCCCCCCCCCGCAACAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct < 3:156183/139‑1 (MQ=255) aCGTTCAATACAAAACAAACCCTCCTGTTTTTTTTGGCAGGCTGCCTTTTACGCAGCCGTTTTTCTCCGCCGCAGCCGGAAATAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct < 7:338314/139‑1 (MQ=255) cGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCCGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAACACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTAGAGCAAAACTCAAGActc > 8:49280/1‑139 (MQ=255) ggcaggcaGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGtga > 3:54012/1‑136 (MQ=255) agcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAAcc > 7:230389/1‑138 (MQ=255) cccGTTTTCTCAGCCG‑AGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAAc > 7:13104/1‑138 (MQ=255) cAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTACtttttt > 7:257115/1‑137 (MQ=255) | CCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTA > BS168_reference_genome/3907260‑3907509 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |