Practical Bioinformatics Mark Voorhies 6/3/2013 Mark Voorhies - PowerPoint PPT Presentation

Feb 11, 2023 •459 likes •569 views

Practical Bioinformatics Mark Voorhies 6/3/2013 Mark Voorhies Practical Bioinformatics Needleman-Wunsch with g=0 A G C G G T A G A G C G G A Mark Voorhies Practical Bioinformatics Needleman-Wunsch with g=0 A G C G G T A G A

Practical Bioinformatics Mark Voorhies 6/3/2013 Mark Voorhies Practical Bioinformatics
Needleman-Wunsch with g=0 A G C G G T A G A G C G G A Mark Voorhies Practical Bioinformatics
Needleman-Wunsch with g=0 A G C G G T A G A def n w f i l l ( seq1 , seq2 , s , e ) : G # m[ i ] [ j ] = best s c o r e f o r subalignment # of seq1 [ : i ] , seq2 [ : j ] C m = [ [ ] ] G G A Mark Voorhies Practical Bioinformatics
Needleman-Wunsch with g=0 A G C G G T A 0 -1 -2 -3 -4 -5 -6 -7 G -1 # Top l e f t c o r n e r A -2 m[ − 1]. append (0) # F i l l f i r s t row as l e a d i n g gaps G -3 f o r j i n range ( l e n ( seq2 ) ) : m[ − 1]. append (m[ 0 ] [ j ]+e ) C -4 # F i l l f i r s t column as l e a d i n g gaps i range ( l e n ( seq1 ) ) : f o r i n G -5 m. append ( [m[ i ] [ 0 ] + e ] ) G -6 A -7 Mark Voorhies Practical Bioinformatics
Needleman-Wunsch with g=0 A G C G G T A 0 -1 -2 -3 -4 -5 -6 -7 -1 -2 -3 -4 -5 G -1 f o r i i n range ( l e n ( seq1 ) ) : -1 0 f o r j i n range ( l e n ( seq2 ) ) : A -2 0 -1 -1 -2 -3 -4 -3 # Score f o r a l i g n i n g seq1 [ i ] with seq2 [ j ] match = m[ i ] [ j ]+ s [ seq1 [ i ] ] [ seq2 [ j ] ] -1 -2 -3 # Score f o r a l i g n i n g seq1 [ i ] with a gap G -3 -1 1 0 0 hgap = m[ i +1][ j ]+e # Score f o r a l i g n i n g seq2 [ i ] with a gap 2 1 0 -1 -2 C -4 -2 0 vgap = m[ i ] [ j +1]+e G -5 -3 -1 1 3 2 1 0 best = max( match , vgap , hgap ) m[ − 1]. append ( best ) -4 -2 2 4 3 2 G -6 0 A -7 -5 -3 -1 1 3 3 4 Mark Voorhies Practical Bioinformatics
Needleman-Wunsch with g=0 # S t a r t at bottom r i g h t c o r n e r curpos = ( l e n ( seq1 ) , l e n ( seq2 ) ) a l i g n e d 1 = ”” a l i g n e d 2 = ”” A G C G G T A e x i t F l a g = F a l s e 0 -1 -2 -3 -4 -5 -6 -7 f o r i i n range ( l e n ( seq1)+ l e n ( seq2 ) ) : p l i s t = p [ curpos [ 0 ] ] [ curpos [ 1 ] ] -1 -2 -3 -4 -5 G -1 -1 0 i f ( p l i s t i s None ) : e x i t F l a g = True -1 -1 -2 -3 -4 -3 A -2 0 break nextpos = p l i s t [ 0 ] -1 -2 -3 # Check f o r vgap G -3 -1 1 0 0 i f ( nextpos [ 0 ] == curpos [ 0 ] ) : a l i g n e d 1 = ” − ”+a l i g n e d 1 C 2 1 0 -1 -2 -4 -2 0 e l s e : a l i g n e d 1 = seq1 [ nextpos [ 0 ] ] + a l i g n e d 1 -3 -1 1 3 2 1 0 G -5 # Check f o r hgap i f ( nextpos [ 1 ] == curpos [ 1 ] ) : G -4 -2 2 4 3 2 -6 0 a l i g n e d 2 = ” − ”+a l i g n e d 2 e l s e : A -5 -3 -1 1 3 3 4 -7 a l i g n e d 2 = seq2 [ nextpos [ 1 ] ] + a l i g n e d 2 curpos = nextpos i f ( e x i t F l a g == F a l s e ) : ”WARNING: Unexpected e x i t from traceback ” p r i n t Mark Voorhies Practical Bioinformatics
Change of Coordinates http://xkcd.com/123/ Mark Voorhies Practical Bioinformatics
Phase-specific gene expression in Histoplasma capsulatum Diane Inglis Mark Voorhies Practical Bioinformatics
Phase-specific gene expression in
Phase-specific gene expression in Histoplasma capsulatum Diane Inglis Mark Voorhies Practical Bioinformatics

