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

chromosome name #170

Open
sotaro-kanematsu opened this issue Oct 12, 2021 · 1 comment
Open

chromosome name #170

sotaro-kanematsu opened this issue Oct 12, 2021 · 1 comment

Comments

@sotaro-kanematsu
Copy link

I am using facets to uncover allele-specific somatic copy number alternations using clinical tumor-normal pairs.
I run the cnv_facets.R script as below.

cnv_facets.R --gbuild hg19 -t Tumor.bam -n Normal.bam -vcf common_all_20170710.vcf.gz -o Nomal_Tumor_cnv_facets_output

I used hs37d5a.fa as a reference sequence in mapping, and so my bam files do not have prefixes (e.g. 1, 2, 3.....Y) as well as the common SNP vcf file.
On the other hand, hg19 has prefixes (e.g. chr1, chr2...chrY), so If I run the above script, an error occurs. (I know the script work well if I change the header of bam files and add "chr" prefix to vcf file, however I thought it is supposed inefficient to fix bam files or change the reference.
How should I fix the cnv_facets.R script?

@veseshan
Copy link
Collaborator

Since cnv_facets is at https://github.com/dariober/cnv_facets have you asked @dariober

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