Thankx for reply

I used this script

convert the sam file to a mapability and gc-content tracks

perl mapAndGc.pl ~/path to dir/outdir/chr1.sam canfam3.1.genome $windowSize $chr $readLength $outdir

gzip $outdir/$chr.gc

gzip $outdir/$chr.map

rm -f $outdir/$chr.sam

Without more information, I'm afraid that I can't help you. The scripts provided are short and straightforward. Your error relates to not being able to find specific files that it expects to be there. Make sure all preceding steps have completed successfully. If you go line by line through the process, verifying that your inputs and outputs from each step look sane, I imagine you'll figure out what's happening.

Actually I used "createCustomAnnotations" scripts from here for non human data:

https://xfer.genome.wustl.edu/gxfer1/project/cancer-genomics/readDepth/index.html

I have bam files and reference genome, so i convert bam file into sam file and run only last command of "runEachchr.sh"

!/gsc/bin/bash

arguments

chr=chromosome37 readLength=100 fastaDir=$3 outdir=outdir windowSize=100

generate all possible reads of the given size from that chromsome

perl allSeq.pl $fastaDir/$chr.fa $readLength > $outdir/$chr.fastq

align those reads using bwa

bwa aln -t 4 $fastaDir/all_sequences.fa $outdir/$chr.fastq >$outdir/$chr.aln.sai

bwa samse -n 1 $fastaDir/all_sequences.fa $outdir/$chr.aln.sai $outdir/$chr.fastq >$outdir/$chr.sam

clean up

rm -f $outdir/$chr.fastq $outdir/$chr.aln.sai

convert the sam file to a mapability and gc-content tracks

perl mapAndGc.pl /home/misbah/InstalledUtilites/Dogsd_bamfiles/outdir/chr37.sam canfam3.1.txt $windowSize $chr $readLength $outdir gzip $outdir/$chr.gc gzip $outdir/$chr.map rm -f $outdir/$chr.sam

but it give error

Use of uninitialized value $len in division (/) at mapAndGc.pl line 83, <gen1> line 2097423.

i am in trouble, please tell me about this issue

thankx