I have a fasta file that is formatted in the following way:
> gene1
atctgtctgct
atcgtc
at
and I want to put it in a format like the following:
> gene1 atctgtctgctatcgtcat
I am struggling to get rid of the whitespace between the lines under the header. Does anyone have any idea how I can do this in python?
•
written
by
a.rex
3 answers
Do you want to convert FASTA format to tab-delimited format?
Writing Python script with Biopython or Perl with BioPerl is very convenient. It can also be achieved by one or more shell commands.
I'd like to introduce the FASTA/Q toolkit SeqKit, which can do this with one command:
seqkit fxtab seqs.fa > formated.txt
Since spaces exited in the sequence of your sample data, a cleaning step was used to remove the spaces.
And seqkit fx2tab outputs 3 columns for compatibility of FASTA and FASTQ, cut was used to remove the empty third column:
$ seqkit seq --remove-gaps seqs.fa | seqkit fx2tab | cut -f 1,2
gene1 atctgtctgctatcgtcat
After manipulations of the tabular format, you can use seqkit tab2fx to convert it back to FASTA format.
$ pip install pyfaidx
$ faidx --transform transposed input.fa | cut -f1,4 > out.tab
If you really want to do this, but I agree it might be better to answer the real question you might have not yet asked...
fold -w 60 input.fa
Log in to answer this question.
-
Global alignment of protein sequences all x all to identify similarities between them
written by Bruna Avelino •
-
How to find conserved enhancers between human cell lines?
written by szp770 •
-
is the sequences data in GISAID updated ?
written by 2001linana •
-
Extracting fasta sequence from CDS matching with sequence ID
written by metallicasster •
-
How to split sequence in fasta file based on long runs of Ns?
written by stacy734 •
-
Extracting subsequence from FASTA file using python
written by shawnt1234 •
-
Extract out fasta sequences using sequence headers.
written by a.rex
-
recon issue for repeatmodeler
written by a.rex
-
Replace tip names in a newick file
written by tlorin
-
Add a seq (the outgroup) to the top of its corresponding gene
written by hosseinv •
Why? Why not leave it in fasta (a format most things already accept)?
As Devon is getting at, what ever method you're looking to use to 'read' this FASTA file is probably not a good one. Is it, by any chance, awk/sed/td/grep?
If you tell us what your bigger problem is that you're looking to solve, we might be able to help you there. But creating a new format, unless for a really good reason/implementation, is generally a bad idea for everyone.
You should look into using Biopython instead of trying to use text parsing. See these pages:
http://www.bioinformatics.org/bradstuff/bp/tut/Tutorial002.html
http://biopython.org/wiki/SeqIO
http://biopython.org/DIST/docs/tutorial/Tutorial.html#htoc11
http://biopython.org/DIST/docs/api/Bio.SeqIO-module.html