Generated a simple qsub script using echo:
echo "zcat input.fq.gz | awk {if (NR % 4 == 0 || NR % 4 == 2) print substr($0, 0, 91); else print} \| gzip > output.fq.gz" > process.sh
In fact, the special variable $0 was first interpreted by echo instead of awk, which makes the actual script looks like this:
zcat input.fq.gz | awk {if (NR % 4 == 0 || NR % 4 == 2) print substr(bash, 0, 91); else print} \| gzip > output.fq.gz
As a result, the script prints nothing between read name and the plus symbol, like this:
Read 1
[blank line]
+
[blank line]
Read 2
[blank line]
+
[blank line]
...
3 Responses
escort bayan
I believe you have mentioned some very interesting points, appreciate it for the post. Edeline Sebastian Jed
erotik
Hey there. I discovered your site by way of Google even as looking for a similar topic, your site came up. It seems great. I have bookmarked it in my google bookmarks to visit then. Viviana Salomon Townshend
auto glass san diego
I like to see that people are expanding their idea of what a basement should look like. Nice job Rob. Maritsa Stanislaus Lorenza