smithlabcode / methpipe

A pipeline for analyzing DNA methylation data from bisulfite sequencing.
http://smithlabresearch.org/methpipe
67 stars 27 forks source link

Getting started with amrtester #209

Closed ctsa closed 2 years ago

ctsa commented 2 years ago

I'm having trouble getting a first test run working with amrtester. I'm working on GRCh38. I have a bed file of test regions and a bam produced by bwa-meth. I am having trouble setting the required -chrom argument correctly. So far I've tried submitting the reference fasta, the fai, the output from fetchChromSizes hg38, and just the string 'hg38', but the output I get in all cases is the following:

NUMBER OF REGIONS: 31448
no chrom file for:
chr1    28735   29737   CPGI    0       +

Do you have a small working demo of this tool to illustrate the expected format?