hello everyone, im new in coding areas and I have a question. I have a bedfile for a complex disease(example obesity, type 2 diabetes) releated snp variants and I have another file for 10 people genome data. I want to calculate my variants 10kb arounds genome coverage for this 10 people genome data. My language is not okey, I know that but I hope you all understand my problem. For this problem I want to write a script for linux but I don't know where to start or what to do. Thanks for helping.
its bam file, actually they are ancient people form different ages. I want tto understand if that person carries variants that increase the likelihood of the disease occurring. my bed file is look like this , ı also look my variants chromosome positions 10kb around. chr start end chr1 219465882 219475882 chr1 65508597 65518597 chr1 177878445 177888445 chr1 72367846 72377846 chr1 72280502 72290502 chr1 72294433 72304433 chr1 177914890 177924890 chr2 626183 636183 . . .
Ok, then as said, run them through a variant caller and intersect them with the BED files. The VCF contains depth for REF and ALT.
thanks al lot for helping !