Edit Seq object with Biopython
I am using SeqIO to load a fasta file and play my sequences. For the sake of my analysis, I want to eliminate a fixed number of bases to the left and right of my sequences. Nevertheless, whenever I try to eliminate the last ten bases I get an error saying
sequences = list (SeqIO.parse(file("myfile.fasta"), "fasta")
for sequence in sequences:
del sequence.seq[-10:]
TypeError: 'Seq' object doesn't support item deletion
I assume trimming borders and doing other edits are very basic functionalities in bioinformatics (the same can be done in Perl in a couple of lines) but I am bit surprised I can't find a method to do such things in Python because the Seq objects are immutable. Is there a way to trim a Seq object?
• 2,585 views
•
link
updated
by
Ram
•
written
by
joreamayarom
•
written
by
joreamayarom
1 answer
Log in to answer this question.
More posts like this
-
Forum: WIG TO BIGWIG - UCSC SOFTWARE - ERROR IN "ITEM ENDS". URGENT
written by matteo.zoia •
-
Delete character from sequence id
written by ericcabarroso •
-
parse of file multifasta in python
written by erick_rc93 •
-
TypeError: 'str' object is not callable
written by S AR
-
Any tools that can score for alignment of two short sequences allowing ambiguous base
written by Zealseeker •
-
Using Bio::SeqIO to write with fasta format
written by mahmood.nt •
-
Problems when parsing VCF file with Python
written by joreamayarom
-
Split fasta to multi fasta
written by felipelira3 •
-
Most efficient way to trim overhanging bases after alignment
written by joreamayarom
-
Why does Seq accept only a string as a sequence?
written by athinkerer •
