danielecook
6/26/2015 - 5:17 PM

metaseq_extract_unmapped reads

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