Databases
StackOverflow2013
Postsdbo
POFinding/Replacing substrings with annotations in an ASCII file in Python
 

Note that there are some explanatory texts on larger screens.

plurals
  1. POFinding/Replacing substrings with annotations in an ASCII file in Python
    primarykey
    data
    text
    <p>I'm having a little coding issue in a bioinformatics project I'm working on. Basically, my task is to extract motif sequences from a database and use the information to annotate a sequence alignment file. The alignment file is plain text, so the annotation will not be anything elaborate, at best simply replacing the extracted sequences with asterisks in the alignment file itself. </p> <p>I have a script which scans the database file, extracts all sequences I need, and writes them to an output file. What I need is, given a query, to read these sequences and match them to their corresponding substrings in the ASCII alignment files. Finally, for every occurrence of a motif sequence (substring of a very large string of characters) I would replace motif sequence XXXXXXX with a sequence of asterisks <strong><em>*</em></strong>. </p> <p>The code I am using goes like this (11SGLOBULIN is the name of the protein entry in the database):</p> <pre><code>motif_file = open('/users/myfolder/final motifs_11SGLOBULIN','r') align_file = open('/Users/myfolder/alignmentfiles/11sglobulin.seqs', 'w+') finalmotifs = motif_file.readlines() seqalign = align_file.readlines() for line in seqalign: if motif[i] in seqalign: # I have stored all motifs in a list called "motif" replace(motif, '*****') </code></pre> <p>But instead of replacing each string with a sequence of asterisks, it deletes the entire file. Can anyone see why this is happening? </p> <p>I suspect that the problem may lie in the fact that my ASCII file is basically just one very long list of amino acids, and Python cannot know how to replace a particular substring hidden within a very long string. </p>
    singulars
    1. PO
      singulars
      1. PTAnswer
    1. This table or related slice is empty.
    1. PTQuestion
    1. This table or related slice is empty.
    1. USSpyros
    plurals
    1. This table or related slice is empty.
    1. This table or related slice is empty.
    1. This table or related slice is empty.
    1. PO
      singulars
      1. PTAnswer
    2. PO
      singulars
      1. PTAnswer
    3. PO
      singulars
      1. PTAnswer
    1. VO
      singulars
      1. POFinding/Replacing substrings with annotations in an ASCII file in Python
      1. This table or related slice is empty.
      1. VTUpMod
    2. VO
      singulars
      1. POFinding/Replacing substrings with annotations in an ASCII file in Python
      1. This table or related slice is empty.
      1. VTUpMod
    3. VO
      singulars
      1. POFinding/Replacing substrings with annotations in an ASCII file in Python
      1. This table or related slice is empty.
      1. VTUpMod
    1. COwhat are `i` and `replace`?
      singulars
      1. POFinding/Replacing substrings with annotations in an ASCII file in Python
      1. USgeoffspear
    2. COI knew a Spyros who studied Biochemistry in London between 1996 and 1999... Anyway. What's this `replace` function you're calling? If you want to change the contents of a file, you'll probably be better off creating a new file and then replacing the old file. That way, if there's a problem in execution, you can check exit state and not lose data.
      singulars
      1. POFinding/Replacing substrings with annotations in an ASCII file in Python
      1. USMattH
    3. CO@Wooble, replace is a string processing command to replace one string with another and motif[i] is an indexing method for the list called motif, which contains a set of strings.
      singulars
      1. POFinding/Replacing substrings with annotations in an ASCII file in Python
      1. USSpyros
 

Querying!

 
Guidance
A row detail

Detail views are divided into sections. All the information in the data section comes from columns in the selected row. The other sections display data from other, related rows.

Related data can be related in a to-one or a to-many fashion. Captions of data related in a to-many fashion link to a list view showing a filtered view of the table.

Try moving around until you find a non-empty to-many entry and click on the label to get to one. You can move back to the root by clicking on the database name in the header.

SQuiL has stopped working due to an internal error.

If you are curious you may find further information in the browser console, which is accessible through the devtools (F12).

Reload