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

beta2bed warning - incomatible (sic) with current genome reference #67

Open
cb4github opened this issue Apr 5, 2024 · 0 comments
Open

Comments

@cb4github
Copy link

cb4github commented Apr 5, 2024

Dear All,

Thanks for your efforts.

Steps follow.
# Clone
git clone https://github.com/nloyfer/wgbs_tools.git
cd wgbs_tools

# compile
python setup.py
module load anaconda3/5.1.0 gcc/8.5.0 samtools/1.16.1 htslib/1.19
export PATH=$PATH:$PWD
cd
wgbstools init_genome --force --no_default --debug --threads 19 hg38
wget --directory-prefix=data ftp://ftp.ncbi.nlm.nih.gov/geo/samples/GSM5652nnn/GSM5652231/suppl/GSM5652231_Cerebellum-Neuron-Z000000TB.hg38.beta
mkdir output
wgbstools beta2bed --outpath output/GSM5652231_Cerebellum-Neuron-Z000000TB.hg38.bed --genome hg38 data/GSM5652231_Cerebellum-Neuron-Z000000TB.hg38.beta

Output follows.
[wt beta] WARNING: beta file size (29,401,795 sites)
incomatible with current genome reference (1,822 sites)

Also, the bed file apparently is successfully created.
$ls -lt output
total 20
-rw-r--r-- 1 18740 Apr 5 15:02 GSM5652231_Cerebellum-Neuron-Z000000TB.hg38.bed

Please advise. Thanks.

Best,
CB

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

1 participant