We usually want to maintain header when we do something for the body of the text.
You can find a body function here which can only accept streams and assumes only one line header.
Use subshells and awk can do the same job and potentially more flexible.
The original credits go to Aaron Quinlan. see a gist below to sort vcf files in natural chromosome order :chr1 chr2 chr3.... rather than chr1 chr10 chr11....
This blog by Tommy Tang is licensed under a Creative Commons Attribution-ShareAlike 4.0 International License.
No comments:
Post a Comment