1#!/bin/bash 2 3paste <(mawk '/Contig.*:/' $1 ) <(mawk '/Over/' $1) | mawk '$4 > 0.95' | cut -f1 | sed 's/://g' > Hap.loc.belowMAF05 4 5LOC=$(mawk '/Contig/' $2 | wc -l) ####This gives you the number of loci, use it in the next line of code 6 7paste <(seq $LOC ) <(mawk '/Contig/' $2) > BS.loci.numbered 8 9grep -wf Hap.loc.belowMAF05 <(sed 's/://g' BS.loci.numbered) | cut -f1 > bs.loci.lowmaf 10 11