Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Add utility methods for adam-shell #1958

Merged
merged 1 commit into from
Mar 23, 2018

Conversation

heuermh
Copy link
Member

@heuermh heuermh commented Mar 21, 2018

E.g.

scala> import org.bdgenomics.adam.rdd.ADAMContext._
import org.bdgenomics.adam.rdd.ADAMContext._

scala> import org.bdgenomics.adam.util.ADAMShell._
import org.bdgenomics.adam.util.ADAMShell._

scala> val variants = sc.loadVariants("Mills_and_1000G_gold_standard.indels.variants.adam")
variants: org.bdgenomics.adam.rdd.variant.VariantRDD = ParquetUnboundVariantRDD with 3366 reference sequences

scala> printVariantFilters(variants)
Filter Header Lines
FILTER=<ID=PASS,Description="Truth sensitivity tranche level at VQS Lod < -381.3592">

Variant Filters
+-------------+--------+--------+-------+--------+-----------------+----------------+----------------+
| Contig Name | Start  |  End   |  Ref  |  Alt   | Filters Applied | Filters Passed | Filters Failed |
+-------------+--------+--------+-------+--------+-----------------+----------------+----------------+
|        chr1 |  55248 |  55249 | C     | CTATGG |            true |           true |             [] |
|        chr1 |  87113 |  87115 | CT    | C      |            true |           true |             [] |
|        chr1 | 263834 | 263837 | CCT   | C      |            true |           true |             [] |
|        chr1 | 756186 | 756191 | CACAG | C      |            true |           true |             [] |
|        chr1 | 779928 | 779930 | AG    | A      |            true |           true |             [] |
|        chr1 | 786693 | 786694 | A     | ATT    |            true |           true |             [] |
|        chr1 | 796753 | 796754 | T     | TA     |            true |           true |             [] |
|        chr1 | 816926 | 816928 | AT    | A      |            true |           true |             [] |
|        chr1 | 822795 | 822798 | CCT   | C      |            true |           true |             [] |
|        chr1 | 824063 | 824068 | ATATT | A      |            true |           true |             [] |
+-------------+--------+--------+-------+--------+-----------------+----------------+----------------+

@AmplabJenkins
Copy link

Test PASSed.
Refer to this link for build results (access rights to CI server needed):
https://amplab.cs.berkeley.edu/jenkins//job/ADAM-prb/2732/
Test PASSed.

Copy link
Member

@fnothaft fnothaft left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Nice! This is a great idea! Did you want to add to this, or would you be OK with me merging this into 0.24.0?

@heuermh
Copy link
Member Author

heuermh commented Mar 21, 2018

Did you want to add to this, or would you be OK with me merging this into 0.24.0?

There is probably more to add, but if you think the pattern is right, then merging into 0.24.0 is ok with me. I'll squash into one commit.

@AmplabJenkins
Copy link

Test PASSed.
Refer to this link for build results (access rights to CI server needed):
https://amplab.cs.berkeley.edu/jenkins//job/ADAM-prb/2733/
Test PASSed.

@fnothaft fnothaft merged commit b86c4ad into bigdatagenomics:master Mar 23, 2018
@fnothaft
Copy link
Member

Merged! Thanks @heuermh!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

3 participants