metaseq_extract_unmapped reads
for i in `ls *.bam`; do # Fetch ummapped reads (unmapped reads and pairs) samtools view -f uU ${i}.bam | awk 'BEGIN {FS="\t"} {print "@" $1 "\n" $10 "\n+\n" $11}' | gzip > ../unmapped/$i.gz