免疫组库

igblast

conda  create -n repe -y  python=2
source activate repe
conda install -y igblast flash 
conda install -c bioconda seqkit 
for i in `seq -w 7 10`;
do prefetch ERR34450${i} \ done
for i in `seq -w 7 10`;
fastq-dump -A ERR34450${i}  -O ../raw --gzip --split-e ../ERR34450${i} | trim_galore -q 25 --phred33 --length 36 --stringency 3 --paired -o clean ERR34450${i}_1.fastq.gz ERR34450${i}_2.fastq.gz 
for i in `seq -w 7 10`;
flash ERR34450${i}_val_1.fastq.gz  ERR34450${i}_val_2.fastq.gz  -M 250 -o ./flash/ERR34450${i}; done
for i in `seq -w 7 10`; do seqkit fq2fa ERR34450${i}.extendedFrags.fastq > ERR34450${i}.fa; done

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。