methylKit icon indicating copy to clipboard operation
methylKit copied to clipboard

Error: The sam file is not sorted properly on positions

Open iranianuser opened this issue 4 months ago • 10 comments

Hello I m using methylkit on bismark alignment files. I used the command below and used the bam output from bismark alignment process: methRaw = processBismarkAln( location = "/media/genome2/SP PHD U3/vahedi/802/SRR4298802_1_val_1_bismark_bt2_pe.bam", sample.id="Cancer_802", assembly = "hg38", treatment = (1), read.context = "CpG") I have 3 cancerous samples and 1 control but first I test just one of my samples. I got this Error :The sam file is not sorted properly on positions: chr: 11 pos: 36283681 is followed by chr: 11 pos: 36283418 You can sort the file in unix-like machines using: grep -v \'^[[:space:]]*\@\' test.sam | sort -k3,3 -k4,4n > test.sorted.sam

what can I do to solve this problem??

thanks alot

iranianuser avatar Oct 03 '24 16:10 iranianuser