I have to do alignments of many sequences to a genome, keep the aligned sequences and align them to other databases and keep unmapped. I have written scripts to do so in loops (
