Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

B-allele frequencies plot #14

Open
thangnx1012 opened this issue Feb 24, 2025 · 2 comments
Open

B-allele frequencies plot #14

thangnx1012 opened this issue Feb 24, 2025 · 2 comments

Comments

@thangnx1012
Copy link

Hi,
Thank you so much for developing the tools.
I have a question related to the Wakhan plot, and I could not see the B-allele frequencies in my final plot. Here is my code to run the plot
python wakhan.py --threads 50
--reference ${hg38_reference}
--target-bam ${longphase_clair3_outdir}/${samplename}_haplotagged.bam
--tumor-vcf ${clair3_outdir}/${samplename}/phased_merge_output.vcf.gz
--genome-name ${samplename}
--reference-name hg38
--contigs chr1-22,X,Y
--loh-enable
--pdf-enable
--out-dir-plots ${wakhan_outdir}/${samplename}
--breakpoints ${severus_outdir}/${samplename}/all_SVs/severus_all.vcf

Here is the plot generated from this code.
https://drive.google.com/drive/folders/1cuhhqe_gdA1r2SdgqQkaaZJpRZkPsHrY?usp=sharing
I do see the B-allele frequencies in each chromosome plot but not in the summary plot. Thank you

@tahashmi
Copy link
Collaborator

Hi @thangnx1012 , actually you have to enable b-allele freq by clicking on BAF legend in main plots (these are disabled by default for visualization as it can take some extra memory in browser). Similarly what you see in red/blue dots are haplotype specific coverage bins. I will update the doc soon.

@thangnx1012
Copy link
Author

Oh, great! I did!
Thanks for the quick reply. I only look at the PDF file, so I don't think about that option in the HTML file.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants