Hello,
I have problem with VCFtools running
I have code like this, but tabix is not recognizing input file as bgzf
(however, i am bgzipping it)
#!bin/bash
module load bio/vcftools/0.1.17
module load bio/samtools/1.10
export inpath=/dir
export outpath=/dir
mkdir -p $outpath
chmod u+rwx $outpath
export reference=/ref
$reference
for file in $inpath/* ;
do
bname=$(basename $file)
echo "base name is $bname"
bfile=$outpath/$bname".bgz"
tabfile=$outpath/$bname".tabix"
outfile=$outpath/$bname".fasta"
bgzip $file $bfile
cd $outpath
tabix -f $bfile > $tabfile
cat $reference | vcf-consensus $tabfile > $outfile
done
--
Best regards
Baiba Alkšere
_______________________________________________
Samtools-help mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/samtools-help
