#!/bin/bash

grep "^>" Dry_Wet.cluster.faa | tr ">" " " |awk '{print $1}' > Dry_Wet.geneids

seqkit grep -f Dry_Wet.geneids Wet2014.ffn > Dry_Wet.ffn
seqkit grep -f Dry_Wet.geneids Dry2014.ffn >> Dry_Wet.ffn

bwa index Dry_Wet.ffn