Kalain Hündin

Ceci est une ancienne révision du document !

for filename in *.fa; do sed -i '/>Recurrent fasta name/c\>'$filename'' $filename; done

all VCF to vcf.gz

ls *.vcf | xargs -n1 -P0 bgzip

AWK

awk '{if(NR%4==2) {count++; bases += length} } END{print bases/count}' <fastq_file>