Recommend

Practical Bioinformatics Mark Voorhies 5/ 24/ 2013 Mark Voorhies Practical Bioinformatics

Practical Bioinformatics Mark Voorhies 5/ 24/ 2013 Mark Voorhies Practical Bioinformatics Clustering exercises { Visualizing the distance matrix Mark Voorhies Practical Bioinformatics Scripting Cluster Running Cluster3 from the command line

401 views • 7 slides

Practical Bioinformatics Mark Voorhies 5/22/2015 Mark Voorhies Practical Bioinformatics PAM

Practical Bioinformatics Mark Voorhies 5/22/2015 Mark Voorhies Practical Bioinformatics PAM (Dayhoff) and BLOSUM matrices PAM1 matrix originally calculated from manual alignments of highly conserved sequences (myoglobin, cytochrome C, etc.)

979 views • 71 slides

Practical Bioinformatics Mark Voorhies 5/15/2015 Mark Voorhies Practical Bioinformatics

Practical Bioinformatics Mark Voorhies 5/15/2015 Mark Voorhies Practical Bioinformatics Gotchas Indentation matters Mark Voorhies Practical Bioinformatics Clustering exercises Visualizing the distance matrix Mark Voorhies Practical

705 views • 46 slides

Practical Bioinformatics Mark Voorhies 4/16/2018 Mark Voorhies Practical Bioinformatics

Practical Bioinformatics Mark Voorhies 4/16/2018 Mark Voorhies Practical Bioinformatics JavaTreeView link-out for ENSEMBL Mouse http://www.ensembl.org/Mus musculus/Gene/Summary?g=HEADER Mark Voorhies Practical Bioinformatics Science! Mark

767 views • 20 slides

Practical Bioinformatics Mark Voorhies 5/23/2019 Mark Voorhies Practical Bioinformatics

Practical Bioinformatics Mark Voorhies 5/23/2019 Mark Voorhies Practical Bioinformatics JavaTreeView link-out for ENSEMBL Mouse http://www.ensembl.org/Mus musculus/Gene/Summary?g=HEADER Mark Voorhies Practical Bioinformatics Science! Mark

330 views • 13 slides

Practical Bioinformatics Mark Voorhies 5/21/2019 Mark Voorhies Practical Bioinformatics Change

Practical Bioinformatics Mark Voorhies 5/21/2019 Mark Voorhies Practical Bioinformatics Change of Coordinates http://xkcd.com/123/ Mark Voorhies Practical Bioinformatics Principal Components Analysis (PCA) Is there a point of view that

264 views • 11 slides

Practical Bioinformatics Mark Voorhies 5/26/2015 Mark Voorhies Practical Bioinformatics Habits

Practical Bioinformatics Mark Voorhies 5/26/2015 Mark Voorhies Practical Bioinformatics Habits are things you get for free, without requiring any special work. Cory Doctorow Advice to Writers, 4/5/2012 Mark Voorhies Practical

521 views • 27 slides

Practical Bioinformatics Mark Voorhies 5/11/2015 Mark Voorhies Practical Bioinformatics

Introduction to Python Practical Bioinformatics Mark Voorhies 5/11/2015 Mark Voorhies Practical Bioinformatics Introduction to Python Good ideas shamelessly lifted from David Erle and Josh Pollack Mark Voorhies Practical Bioinformatics

718 views • 46 slides

Practical Bioinformatics Mark Voorhies 4/9/2018 Mark Voorhies Practical Bioinformatics

Practical Bioinformatics Mark Voorhies 4/9/2018 Mark Voorhies Practical Bioinformatics Clustering exercises Visualizing the distance matrix Mark Voorhies Practical Bioinformatics Dictionaries d i c t i o n a r y = { A : T ,

77 views • 5 slides

Practical Bioinformatics Mark Voorhies 5/14/2019 Mark Voorhies Practical Bioinformatics Course

Practical Bioinformatics Mark Voorhies 5/14/2019 Mark Voorhies Practical Bioinformatics Course platform: VirtualBox Host operating system (e.g., OS X) VirtualBox Debian Linux Web 8888 8088 Jupyter Browser launches Python3 Bash Bash

756 views • 32 slides

Practical Bioinformatics Mark Voorhies 4/19/2011 Mark Voorhies Practical Bioinformatics Review

Practical Bioinformatics Mark Voorhies 4/19/2011 Mark Voorhies Practical Bioinformatics Review Topic Learning Python Dive Into Python Types (string, int, float, list) 2 3 Assignment (=, +=) 2 3 Slicing and indexing 4 3.2 References

443 views • 6 slides

Practical Bioinformatics Mark Voorhies 4/29/2011 Mark Voorhies Practical Bioinformatics Our

Practical Bioinformatics Mark Voorhies 4/29/2011 Mark Voorhies Practical Bioinformatics Our current tool set data (strings, floats, lists, nested lists) logic (if/then, try/except, for, while) functions (def, import, reload) Mark Voorhies

615 views • 25 slides

Practical Bioinformatics Mark Voorhies 5/29/2019 Mark Voorhies Practical Bioinformatics

Practical Bioinformatics Mark Voorhies 5/29/2019 Mark Voorhies Practical Bioinformatics Example Pipeline: Overview Mark Voorhies Practical Bioinformatics Example Pipeline: Overview Generate Genome Transfer Pre-process Samples Coverage

521 views • 16 slides

Practical Bioinformatics Mark Voorhies 5/19/2015 Mark Voorhies Practical Bioinformatics Review

Practical Bioinformatics Mark Voorhies 5/19/2015 Mark Voorhies Practical Bioinformatics Review Documentation A program should communicate its intent to both the computer and human programmers. Comments Docstrings Mark Voorhies

495 views • 26 slides

Practical Bioinformatics Mark Voorhies 5/31/2013 Mark Voorhies Practical Bioinformatics

Practical Bioinformatics Mark Voorhies 5/31/2013 Mark Voorhies Practical Bioinformatics Exercise: Scoring a gapped alignment 1 Given two equal length gapped sequences (where - represents a gap) and a scoring matrix, calculate an

335 views • 5 slides

Practical Bioinformatics Mark Voorhies 4/2/2018 Mark Voorhies Practical Bioinformatics

Introduction to Python Practical Bioinformatics Mark Voorhies 4/2/2018 Mark Voorhies Practical Bioinformatics Introduction to Python Resources Course website: http://histo.ucsf.edu/BMS270/ Resources on the course website: Syllabus Papers

660 views • 37 slides

Practical Bioinformatics Mark Voorhies 4/24/2017 Mark Voorhies Practical Bioinformatics

Introduction to Python Practical Bioinformatics Mark Voorhies 4/24/2017 Mark Voorhies Practical Bioinformatics Introduction to Python Resources Course website: http://histo.ucsf.edu/BMS270/ Resources on the course website: Syllabus Papers

701 views • 43 slides

Practical Bioinformatics Mark Voorhies 5/13/2019 Mark Voorhies Practical Bioinformatics

Introduction to Python Practical Bioinformatics Mark Voorhies 5/13/2019 Mark Voorhies Practical Bioinformatics Introduction to Python Resources Course website: http://histo.ucsf.edu/BMS270/ Resources on the course website: Syllabus Papers

594 views • 42 slides

Practical Bioinformatics Mark Voorhies 4/6/2017 Mark Voorhies Practical Bioinformatics Loading

Practical Bioinformatics Mark Voorhies 4/6/2017 Mark Voorhies Practical Bioinformatics Loading and re-loading your functions # Use import the f i r s t time you load a module # (And keep using import u n t i l i t loads # s u c c

775 views • 32 slides

Practical Bioinformatics Mark Voorhies 5/12/2015 Mark Voorhies Practical Bioinformatics

Practical Bioinformatics Mark Voorhies 5/12/2015 Mark Voorhies Practical Bioinformatics Gotchas Strings are quoted, names of things are not. mystring = mystring Mark Voorhies Practical Bioinformatics Gotchas Strings are quoted,

554 views • 27 slides

Practical Bioinformatics Mark Voorhies 5/21/2013 Mark Voorhies Practical Bioinformatics

Practical Bioinformatics Mark Voorhies 5/21/2013 Mark Voorhies Practical Bioinformatics Gotchas Enclosing symbols (,,(),[],...) need to be matched Mark Voorhies Practical Bioinformatics Gotchas Enclosing symbols

99 views • 7 slides

Practical Bioinformatics Mark Voorhies 4/20/2011 Mark Voorhies Practical Bioinformatics Review

Practical Bioinformatics Mark Voorhies 4/20/2011 Mark Voorhies Practical Bioinformatics Review Corrected slides for days 1 and 2 are on the website Mark Voorhies Practical Bioinformatics Review Corrected slides for days 1 and 2 are on the

248 views • 10 slides

Practical Bioinformatics Mark Voorhies 4/28/2017 Mark Voorhies Practical Bioinformatics

Practical Bioinformatics Mark Voorhies 4/28/2017 Mark Voorhies Practical Bioinformatics Pearson distances Pearson similarity N i ( x i x offset )( y i y offset ) s ( x , y ) = N N i ( x i x offset ) 2 i ( y i

502 views • 28 slides

Intro to NGS Theoretical and Practical HiC Workshop: Wet-lab and Bioinformatics November 4th,

Intro to NGS Theoretical and Practical HiC Workshop: Wet-lab and Bioinformatics November 4th, 2019 Selene L. Fernndez-Valverde regRNAlab.github.io @SelFdz 1 Learning objectives In this class we will learn How high-throughput (NGS)

728 views • 30 slides