gunzip *
for shname in `ls *R1.fastq`;
do samplename=`echo "$shname" | awk -F_ '{print $1}'`;
less $shname|awk '{ if(NR%4==1) gsub(/@/,"@'$samplename'_"(NR+3)/4" ");
print $1}'>$samplename.R1.fq ;
done
for shname in `ls *R2.fastq`;
do samplename=`echo "$shname" | awk -F_ '{print $1}'`;
less $shname|awk '{ if(NR%4==1) gsub(/@/,"@'$samplename'_"(NR+3)/4" ");
print $1}'>$samplename.R2.fq ;
done
gzip *.fq
ls *.R1.fq.gz | awk -F. '{print "mv "$0" "$1"_R1.fastq.gz"}' |bash
ls *.R2.fq.gz | awk -F. '{print "mv "$0" "$1"_R2.fastq.gz"}' |bash
rm *.fastq
#基因组
Instrument: Illumina HiSeq 1000
Strategy: WGS
Source: GENOMIC
Selection: RANDOM
Layout: PAIRED
#16s
Instrument: Illumina MiSeq
Strategy: AMPLICON
Source: METAGENOMIC
Selection: PCR
Layout: PAIRED
#宏基因组
Instrument: Illumina HiSeq 4000
Strategy: WGS
Source: METAGENOMIC
Selection: cDNA
Layout: PAIRED
# Further information and requests for resources and reagents should be directed to Guopeng Lin ([email protected]).
暂无评论