# done on bufo using cdhit version 4.6

cdhit-est \
-i ../../_A_01_renameFASTA-Bioinfo/output/soft_filtered_transcripts_tr1tr2.fasta \
-o ../output/soft_filtered_transcripts_tr1tr2.cdhitest.fasta \
-c 1 \
-M 100000 \
-t 30


#input tr1tr2 sequences
grep -c ">" ../../_A_01_renameFASTA-Bioinfo/output/soft_filtered_transcripts_tr1tr2.fasta
#451770

#output of unique sequences
grep -c ">" ../output/soft_filtered_transcripts_tr1tr2.cdhitest.fasta
#389962 (68.3% of input)