...
sed -n '1~4s/^@/>/p;2~4p' ./*.fq > ./TrnlTest16S.fa
vsearch -derep_fulllength --input TrnlTest16S.fa --output uniqueSequences.fa --sizeout --sizein
Code Block |
---|
vsearch --derep_fulllength $s TrnlTest16S.fa \ --strand plus \ --output $s derep.fa \ --sizeout \ --uc $s.derep.uc \ --relabel $s. \ --fasta_width 0 |
...