Difference between revisions of "User:Timothee Flutre/Notebook/Postdoc/2012/05/25"

From OpenWetWare
Jump to: navigation, search
(Autocreate 2012/05/25 Entry for User:Timothee_Flutre/Notebook/Postdoc)
 
(Entry title)
Line 6: Line 6:
 
| colspan="2"|
 
| colspan="2"|
 
<!-- ##### DO NOT edit above this line unless you know what you are doing. ##### -->
 
<!-- ##### DO NOT edit above this line unless you know what you are doing. ##### -->
==Entry title==
+
==Awk one-liner to convert fasta file into fastq format==
* Insert content here...
 
  
 +
* We can use the built-in variable "RS" ([http://www.gnu.org/software/gawk/manual/gawk.html#Records split records]) and use "split" ([http://www.gnu.org/software/gawk/manual/gawk.html#String-Functions string function]):
 +
 +
<nowiki>
 +
awk 'BEGIN{RS=">"} {if(NF==0)next; split($0,a,"\n"); printf "@"a[1]"\n"a[2]"\n+\n"; \
 +
for(i=1;i<=length(a[2]);i++)printf "}"; printf"\n"}' probes.fa > probes.fq
 +
</nowiki>
  
 
<!-- ##### DO NOT edit below this line unless you know what you are doing. ##### -->
 
<!-- ##### DO NOT edit below this line unless you know what you are doing. ##### -->

Revision as of 13:09, 25 May 2012

Owwnotebook icon.png Project name <html><img src="/images/9/94/Report.png" border="0" /></html> Main project page
<html><img src="/images/c/c3/Resultset_previous.png" border="0" /></html>Previous entry<html>      </html>Next entry<html><img src="/images/5/5c/Resultset_next.png" border="0" /></html>

Awk one-liner to convert fasta file into fastq format

awk 'BEGIN{RS=">"} {if(NF==0)next; split($0,a,"\n"); printf "@"a[1]"\n"a[2]"\n+\n"; \
for(i=1;i<=length(a[2]);i++)printf "}"; printf"\n"}' probes.fa > probes.fq