From 7c4c07abb60b9a194f0ae0e731b1f9cd80ea55d1 Mon Sep 17 00:00:00 2001 From: mizraelson <18702359+mizraelson@users.noreply.github.com> Date: Thu, 26 Sep 2024 01:04:28 +0000 Subject: [PATCH] regression tests automated change --- ...gen-human-rna-tcr-umi-targeted-qiaseq.yaml | 477 ++++++++++++++++++ ...gen-mouse-rna-tcr-umi-targeted-qiaseq.yaml | 477 ++++++++++++++++++ regression/presets/list.txt | 4 + 3 files changed, 958 insertions(+) create mode 100644 regression/presets/analyze/qiagen-human-rna-tcr-umi-targeted-qiaseq.yaml create mode 100644 regression/presets/analyze/qiagen-mouse-rna-tcr-umi-targeted-qiaseq.yaml diff --git a/regression/presets/analyze/qiagen-human-rna-tcr-umi-targeted-qiaseq.yaml b/regression/presets/analyze/qiagen-human-rna-tcr-umi-targeted-qiaseq.yaml new file mode 100644 index 000000000..5c9eef904 --- /dev/null +++ b/regression/presets/analyze/qiagen-human-rna-tcr-umi-targeted-qiaseq.yaml @@ -0,0 +1,477 @@ +flags: [] +pipeline: + - align + - refineTagsAndSort + - assemble + - exportClones +align: + species: hsa + libraryName: default + trimmingQualityThreshold: 10 + trimmingWindowSize: 6 + chains: ALL + replaceWildcards: true + overlapPairedReads: true + bamDropNonVDJ: false + writeFailedAlignments: false + tagPattern: "^(R1:*) \\ ^(UMI:N{14:18})atgcattggagtcct" + tagUnstranded: false + tagMaxBudget: 10.0 + headerExtractors: [] + readIdAsCellTag: false + sampleTable: null + tagsValidations: + - type: MustContainTagType + tagType: Molecule + splitBySample: true + limit: null + parameters: + vParameters: + geneFeatureToAlign: VTranscriptWithP + minSumScore: 150 + relativeMinScore: 0.97 + maxHits: 5 + parameters: + type: kaligner2 + mapperNValue: 8 + mapperKValue: 1 + floatingLeftBound: false + floatingRightBound: true + mapperAbsoluteMinClusterScore: 102 + mapperExtraClusterScore: -38 + mapperMatchScore: 98 + mapperMismatchScore: -14 + mapperOffsetShiftScore: -82 + mapperSlotCount: 8 + mapperMaxClusters: 11 + mapperMaxClusterIndels: 4 + mapperKMersPerPosition: 4 + mapperAbsoluteMinScore: 103 + mapperRelativeMinScore: 0.81 + mapperMinSeedsDistance: 8 + mapperMaxSeedsDistance: 8 + alignmentStopPenalty: 0 + absoluteMinScore: 150 + relativeMinScore: 0.8 + maxHits: 2 + scoring: + type: affine + alphabet: nucleotide + subsMatrix: "simple(match = 10, mismatch = -19)" + gapOpenPenalty: -62 + gapExtensionPenalty: -11 + dParameters: + geneFeatureToAlign: DRegionWithP + relativeMinScore: 0.85 + absoluteMinScore: 25.0 + maxHits: 3 + scoring: + type: linear + alphabet: nucleotide + subsMatrix: "simple(match = 5, mismatch = -9)" + gapPenalty: -12 + jParameters: + geneFeatureToAlign: JRegionWithP + minSumScore: 150 + relativeMinScore: 0.97 + maxHits: 5 + parameters: + type: kaligner2 + mapperNValue: 8 + mapperKValue: 1 + floatingLeftBound: true + floatingRightBound: false + mapperAbsoluteMinClusterScore: 102 + mapperExtraClusterScore: -38 + mapperMatchScore: 95 + mapperMismatchScore: -14 + mapperOffsetShiftScore: -82 + mapperSlotCount: 6 + mapperMaxClusters: 4 + mapperMaxClusterIndels: 4 + mapperKMersPerPosition: 4 + mapperAbsoluteMinScore: 100 + mapperRelativeMinScore: 0.8 + mapperMinSeedsDistance: 5 + mapperMaxSeedsDistance: 5 + alignmentStopPenalty: 0 + absoluteMinScore: 150 + relativeMinScore: 0.8 + maxHits: 3 + scoring: + type: affine + alphabet: nucleotide + subsMatrix: "simple(match = 10, mismatch = -19)" + gapOpenPenalty: -62 + gapExtensionPenalty: -11 + cParameters: + geneFeatureToAlign: CExon1 + minSumScore: 40 + relativeMinScore: 0.97 + maxHits: 5 + parameters: + type: kaligner + mapperKValue: 5 + floatingLeftBound: false + floatingRightBound: true + mapperAbsoluteMinScore: 1.5 + mapperRelativeMinScore: 0.75 + mapperMatchScore: 1.0 + mapperMismatchPenalty: -0.1 + mapperOffsetShiftPenalty: -0.3 + mapperMinSeedsDistance: 4 + mapperMaxSeedsDistance: 10 + minAlignmentLength: 15 + maxAdjacentIndels: 2 + alignmentStopPenalty: -1000 + absoluteMinScore: 40.0 + relativeMinScore: 0.87 + maxHits: 7 + scoring: + type: linear + alphabet: nucleotide + subsMatrix: "simple(match = 5, mismatch = -9)" + gapPenalty: -12 + vjAlignmentOrder: VThenJ + libraryStructure: Unknown + includeDScore: false + includeCScore: false + minSumScore: 120.0 + relativeMinVFR3CDR3Score: 0.7 + allowPartialAlignments: false + allowNoCDR3PartAlignments: false + allowChimeras: false + readsLayout: ReverseOnly + mergerParameters: + qualityMergingAlgorithm: MaxSubtraction + partsLayout: null + minimalOverlap: 13 + minimalMatchQualitySum: 364 + maxQuality: 50 + minimalIdentity: 0.7 + identityType: MinimalQualityWeighted + fixSeed: true + alignmentBoundaryTolerance: 5 + minChimeraDetectionScore: 120 + vjOverlapWindow: 3 + saveOriginalSequence: false + saveOriginalReads: false + smartForceEdgeAlignments: true + retriesForPartialAlignments: 3 + retriesAlignmentOfNotCoveredPart: 2 +refineTagsAndSort: + whitelists: {} + runCorrection: true + parameters: + correctionPower: 0.001 + backgroundSubstitutionRate: 0.001 + backgroundIndelRate: 1.0E-5 + minQuality: 12 + maxSubstitutions: 2 + maxIndels: 2 + maxTotalErrors: 3 + postFilter: + type: filter_groups + groupingKeys: + - allTags:Molecule + predicates: + - metrics: + - type: group_metric_sum_weight + reportHist: + log: true + binNumber: 0 + minBinWidth: 0.2 + multiplyWeightByKey: false + operator: + type: group_operator_lowest_threshold + operators: + - type: group_operator_advanced_thresholding + algo: + type: otsu + histTransformers: + - type: inflate + fraction: 0.15 + value: 1.0 + logX: true + minimalSample: 20 + fallbackThreshold: 1.0 + accept: High + - type: group_operator_cumtop + share: 0.85 + round: Down + accept: High + accept: High + expectedSorting: [] + requiredSequences: [] + expectedSorting: + - allTags:Molecule + requiredSequences: [] +assemblePartial: + overlappedOnly: false + dropPartial: false + cellLevel: false + parameters: + kValue: 12 + kOffset: -7 + minimalAssembleOverlap: 12 + minimalNOverlap: 7 + minimalNOverlapShare: 0.65 + minimalAlignmentMergeIdentity: 0.85 + mergerParameters: + qualityMergingAlgorithm: SumSubtraction + partsLayout: CollinearDirect + minimalOverlap: 20 + minimalMatchQualitySum: 0 + maxQuality: 45 + minimalIdentity: 0.95 + identityType: Unweighted + maxLeftParts: 256000 + maxLeftMatches: 6144 +extend: + vAnchor: CDR3Begin + jAnchor: FR4Begin + minimalVScore: 50 + minimalJScore: 50 +assemble: + sortBySequence: false + clnaOutput: false + cellLevel: false + consensusAssemblerParameters: + assembler: + aAssemblerParameters: + bandWidth: 4 + scoring: + type: linear + alphabet: nucleotide + subsMatrix: "simple(match = 5, mismatch = -4)" + gapPenalty: -14 + minAlignmentScore: 40 + maxNormalizedAlignmentPenalty: 0.2 + trimMinimalSumQuality: 0 + trimReferenceRegion: false + maxQuality: 45 + maxIterations: 6 + minAltSeedQualityScore: 11 + minAltSeedNormalizedPenalty: 0.35 + altSeedPenaltyTolerance: 0.5 + minRecordSharePerConsensus: 0.2 + minRecordsPerConsensus: 0 + minRecursiveRecordShare: 0.4 + minQualityScore: 0 + maxConsensuses: 3 + minTagSuffixShare: 0.0 + isolateChains: false + empiricalAlignmentAssignment: + uniqGeneOccurrenceWithinTag: + - V + - J + uniqGeneOccurrenceWithinCluster: + - V + - J + matchAlignmentsWithoutVJHitsByTags: true + cloneAssemblerParameters: + assemblingFeatures: + - CDR3 + minimalClonalSequenceLength: 12 + qualityAggregationType: BetaScore + cloneClusteringParameters: + searchDepth: 2 + allowedMutationsInNRegions: 1 + searchParameters: twoMismatchesOrIndels + clusteringFilter: + type: advanced + correctionPower: 0.001 + backgroundSubstitutionRate: 5.0E-4 + backgroundIndelRate: 2.0E-4 + cloneFactoryParameters: + vParameters: + maxAlignmentWidthLinear: 5 + maxAlignmentWidthAffine: 500 + jParameters: + maxAlignmentWidthLinear: 5 + maxAlignmentWidthAffine: 500 + cParameters: + maxAlignmentWidthLinear: 5 + maxAlignmentWidthAffine: 500 + dParameters: + relativeMinScore: null + absoluteMinScore: null + maxHits: null + scoring: null + separateByV: false + separateByJ: false + separateByC: false + maximalPreClusteringRatio: 1.0 + preClusteringScoreFilteringRatio: 2.0 + preClusteringCountFilteringRatio: 2.0 + addReadsCountOnClustering: false + badQualityThreshold: 15 + maxBadPointsPercent: 0.7 + mappingThreshold: 2of5 + minimalQuality: 20 + postFilters: null + inferMinRecordsPerConsensus: true +assembleContigs: + ignoreTags: false + parameters: + branchingMinimalQualityShare: 0.1 + branchingMinimalSumQuality: 60 + decisiveBranchingSumQualityThreshold: 120 + alignedSequenceEdgeDelta: 3 + alignmentEdgeRegionSize: 7 + minimalNonEdgePointsFraction: 0.25 + minimalMeanNormalizedQuality: 5.0 + outputMinimalQualityShare: 0.75 + outputMinimalSumQuality: 0 + subCloningRegions: null + assemblingRegions: null + postFiltering: + type: NoFiltering + trimmingParameters: + averageQualityThreshold: 10.0 + windowSize: 8 + minimalContigLength: 20 + alignedRegionsOnly: false + discardAmbiguousNucleotideCalls: false + cloneFactoryParameters: + vParameters: + relativeMinScore: 0.97 + maxAlignmentWidthLinear: 5 + maxAlignmentWidthAffine: 500 + jParameters: + maxAlignmentWidthLinear: 5 + maxAlignmentWidthAffine: 500 + cParameters: + maxAlignmentWidthLinear: 5 + maxAlignmentWidthAffine: 500 + dParameters: + relativeMinScore: null + absoluteMinScore: null + maxHits: null + scoring: null +exportAlignments: + chains: ALL + noHeader: false + fields: + - field: -readIds + - field: -tags + args: + - Molecule + - field: -targetSequences + - field: -targetQualities + - field: -vHitsWithScore + - field: -dHitsWithScore + - field: -jHitsWithScore + - field: -cHitsWithScore + - field: -vAlignments + - field: -dAlignments + - field: -jAlignments + - field: -cAlignments + - field: -allNFeaturesWithMinQuality + - field: -allAAFeatures + - field: -defaultAnchorPoints + - field: -topChains +exportClones: + splitByTagType: null + filterOutOfFrames: false + filterStops: false + chains: ALL + noHeader: false + fields: + - field: -cloneId + - field: -readCount + - field: -readFraction + - field: -uniqueTagCount + args: + - Molecule + - field: -uniqueTagFraction + args: + - Molecule + - field: -targetSequences + - field: -targetQualities + - field: -vHitsWithScore + - field: -dHitsWithScore + - field: -jHitsWithScore + - field: -cHitsWithScore + - field: -vAlignments + - field: -dAlignments + - field: -jAlignments + - field: -cAlignments + - field: -allNFeaturesWithMinQuality + - field: -allAAFeatures + - field: -defaultAnchorPoints + splitFilesBy: + - chain + groupClonesBy: [] +qc: + checks: + - type: SuccessfullyAlignedReads + upper: 0.85 + middle: 0.7 + label: Successfully aligned reads + - type: OffTargetReads + upper: 0.2 + middle: 0.1 + label: Off target (non TCR/IG) reads + - type: ReadsWithNoVOrJHits + upper: 0.2 + middle: 0.1 + label: Reads with no V or J hits + - type: AlignmentsWithNoAssemblingFeature + upper: 0.15 + middle: 0.05 + label: Alignments without assembling feature + - type: ReadsWithNoBarcode + upper: 0.2 + middle: 0.1 + label: Reads with no barcode + - type: ReadsDroppedInTagRefinement + upper: 0.1 + middle: 0.05 + label: Reads dropped in tags error correction and filtering + - type: TagArtificialDiversityEliminated + tag: UMI + upper: 0.5 + middle: 0.3 + label: UMIs artificial diversity eliminated + - type: ReadsDroppedInTagCorrection + tag: UMI + upper: 0.1 + middle: 0.05 + label: Reads dropped in UMI error correction and whitelist + - type: ReadsDroppedInTagFiltering + upper: 0.1 + middle: 0.05 + label: Reads dropped in tags filtering + - type: TagGroupsWithNoAssemblingFeature + upper: 0.2 + middle: 0.1 + label: Tag groups with no assembling feature + - type: TagGroupsWithMultipleConsensuses + upper: 0.1 + middle: 0.01 + label: Barcode collisions in clonotype assembly + - type: UnassignedAlignments + upper: 0.1 + middle: 0.05 + label: Unassigned alignments in clonotype assembly + - type: ReadsUsedInClonotypes + upper: 0.9 + middle: 0.7 + label: Reads used in clonotypes + - type: AlignmentsDroppedLowQuality + upper: 0.05 + middle: 0.01 + label: Alignments dropped due to low sequence quality + - type: ClonesDroppedInPostFiltering + upper: 0.05 + middle: 0.01 + label: Clones dropped in post-filtering + - type: AlignmentsDroppedInPostFiltering + upper: 0.05 + middle: 0.01 + label: Alignments dropped in clones post-filtering + - type: OverlappedReadsLessBetter + upper: 0.1 + middle: 0.05 + label: Overlapped paired-end reads diff --git a/regression/presets/analyze/qiagen-mouse-rna-tcr-umi-targeted-qiaseq.yaml b/regression/presets/analyze/qiagen-mouse-rna-tcr-umi-targeted-qiaseq.yaml new file mode 100644 index 000000000..ba22426e4 --- /dev/null +++ b/regression/presets/analyze/qiagen-mouse-rna-tcr-umi-targeted-qiaseq.yaml @@ -0,0 +1,477 @@ +flags: [] +pipeline: + - align + - refineTagsAndSort + - assemble + - exportClones +align: + species: mmu + libraryName: default + trimmingQualityThreshold: 10 + trimmingWindowSize: 6 + chains: ALL + replaceWildcards: true + overlapPairedReads: true + bamDropNonVDJ: false + writeFailedAlignments: false + tagPattern: "^(R1:*) \\ ^(UMI:N{14:18})atgcattggagtcct" + tagUnstranded: false + tagMaxBudget: 10.0 + headerExtractors: [] + readIdAsCellTag: false + sampleTable: null + tagsValidations: + - type: MustContainTagType + tagType: Molecule + splitBySample: true + limit: null + parameters: + vParameters: + geneFeatureToAlign: VTranscriptWithP + minSumScore: 150 + relativeMinScore: 0.97 + maxHits: 5 + parameters: + type: kaligner2 + mapperNValue: 8 + mapperKValue: 1 + floatingLeftBound: false + floatingRightBound: true + mapperAbsoluteMinClusterScore: 102 + mapperExtraClusterScore: -38 + mapperMatchScore: 98 + mapperMismatchScore: -14 + mapperOffsetShiftScore: -82 + mapperSlotCount: 8 + mapperMaxClusters: 11 + mapperMaxClusterIndels: 4 + mapperKMersPerPosition: 4 + mapperAbsoluteMinScore: 103 + mapperRelativeMinScore: 0.81 + mapperMinSeedsDistance: 8 + mapperMaxSeedsDistance: 8 + alignmentStopPenalty: 0 + absoluteMinScore: 150 + relativeMinScore: 0.8 + maxHits: 2 + scoring: + type: affine + alphabet: nucleotide + subsMatrix: "simple(match = 10, mismatch = -19)" + gapOpenPenalty: -62 + gapExtensionPenalty: -11 + dParameters: + geneFeatureToAlign: DRegionWithP + relativeMinScore: 0.85 + absoluteMinScore: 25.0 + maxHits: 3 + scoring: + type: linear + alphabet: nucleotide + subsMatrix: "simple(match = 5, mismatch = -9)" + gapPenalty: -12 + jParameters: + geneFeatureToAlign: JRegionWithP + minSumScore: 150 + relativeMinScore: 0.97 + maxHits: 5 + parameters: + type: kaligner2 + mapperNValue: 8 + mapperKValue: 1 + floatingLeftBound: true + floatingRightBound: false + mapperAbsoluteMinClusterScore: 102 + mapperExtraClusterScore: -38 + mapperMatchScore: 95 + mapperMismatchScore: -14 + mapperOffsetShiftScore: -82 + mapperSlotCount: 6 + mapperMaxClusters: 4 + mapperMaxClusterIndels: 4 + mapperKMersPerPosition: 4 + mapperAbsoluteMinScore: 100 + mapperRelativeMinScore: 0.8 + mapperMinSeedsDistance: 5 + mapperMaxSeedsDistance: 5 + alignmentStopPenalty: 0 + absoluteMinScore: 150 + relativeMinScore: 0.8 + maxHits: 3 + scoring: + type: affine + alphabet: nucleotide + subsMatrix: "simple(match = 10, mismatch = -19)" + gapOpenPenalty: -62 + gapExtensionPenalty: -11 + cParameters: + geneFeatureToAlign: CExon1 + minSumScore: 40 + relativeMinScore: 0.97 + maxHits: 5 + parameters: + type: kaligner + mapperKValue: 5 + floatingLeftBound: false + floatingRightBound: true + mapperAbsoluteMinScore: 1.5 + mapperRelativeMinScore: 0.75 + mapperMatchScore: 1.0 + mapperMismatchPenalty: -0.1 + mapperOffsetShiftPenalty: -0.3 + mapperMinSeedsDistance: 4 + mapperMaxSeedsDistance: 10 + minAlignmentLength: 15 + maxAdjacentIndels: 2 + alignmentStopPenalty: -1000 + absoluteMinScore: 40.0 + relativeMinScore: 0.87 + maxHits: 7 + scoring: + type: linear + alphabet: nucleotide + subsMatrix: "simple(match = 5, mismatch = -9)" + gapPenalty: -12 + vjAlignmentOrder: VThenJ + libraryStructure: Unknown + includeDScore: false + includeCScore: false + minSumScore: 120.0 + relativeMinVFR3CDR3Score: 0.7 + allowPartialAlignments: false + allowNoCDR3PartAlignments: false + allowChimeras: false + readsLayout: ReverseOnly + mergerParameters: + qualityMergingAlgorithm: MaxSubtraction + partsLayout: null + minimalOverlap: 13 + minimalMatchQualitySum: 364 + maxQuality: 50 + minimalIdentity: 0.7 + identityType: MinimalQualityWeighted + fixSeed: true + alignmentBoundaryTolerance: 5 + minChimeraDetectionScore: 120 + vjOverlapWindow: 3 + saveOriginalSequence: false + saveOriginalReads: false + smartForceEdgeAlignments: true + retriesForPartialAlignments: 3 + retriesAlignmentOfNotCoveredPart: 2 +refineTagsAndSort: + whitelists: {} + runCorrection: true + parameters: + correctionPower: 0.001 + backgroundSubstitutionRate: 0.001 + backgroundIndelRate: 1.0E-5 + minQuality: 12 + maxSubstitutions: 2 + maxIndels: 2 + maxTotalErrors: 3 + postFilter: + type: filter_groups + groupingKeys: + - allTags:Molecule + predicates: + - metrics: + - type: group_metric_sum_weight + reportHist: + log: true + binNumber: 0 + minBinWidth: 0.2 + multiplyWeightByKey: false + operator: + type: group_operator_lowest_threshold + operators: + - type: group_operator_advanced_thresholding + algo: + type: otsu + histTransformers: + - type: inflate + fraction: 0.15 + value: 1.0 + logX: true + minimalSample: 20 + fallbackThreshold: 1.0 + accept: High + - type: group_operator_cumtop + share: 0.85 + round: Down + accept: High + accept: High + expectedSorting: [] + requiredSequences: [] + expectedSorting: + - allTags:Molecule + requiredSequences: [] +assemblePartial: + overlappedOnly: false + dropPartial: false + cellLevel: false + parameters: + kValue: 12 + kOffset: -7 + minimalAssembleOverlap: 12 + minimalNOverlap: 7 + minimalNOverlapShare: 0.65 + minimalAlignmentMergeIdentity: 0.85 + mergerParameters: + qualityMergingAlgorithm: SumSubtraction + partsLayout: CollinearDirect + minimalOverlap: 20 + minimalMatchQualitySum: 0 + maxQuality: 45 + minimalIdentity: 0.95 + identityType: Unweighted + maxLeftParts: 256000 + maxLeftMatches: 6144 +extend: + vAnchor: CDR3Begin + jAnchor: FR4Begin + minimalVScore: 50 + minimalJScore: 50 +assemble: + sortBySequence: false + clnaOutput: false + cellLevel: false + consensusAssemblerParameters: + assembler: + aAssemblerParameters: + bandWidth: 4 + scoring: + type: linear + alphabet: nucleotide + subsMatrix: "simple(match = 5, mismatch = -4)" + gapPenalty: -14 + minAlignmentScore: 40 + maxNormalizedAlignmentPenalty: 0.2 + trimMinimalSumQuality: 0 + trimReferenceRegion: false + maxQuality: 45 + maxIterations: 6 + minAltSeedQualityScore: 11 + minAltSeedNormalizedPenalty: 0.35 + altSeedPenaltyTolerance: 0.5 + minRecordSharePerConsensus: 0.2 + minRecordsPerConsensus: 0 + minRecursiveRecordShare: 0.4 + minQualityScore: 0 + maxConsensuses: 3 + minTagSuffixShare: 0.0 + isolateChains: false + empiricalAlignmentAssignment: + uniqGeneOccurrenceWithinTag: + - V + - J + uniqGeneOccurrenceWithinCluster: + - V + - J + matchAlignmentsWithoutVJHitsByTags: true + cloneAssemblerParameters: + assemblingFeatures: + - CDR3 + minimalClonalSequenceLength: 12 + qualityAggregationType: BetaScore + cloneClusteringParameters: + searchDepth: 2 + allowedMutationsInNRegions: 1 + searchParameters: twoMismatchesOrIndels + clusteringFilter: + type: advanced + correctionPower: 0.001 + backgroundSubstitutionRate: 5.0E-4 + backgroundIndelRate: 2.0E-4 + cloneFactoryParameters: + vParameters: + maxAlignmentWidthLinear: 5 + maxAlignmentWidthAffine: 500 + jParameters: + maxAlignmentWidthLinear: 5 + maxAlignmentWidthAffine: 500 + cParameters: + maxAlignmentWidthLinear: 5 + maxAlignmentWidthAffine: 500 + dParameters: + relativeMinScore: null + absoluteMinScore: null + maxHits: null + scoring: null + separateByV: false + separateByJ: false + separateByC: false + maximalPreClusteringRatio: 1.0 + preClusteringScoreFilteringRatio: 2.0 + preClusteringCountFilteringRatio: 2.0 + addReadsCountOnClustering: false + badQualityThreshold: 15 + maxBadPointsPercent: 0.7 + mappingThreshold: 2of5 + minimalQuality: 20 + postFilters: null + inferMinRecordsPerConsensus: true +assembleContigs: + ignoreTags: false + parameters: + branchingMinimalQualityShare: 0.1 + branchingMinimalSumQuality: 60 + decisiveBranchingSumQualityThreshold: 120 + alignedSequenceEdgeDelta: 3 + alignmentEdgeRegionSize: 7 + minimalNonEdgePointsFraction: 0.25 + minimalMeanNormalizedQuality: 5.0 + outputMinimalQualityShare: 0.75 + outputMinimalSumQuality: 0 + subCloningRegions: null + assemblingRegions: null + postFiltering: + type: NoFiltering + trimmingParameters: + averageQualityThreshold: 10.0 + windowSize: 8 + minimalContigLength: 20 + alignedRegionsOnly: false + discardAmbiguousNucleotideCalls: false + cloneFactoryParameters: + vParameters: + relativeMinScore: 0.97 + maxAlignmentWidthLinear: 5 + maxAlignmentWidthAffine: 500 + jParameters: + maxAlignmentWidthLinear: 5 + maxAlignmentWidthAffine: 500 + cParameters: + maxAlignmentWidthLinear: 5 + maxAlignmentWidthAffine: 500 + dParameters: + relativeMinScore: null + absoluteMinScore: null + maxHits: null + scoring: null +exportAlignments: + chains: ALL + noHeader: false + fields: + - field: -readIds + - field: -tags + args: + - Molecule + - field: -targetSequences + - field: -targetQualities + - field: -vHitsWithScore + - field: -dHitsWithScore + - field: -jHitsWithScore + - field: -cHitsWithScore + - field: -vAlignments + - field: -dAlignments + - field: -jAlignments + - field: -cAlignments + - field: -allNFeaturesWithMinQuality + - field: -allAAFeatures + - field: -defaultAnchorPoints + - field: -topChains +exportClones: + splitByTagType: null + filterOutOfFrames: false + filterStops: false + chains: ALL + noHeader: false + fields: + - field: -cloneId + - field: -readCount + - field: -readFraction + - field: -uniqueTagCount + args: + - Molecule + - field: -uniqueTagFraction + args: + - Molecule + - field: -targetSequences + - field: -targetQualities + - field: -vHitsWithScore + - field: -dHitsWithScore + - field: -jHitsWithScore + - field: -cHitsWithScore + - field: -vAlignments + - field: -dAlignments + - field: -jAlignments + - field: -cAlignments + - field: -allNFeaturesWithMinQuality + - field: -allAAFeatures + - field: -defaultAnchorPoints + splitFilesBy: + - chain + groupClonesBy: [] +qc: + checks: + - type: SuccessfullyAlignedReads + upper: 0.85 + middle: 0.7 + label: Successfully aligned reads + - type: OffTargetReads + upper: 0.2 + middle: 0.1 + label: Off target (non TCR/IG) reads + - type: ReadsWithNoVOrJHits + upper: 0.2 + middle: 0.1 + label: Reads with no V or J hits + - type: AlignmentsWithNoAssemblingFeature + upper: 0.15 + middle: 0.05 + label: Alignments without assembling feature + - type: ReadsWithNoBarcode + upper: 0.2 + middle: 0.1 + label: Reads with no barcode + - type: ReadsDroppedInTagRefinement + upper: 0.1 + middle: 0.05 + label: Reads dropped in tags error correction and filtering + - type: TagArtificialDiversityEliminated + tag: UMI + upper: 0.5 + middle: 0.3 + label: UMIs artificial diversity eliminated + - type: ReadsDroppedInTagCorrection + tag: UMI + upper: 0.1 + middle: 0.05 + label: Reads dropped in UMI error correction and whitelist + - type: ReadsDroppedInTagFiltering + upper: 0.1 + middle: 0.05 + label: Reads dropped in tags filtering + - type: TagGroupsWithNoAssemblingFeature + upper: 0.2 + middle: 0.1 + label: Tag groups with no assembling feature + - type: TagGroupsWithMultipleConsensuses + upper: 0.1 + middle: 0.01 + label: Barcode collisions in clonotype assembly + - type: UnassignedAlignments + upper: 0.1 + middle: 0.05 + label: Unassigned alignments in clonotype assembly + - type: ReadsUsedInClonotypes + upper: 0.9 + middle: 0.7 + label: Reads used in clonotypes + - type: AlignmentsDroppedLowQuality + upper: 0.05 + middle: 0.01 + label: Alignments dropped due to low sequence quality + - type: ClonesDroppedInPostFiltering + upper: 0.05 + middle: 0.01 + label: Clones dropped in post-filtering + - type: AlignmentsDroppedInPostFiltering + upper: 0.05 + middle: 0.01 + label: Alignments dropped in clones post-filtering + - type: OverlappedReadsLessBetter + upper: 0.1 + middle: 0.05 + label: Overlapped paired-end reads diff --git a/regression/presets/list.txt b/regression/presets/list.txt index 33982fee3..1881cf814 100644 --- a/regression/presets/list.txt +++ b/regression/presets/list.txt @@ -302,8 +302,12 @@ bd-human-sc-xcr-rhapsody-cdr3 (Rhapsody Human VDJ CDR3 Assay) -----QIAGEN----- qiagen-human-rna-tcr-umi-qiaseq (QIAseq Human TCR Panel Immune Repertoire RNA Library Kit) +qiagen-mouse-rna-tcr-umi-targeted-qiaseq (QIAseq Targeted RNA Panel Mouse TCR) + qiagen-mouse-rna-tcr-umi-qiaseq (QIAseq Mouse TCR Panel Immune Repertoire RNA Library Kit) +qiagen-human-rna-tcr-umi-targeted-qiaseq (QIAseq Targeted RNA Panel Human TCR) + -----New England BioLabs----- neb-human-rna-xcr-umi-nebnext (NEBNext® Immune Sequencing Kit (Human) BCR & TCR)