cat yourname_vs_hg19.align | grep NM | head -500 | cut -f 3 | awk '{print ">" $1 "\n" $1}'
As part of our genotyping of populations we always check 1000 reads from each sample. It often explains some discordant results (lots of reads but low depth at the loci because most the sample is something else!).
Leave a comment: