beginning script to make alignmentsof nterm and cterm of tyrosinase seqs.
can now pull GID.
can get gid, start position, and end position
working on multople domain issue figure out the start part
Now output beginning of domains, and end of domains, and how many tyrosinase domains present for each sequence.
commit after merge
now extract full seqs from blastcmd ready to get nterm and cterm pieces.
get nterm and cterm and full seqs in dict for md and sd sequences. Now can begin making trees.
output four files with seqs to be aligned from cterm and nterm md and sd.
added script to run clustalo on nterm_md file.
issue with nterm being too short will have to filter out things with a short nterm.
input seqs and get their size distribution
script can get each seq and its length.
lengths now in dict
just commiting things
at the stage where I will now extract seqs from the original seqs that have the length I am looking for.
at the stage where I will now extract seqs from the original seqs that have the length I am looking for.
opening graph does not stop flow of the script.
opening graph does not stop flow of the script.
Have sequence and GID for seqs above my percentile threshold.