Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- Thu Aug 5 10:25:13 CDT 2021
- /tassel-5-standalone/lib/ahocorasick-0.2.4.jar:/tassel-5-standalone/lib/biojava-alignment-4.0.0.jar:/tassel-5-standalone/lib/biojava-core-4.0.0.jar:/tassel-5-standalone/lib/biojava-phylo-4.0.0.jar:/tassel-5-standalone/lib/colt-1.2.0.jar:/tassel-5-standalone/lib/commons-codec-1.10.jar:/tassel-5-standalone/lib/commons-math3-3.4.1.jar:/tassel-5-standalone/lib/ejml-0.23.jar:/tassel-5-standalone/lib/fastutil-8.2.2.jar:/tassel-5-standalone/lib/forester-1.038.jar:/tassel-5-standalone/lib/gs-core-1.3.jar:/tassel-5-standalone/lib/gs-ui-1.3.jar:/tassel-5-standalone/lib/guava-22.0.jar:/tassel-5-standalone/lib/htsjdk-2.23.0.jar:/tassel-5-standalone/lib/itextpdf-5.1.0.jar:/tassel-5-standalone/lib/javax.json-1.0.4.jar:/tassel-5-standalone/lib/jcommon-1.0.23.jar:/tassel-5-standalone/lib/jfreechart-1.0.19.jar:/tassel-5-standalone/lib/jfreesvg-3.2.jar:/tassel-5-standalone/lib/jhdf5-14.12.5.jar:/tassel-5-standalone/lib/json-simple-1.1.1.jar:/tassel-5-standalone/lib/junit-4.10.jar:/tassel-5-standalone/lib/kotlin-stdlib-1.3.50.jar:/tassel-5-standalone/lib/kotlin-stdlib-jdk7-1.3.50.jar:/tassel-5-standalone/lib/kotlin-stdlib-jdk8-1.3.50.jar:/tassel-5-standalone/lib/kotlinx-coroutines-core-1.3.0.jar:/tassel-5-standalone/lib/log4j-1.2.13.jar:/tassel-5-standalone/lib/mail-1.4.jar:/tassel-5-standalone/lib/phg.jar:/tassel-5-standalone/lib/postgresql-9.4-1201.jdbc41.jar:/tassel-5-standalone/lib/scala-library-2.10.1.jar:/tassel-5-standalone/lib/slf4j-api-1.7.10.jar:/tassel-5-standalone/lib/slf4j-simple-1.7.10.jar:/tassel-5-standalone/lib/snappy-java-1.1.1.6.jar:/tassel-5-standalone/lib/sqlite-jdbc-3.8.5-pre1.jar:/tassel-5-standalone/lib/trove-3.0.3.jar:/tassel-5-standalone/sTASSEL.jar
- Memory Settings: -Xms512m -Xmx215040m
- Tassel Pipeline Arguments: -debug -configParameters /phg/config_ImputePipelinePlugin_fastq.txt -ImputePipelinePlugin -imputeTarget pathToVCF -endPlugin
- [main] INFO net.maizegenetics.plugindef.ParameterCache - load: loading parameter cache with: /phg/config_ImputePipelinePlugin_fastq.txt
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: DBtype value: sqlite
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: maxSecondary value: 20
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: host value: localHost
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: minTransitionProb value: 0.001
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: DB value: /phg/95HRSW_attempt4.db
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: inputType value: fastq
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: pangenomeDir value: /phg/outputDir/pangenome
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: outputSecondaryStats value: false
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: maxNodes value: 1000
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: samDir value: /phg/inputDir/imputation/sam/
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: fParameter value: f1000,5000
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: minimapLocation value: minimap2
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: numThreads value: 3
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: pathMethod value: pathMethodwithnewdb_fq_GATK_PIPELINE
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: probCorrect value: 0.99
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: indexNumberBases value: 90G
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: pangenomeHaplotypeMethod value: GATK_PIPELINE
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: minTaxa value: 1
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: indexKmerLength value: 21
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: usebf value: false
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: removeEqual value: true
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: pathHaplotypeMethod value: GATK_PIPELINE
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: minReads value: 1
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: keyFile value: /phg/readMapping_key_file_17C23-2fq.txt
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: user value: sqlite
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: lowMemMode value: false
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: password value: sqlite
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: splitNodes value: true
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: outVcfFile value: 95HRSW_imputedwith17C23-2fastq_VCF
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: maxReads value: 10000
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: maxRefRangeErr value: 0.25
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: indexWindowSize value: 11
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: readMethodDescription value: readMethod+newdb_filetype
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: minP value: 0.8
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: splitProb value: 0.99
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: fastqDir value: /phg/inputDir/imputation/fastq/
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: pathMethodDescription value: pathMethod+newdb_filetype_HaplotypeMethod
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: readMethod value: GATK_PIPELINE
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: configFile value: /phg/config_ImputePipelinePlugin_fastq.txt
- [main] INFO net.maizegenetics.plugindef.ParameterCache - ParameterCache: key: liquibaseOutdir value: /phg/outputDir
- [main] INFO net.maizegenetics.tassel.TasselLogging - Tassel Version: 5.2.71 Date: March 26, 2021
- [main] INFO net.maizegenetics.tassel.TasselLogging - Max Available Memory Reported by JVM: 191147 MB
- [main] INFO net.maizegenetics.tassel.TasselLogging - Java Version: 1.8.0_242
- [main] INFO net.maizegenetics.tassel.TasselLogging - OS: Linux
- [main] INFO net.maizegenetics.tassel.TasselLogging - Number of Processors: 72
- [main] INFO net.maizegenetics.pipeline.TasselPipeline - Tassel Pipeline Arguments: [-fork1, -ImputePipelinePlugin, -imputeTarget, pathToVCF, -endPlugin, -runfork1]
- net.maizegenetics.pangenome.pipeline.ImputePipelinePlugin
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin - Starting net.maizegenetics.pangenome.pipeline.ImputePipelinePlugin: time: Aug 5, 2021 15:25:32
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin -
- ImputePipelinePlugin Parameters
- imputeTarget: pathToVCF
- inputType: fastq
- configFile: /phg/config_ImputePipelinePlugin_fastq.txt
- pangenomeHaplotypeMethod: GATK_PIPELINE
- pathHaplotypeMethod: GATK_PIPELINE
- pangenomeDir: /phg/outputDir/pangenome
- pangenomeIndexName: null
- indexKmerLength: 21
- indexWindowSize: 11
- indexNumberBases: 90G
- minimapLocation: minimap2
- readMethod: GATK_PIPELINE
- readMethodDescription: readMethod+newdb_filetype
- outVcfFile: 95HRSW_imputedwith17C23-2fastq_VCF
- forceDBUpdate: false
- liquibaseOutdir: /phg/outputDir
- skipLiquibaseCheck: false
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.pipeline.ImputePipelinePlugin - Checking if Liquibase can be run.
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin - Starting net.maizegenetics.pangenome.liquibase.CheckDBVersionPlugin: time: Aug 5, 2021 15:25:32
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin -
- CheckDBVersionPlugin Parameters
- outputDir: /phg/outputDir
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.liquibase.CheckDBVersionPlugin - Deleting yesFile /phg/outputDir/run_yes.txt if it exists
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.liquibase.CheckDBVersionPlugin - Deleting noFile /phg/outputDirrun_no.txt if it exists
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.DBLoadingUtils - first connection: dbName from config file = /phg/95HRSW_attempt4.db host: localHost user: sqlite type: sqlite
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.DBLoadingUtils - Database URL: jdbc:sqlite:/phg/95HRSW_attempt4.db
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.DBLoadingUtils - Connected to database:
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.liquibase.CheckDBVersionPlugin - queueHaplotypeNodesByRange: query: select name FROM sqlite_master where type='table' and name='variants';
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin - Finished net.maizegenetics.pangenome.liquibase.CheckDBVersionPlugin: time: Aug 5, 2021 15:25:32
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin - net.maizegenetics.pangenome.liquibase.CheckDBVersionPlugin Citation: Bradbury PJ, Zhang Z, Kroon DE, Casstevens TM, Ramdoss Y, Buckler ES. (2007) TASSEL: Software for association mapping of complex traits in diverse samples. Bioinformatics 23:2633-2635.
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin - Starting net.maizegenetics.pangenome.liquibase.LiquibaseUpdatePlugin: time: Aug 5, 2021 15:25:32
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin -
- LiquibaseUpdatePlugin Parameters
- outputDir: /phg/outputDir
- command: status
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.liquibase.LiquibaseUpdatePlugin - Please wait, begin Command:liquibase --driver=org.sqlite.JDBC --url=jdbc:sqlite:/phg/95HRSW_attempt4.db --username=sqlite --password=sqlite --changeLogFile=/liquibase/changelogs/db.changelog-master.xml status --verbose
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin - Finished net.maizegenetics.pangenome.liquibase.LiquibaseUpdatePlugin: time: Aug 5, 2021 15:25:34
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.pipeline.ImputePipelinePlugin - PHG DB is up to date. Proceeding with Populating the PHG DB.
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin - Starting net.maizegenetics.pangenome.api.HaplotypeGraphBuilderPlugin: time: Aug 5, 2021 15:25:34
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin -
- HaplotypeGraphBuilderPlugin Parameters
- configFile: /phg/config_ImputePipelinePlugin_fastq.txt
- methods: GATK_PIPELINE
- includeSequences: true
- includeVariantContexts: false
- haplotypeIds: null
- chromosomes: null
- taxa: null
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.DBLoadingUtils - first connection: dbName from config file = /phg/95HRSW_attempt4.db host: localHost user: sqlite type: sqlite
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.DBLoadingUtils - Database URL: jdbc:sqlite:/phg/95HRSW_attempt4.db
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.DBLoadingUtils - Connected to database:
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - referenceRangesAsMap: query statement: select reference_ranges.ref_range_id, chrom, range_start, range_end, methods.name from reference_ranges INNER JOIN ref_range_ref_range_method on ref_range_ref_range_method.ref_range_id=reference_ranges.ref_range_id INNER JOIN methods on ref_range_ref_range_method.method_id = methods.method_id AND methods.method_type = 7 ORDER BY reference_ranges.ref_range_id
- methods size: 1
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - referenceRangesAsMap: number of reference ranges: 94229
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - referenceRangesAsMap: time: 0.624981529 secs.
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - taxaListMap: query statement: SELECT gamete_haplotypes.gamete_grp_id, genotypes.line_name FROM gamete_haplotypes INNER JOIN gametes ON gamete_haplotypes.gameteid = gametes.gameteid INNER JOIN genotypes on gametes.genoid = genotypes.genoid ORDER BY gamete_haplotypes.gamete_grp_id;
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - taxaListMap: number of taxa lists: 136337
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - taxaListMap: time: 6.726373845 secs.
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - createHaplotypeNodes: haplotype method: GATK_PIPELINE range group method: null
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - createHaplotypeNodes: query statement: SELECT haplotypes_id, gamete_grp_id, haplotypes.ref_range_id, asm_contig, asm_start_coordinate, asm_end_coordinate, genome_file_id, sequence, seq_hash, seq_len FROM haplotypes WHERE method_id = 3;
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - addNodes: number of nodes: 8857526
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - addNodes: number of reference ranges: 94229
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - createHaplotypeNodes: time: 440.934323876 secs.
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.HaplotypeGraph - Created graph edges: created when requested number of nodes: 8857526 number of reference ranges: 94229
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin - Finished net.maizegenetics.pangenome.api.HaplotypeGraphBuilderPlugin: time: Aug 5, 2021 15:33:9
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin - Starting net.maizegenetics.pangenome.hapCalling.FastqToMappingPlugin: time: Aug 5, 2021 15:33:9
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin -
- FastqToMappingPlugin Parameters
- minimap2IndexFile: /phg/outputDir/pangenome/pangenome_GATK_PIPELINE_k21w11I90G.mmi
- keyFile: /phg/readMapping_key_file_17C23-2fq.txt
- fastqDir: /phg/inputDir/imputation/fastq/
- maxRefRangeErr: 0.25
- lowMemMode: false
- maxSecondary: 20
- fParameter: f1000,5000
- minimapLocation: minimap2
- methodName: GATK_PIPELINE
- methodDescription: readMethod+newdb_filetype
- debugDir:
- outputSecondaryStats: false
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.DBLoadingUtils - first connection: dbName from config file = /phg/95HRSW_attempt4.db host: localHost user: sqlite type: sqlite
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.DBLoadingUtils - Database URL: jdbc:sqlite:/phg/95HRSW_attempt4.db
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.DBLoadingUtils - Connected to database:
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess - PHGdbAccess - db is setup, init prepared statements, load hash table
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess -
- beginning - isSqlite is true
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess - before loading hash, size of all geneotypes in genotype table=95
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess - refRangeRefRangeIDMap is null, creating new one with size : 94229
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess - loadAnchorHash: at end, size of refRangeRefRangeIDMap: 94229, number of rs.next processed: 94229
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess - before loading hash, size of all methods in method table=10
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess - before loading hash, size of all groups in taxa_groups table=0
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess - before loading hash, size of all groups in gamete_groups table=136337
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess - before loading hash, size of all gametes in gametes table=95
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess - before loading readMappingHash, size of all read_mappings in read_mapping table=1
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.hapCalling.Minimap2Utils - Skipping Keyfile entry: cultivar 17C23-2, flowcell_lane wgs has already been processed and loaded into the DB.
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess - Closing DB
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin - Finished net.maizegenetics.pangenome.hapCalling.FastqToMappingPlugin: time: Aug 5, 2021 15:33:13
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin - Starting net.maizegenetics.pangenome.api.HaplotypeGraphBuilderPlugin: time: Aug 5, 2021 15:33:13
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin -
- HaplotypeGraphBuilderPlugin Parameters
- configFile: /phg/config_ImputePipelinePlugin_fastq.txt
- methods: GATK_PIPELINE
- includeSequences: false
- includeVariantContexts: false
- haplotypeIds: null
- chromosomes: null
- taxa: null
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.DBLoadingUtils - first connection: dbName from config file = /phg/95HRSW_attempt4.db host: localHost user: sqlite type: sqlite
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.DBLoadingUtils - Database URL: jdbc:sqlite:/phg/95HRSW_attempt4.db
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.DBLoadingUtils - Connected to database:
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - referenceRangesAsMap: query statement: select reference_ranges.ref_range_id, chrom, range_start, range_end, methods.name from reference_ranges INNER JOIN ref_range_ref_range_method on ref_range_ref_range_method.ref_range_id=reference_ranges.ref_range_id INNER JOIN methods on ref_range_ref_range_method.method_id = methods.method_id AND methods.method_type = 7 ORDER BY reference_ranges.ref_range_id
- methods size: 1
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - referenceRangesAsMap: number of reference ranges: 94229
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - referenceRangesAsMap: time: 0.421060762 secs.
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - taxaListMap: query statement: SELECT gamete_haplotypes.gamete_grp_id, genotypes.line_name FROM gamete_haplotypes INNER JOIN gametes ON gamete_haplotypes.gameteid = gametes.gameteid INNER JOIN genotypes on gametes.genoid = genotypes.genoid ORDER BY gamete_haplotypes.gamete_grp_id;
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - taxaListMap: number of taxa lists: 136337
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - taxaListMap: time: 4.490692238 secs.
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - createHaplotypeNodes: haplotype method: GATK_PIPELINE range group method: null
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - createHaplotypeNodes: query statement: SELECT haplotypes_id, gamete_grp_id, haplotypes.ref_range_id, asm_contig, asm_start_coordinate, asm_end_coordinate, genome_file_id, seq_hash, seq_len FROM haplotypes WHERE method_id = 3;
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - addNodes: number of nodes: 8857526
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - addNodes: number of reference ranges: 94229
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - createHaplotypeNodes: time: 117.803670678 secs.
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.HaplotypeGraph - Created graph edges: created when requested number of nodes: 8857526 number of reference ranges: 94229
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin - Finished net.maizegenetics.pangenome.api.HaplotypeGraphBuilderPlugin: time: Aug 5, 2021 15:35:20
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin - Starting net.maizegenetics.pangenome.hapCalling.BestHaplotypePathPlugin: time: Aug 5, 2021 15:35:20
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin -
- BestHaplotypePathPlugin Parameters
- keyFile: /phg/readMapping_key_file_17C23-2fq_pathKeyFile.txt
- readFile: null
- readDir: null
- outDir: null
- readMethod: GATK_PIPELINE
- pathMethod: pathMethodwithnewdb_fq_GATK_PIPELINE
- pathMethodDescription: pathMethod+newdb_filetype_HaplotypeMethod
- overwrite: false
- minTaxa: 1
- minReads: 1
- maxReads: 10000
- maxNodes: 1000
- minTransitionProb: 0.001
- probCorrect: 0.99
- splitNodes: true
- splitProb: 0.99
- usebf: false
- minP: 0.8
- bfInfoFile: null
- removeEqual: true
- numThreads: 3
- requiredTaxa: null
- algorithmType: classic
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.hapCalling.BestHaplotypePathPlugin - Classic filter found 0 ranges missing one of the required taxa, 0 ranges with too many nodes, and 0 with too few taxa}
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.hapCalling.BestHaplotypePathPlugin - Classic filter removed 0 ranges
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - createEdges: creating edges from nodes.
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - createEdges: time: 6.264702469 secs.
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.HaplotypeGraph - Created graph number of edges: 8855552 number of nodes: 8857526 number of reference ranges: 94229
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - createEdgesFullyConnected: creating edges from nodes.
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - createEdgesFullyConnected: time: 1473.881245856 secs.
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.HaplotypeGraph - Created graph number of edges: 832421888 number of nodes: 8857526 number of reference ranges: 94229
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.DBLoadingUtils - first connection: dbName from config file = /phg/95HRSW_attempt4.db host: localHost user: sqlite type: sqlite
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.DBLoadingUtils - Database URL: jdbc:sqlite:/phg/95HRSW_attempt4.db
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.DBLoadingUtils - Connected to database:
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess - PHGdbAccess - db is setup, init prepared statements, load hash table
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess -
- beginning - isSqlite is true
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess - before loading hash, size of all geneotypes in genotype table=95
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess - refRangeRefRangeIDMap is null, creating new one with size : 94229
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess - loadAnchorHash: at end, size of refRangeRefRangeIDMap: 94229, number of rs.next processed: 94229
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess - before loading hash, size of all methods in method table=10
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess - before loading hash, size of all groups in taxa_groups table=0
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess - before loading hash, size of all groups in gamete_groups table=136337
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess - before loading hash, size of all gametes in gametes table=95
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess - before loading readMappingHash, size of all read_mappings in read_mapping table=1
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.DBLoadingUtils - first connection: dbName from config file = /phg/95HRSW_attempt4.db host: localHost user: sqlite type: sqlite
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.DBLoadingUtils - Database URL: jdbc:sqlite:/phg/95HRSW_attempt4.db
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.DBLoadingUtils - Connected to database:
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess - PHGdbAccess - db is setup, init prepared statements, load hash table
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess -
- beginning - isSqlite is true
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess - before loading hash, size of all geneotypes in genotype table=95
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess - refRangeRefRangeIDMap is null, creating new one with size : 94229
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess - loadAnchorHash: at end, size of refRangeRefRangeIDMap: 94229, number of rs.next processed: 94229
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess - before loading hash, size of all methods in method table=10
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess - before loading hash, size of all groups in taxa_groups table=0
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess - before loading hash, size of all groups in gamete_groups table=136337
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess - before loading hash, size of all gametes in gametes table=95
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess - before loading readMappingHash, size of all read_mappings in read_mapping table=1
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess - Closing DB
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.hapCalling.BestHaplotypePathPlugin - The pathKeyFile has 1 rows. Paths will be found for 0. The others already have paths in the db.
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.PHGdbAccess - Closing DB
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin - Finished net.maizegenetics.pangenome.hapCalling.BestHaplotypePathPlugin: time: Aug 5, 2021 16:31:23
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin - Starting net.maizegenetics.pangenome.api.HaplotypeGraphBuilderPlugin: time: Aug 5, 2021 16:31:23
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin -
- HaplotypeGraphBuilderPlugin Parameters
- configFile: /phg/config_ImputePipelinePlugin_fastq.txt
- methods: GATK_PIPELINE
- includeSequences: false
- includeVariantContexts: true
- haplotypeIds: null
- chromosomes: null
- taxa: null
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.DBLoadingUtils - first connection: dbName from config file = /phg/95HRSW_attempt4.db host: localHost user: sqlite type: sqlite
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.DBLoadingUtils - Database URL: jdbc:sqlite:/phg/95HRSW_attempt4.db
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.DBLoadingUtils - Connected to database:
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - referenceRangesAsMap: query statement: select reference_ranges.ref_range_id, chrom, range_start, range_end, methods.name from reference_ranges INNER JOIN ref_range_ref_range_method on ref_range_ref_range_method.ref_range_id=reference_ranges.ref_range_id INNER JOIN methods on ref_range_ref_range_method.method_id = methods.method_id AND methods.method_type = 7 ORDER BY reference_ranges.ref_range_id
- methods size: 1
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - referenceRangesAsMap: number of reference ranges: 94229
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - referenceRangesAsMap: time: 0.403850626 secs.
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - taxaListMap: query statement: SELECT gamete_haplotypes.gamete_grp_id, genotypes.line_name FROM gamete_haplotypes INNER JOIN gametes ON gamete_haplotypes.gameteid = gametes.gameteid INNER JOIN genotypes on gametes.genoid = genotypes.genoid ORDER BY gamete_haplotypes.gamete_grp_id;
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - taxaListMap: number of taxa lists: 136337
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - taxaListMap: time: 4.734796851 secs.
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.VariantUtils - variantIdsToVariantMap: query statement: SELECT variant_id, chrom, position, ref_allele_id, alt_allele_id FROM variants;
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - createHaplotypeNodes: haplotype method: GATK_PIPELINE range group method: null
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - createHaplotypeNodes: query statement: SELECT haplotypes_id, gamete_grp_id, haplotypes.ref_range_id, asm_contig, asm_start_coordinate, asm_end_coordinate, genome_file_id, seq_hash, seq_len, variant_list FROM haplotypes WHERE method_id = 3;
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - addNodes: number of nodes: 8857526
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - addNodes: number of reference ranges: 94229
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.CreateGraphUtils - createHaplotypeNodes: time: 202.517586763 secs.
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.api.HaplotypeGraph - Created graph edges: created when requested number of nodes: 8857526 number of reference ranges: 94229
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin - Finished net.maizegenetics.pangenome.api.HaplotypeGraphBuilderPlugin: time: Aug 5, 2021 16:34:56
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin - Starting net.maizegenetics.pangenome.hapCalling.ImportDiploidPathPlugin: time: Aug 5, 2021 16:34:56
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin -
- ImportDiploidPathPlugin Parameters
- pathMethodName: pathMethodwithnewdb_fq_GATK_PIPELINE
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.DBLoadingUtils - first connection: dbName from config file = /phg/95HRSW_attempt4.db host: localHost user: sqlite type: sqlite
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.DBLoadingUtils - Database URL: jdbc:sqlite:/phg/95HRSW_attempt4.db
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.db_loading.DBLoadingUtils - Connected to database:
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.hapCalling.ImportDiploidPathPlugin - importPathsFromDB: query: SELECT line_name, paths_data FROM paths, genotypes, methods WHERE paths.genoid=genotypes.genoid AND methods.method_id=paths.method_id AND methods.name='pathMethodwithnewdb_fq_GATK_PIPELINE'
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.hapCalling.ImportDiploidPathPlugin - importPathsFromDB: number of path list: 1
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin - Finished net.maizegenetics.pangenome.hapCalling.ImportDiploidPathPlugin: time: Aug 5, 2021 16:34:57
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin - Starting net.maizegenetics.pangenome.hapCalling.PathsToVCFPlugin: time: Aug 5, 2021 16:34:57
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin -
- PathsToVCFPlugin Parameters
- outputFile: 95HRSW_imputedwith17C23-2fastq_VCF.vcf
- refRangeFileVCF: null
- referenceFasta: null
- makeDiploid: true
- positions: null
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.hapCalling.PathsToVCFPlugin - PathsToVCFPlugin: processData: number of ranges: 94229
- [pool-1-thread-1] INFO net.maizegenetics.pangenome.hapCalling.PathsToVCFPlugin - PathsToVCFPlugin: processData: number of taxa: 1
- [pool-1-thread-1] DEBUG net.maizegenetics.plugindef.AbstractPlugin - Parent job is Cancelling
- kotlinx.coroutines.JobCancellationException: Parent job is Cancelling; job=UndispatchedCoroutine{Cancelled}@7dbcd962
- Caused by: java.lang.IllegalStateException: Allele in genotype C not in the variant context [C, C]
- at htsjdk.variant.variantcontext.VariantContext$Validation.validateGenotypes(VariantContext.java:382)
- at htsjdk.variant.variantcontext.VariantContext$Validation.access$200(VariantContext.java:323)
- at htsjdk.variant.variantcontext.VariantContext$Validation$2.validate(VariantContext.java:331)
- at htsjdk.variant.variantcontext.VariantContext.lambda$validate$0(VariantContext.java:1384)
- at java.lang.Iterable.forEach(Iterable.java:75)
- at htsjdk.variant.variantcontext.VariantContext.validate(VariantContext.java:1384)
- at htsjdk.variant.variantcontext.VariantContext.<init>(VariantContext.java:489)
- at htsjdk.variant.variantcontext.VariantContextBuilder.make(VariantContextBuilder.java:647)
- at htsjdk.variant.variantcontext.VariantContextBuilder.make(VariantContextBuilder.java:638)
- at net.maizegenetics.pangenome.hapCalling.PathsToVCFPlugin.createVariantContext(PathsToVCFPlugin.kt:342)
- at net.maizegenetics.pangenome.hapCalling.PathsToVCFPlugin.variantContexts(PathsToVCFPlugin.kt:475)
- at net.maizegenetics.pangenome.hapCalling.PathsToVCFPlugin.access$variantContexts(PathsToVCFPlugin.kt:53)
- at net.maizegenetics.pangenome.hapCalling.PathsToVCFPlugin$infosByRange$2$invokeSuspend$$inlined$forEach$lambda$1.invokeSuspend(PathsToVCFPlugin.kt:216)
- at kotlin.coroutines.jvm.internal.BaseContinuationImpl.resumeWith(ContinuationImpl.kt:33)
- at kotlinx.coroutines.DispatchedTask.run(Dispatched.kt:241)
- at kotlinx.coroutines.scheduling.CoroutineScheduler.runSafely(CoroutineScheduler.kt:594)
- at kotlinx.coroutines.scheduling.CoroutineScheduler.access$runSafely(CoroutineScheduler.kt:60)
- at kotlinx.coroutines.scheduling.CoroutineScheduler$Worker.run(CoroutineScheduler.kt:740)
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin -
- Usage:
- PathsToVCFPlugin <options>
- -outputFile <Output VCF File Name> : Output file name (required)
- -refRangeFileVCF <Reference Range File> : Reference Range file used to subset the paths for only specified regions of the genome.
- -referenceFasta <Reference Genome> : Reference Genome.
- -makeDiploid <true | false> : Whether to report haploid paths as homozygousdiploid (Default: true)
- -positions <Position List> : Positions to include in VCF. Can be specified by Genotype file (i.e. VCF, Hapmap, etc.), bed file, or json file containing the requested positions.
- [pool-1-thread-1] ERROR net.maizegenetics.plugindef.AbstractPlugin - Parent job is Cancelling
- [pool-1-thread-1] INFO net.maizegenetics.plugindef.AbstractPlugin - Finished net.maizegenetics.pangenome.pipeline.ImputePipelinePlugin: time: Aug 5, 2021 16:34:58
- [pool-1-thread-1] INFO net.maizegenetics.pipeline.TasselPipeline - net.maizegenetics.pangenome.pipeline.ImputePipelinePlugin: time: Aug 5, 2021 16:34:58: progress: 100%
Add Comment
Please, Sign In to add comment