bio-samtools 0.4.2 → 0.5.0
Sign up to get free protection for your applications and to get access to all the features.
- data/Gemfile.lock +1 -0
- data/README.rdoc +1 -0
- data/VERSION +1 -1
- data/bio-samtools.gemspec +19 -19
- data/doc/Bio.html +71 -4
- data/doc/Bio/DB.html +6 -4
- data/doc/Bio/DB/Alignment.html +8 -8
- data/doc/Bio/DB/Pileup.html +515 -0
- data/doc/Bio/DB/SAM.html +37 -37
- data/doc/Bio/DB/SAM/Library.html +4 -4
- data/doc/Bio/DB/SAM/Tools.html +4 -4
- data/doc/Bio/DB/SAM/Tools/Bam1CoreT.html +4 -4
- data/doc/Bio/DB/SAM/Tools/Bam1T.html +4 -4
- data/doc/Bio/DB/SAM/Tools/BamHeaderT.html +4 -4
- data/doc/Bio/DB/SAM/Tools/BamPileup1T.html +4 -4
- data/doc/Bio/DB/SAM/Tools/SamfileT.html +4 -4
- data/doc/Bio/DB/SAM/Tools/SamfileTX.html +4 -4
- data/doc/Bio/DB/SAMException.html +5 -5
- data/doc/Bio/DB/Tag.html +5 -5
- data/doc/Bio/DB/Vcf.html +467 -0
- data/doc/LICENSE_txt.html +57 -60
- data/doc/LibC.html +6 -4
- data/doc/created.rid +7 -7
- data/doc/index.html +4 -4
- data/doc/js/search_index.js +1 -1
- data/doc/lib/bio-samtools_rb.html +4 -4
- data/doc/lib/bio/db/pileup_rb.html +171 -0
- data/doc/lib/bio/db/sam/bam_rb.html +4 -4
- data/doc/lib/bio/db/sam/faidx_rb.html +4 -4
- data/doc/lib/bio/db/sam/library_rb.html +4 -4
- data/doc/lib/bio/db/sam/sam_rb.html +4 -4
- data/doc/lib/bio/db/sam_rb.html +4 -4
- data/doc/lib/bio/db/vcf_rb.html +124 -0
- data/doc/table_of_contents.html +16 -26
- data/doc/tutorial.html +26 -1
- data/doc/tutorial.pdf +0 -0
- data/lib/bio/db/{sam/pileup.rb → pileup.rb} +6 -4
- data/lib/bio/db/sam.rb +13 -12
- data/lib/bio/db/{sam/vcf.rb → vcf.rb} +6 -2
- data/test/samples/small/sorted.bam.bai +0 -0
- data/test/test_basic.rb +7 -5
- data/test/test_vcf.rb +5 -5
- metadata +29 -26
- data/test/basictest.rb +0 -322
- data/test/pileup.rb +0 -68
data/doc/LICENSE_txt.html
CHANGED
@@ -65,6 +65,8 @@
|
|
65
65
|
|
66
66
|
<li><a href="./Bio/DB/Alignment.html">Bio::DB::Alignment</a>
|
67
67
|
|
68
|
+
<li><a href="./Bio/DB/Pileup.html">Bio::DB::Pileup</a>
|
69
|
+
|
68
70
|
<li><a href="./Bio/DB/SAM.html">Bio::DB::SAM</a>
|
69
71
|
|
70
72
|
<li><a href="./Bio/DB/SAM/Library.html">Bio::DB::SAM::Library</a>
|
@@ -89,11 +91,9 @@
|
|
89
91
|
|
90
92
|
<li><a href="./Bio/DB/Tag.html">Bio::DB::Tag</a>
|
91
93
|
|
92
|
-
<li><a href="./
|
93
|
-
|
94
|
-
<li><a href="./Pileup.html">Pileup</a>
|
94
|
+
<li><a href="./Bio/DB/Vcf.html">Bio::DB::Vcf</a>
|
95
95
|
|
96
|
-
<li><a href="./
|
96
|
+
<li><a href="./LibC.html">LibC</a>
|
97
97
|
|
98
98
|
</ul>
|
99
99
|
</nav>
|
@@ -106,18 +106,17 @@
|
|
106
106
|
<p>Copyright © 2011 Raoul J.P. Bonnal</p>
|
107
107
|
|
108
108
|
<p>Permission is hereby granted, free of charge, to any person obtaining a
|
109
|
-
copy of this software and associated documentation files (the
|
110
|
-
|
111
|
-
|
112
|
-
|
113
|
-
|
114
|
-
following conditions:</p>
|
109
|
+
copy of this software and associated documentation files (the “Software”),
|
110
|
+
to deal in the Software without restriction, including without limitation
|
111
|
+
the rights to use, copy, modify, merge, publish, distribute, sublicense,
|
112
|
+
and/or sell copies of the Software, and to permit persons to whom the
|
113
|
+
Software is furnished to do so, subject to the following conditions:</p>
|
115
114
|
|
116
115
|
<p>The above copyright notice and this permission notice shall be included in
|
117
116
|
all copies or substantial portions of the Software.</p>
|
118
117
|
|
119
|
-
<p>THE SOFTWARE IS PROVIDED “AS IS”, WITHOUT WARRANTY OF ANY KIND, EXPRESS
|
120
|
-
|
118
|
+
<p>THE SOFTWARE IS PROVIDED “AS IS”, WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
|
119
|
+
IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
|
121
120
|
FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
|
122
121
|
AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
|
123
122
|
LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
|
@@ -143,11 +142,10 @@ of this license document, but changing it is not allowed.
|
|
143
142
|
|
144
143
|
<p>to take away your freedom to share and change the works. By contrast, the
|
145
144
|
GNU General Public License is intended to guarantee your freedom to share
|
146
|
-
and change all versions of a program–to make sure it remains free
|
147
|
-
|
148
|
-
|
149
|
-
|
150
|
-
programs, too.</p>
|
145
|
+
and change all versions of a program–to make sure it remains free software
|
146
|
+
for all its users. We, the Free Software Foundation, use the GNU General
|
147
|
+
Public License for most of our software; it applies also to any other work
|
148
|
+
released this way by its authors. You can apply it to your programs, too.</p>
|
151
149
|
|
152
150
|
<pre>When we speak of free software, we are referring to freedom, not</pre>
|
153
151
|
|
@@ -258,8 +256,8 @@ list meets this criterion.</p>
|
|
258
256
|
|
259
257
|
The "source code" for a work means the preferred form of the work</pre>
|
260
258
|
|
261
|
-
<p>for making modifications to it. “Object code” means any non-source
|
262
|
-
|
259
|
+
<p>for making modifications to it. “Object code” means any non-source form of
|
260
|
+
a work.</p>
|
263
261
|
|
264
262
|
<pre>A "Standard Interface" means an interface that either is an official</pre>
|
265
263
|
|
@@ -273,11 +271,11 @@ widely used among developers working in that language.</p>
|
|
273
271
|
packaging a Major Component, but which is not part of that Major Component,
|
274
272
|
and (b) serves only to enable use of the work with that Major Component, or
|
275
273
|
to implement a Standard Interface for which an implementation is available
|
276
|
-
to the public in source code form. A “Major Component”, in this
|
277
|
-
|
278
|
-
|
279
|
-
|
280
|
-
|
274
|
+
to the public in source code form. A “Major Component”, in this context,
|
275
|
+
means a major essential component (kernel, window system, and so on) of the
|
276
|
+
specific operating system (if any) on which the executable work runs, or a
|
277
|
+
compiler used to produce the work, or an object code interpreter used to
|
278
|
+
run it.</p>
|
281
279
|
|
282
280
|
<pre>The "Corresponding Source" for a work in object code form means all</pre>
|
283
281
|
|
@@ -341,8 +339,8 @@ prohibiting or restricting circumvention of such measures.</p>
|
|
341
339
|
<p>circumvention of technological measures to the extent such circumvention is
|
342
340
|
effected by exercising rights under this License with respect to the
|
343
341
|
covered work, and you disclaim any intention to limit operation or
|
344
|
-
modification of the work as a means of enforcing, against the work’s
|
345
|
-
|
342
|
+
modification of the work as a means of enforcing, against the work’s users,
|
343
|
+
your or third parties’ legal rights to forbid circumvention of
|
346
344
|
technological measures.</p>
|
347
345
|
|
348
346
|
<pre>4. Conveying Verbatim Copies.
|
@@ -392,11 +390,11 @@ A compilation of a covered work with other separate and independent</pre>
|
|
392
390
|
|
393
391
|
<p>works, which are not by their nature extensions of the covered work, and
|
394
392
|
which are not combined with it such as to form a larger program, in or on a
|
395
|
-
volume of a storage or distribution medium, is called an “aggregate” if
|
396
|
-
|
397
|
-
|
398
|
-
|
399
|
-
|
393
|
+
volume of a storage or distribution medium, is called an “aggregate” if the
|
394
|
+
compilation and its resulting copyright are not used to limit the access or
|
395
|
+
legal rights of the compilation’s users beyond what the individual works
|
396
|
+
permit. Inclusion of a covered work in an aggregate does not cause this
|
397
|
+
License to apply to the other parts of the aggregate.</p>
|
400
398
|
|
401
399
|
<pre>6. Conveying Non-Source Forms.
|
402
400
|
|
@@ -458,13 +456,13 @@ conveying the object code work.</p>
|
|
458
456
|
household purposes, or (2) anything designed or sold for incorporation into
|
459
457
|
a dwelling. In determining whether a product is a consumer product,
|
460
458
|
doubtful cases shall be resolved in favor of coverage. For a particular
|
461
|
-
product received by a particular user, “normally used” refers to a
|
462
|
-
|
463
|
-
|
464
|
-
|
465
|
-
|
466
|
-
|
467
|
-
|
459
|
+
product received by a particular user, “normally used” refers to a typical
|
460
|
+
or common use of that class of product, regardless of the status of the
|
461
|
+
particular user or of the way in which the particular user actually uses,
|
462
|
+
or expects or is expected to use, the product. A product is a consumer
|
463
|
+
product regardless of whether the product has substantial commercial,
|
464
|
+
industrial or non-consumer uses, unless such uses represent the only
|
465
|
+
significant mode of use of the product.</p>
|
468
466
|
|
469
467
|
<pre>"Installation Information" for a User Product means any methods,</pre>
|
470
468
|
|
@@ -671,9 +669,9 @@ the contents of its contributor version.</p>
|
|
671
669
|
|
672
670
|
<p>agreement or commitment, however denominated, not to enforce a patent (such
|
673
671
|
as an express permission to practice a patent or covenant not to sue for
|
674
|
-
patent infringement). To “grant” such a patent license to a party
|
675
|
-
|
676
|
-
|
672
|
+
patent infringement). To “grant” such a patent license to a party means to
|
673
|
+
make such an agreement or commitment not to enforce a patent against the
|
674
|
+
party.</p>
|
677
675
|
|
678
676
|
<pre>If you convey a covered work, knowingly relying on a patent license,</pre>
|
679
677
|
|
@@ -684,10 +682,10 @@ you must either (1) cause the Corresponding Source to be so available, or
|
|
684
682
|
(2) arrange to deprive yourself of the benefit of the patent license for
|
685
683
|
this particular work, or (3) arrange, in a manner consistent with the
|
686
684
|
requirements of this License, to extend the patent license to downstream
|
687
|
-
recipients. “Knowingly relying” means you have actual knowledge that,
|
688
|
-
|
689
|
-
|
690
|
-
|
685
|
+
recipients. “Knowingly relying” means you have actual knowledge that, but
|
686
|
+
for the patent license, your conveying the covered work in a country, or
|
687
|
+
your recipient’s use of the covered work in a country, would infringe one
|
688
|
+
or more identifiable patents in that country that you have reason to
|
691
689
|
believe are valid.</p>
|
692
690
|
|
693
691
|
<pre>If, pursuant to or in connection with a single transaction or</pre>
|
@@ -764,9 +762,9 @@ Foundation.</p>
|
|
764
762
|
|
765
763
|
<pre>If the Program specifies that a proxy can decide which future</pre>
|
766
764
|
|
767
|
-
<p>versions of the GNU General Public License can be used, that proxy’s
|
768
|
-
|
769
|
-
|
765
|
+
<p>versions of the GNU General Public License can be used, that proxy’s public
|
766
|
+
statement of acceptance of a version permanently authorizes you to choose
|
767
|
+
that version for the Program.</p>
|
770
768
|
|
771
769
|
<pre>Later license versions may give you additional or different</pre>
|
772
770
|
|
@@ -778,12 +776,12 @@ or copyright holder as a result of your choosing to follow a later version.</p>
|
|
778
776
|
THERE IS NO WARRANTY FOR THE PROGRAM, TO THE EXTENT PERMITTED BY</pre>
|
779
777
|
|
780
778
|
<p>APPLICABLE LAW. EXCEPT WHEN OTHERWISE STATED IN WRITING THE COPYRIGHT
|
781
|
-
HOLDERS AND/OR OTHER PARTIES PROVIDE THE PROGRAM “AS IS” WITHOUT
|
782
|
-
|
783
|
-
|
784
|
-
|
785
|
-
|
786
|
-
|
779
|
+
HOLDERS AND/OR OTHER PARTIES PROVIDE THE PROGRAM “AS IS” WITHOUT WARRANTY
|
780
|
+
OF ANY KIND, EITHER EXPRESSED OR IMPLIED, INCLUDING, BUT NOT LIMITED TO,
|
781
|
+
THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
|
782
|
+
PURPOSE. THE ENTIRE RISK AS TO THE QUALITY AND PERFORMANCE OF THE PROGRAM
|
783
|
+
IS WITH YOU. SHOULD THE PROGRAM PROVE DEFECTIVE, YOU ASSUME THE COST OF
|
784
|
+
ALL NECESSARY SERVICING, REPAIR OR CORRECTION.</p>
|
787
785
|
|
788
786
|
<pre>16. Limitation of Liability.
|
789
787
|
|
@@ -850,15 +848,14 @@ This program comes with ABSOLUTELY NO WARRANTY; for details type `show w'.
|
|
850
848
|
This is free software, and you are welcome to redistribute it
|
851
849
|
under certain conditions; type `show c' for details.</pre>
|
852
850
|
|
853
|
-
<p>The hypothetical commands `show w’ and `show c’ should show the
|
854
|
-
|
855
|
-
|
856
|
-
an “about box”.</p>
|
851
|
+
<p>The hypothetical commands `show w’ and `show c’ should show the appropriate
|
852
|
+
parts of the General Public License. Of course, your program’s commands
|
853
|
+
might be different; for a GUI interface, you would use an “about box”.</p>
|
857
854
|
|
858
855
|
<pre>You should also get your employer (if you work as a programmer) or school,</pre>
|
859
856
|
|
860
|
-
<p>if any, to sign a “copyright disclaimer” for the program, if necessary.
|
861
|
-
|
857
|
+
<p>if any, to sign a “copyright disclaimer” for the program, if necessary. For
|
858
|
+
more information on this, and how to apply and follow the GNU GPL, see
|
862
859
|
<<a href="http://www.gnu.org/licenses/">www.gnu.org/licenses/</a>>.</p>
|
863
860
|
|
864
861
|
<pre>The GNU General Public License does not permit incorporating your program</pre>
|
data/doc/LibC.html
CHANGED
@@ -83,6 +83,8 @@
|
|
83
83
|
|
84
84
|
<li><a href="./Bio/DB/Alignment.html">Bio::DB::Alignment</a>
|
85
85
|
|
86
|
+
<li><a href="./Bio/DB/Pileup.html">Bio::DB::Pileup</a>
|
87
|
+
|
86
88
|
<li><a href="./Bio/DB/SAM.html">Bio::DB::SAM</a>
|
87
89
|
|
88
90
|
<li><a href="./Bio/DB/SAM/Library.html">Bio::DB::SAM::Library</a>
|
@@ -107,11 +109,9 @@
|
|
107
109
|
|
108
110
|
<li><a href="./Bio/DB/Tag.html">Bio::DB::Tag</a>
|
109
111
|
|
110
|
-
<li><a href="./
|
111
|
-
|
112
|
-
<li><a href="./Pileup.html">Pileup</a>
|
112
|
+
<li><a href="./Bio/DB/Vcf.html">Bio::DB::Vcf</a>
|
113
113
|
|
114
|
-
<li><a href="./
|
114
|
+
<li><a href="./LibC.html">LibC</a>
|
115
115
|
|
116
116
|
</ul>
|
117
117
|
</nav>
|
@@ -124,6 +124,8 @@
|
|
124
124
|
|
125
125
|
<div id="description" class="description">
|
126
126
|
|
127
|
+
<p>require ‘bio/db/pileup’ require ‘bio/db/vcf’</p>
|
128
|
+
|
127
129
|
</div><!-- description -->
|
128
130
|
|
129
131
|
|
data/doc/created.rid
CHANGED
@@ -1,10 +1,10 @@
|
|
1
|
-
|
1
|
+
Fri, 13 Jan 2012 10:07:15 +0000
|
2
|
+
./lib/bio/db/pileup.rb Thu, 12 Jan 2012 15:37:18 +0000
|
3
|
+
./lib/bio/db/sam/bam.rb Thu, 12 Jan 2012 15:31:33 +0000
|
2
4
|
./lib/bio/db/sam/faidx.rb Wed, 02 Nov 2011 11:50:27 +0000
|
3
|
-
./LICENSE.txt Wed, 02 Nov 2011 11:50:27 +0000
|
4
|
-
./lib/bio/db/sam/pileup.rb Wed, 02 Nov 2011 11:50:27 +0000
|
5
|
-
./lib/bio/db/sam/sam.rb Mon, 07 Nov 2011 10:04:59 +0000
|
6
5
|
./lib/bio/db/sam/library.rb Wed, 02 Nov 2011 13:28:41 +0000
|
6
|
+
./lib/bio/db/sam/sam.rb Mon, 07 Nov 2011 10:04:59 +0000
|
7
|
+
./lib/bio/db/sam.rb Thu, 12 Jan 2012 15:38:07 +0000
|
8
|
+
./lib/bio/db/vcf.rb Thu, 12 Jan 2012 15:52:37 +0000
|
7
9
|
./lib/bio-samtools.rb Wed, 02 Nov 2011 16:15:52 +0000
|
8
|
-
./
|
9
|
-
./lib/bio/db/sam/bam.rb Wed, 02 Nov 2011 13:43:07 +0000
|
10
|
-
./lib/bio/db/sam.rb Wed, 09 Nov 2011 21:03:13 +0000
|
10
|
+
./LICENSE.txt Wed, 02 Nov 2011 11:50:27 +0000
|
data/doc/index.html
CHANGED
@@ -65,6 +65,8 @@
|
|
65
65
|
|
66
66
|
<li><a href="./Bio/DB/Alignment.html">Bio::DB::Alignment</a>
|
67
67
|
|
68
|
+
<li><a href="./Bio/DB/Pileup.html">Bio::DB::Pileup</a>
|
69
|
+
|
68
70
|
<li><a href="./Bio/DB/SAM.html">Bio::DB::SAM</a>
|
69
71
|
|
70
72
|
<li><a href="./Bio/DB/SAM/Library.html">Bio::DB::SAM::Library</a>
|
@@ -89,11 +91,9 @@
|
|
89
91
|
|
90
92
|
<li><a href="./Bio/DB/Tag.html">Bio::DB::Tag</a>
|
91
93
|
|
92
|
-
<li><a href="./
|
93
|
-
|
94
|
-
<li><a href="./Pileup.html">Pileup</a>
|
94
|
+
<li><a href="./Bio/DB/Vcf.html">Bio::DB::Vcf</a>
|
95
95
|
|
96
|
-
<li><a href="./
|
96
|
+
<li><a href="./LibC.html">LibC</a>
|
97
97
|
|
98
98
|
</ul>
|
99
99
|
</nav>
|
data/doc/js/search_index.js
CHANGED
@@ -1 +1 @@
|
|
1
|
-
var search_data = {"index":{"info":[["Bio","","Bio.html","","<p>
|
1
|
+
var search_data = {"index":{"searchIndex":["bio","db","alignment","pileup","sam","library","tools","bam1coret","bam1t","bamheadert","bampileup1t","samfilet","samfiletx","samexception","sam","tag","vcf","libc","average_coverage()","bam_sort()","chromosome_coverage()","close()","deprecated_pileup()","fetch()","fetch_reference()","fetch_with_function()","filename()","finalize()","finalize()","index_stats()","int_or_raw()","load_index()","load_reference()","merge()","mpileup()","mpileup_plus()","new()","new()","new()","new()","new()","non_ref_count()","non_refs()","open()","parse_line()","qname()","query_string()","ref_count()","sam=()","set()","set()","text()","text=()","to_s()","license"],"longSearchIndex":["bio","bio::db","bio::db::alignment","bio::db::pileup","bio::db::sam","bio::db::sam::library","bio::db::sam::tools","bio::db::sam::tools::bam1coret","bio::db::sam::tools::bam1t","bio::db::sam::tools::bamheadert","bio::db::sam::tools::bampileup1t","bio::db::sam::tools::samfilet","bio::db::sam::tools::samfiletx","bio::db::samexception","bio::db::sam","bio::db::tag","bio::db::vcf","libc","bio::db::sam#average_coverage()","bio::db::sam::tools::bam_sort()","bio::db::sam#chromosome_coverage()","bio::db::sam#close()","bio::db::sam#deprecated_pileup()","bio::db::sam#fetch()","bio::db::sam#fetch_reference()","bio::db::sam#fetch_with_function()","bio::db::sam::library::filename()","bio::db::alignment::finalize()","bio::db::sam::finalize()","bio::db::sam#index_stats()","bio::db::vcf#int_or_raw()","bio::db::sam#load_index()","bio::db::sam#load_reference()","bio::db::sam::merge()","bio::db::sam#mpileup()","bio::db::sam#mpileup_plus()","bio::db::alignment::new()","bio::db::pileup::new()","bio::db::samexception::new()","bio::db::sam::new()","bio::db::vcf::new()","bio::db::pileup#non_ref_count()","bio::db::pileup#non_refs()","bio::db::sam#open()","bio::db::vcf#parse_line()","bio::db::sam::tools::bam1t#qname()","bio::db::sam#query_string()","bio::db::pileup#ref_count()","bio::db::alignment#sam=()","bio::db::alignment#set()","bio::db::tag#set()","bio::db::sam::tools::bamheadert#text()","bio::db::sam::tools::bamheadert#text=()","bio::db::sam#to_s()",""],"info":[["Bio","","Bio.html","","<p>Bio::DB::Pileup \n<p>A class representing information in SAMTools pileup format\n<p>Author — Dan MacLean (dan.maclean@tsl …\n"],["Bio::DB","","Bio/DB.html","",""],["Bio::DB::Alignment","","Bio/DB/Alignment.html","",""],["Bio::DB::Pileup","","Bio/DB/Pileup.html","",""],["Bio::DB::SAM","","Bio/DB/SAM.html","",""],["Bio::DB::SAM::Library","","Bio/DB/SAM/Library.html","",""],["Bio::DB::SAM::Tools","","Bio/DB/SAM/Tools.html","",""],["Bio::DB::SAM::Tools::Bam1CoreT","","Bio/DB/SAM/Tools/Bam1CoreT.html","",""],["Bio::DB::SAM::Tools::Bam1T","","Bio/DB/SAM/Tools/Bam1T.html","",""],["Bio::DB::SAM::Tools::BamHeaderT","","Bio/DB/SAM/Tools/BamHeaderT.html","","<p>typedef struct {\n\n<pre>int32_t n_targets;\nchar **target_name;\nuint32_t *target_len;\nvoid *dict, *hash, *rg2lib; ...</pre>\n"],["Bio::DB::SAM::Tools::BamPileup1T","","Bio/DB/SAM/Tools/BamPileup1T.html","",""],["Bio::DB::SAM::Tools::SamfileT","","Bio/DB/SAM/Tools/SamfileT.html","","<p>typedef struct {\n\n<pre>int type;\nunion {\n tamFile tamr;\n bamFile bam;\n FILE *tamw;\n} x; ...</pre>\n"],["Bio::DB::SAM::Tools::SamfileTX","","Bio/DB/SAM/Tools/SamfileTX.html","",""],["Bio::DB::SAMException","","Bio/DB/SAMException.html","",""],["Bio::DB::Sam","","Bio/DB/Sam.html","",""],["Bio::DB::Tag","","Bio/DB/Tag.html","",""],["Bio::DB::Vcf","","Bio/DB/Vcf.html","",""],["LibC","","LibC.html","","<p>require ‘bio/db/pileup’ require ‘bio/db/vcf’\n"],["average_coverage","Bio::DB::Sam","Bio/DB/Sam.html#method-i-average_coverage","(chromosome, qstart, len)","<p>Returns the average coverage of a region in a bam file.\n"],["bam_sort","Bio::DB::SAM::Tools","Bio/DB/SAM/Tools.html#method-c-bam_sort","(bam_filename, bam_output_prefix)",""],["chromosome_coverage","Bio::DB::Sam","Bio/DB/Sam.html#method-i-chromosome_coverage","(chromosome, qstart, len)","<p>Returns an array with the coverage at each possition in the queried region\nThis is a simple average coverage …\n"],["close","Bio::DB::Sam","Bio/DB/Sam.html#method-i-close","()","<p>Closes the sam file and destroys the C pointers using the functions\nprovided by libbam\n"],["deprecated_pileup","Bio::DB::Sam","Bio/DB/Sam.html#method-i-deprecated_pileup","( cmd )","<p>utility method that does not use the samtools API, it calls samtools\ndirectly as if on the command line …\n"],["fetch","Bio::DB::Sam","Bio/DB/Sam.html#method-i-fetch","(chromosome, qstart, qend)","<p>Returns an array of Alignments on a given region.\n"],["fetch_reference","Bio::DB::Sam","Bio/DB/Sam.html#method-i-fetch_reference","(chromosome, qstart,qend)","<p>Returns the sequence for a given region.\n"],["fetch_with_function","Bio::DB::Sam","Bio/DB/Sam.html#method-i-fetch_with_function","(chromosome, qstart, qend, function)","<p>Executes a function on each Alignment inside the queried region of the\nchromosome. The chromosome can …\n"],["filename","Bio::DB::SAM::Library","Bio/DB/SAM/Library.html#method-c-filename","()","<p>Return the path with the file name of the library for the specific\noperating system\n"],["finalize","Bio::DB::Alignment","Bio/DB/Alignment.html#method-c-finalize","(object_id)",""],["finalize","Bio::DB::Sam","Bio/DB/Sam.html#method-c-finalize","(id)","<p>Destructor method that closes the file before letting the object be garbage\ncollected.\n"],["index_stats","Bio::DB::Sam","Bio/DB/Sam.html#method-i-index_stats","()",""],["int_or_raw","Bio::DB::Vcf","Bio/DB/Vcf.html#method-i-int_or_raw","(x)",""],["load_index","Bio::DB::Sam","Bio/DB/Sam.html#method-i-load_index","()","<p>Loads the bam index to be used for fetching. If the index doesn’t exists\nthe index is built provided …\n"],["load_reference","Bio::DB::Sam","Bio/DB/Sam.html#method-i-load_reference","()","<p>Loads the reference file to be able to query regions of it. This requires\nthe fai index to exist in the …\n"],["merge","Bio::DB::Sam","Bio/DB/Sam.html#method-c-merge","(files, merged_file, headers, add_RG, by_qname)","<p>Merges n BAM files. This doesn’t require to create a SAM object\n<p>files — An array with the paths to the files. …\n"],["mpileup","Bio::DB::Sam","Bio/DB/Sam.html#method-i-mpileup","( opts )","<p>calls the mpileup function, opts is a hash of options identical to the\ncommand line options for mpileup. …\n"],["mpileup_plus","Bio::DB::Sam","Bio/DB/Sam.html#method-i-mpileup_plus","( opts )","<p>experimental method that spawns a samtools mpileup | bcftools view process\nand supports returning of …\n"],["new","Bio::DB::Alignment","Bio/DB/Alignment.html#method-c-new","()",""],["new","Bio::DB::Pileup","Bio/DB/Pileup.html#method-c-new","(pile_up_line)","<p>creates the Pileup object\n\n<pre>pile_up_line = "seq2\\t151\\tG\\tG\\t36\\t0\\t99\\t12\\t...........A\\t:9<;;7=<<<<<" ...</pre>\n"],["new","Bio::DB::SAMException","Bio/DB/SAMException.html#method-c-new","()","<p>we can add further variables to give information of the excpetion\n"],["new","Bio::DB::Sam","Bio/DB/Sam.html#method-c-new","(optsa={})","<p>To make a new sam object. Initialize expects a hash optsa with the\nfollowing elemets:\n<p>fasta — The fasta …\n"],["new","Bio::DB::Vcf","Bio/DB/Vcf.html#method-c-new","(line=nil, sample_names=nil)","<p>create the vcf object, use the ordered list of sample names to label\nsamples if provided [‘A’, ‘B’, ‘C’], …\n"],["non_ref_count","Bio::DB::Pileup","Bio/DB/Pileup.html#method-i-non_ref_count","()","<p>returns the total non-reference bases in the reads at this position\n"],["non_refs","Bio::DB::Pileup","Bio/DB/Pileup.html#method-i-non_refs","()","<p>Calculate the total count of each non-reference nucleotide and return a\nhash of all 4 nt counts, returns …\n"],["open","Bio::DB::Sam","Bio/DB/Sam.html#method-i-open","()","<p>Function that actually opens the sam file Throws a SAMException if the file\ncan’t be open.\n"],["parse_line","Bio::DB::Vcf","Bio/DB/Vcf.html#method-i-parse_line","(line, sample_names=nil)","<p>gets the info in the Vcf lines and parses it, setting the attributes\n"],["qname","Bio::DB::SAM::Tools::Bam1T","Bio/DB/SAM/Tools/Bam1T.html#method-i-qname","()",""],["query_string","Bio::DB::Sam","Bio/DB/Sam.html#method-i-query_string","(chromosome, qstart,qend)","<p>Generates a query sting to be used by the region parser in samtools. In\nprinciple, you shouldn’t need …\n"],["ref_count","Bio::DB::Pileup","Bio/DB/Pileup.html#method-i-ref_count","()","<p>returns the count of reference-bases in the reads at this position\n"],["sam=","Bio::DB::Alignment","Bio/DB/Alignment.html#method-i-sam-3D","(sam)",""],["set","Bio::DB::Alignment","Bio/DB/Alignment.html#method-i-set","(bam_alignment, header)",""],["set","Bio::DB::Tag","Bio/DB/Tag.html#method-i-set","(str)",""],["text","Bio::DB::SAM::Tools::BamHeaderT","Bio/DB/SAM/Tools/BamHeaderT.html#method-i-text","()",""],["text=","Bio::DB::SAM::Tools::BamHeaderT","Bio/DB/SAM/Tools/BamHeaderT.html#method-i-text-3D","(str)",""],["to_s","Bio::DB::Sam","Bio/DB/Sam.html#method-i-to_s","()","<p>Prints a description of the sam file in a text format containg if it is\nbinary or text, the path and …\n"],["LICENSE","","LICENSE_txt.html","","<p>Copyright © 2011 Raoul J.P. Bonnal\n<p>Permission is hereby granted, free of charge, to any person obtaining …\n"]]}}
|
@@ -65,6 +65,8 @@
|
|
65
65
|
|
66
66
|
<li><a href="../Bio/DB/Alignment.html">Bio::DB::Alignment</a>
|
67
67
|
|
68
|
+
<li><a href="../Bio/DB/Pileup.html">Bio::DB::Pileup</a>
|
69
|
+
|
68
70
|
<li><a href="../Bio/DB/SAM.html">Bio::DB::SAM</a>
|
69
71
|
|
70
72
|
<li><a href="../Bio/DB/SAM/Library.html">Bio::DB::SAM::Library</a>
|
@@ -89,11 +91,9 @@
|
|
89
91
|
|
90
92
|
<li><a href="../Bio/DB/Tag.html">Bio::DB::Tag</a>
|
91
93
|
|
92
|
-
<li><a href="../
|
93
|
-
|
94
|
-
<li><a href="../Pileup.html">Pileup</a>
|
94
|
+
<li><a href="../Bio/DB/Vcf.html">Bio::DB::Vcf</a>
|
95
95
|
|
96
|
-
<li><a href="../
|
96
|
+
<li><a href="../LibC.html">LibC</a>
|
97
97
|
|
98
98
|
</ul>
|
99
99
|
</nav>
|
@@ -0,0 +1,171 @@
|
|
1
|
+
<!DOCTYPE html>
|
2
|
+
|
3
|
+
<html>
|
4
|
+
<head>
|
5
|
+
<meta content="text/html; charset=UTF-8" http-equiv="Content-Type">
|
6
|
+
|
7
|
+
<title>LICENSE</title>
|
8
|
+
|
9
|
+
<link type="text/css" media="screen" href="../../../rdoc.css" rel="stylesheet">
|
10
|
+
|
11
|
+
<script type="text/javascript">
|
12
|
+
var rdoc_rel_prefix = "../../../";
|
13
|
+
</script>
|
14
|
+
|
15
|
+
<script type="text/javascript" charset="utf-8" src="../../../js/jquery.js"></script>
|
16
|
+
<script type="text/javascript" charset="utf-8" src="../../../js/navigation.js"></script>
|
17
|
+
<script type="text/javascript" charset="utf-8" src="../../../js/search_index.js"></script>
|
18
|
+
<script type="text/javascript" charset="utf-8" src="../../../js/search.js"></script>
|
19
|
+
<script type="text/javascript" charset="utf-8" src="../../../js/searcher.js"></script>
|
20
|
+
<script type="text/javascript" charset="utf-8" src="../../../js/darkfish.js"></script>
|
21
|
+
|
22
|
+
|
23
|
+
<body class="file">
|
24
|
+
<nav id="metadata">
|
25
|
+
<nav id="home-section" class="section">
|
26
|
+
<h3 class="section-header">
|
27
|
+
<a href="../../../index.html">Home</a>
|
28
|
+
<a href="../../../table_of_contents.html#classes">Classes</a>
|
29
|
+
<a href="../../../table_of_contents.html#methods">Methods</a>
|
30
|
+
</h3>
|
31
|
+
</nav>
|
32
|
+
|
33
|
+
|
34
|
+
<nav id="search-section" class="section project-section" class="initially-hidden">
|
35
|
+
<form action="#" method="get" accept-charset="utf-8">
|
36
|
+
<h3 class="section-header">
|
37
|
+
<input type="text" name="search" placeholder="Search" id="search-field"
|
38
|
+
title="Type to search, Up and Down to navigate, Enter to load">
|
39
|
+
</h3>
|
40
|
+
</form>
|
41
|
+
|
42
|
+
<ul id="search-results" class="initially-hidden"></ul>
|
43
|
+
</nav>
|
44
|
+
|
45
|
+
|
46
|
+
<div id="project-metadata">
|
47
|
+
<nav id="fileindex-section" class="section project-section">
|
48
|
+
<h3 class="section-header">Pages</h3>
|
49
|
+
|
50
|
+
<ul>
|
51
|
+
|
52
|
+
<li class="file"><a href="../../../LICENSE_txt.html">LICENSE</a>
|
53
|
+
|
54
|
+
</ul>
|
55
|
+
</nav>
|
56
|
+
|
57
|
+
<nav id="classindex-section" class="section project-section">
|
58
|
+
<h3 class="section-header">Class and Module Index</h3>
|
59
|
+
|
60
|
+
<ul class="link-list">
|
61
|
+
|
62
|
+
<li><a href="../../../Bio.html">Bio</a>
|
63
|
+
|
64
|
+
<li><a href="../../../Bio/DB.html">Bio::DB</a>
|
65
|
+
|
66
|
+
<li><a href="../../../Bio/DB/Alignment.html">Bio::DB::Alignment</a>
|
67
|
+
|
68
|
+
<li><a href="../../../Bio/DB/Pileup.html">Bio::DB::Pileup</a>
|
69
|
+
|
70
|
+
<li><a href="../../../Bio/DB/SAM.html">Bio::DB::SAM</a>
|
71
|
+
|
72
|
+
<li><a href="../../../Bio/DB/SAM/Library.html">Bio::DB::SAM::Library</a>
|
73
|
+
|
74
|
+
<li><a href="../../../Bio/DB/SAM/Tools.html">Bio::DB::SAM::Tools</a>
|
75
|
+
|
76
|
+
<li><a href="../../../Bio/DB/SAM/Tools/Bam1CoreT.html">Bio::DB::SAM::Tools::Bam1CoreT</a>
|
77
|
+
|
78
|
+
<li><a href="../../../Bio/DB/SAM/Tools/Bam1T.html">Bio::DB::SAM::Tools::Bam1T</a>
|
79
|
+
|
80
|
+
<li><a href="../../../Bio/DB/SAM/Tools/BamHeaderT.html">Bio::DB::SAM::Tools::BamHeaderT</a>
|
81
|
+
|
82
|
+
<li><a href="../../../Bio/DB/SAM/Tools/BamPileup1T.html">Bio::DB::SAM::Tools::BamPileup1T</a>
|
83
|
+
|
84
|
+
<li><a href="../../../Bio/DB/SAM/Tools/SamfileT.html">Bio::DB::SAM::Tools::SamfileT</a>
|
85
|
+
|
86
|
+
<li><a href="../../../Bio/DB/SAM/Tools/SamfileTX.html">Bio::DB::SAM::Tools::SamfileTX</a>
|
87
|
+
|
88
|
+
<li><a href="../../../Bio/DB/SAMException.html">Bio::DB::SAMException</a>
|
89
|
+
|
90
|
+
<li><a href="../../../Bio/DB/Sam.html">Bio::DB::Sam</a>
|
91
|
+
|
92
|
+
<li><a href="../../../Bio/DB/Tag.html">Bio::DB::Tag</a>
|
93
|
+
|
94
|
+
<li><a href="../../../Bio/DB/Vcf.html">Bio::DB::Vcf</a>
|
95
|
+
|
96
|
+
<li><a href="../../../LibC.html">LibC</a>
|
97
|
+
|
98
|
+
</ul>
|
99
|
+
</nav>
|
100
|
+
|
101
|
+
</div>
|
102
|
+
</nav>
|
103
|
+
|
104
|
+
<div id="documentation" class="description">
|
105
|
+
|
106
|
+
<h1 id="label-Bio%3A%3ADB%3A%3APileup+"><a href="../../../Bio/DB/Pileup.html">Bio::DB::Pileup</a> </h1>
|
107
|
+
|
108
|
+
<p>A class representing information in SAMTools pileup format</p>
|
109
|
+
<dl class="rdoc-list note-list"><dt>Author
|
110
|
+
<dd>
|
111
|
+
<p>Dan MacLean (tsl at dan.maclean.ac.uk)</p>
|
112
|
+
</dd></dl>
|
113
|
+
|
114
|
+
<p>Pileup is described at <a
|
115
|
+
href="http://sourceforge.net/apps/mediawiki/samtools/index.php?title=SAM_FAQ#I_do_not_understand_the_columns_in_the_pileup_output">sourceforge.net/apps/mediawiki/samtools/index.php?title=SAM_FAQ#I_do_not_understand_the_columns_in_the_pileup_output</a>.
|
116
|
+
Briefly (when you invoke pileup with the -c option):</p>
|
117
|
+
<ul><li>
|
118
|
+
<p>1 reference sequence name</p>
|
119
|
+
</li><li>
|
120
|
+
<p>2 reference coordinate</p>
|
121
|
+
</li><li>
|
122
|
+
<p>(3) reference base, or `*‘ for an indel line</p>
|
123
|
+
</li><li>
|
124
|
+
<p>(4) genotype where heterozygotes are encoded in the IUB code: M=A/C, R=A/G,
|
125
|
+
W=A/T, S=C/G, Y=C/T and K=G/T; indels are indicated by, for example, +A,
|
126
|
+
-A or +CC/-C. There is no difference between */+A or +A/*.</p>
|
127
|
+
</li><li>
|
128
|
+
<p>(5) Phred-scaled likelihood that the genotype is wrong, which is also
|
129
|
+
called `consensus quality’.</p>
|
130
|
+
</li><li>
|
131
|
+
<p>(6) Phred-scaled likelihood that the genotype is identical to the
|
132
|
+
reference, which is also called `SNP quality’. Suppose the reference base
|
133
|
+
is A and in alignment we see 17 G and 3 A. We will get a low consensus
|
134
|
+
quality because it is difficult to distinguish an A/G heterozygote from a
|
135
|
+
G/G homozygote. We will get a high SNP quality, though, because the
|
136
|
+
evidence of a SNP is very strong.</p>
|
137
|
+
</li><li>
|
138
|
+
<p>(7) root mean square (RMS) mapping quality</p>
|
139
|
+
</li><li>
|
140
|
+
<p>8 # reads covering the position</p>
|
141
|
+
</li><li>
|
142
|
+
<p>9 read bases at a SNP line (check the manual page for more information);
|
143
|
+
the 1st indel allele otherwise</p>
|
144
|
+
</li><li>
|
145
|
+
<p>10 base quality at a SNP line; the 2nd indel allele otherwise</p>
|
146
|
+
</li><li>
|
147
|
+
<p>(11) indel line only: # reads directly supporting the 1st indel allele</p>
|
148
|
+
</li><li>
|
149
|
+
<p>(12) indel line only: # reads directly supporting the 2nd indel allele</p>
|
150
|
+
</li><li>
|
151
|
+
<p>(13) indel line only: # reads supporting a third indel allele</p>
|
152
|
+
</li></ul>
|
153
|
+
|
154
|
+
<p>If pileup is invoked without `-c’, indel lines and columns between 3 and 7
|
155
|
+
inclusive will not be outputted.</p>
|
156
|
+
|
157
|
+
<p>NB mpileup uses the 6 column output format eg
|
158
|
+
“seq2t151tGtGt36t0t99t12t.….……At:9<;;7=<<<<<” Pileup
|
159
|
+
provides accessors for all columns (6 or 10 column format) and a few other
|
160
|
+
useful methods</p>
|
161
|
+
|
162
|
+
</div>
|
163
|
+
|
164
|
+
|
165
|
+
|
166
|
+
<footer id="validator-badges">
|
167
|
+
<p><a href="http://validator.w3.org/check/referer">[Validate]</a>
|
168
|
+
<p>Generated by <a href="https://github.com/rdoc/rdoc">RDoc</a> 3.11.
|
169
|
+
<p>Generated with the <a href="http://deveiate.org/projects/Darkfish-Rdoc/">Darkfish Rdoc Generator</a> 3.
|
170
|
+
</footer>
|
171
|
+
|