awk

最后发布时间:2023-12-14 14:22:02 浏览量:
/data/MGI2000/V350200125/L04/V350200125_L04_undecoded_1.fq.gz

TAACACGACGCGTCGATGAC
TGTTCTCTTCATATAAGGCG
CAGATCATTTATTTTCTCGC

barcodes.txt
CTAATAACTTAT
paste <(zcat /data/MGI2000/V350200125/L04/V350200125_L04_undecoded_1.fq.gz|paste - - - -) <(zcat /data/MGI2000/V350200125/L04/V350200125_L04_undecoded_2.fq.gz|paste - - - - )| \
 awk -v FS="\t" -v OFS="\n" 'FNR==NR {samples[$2]=$1; next} {barcode = substr($6,0,8); if(samples[barcode]) { print $1,$2,$3,$4>>samples[barcode]"_1.fastq"; print $5,$6,$7,$8>>samples[barcode]"_2.fastq"}}' samples.txt -

paste <(zcat /data/MGI2000/V350200125/L04/V350200125_L04_undecoded_1.fq.gz|paste - - - -) <(zcat /data/MGI2000/V350200125/L04/V350200125_L04_undecoded_2.fq.gz|paste - - - - )| \
 awk -v FS="\t" -v OFS="\n" 'BEGIN { myArray["TAACACGACGCGTCGATGAC"] = 385 ; myArray["TGTTCTCTTCATATAAGGCG"] = 386;}{barcode=substr($6,151,170); if(myArray[barcode]) print $1,$2,$3,$4>>myArray[barcode]"_1.fastq"; print $5,$6,$7,$8>>myArray[barcode]"_2.fastq" }'  -



awk -v FS="\t" -v OFS="\n" 'BEGIN { myArray["TAACACGACGCGTCGATGAC"] = 385 ; myArray["TGTTCTCTTCATATAAGGCG"] = 386;barcode="TGTTCTCTTCATATAAGGCG" ;if(myArray[barcode]) print myArray[barcode]; }' samples.txt -

paste <(zcat /data/MGI2000/V350200125/L04/V350200125_L04_undecoded_1.fq.gz|paste - - - -) <(zcat /data/MGI2000/V350200125/L04/V350200125_L04_undecoded_2.fq.gz|paste - - - - )| \
     awk -v FS="\t" -v OFS="\n"  '{ barcode = substr($6,0,8); print $barcode }' samples.txt   -


paste <(zcat /data/MGI2000/V350200125/L04/V350200125_L04_undecoded_1.fq.gz|paste - - - -) <(zcat /data/MGI2000/V350200125/L04/V350200125_L04_undecoded_2.fq.gz|paste - - - - )| \
 awk 'BEGIN { myArray["TAACACGACGCGTCGATGAC"] = 385 ; myArray["TGTTCTCTTCATATAAGGCG"] = 386;barcode="TGTTCTCTTCATATAAGGCG" ;if(myArray[barcode]) print myArray[barcode]; }' 

快捷入口
Linux 思维导图 浏览PDF 下载PDF
分享到:
标签