grep \ACAGGTTGGCTGATAAGTCCCCGGTCTGACACATAGATGG ZY392.clip.1.fq |awk '{s++}{print ">r1_s"s"\n"$0}' >r1.IR-L.fa
grep \TGTTTTTCTTTGTTAGACCGGGGACTTATCAGCCAACCTG ZY392.clip.1.fq |awk '{s++}{print ">r1_s"s"\n"$0}' >r1.IR-R.fa
grep \ACAGGTTGGCTGATAAGTCCCCGGTCTGACACATAGATGG ZY392.clip.2.fq |awk '{s++}{print ">r2_s"s"\n"$0}' >r2.IR-L.fa
grep \TGTTTTTCTTTGTTAGACCGGGGACTTATCAGCCAACCTG ZY392.clip.2.fq |awk '{s++}{print ">r2_s"s"\n"$0}' >r2.IR-R.fa
formatdb -i OE1-1.fa -p F
blastn -db OE1-1.fa -query r1.IR-L.fa -out r1.IR-L.m8 -max_target_seqs 5 -evalue 1e-5 -outfmt 6 -num_threads 8 &
blastn -db OE1-1.fa -query r1.IR-R.fa -out r1.IR-R.m8 -max_target_seqs 5 -evalue 1e-5 -outfmt 6 -num_threads 8 &
blastn -db OE1-1.fa -query r2.IR-L.fa -out r2.IR-L.m8 -max_target_seqs 5 -evalue 1e-5 -outfmt 6 -num_threads 8 &
blastn -db OE1-1.fa -query r2.IR-R.fa -out r2.IR-R.m8 -max_target_seqs 5 -evalue 1e-5 -outfmt 6 -num_threads 8 &
ZY系列,提取barcode,比对到参考基因组
最新推荐文章于 2024-09-04 10:07:44 发布