Jump to content

"GC" A script to rapidly manipulate DNA sequences


Recommended Posts

Hi all,

This is a short and again messy script (as far as the source code is concerned), but its been very very usefull to me over the past 3-4 months so I thought id post it up.

simply if you have some DNA sequences with other characters or spaces or /n in it, It will strip it down to the ATGC.

you simply highlight and select whatever you want to manipluate, copy it, then press the past from clipboard button to put it into the GUI. (to copy anything out of the GUI simply press the copy button)

once in the edit you can do usefull things such as reverse complement, switch all to upper or lower case, and analyse the GC content and possible methylation sites.

Included in download is the source and the icon, so it can be compiled whenever.

Any suggestions as to what else would be usefull to add would be greatly appreciated.


Sebastian Eves-van den Akker.

GC autoit.zip

GC - Program to rapidly manipulate DNA SequencesRotaMol - Program to measure Protein Size
Link to comment
Share on other sites

  • 2 months later...

You can simplify and speed up things a lot.

For instance:

if $msg = $getclip Then
  $promoter = clipget()
  $promoterr = StringStripWS($promoter,8)
  $sStripped_1 = StringRegExpReplace($promoterr,"[1,2,3,4,5,6,7,8,9,0,-,=,+,_,q,w,e,r,y,u,i,o,p,},{,s,d,f,h,j,k,l,;,:,',@,#,~,,|,z,x,v,b,n,m,<,>,.,?,/,Q,W,E,R,Y,U,I,O,P,S,D,F,H,J,K,L,Z,X,V,B,N,M,`,!]","")
$sStripped_2 = StringReplace($sStripped_1, "[", "")
$sStripped_3 = StringReplace($sStripped_2, "]", "")
$sStripped_4 = StringReplace($sStripped_3, "$", "")
$sStripped_5 = StringReplace($sStripped_4, "â, "")
$sStripped_6 = StringReplace($sStripped_5, "%", "")
$sStripped_7 = StringReplace($sStripped_6, "*", "")
$sStripped_8 = StringReplace($sStripped_7, "&", "")
$sStripped_9 = StringReplace($sStripped_8, "", "")
$sStripped_10 = StringReplace($sStripped_9, "-", "")
$sStripped_11 = StringReplace($sStripped_10, ")", "")
$sStripped_12= StringReplace($sStripped_11, "(", "")
$sStripped_13= StringReplace($sStripped_12, "]", "")
$sStripped_14= StringReplace($sStripped_13, "Ң, "")
;$sStripped_15= StringReplace($sStripped_14, "C", "T") ;for methylation if you want
$string7 = StringUpper($sStripped_14)
make that:

If $msg = $getclip Then
  $promoter = StringUpper(ClipGet())
  $promoter = StringRegExpReplace($promoter, "(?m)[^AGCT]", "")
  ;$promoter_methyl = StringReplace($promoter, "C", "T") ;for methylation if you want
  GUICtrlSetData($string1, $promoter)

This wonderful site allows debugging and testing regular expressions (many flavors available). An absolute must have in your bookmarks.
Another excellent RegExp tutorial. Don't forget downloading your copy of up-to-date pcretest.exe and pcregrep.exe here
RegExp tutorial: enough to get started
PCRE v8.33 regexp documentation latest available release and currently implemented in AutoIt beta.

SQLitespeed is another feature-rich premier SQLite manager (includes import/export). Well worth a try.
SQLite Expert (freeware Personal Edition or payware Pro version) is a very useful SQLite database manager.
An excellent eBook covering almost every aspect of SQLite3: a must-read for anyone doing serious work.
SQL tutorial (covers "generic" SQL, but most of it applies to SQLite as well)
A work-in-progress SQLite3 tutorial. Don't miss other LxyzTHW pages!
SQLite official website with full documentation (may be newer than the SQLite library that comes standard with AutoIt)

Link to comment
Share on other sites

Create an account or sign in to comment

You need to be a member in order to leave a comment

Create an account

Sign up for a new account in our community. It's easy!

Register a new account

Sign in

Already have an account? Sign in here.

Sign In Now

  • Create New...