-
Notifications
You must be signed in to change notification settings - Fork 587
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
Calculate and Store site-level QCs #7197
Conversation
CreateFilteringFiles \ | ||
--ref-version 38 \ | ||
--filter-set-name ~{filter_set_name} \ | ||
-mode SNP \ | ||
-V ~{snp_recal_file} \ | ||
-O ~{filter_set_name}.snps.recal.tsv | ||
|
||
gatk --java-options "-Xmx4g" \ | ||
gatk --java-options "-Xmx1g" \ |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
nice catch!
omitFromCommandLine = true | ||
) | ||
public final class CreateSiteFilteringFiles extends VariantWalker { | ||
static final Logger logger = LogManager.getLogger(CreateVariantIngestFiles.class); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
should this use CreateSiteFilteringFiles.class instead?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
yup
|
||
private SimpleXSVWriter writer; | ||
|
||
private List<String> HEADER = |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
for consistency with other tools, move this to SchemaUtils ?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
my future self thanks you
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
looks great, only minor comments/questions
@@ -234,6 +252,14 @@ private void processVQSRRecordForPosition(ExtractFeaturesRecord rec) { | |||
} | |||
builder.attribute(GATKVCFConstants.EXCESS_HET_KEY, String.format("%.3f", excessHetApprox)); | |||
|
|||
if (rec.getDistinctAlleles() > 6) { |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
is this a constant defined somewhere that we can use? (or should we make it one?)
@@ -75,6 +75,9 @@ public static VCFFormatHeaderLine getEquivalentFormatHeaderLine(final String inf | |||
addFilterLine(new VCFFilterHeaderLine(VCFConstants.PASSES_FILTERS_v4, "Site contains at least one allele that passes filters")); | |||
|
|||
addFilterLine(new VCFFilterHeaderLine(NAY_FROM_YNG, "Considered a NAY in the Yay, Nay, Grey table")); | |||
addFilterLine(new VCFFilterHeaderLine(EXCESS_ALLELES, "Site has an excess of alternate alleles")); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
if the excess alleles number (6) was a constant somewhere, you could put that into the header. (maybe users don't care, but i would personally find it helpful!)
Addresses 219
Major changes
feature_extract.sql
location
fieldsMinor changes