31
31
< link href ="_static/styles/theme.css?digest=8878045cc6db502f8baf " rel ="stylesheet " />
32
32
< link href ="_static/styles/pydata-sphinx-theme.css?digest=8878045cc6db502f8baf " rel ="stylesheet " />
33
33
34
- < link rel ="stylesheet " type ="text/css " href ="_static/pygments.css?v=a746c00c " />
34
+ < link rel ="stylesheet " type ="text/css " href ="_static/pygments.css?v=8f2a1f02 " />
35
35
< link rel ="stylesheet " type ="text/css " href ="_static/copybutton.css?v=76b2166b " />
36
36
< link rel ="stylesheet " type ="text/css " href ="_static/plot_directive.css " />
37
37
< link rel ="stylesheet " type ="text/css " href ="_static/sphinx-design.min.css?v=95c83b7e " />
@@ -754,7 +754,7 @@ <h2>Tutorial<a class="headerlink" href="#tutorial" title="Link to this heading">
754
754
< h3 > Alignment data structure< a class ="headerlink " href ="#alignment-data-structure " title ="Link to this heading "> #</ a > </ h3 >
755
755
< p > Load two DNA sequences that have been previously aligned into a < code class ="docutils literal notranslate "> < span class ="pre "> TabularMSA</ span > </ code >
756
756
object, using sequence IDs as the MSA’s index:</ p >
757
- < div class ="doctest highlight-default notranslate "> < div class ="highlight "> < pre > < span > </ span > < span class ="gp "> >>> </ span > < span class ="kn "> from</ span > < span class ="nn "> skbio</ span > < span class ="kn "> import</ span > < span class ="n "> TabularMSA</ span > < span class ="p "> ,</ span > < span class ="n "> DNA</ span >
757
+ < div class ="doctest highlight-default notranslate "> < div class ="highlight "> < pre > < span > </ span > < span class ="gp "> >>> </ span > < span class ="kn "> from</ span > < span class =" w " > </ span > < span class ="nn "> skbio</ span > < span class =" w " > </ span > < span class ="kn "> import</ span > < span class ="n "> TabularMSA</ span > < span class ="p "> ,</ span > < span class ="n "> DNA</ span >
758
758
< span class ="gp "> >>> </ span > < span class ="n "> seqs</ span > < span class ="o "> =</ span > < span class ="p "> [</ span > < span class ="n "> DNA</ span > < span class ="p "> (</ span > < span class ="s2 "> "ACC--G-GGTA.."</ span > < span class ="p "> ,</ span > < span class ="n "> metadata</ span > < span class ="o "> =</ span > < span class ="p "> {</ span > < span class ="s1 "> 'id'</ span > < span class ="p "> :</ span > < span class ="s2 "> "seq1"</ span > < span class ="p "> }),</ span >
759
759
< span class ="gp "> ... </ span > < span class ="n "> DNA</ span > < span class ="p "> (</ span > < span class ="s2 "> "TCC--G-GGCA.."</ span > < span class ="p "> ,</ span > < span class ="n "> metadata</ span > < span class ="o "> =</ span > < span class ="p "> {</ span > < span class ="s1 "> 'id'</ span > < span class ="p "> :</ span > < span class ="s2 "> "seq2"</ span > < span class ="p "> })]</ span >
760
760
< span class ="gp "> >>> </ span > < span class ="n "> msa</ span > < span class ="o "> =</ span > < span class ="n "> TabularMSA</ span > < span class ="p "> (</ span > < span class ="n "> seqs</ span > < span class ="p "> ,</ span > < span class ="n "> minter</ span > < span class ="o "> =</ span > < span class ="s1 "> 'id'</ span > < span class ="p "> )</ span >
@@ -775,7 +775,7 @@ <h3>Alignment data structure<a class="headerlink" href="#alignment-data-structur
775
775
< section id ="using-the-optimized-alignment-algorithm ">
776
776
< h3 > Using the optimized alignment algorithm< a class ="headerlink " href ="#using-the-optimized-alignment-algorithm " title ="Link to this heading "> #</ a > </ h3 >
777
777
< p > Using the convenient < code class ="docutils literal notranslate "> < span class ="pre "> local_pairwise_align_ssw</ span > </ code > function:</ p >
778
- < div class ="doctest highlight-default notranslate "> < div class ="highlight "> < pre > < span > </ span > < span class ="gp "> >>> </ span > < span class ="kn "> from</ span > < span class ="nn "> skbio.alignment</ span > < span class ="kn "> import</ span > < span class ="n "> local_pairwise_align_ssw</ span >
778
+ < div class ="doctest highlight-default notranslate "> < div class ="highlight "> < pre > < span > </ span > < span class ="gp "> >>> </ span > < span class ="kn "> from</ span > < span class =" w " > </ span > < span class ="nn "> skbio.alignment</ span > < span class =" w " > </ span > < span class ="kn "> import</ span > < span class ="n "> local_pairwise_align_ssw</ span >
779
779
< span class ="gp "> >>> </ span > < span class ="n "> alignment</ span > < span class ="p "> ,</ span > < span class ="n "> score</ span > < span class ="p "> ,</ span > < span class ="n "> start_end_positions</ span > < span class ="o "> =</ span > < span class ="n "> local_pairwise_align_ssw</ span > < span class ="p "> (</ span >
780
780
< span class ="gp "> ... </ span > < span class ="n "> DNA</ span > < span class ="p "> (</ span > < span class ="s2 "> "ACTAAGGCTCTCTACCCCTCTCAGAGA"</ span > < span class ="p "> ),</ span >
781
781
< span class ="gp "> ... </ span > < span class ="n "> DNA</ span > < span class ="p "> (</ span > < span class ="s2 "> "ACTAAGGCTCCTAACCCCCTTTTCTCAGA"</ span > < span class ="p "> )</ span >
@@ -796,7 +796,7 @@ <h3>Using the optimized alignment algorithm<a class="headerlink" href="#using-th
796
796
</ pre > </ div >
797
797
</ div >
798
798
< p > Using the < code class ="docutils literal notranslate "> < span class ="pre "> StripedSmithWaterman</ span > </ code > object:</ p >
799
- < div class ="doctest highlight-default notranslate "> < div class ="highlight "> < pre > < span > </ span > < span class ="gp "> >>> </ span > < span class ="kn "> from</ span > < span class ="nn "> skbio.alignment</ span > < span class ="kn "> import</ span > < span class ="n "> StripedSmithWaterman</ span >
799
+ < div class ="doctest highlight-default notranslate "> < div class ="highlight "> < pre > < span > </ span > < span class ="gp "> >>> </ span > < span class ="kn "> from</ span > < span class =" w " > </ span > < span class ="nn "> skbio.alignment</ span > < span class =" w " > </ span > < span class ="kn "> import</ span > < span class ="n "> StripedSmithWaterman</ span >
800
800
< span class ="gp "> >>> </ span > < span class ="n "> query</ span > < span class ="o "> =</ span > < span class ="n "> StripedSmithWaterman</ span > < span class ="p "> (</ span > < span class ="s2 "> "ACTAAGGCTCTCTACCCCTCTCAGAGA"</ span > < span class ="p "> )</ span >
801
801
< span class ="gp "> >>> </ span > < span class ="n "> alignment</ span > < span class ="o "> =</ span > < span class ="n "> query</ span > < span class ="p "> (</ span > < span class ="s2 "> "AAAAAACTCTCTAAACTCACTAAGGCTCTCTACCCCTCTTCAGAGAAGTCGA"</ span > < span class ="p "> )</ span >
802
802
< span class ="gp "> >>> </ span > < span class ="nb "> print</ span > < span class ="p "> (</ span > < span class ="n "> alignment</ span > < span class ="p "> )</ span >
@@ -808,7 +808,7 @@ <h3>Using the optimized alignment algorithm<a class="headerlink" href="#using-th
808
808
</ div >
809
809
< p > Using the < code class ="docutils literal notranslate "> < span class ="pre "> StripedSmithWaterman</ span > </ code > object for multiple targets in an efficient
810
810
way and finding the aligned sequence representations:</ p >
811
- < div class ="doctest highlight-default notranslate "> < div class ="highlight "> < pre > < span > </ span > < span class ="gp "> >>> </ span > < span class ="kn "> from</ span > < span class ="nn "> skbio.alignment</ span > < span class ="kn "> import</ span > < span class ="n "> StripedSmithWaterman</ span >
811
+ < div class ="doctest highlight-default notranslate "> < div class ="highlight "> < pre > < span > </ span > < span class ="gp "> >>> </ span > < span class ="kn "> from</ span > < span class =" w " > </ span > < span class ="nn "> skbio.alignment</ span > < span class =" w " > </ span > < span class ="kn "> import</ span > < span class ="n "> StripedSmithWaterman</ span >
812
812
< span class ="gp "> >>> </ span > < span class ="n "> alignments</ span > < span class ="o "> =</ span > < span class ="p "> []</ span >
813
813
< span class ="gp "> >>> </ span > < span class ="n "> target_sequences</ span > < span class ="o "> =</ span > < span class ="p "> [</ span >
814
814
< span class ="gp "> ... </ span > < span class ="s2 "> "GCTAACTAGGCTCCCTTCTACCCCTCTCAGAGA"</ span > < span class ="p "> ,</ span >
@@ -848,8 +848,8 @@ <h3>Using the slow alignment algorithm<a class="headerlink" href="#using-the-slo
848
848
< p > Here we locally align a pair of protein sequences using gap open penalty
849
849
of 11 and a gap extend penalty of 1 (in other words, it is much more
850
850
costly to open a new gap than extend an existing one).</ p >
851
- < div class ="doctest highlight-default notranslate "> < div class ="highlight "> < pre > < span > </ span > < span class ="gp "> >>> </ span > < span class ="kn "> from</ span > < span class ="nn "> skbio</ span > < span class ="kn "> import</ span > < span class ="n "> Protein</ span >
852
- < span class ="gp "> >>> </ span > < span class ="kn "> from</ span > < span class ="nn "> skbio.alignment</ span > < span class ="kn "> import</ span > < span class ="n "> local_pairwise_align_protein</ span >
851
+ < div class ="doctest highlight-default notranslate "> < div class ="highlight "> < pre > < span > </ span > < span class ="gp "> >>> </ span > < span class ="kn "> from</ span > < span class =" w " > </ span > < span class ="nn "> skbio</ span > < span class =" w " > </ span > < span class ="kn "> import</ span > < span class ="n "> Protein</ span >
852
+ < span class ="gp "> >>> </ span > < span class ="kn "> from</ span > < span class =" w " > </ span > < span class ="nn "> skbio.alignment</ span > < span class =" w " > </ span > < span class ="kn "> import</ span > < span class ="n "> local_pairwise_align_protein</ span >
853
853
< span class ="gp "> >>> </ span > < span class ="n "> s1</ span > < span class ="o "> =</ span > < span class ="n "> Protein</ span > < span class ="p "> (</ span > < span class ="s2 "> "HEAGAWGHEE"</ span > < span class ="p "> )</ span >
854
854
< span class ="gp "> >>> </ span > < span class ="n "> s2</ span > < span class ="o "> =</ span > < span class ="n "> Protein</ span > < span class ="p "> (</ span > < span class ="s2 "> "PAWHEAE"</ span > < span class ="p "> )</ span >
855
855
< span class ="gp "> >>> </ span > < span class ="n "> alignment</ span > < span class ="p "> ,</ span > < span class ="n "> score</ span > < span class ="p "> ,</ span > < span class ="n "> start_end_positions</ span > < span class ="o "> =</ span > < span class ="n "> local_pairwise_align_protein</ span > < span class ="p "> (</ span >
@@ -874,8 +874,8 @@ <h3>Using the slow alignment algorithm<a class="headerlink" href="#using-the-slo
874
874
</ pre > </ div >
875
875
</ div >
876
876
< p > Similarly, we can perform global alignment of nucleotide sequences:</ p >
877
- < div class ="doctest highlight-default notranslate "> < div class ="highlight "> < pre > < span > </ span > < span class ="gp "> >>> </ span > < span class ="kn "> from</ span > < span class ="nn "> skbio</ span > < span class ="kn "> import</ span > < span class ="n "> DNA</ span >
878
- < span class ="gp "> >>> </ span > < span class ="kn "> from</ span > < span class ="nn "> skbio.alignment</ span > < span class ="kn "> import</ span > < span class ="n "> global_pairwise_align_nucleotide</ span >
877
+ < div class ="doctest highlight-default notranslate "> < div class ="highlight "> < pre > < span > </ span > < span class ="gp "> >>> </ span > < span class ="kn "> from</ span > < span class =" w " > </ span > < span class ="nn "> skbio</ span > < span class =" w " > </ span > < span class ="kn "> import</ span > < span class ="n "> DNA</ span >
878
+ < span class ="gp "> >>> </ span > < span class ="kn "> from</ span > < span class =" w " > </ span > < span class ="nn "> skbio.alignment</ span > < span class =" w " > </ span > < span class ="kn "> import</ span > < span class ="n "> global_pairwise_align_nucleotide</ span >
879
879
< span class ="gp "> >>> </ span > < span class ="n "> s1</ span > < span class ="o "> =</ span > < span class ="n "> DNA</ span > < span class ="p "> (</ span > < span class ="s2 "> "GCGTGCCTAAGGTATGCAAG"</ span > < span class ="p "> )</ span >
880
880
< span class ="gp "> >>> </ span > < span class ="n "> s2</ span > < span class ="o "> =</ span > < span class ="n "> DNA</ span > < span class ="p "> (</ span > < span class ="s2 "> "ACGTGCCTAGGTACGCAAG"</ span > < span class ="p "> )</ span >
881
881
< span class ="gp "> >>> </ span > < span class ="n "> alignment</ span > < span class ="p "> ,</ span > < span class ="n "> score</ span > < span class ="p "> ,</ span > < span class ="n "> start_end_positions</ span > < span class ="o "> =</ span > < span class ="n "> global_pairwise_align_nucleotide</ span > < span class ="p "> (</ span >
@@ -981,7 +981,7 @@ <h3>Using the slow alignment algorithm<a class="headerlink" href="#using-the-slo
981
981
982
982
< p class ="copyright ">
983
983
984
- © Copyright 2014-2024 , scikit-bio development team.
984
+ © Copyright 2014-2025 , scikit-bio development team.
985
985
< br />
986
986
987
987
</ p >
0 commit comments