...
cd /project/microbiome/data_queue/seq/TRNLloc_Testad1/tfmergedreads/16S/TRNL1loc_ad1/trimmed
cp *./*R1.fq /project/microbiome/data_queue/seq/TRNLloc_Testad1/R1only/
cd /project/microbiome/data_queue/seq/TRNLloc_Testad1/R1only
sed -n '1~4s/^@/>/p;2~4p' ./*.fq > ./TRNL_Test_16SLocAdR1.fa
vsearch -derep_fulllength --input TrnlTest16S.fa --output uniqueSequences.fa --sizeout --sizein
Code Block |
---|
vsearch --derep_fulllength $s |
...
LocAdR1.fa \ --strand plus \ --output $s derep.fa \ --sizeout \ --uc $s.derep.uc \ --relabel $s. \ --fasta_width 0 |
...
vsearch --cluster_unoise derep.fa --centroids zotus_vsearch.fa --sizein --sizeout