nloyfer / wgbs_tools

tools for working with Bisulfite Sequencing data while preserving reads intrinsic dependencies
Other
125 stars 33 forks source link

segment sites mismatch #64

Open chang-bwl opened 3 months ago

chang-bwl commented 3 months ago

Hello, I tried to split the beta (hg19) of the cell_WGBS in the paper : $ wgbstools segment --betas *beta --min_cpg 3 --max_bp 2000 -o blocks.small.bed

The following error occurred: [wt beta] WARNING: beta file size (28,587,300 sites) incomatible with current genome reference (28,217,448 sites) Invalid input argument [wt segment] ERROR: current genome reference (hg19) does not match the input beta file (GSM5652208_Coronary-Artery-Smooth-Muscle-Z00000420.beta).

How can I solve it? Thanks

chang-bwl commented 3 months ago

I've finally solved this problem.It's because wget -c downloads have some extra lines in the download file, which causes the CSI to lag