Note that there are some explanatory texts on larger screens.

plurals
  1. POloop fails when increasing string from elements in list
    primarykey
    data
    text
    <p>I have some problem with sequence generator. I have a file where each line contain one fragment (8 letters). I load it from file in to list, where each element is one fragment. It is DNA so it should go that way:<br><br> 1. Takes first 8-letter element<br> 2. Check for element in which first 7 letters is the same as last 7 letters in first.<br> 3. Add 8th letter from second element in to sequence.<br></p> <p>It should look like this:</p> <pre><code>ATTGCCAT TTGCCATA TGCAATAC </code></pre> <p>So sequence: ATTGCCATAC</p> <p>Unfortunately it only add one element. :( First element is given (we knew it). I do it that way its first in file (first line).</p> <p>Here is the code:</p> <pre><code>from os import sys import random def frag_get(seqfile): frags = [] f_in = open(seqfile, "r") for i in f_in.readlines(): frags.append(i.strip()) f_in.close() return frags def frag_list_shuffle(frags): random.shuffle(frags) return frags def seq_build(first, frags): seq = first for f in frags: if seq[-7:] == f[:7]: seq += f[-1:] return seq def errors(): pass if __name__ == "__main__": frags = frag_get(sys.argv[1]) first = frags[0] frags.remove(first) frags = frag_list_shuffle(frags) seq = seq_build(first, frags) check(sys.argv[2], seq) spectrum(sys.argv[2], sys.argv[3]) </code></pre> <p>I have deleted check and spectrum functions because it's simple calculations e.g. length comparison, so it is not what cause a problem as I think.</p> <p>I will be very thankfully for help!</p> <p>Regards, Mateusz</p>
    singulars
    1. This table or related slice is empty.
    1. This table or related slice is empty.
    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.
 

Querying!

 
Guidance

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