Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

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