![]() >sequence1 ATGCGGAGCTTAGATTCTCGAGATCTCGATATCGCGCTTATAAAAGGCCCGGATTAGGGCTAGCTAGATATCGCGATAGCTAGGGATATCGAGATGCGATAC If you really want to keep the leading > of the header (which doesn't seem very useful), you could do something like this: $ perl -0pe 's/\n//g s/.>/\n>/g s/$/\n/ ' file.fa This can be a very useful trick when searching a fasta file for a string: TblToFasta file.fa | grep 'foo' | FastaToTbl Sequence2 GTACTCGATACGCTACGCGATATTGCGCGATACGCATAGCTAACGATCGACTAGTGATGCATAGAGCTAGATCAGCTACGATAGCATCGATCGACTACGATCAGCATCACĪnd, to get the Fasta back: $ FastaToTbl file.fa | TblToFasta Sequence1 ATGCGGAGCTTAGATTCTCGAGATCTCGATATCGCGCTTATAAAAGGCCCGGATTAGGGCTAGCTAGATATCGCGATAGCTAGGGATATCGAGATGCGATACG GTACTCGATACGCTACGCGATATTGCGCGATACGCATAGCTAACGATCGACTAGTGATGCĪTAGAGCTAGATCAGCTACGATAGCATCGATCGACTACGATCAGCATCAC ![]() TAGCTAGATATCGCGATAGCTAGGGATATCGAGATGCGATACG Printf "\n%s\t", substr($0,2,length($0)-1)ĪTGCGGAGCTTAGATTCTCGAGATCTCGATATCGCGCTTATAAAAGGCCCGGATTAGGGC If you have multi-line fasta files, as is very common, you can use these scripts 1 to convert between fasta and tbl (sequence_name sequence) format:
0 Comments
Leave a Reply. |