This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Primer3_core

How to make fasta file as input to primer3_core using bash command line

line command primer3_core
• 1,263 views • link
updated by GenoMax GenoMax • written by Mohamed • Mohamed

1 answer 

awk '/^>/ {printf("%s%s\t",(N>0?"\n":""),$0);N++;next;} {printf("%s",$0);} END {printf("\n");}' < input.fa  |\
awk '{printf("SEQUENCE_ID=%s\nSEQUENCE_TEMPLATE=%s\n=\n",substr($1,2),$2);}' |\
primer3_core
• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Thank you so much, this command works perfectly,

awk '/^>/ {printf("%s%s\t",(N>0?"\n":""),$0);N++;next;} {printf("%s",$0);} END {printf("\n");}' < M_gen_seq.fa  |awk '{printf("SEQUENCE_TEMPLATE=%s\nSSEQUENCE_ID=%s\n=\n",substr($1,2),$2);}' | primer3_core
• 0 views • link
updated by GenoMax GenoMax • written by Mohamed • Mohamed

Log in to answer this question.

More posts like this