picard报错的问题
picard报错的问题
parallel --will-cite --jobs 2 --delay 0.5 'picard MarkDuplicates
INPUT={}
OUTPUT={.}.rmdup
METRICS_FILE={.}.dups.log
REMOVE_DUPLICATES=true
BARCODE_TAG=RG
TMP_DIR=/mnt/f/temp/ ’ ::: SubBarcodes.Name.Name.Name.Species.*.bam
#####失败报错 Name的问题
Name=SubBarcodes.sp.atac.first
Species=hg19
parallel --will-cite --jobs 2 --delay 0.5 'picard MarkDuplicates INPUT={} OUTPUT={.}.rmdup METRICS_FILE={.}.dups.log ’ ::: SubBarcodes.Name.Name.Name.Species.*.bam
###成功的话
picard MarkDuplicates INPUT=“SubBarcodes.sp.atac.first.hg19.0000.bam” OUTPUT=“test.bam” METRICS_FILE=“test.dups.log”
picard MarkDuplicates INPUT=“SubBarcodes.sp.atac.first.hg19.0000.bam” OUTPUT=“test2.bam” METRICS_FILE=“test2.dups.log” REMOVE_DUPLICATES=“true”
#####这个又成功
Name=SubBarcodes.sp.atac.first
Species=hg19
parallel --will-cite --jobs 2 --delay 0.5 'picard MarkDuplicates INPUT={} OUTPUT={.}.rmdup METRICS_FILE={.}.dups.log ’ ::: SubBarcodes.sp.atac.first.hg19.0000.bam
#####失败 BARCODE_TAG=RG 问题
Name=sp.atac.first
Species=hg19
parallel --will-cite --jobs 2 --delay 0.5 'picard MarkDuplicates INPUT={} OUTPUT={.}.rmdup METRICS_FILE={.}.dups.log REMOVE_DUPLICATES=true BARCODE_TAG=RG ’ ::: SubBarcodes.sp.atac.first.hg19.0000.bam
picard MarkDuplicates INPUT=“SubBarcodes.sp.atac.first.hg19.0000.bam” OUTPUT=“test2.bam” METRICS_FILE=“test2.dups.log” REMOVE_DUPLICATES=“true” BARCODE_TAG=“RG”
picard MarkDuplicates INPUT=“SubBarcodes.sp.atac.first.hg19.0000.bam” OUTPUT=“test2.bam” METRICS_FILE=“test2.dups.log” REMOVE_DUPLICATES=“true” BARCODE_TAG=RG
(rna) yyp@DESKTOP-LRUCLJJ:~/picard-master$ parallel --will-cite --jobs 2 --delay 0.5 'picard MarkDuplicates INPUT={} OUTPUT={.}.rmdup METRICS_FILE={.}.dups.log REMOVE_DUPLICATES=true VALIDATION_STRINGENCY=SILENT BARCODE_TAG=RG TMP_DIR='dir′/temp/′:::SubBarcodes.dir'/temp/ ' ::: SubBarcodes.dir′/temp/′:::SubBarcodes.Name.$Species.*.bam
INFO 2021-10-31 19:35:12 MarkDuplicates
********** NOTE: Picard’s command line syntax is changing.
********** For more information, please see:
********** https://github.com/broadinstitute/picard/wiki/Command-Line-Syntax-Transition-For-Users-(Pre-Transition)
********** The command line looks like this in the new syntax:
********** MarkDuplicates -INPUT SubBarcodes….bam -OUTPUT SubBarcodes….rmdup -METRICS_FILE SubBarcodes…*.dups.log -REMOVE_DUPLICATES true -VALIDATION_STRINGENCY SILENT -BARCODE_TAG RG -TMP_DIR /temp/
19:35:12.981 INFO NativeLibraryLoader - Loading libgkl_compression.so from jar:file:/home/yyp/miniconda3/envs/rna/share/picard-2.18.29-0/picard.jar!/com/intel/gkl/native/libgkl_compression.so
19:35:12.988 WARN NativeLibraryLoader - Unable to load libgkl_compression.so from native/libgkl_compression.so (No such file or directory)
19:35:12.989 INFO NativeLibraryLoader - Loading libgkl_compression.so from jar:file:/home/yyp/miniconda3/envs/rna/share/picard-2.18.29-0/picard.jar!/com/intel/gkl/native/libgkl_compression.so
19:35:12.990 WARN NativeLibraryLoader - Unable to load libgkl_compression.so from native/libgkl_compression.so (No such file or directory)
[Sun Oct 31 19:35:12 CST 2021] MarkDuplicates BARCODE_TAG=RG INPUT=[SubBarcodes….bam] OUTPUT=SubBarcodes….rmdup METRICS_FILE=SubBarcodes….dups.log REMOVE_DUPLICATES=true TMP_DIR=[/temp] VALIDATION_STRINGENCY=SILENT MAX_SEQUENCES_FOR_DISK_READ_ENDS_MAP=50000 MAX_FILE_HANDLES_FOR_READ_ENDS_MAP=8000 SORTING_COLLECTION_SIZE_RATIO=0.25 TAG_DUPLICATE_SET_MEMBERS=false REMOVE_SEQUENCING_DUPLICATES=false TAGGING_POLICY=DontTag CLEAR_DT=true DUPLEX_UMI=false ADD_PG_TAG_TO_READS=true ASSUME_SORTED=false DUPLICATE_SCORING_STRATEGY=SUM_OF_BASE_QUALITIES PROGRAM_RECORD_ID=MarkDuplicates PROGRAM_GROUP_NAME=MarkDuplicates READ_NAME_REGEX=<optimized capture of last three ‘:’ separated fields as numeric values> OPTICAL_DUPLICATE_PIXEL_DISTANCE=100 MAX_OPTICAL_DUPLICATE_SET_SIZE=300000 VERBOSITY=INFO QUIET=false COMPRESSION_LEVEL=5 MAX_RECORDS_IN_RAM=500000 CREATE_INDEX=false CREATE_MD5_FILE=false GA4GH_CLIENT_SECRETS=client_secrets.json USE_JDK_DEFLATER=false USE_JDK_INFLATER=false
[Sun Oct 31 19:35:12 CST 2021] Executing as yyp@DESKTOP-LRUCLJJ on Linux 4.19.128-microsoft-standard amd64; OpenJDK 64-Bit Server VM 10.0.2+13; Deflater: Jdk; Inflater: Jdk; Provider GCS is not available; Picard version: 2.18.29-SNAPSHOT
[Sun Oct 31 19:35:12 CST 2021] picard.sam.markduplicates.MarkDuplicates done. Elapsed time: 0.00 minutes.
Runtime.totalMemory()=536870912
To get help, see http://broadinstitute.github.io/picard/index.html#GettingHelp
Exception in thread “main” htsjdk.samtools.SAMException: Cannot read non-existent file: file:///home/yyp/picard-master/SubBarcodes….bam
at htsjdk.samtools.util.IOUtil.assertFileIsReadable(IOUtil.java:432)
at htsjdk.samtools.util.IOUtil.assertFileIsReadable(IOUtil.java:419)
at htsjdk.samtools.util.IOUtil.assertInputIsValid(IOUtil.java:395)
at htsjdk.samtools.util.IOUtil.assertInputsAreValid(IOUtil.java:471)
at picard.sam.markduplicates.MarkDuplicates.doWork(MarkDuplicates.java:256)
at picard.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:295)
at picard.cmdline.PicardCommandLine.instanceMain(PicardCommandLine.java:103)
at picard.cmdline.PicardCommandLine.main(PicardCommandLine.java:113)
(rna) yyp@DESKTOP-LRUCLJJ:~/picard-master$
picard报错的问题相关推荐
- 2021.11.23【bug笔记】丨picard运行报错:Exception in thread “main“ java.lang.UnsupportedClassVersionError
项目场景: RNA-seq对比对后bam文件绘制insert图片 问题描述: 执行picard命令时发生报错: JAVA报错代码: at java.lang.ClassLoader.defineCla ...
- Linux运行Java出现“Exception in thread main java.lang.OutOfMemoryError: Java heap space”报错...
在运行如下程序时出现"Exception in thread "main" java.lang.OutOfMemoryError: Java heap space&quo ...
- 2022.11.21【bug笔记】|bam文件报错:Cannot add sequence that already exists in SAMSequenceDictionary
项目场景: sam文件是通过hisat2,bowtie2或者bwa将rawdata进行比对后得到的包含比对信息的数据格式.经过samtools处理后得到的bam文件经常用于后续分析,比如RNA-seq ...
- GATK BaseRecalibrator报错”The contig order in dbsnp_138.hg38.vcf and reference is not the same“
GATK BaseRecalibrator报错"The contig order in dbsnp_138.hg38.vcf and reference is not the same&qu ...
- 继承WebMvcConfigurer 和 WebMvcConfigurerAdapter类依然CORS报错? springboot 两种方式稳定解决跨域问题
继承WebMvcConfigurer 和 WebMvcConfigurerAdapter类依然CORS报错???springboot 两种方式稳定解决跨域问题! 之前我写了一篇文章,来解决CORS报错 ...
- Access to XMLHttpRequest at file from origin ‘null‘ has been blocked by CORS policy谷歌浏览器本地打开项目js文件报错
Access to XMLHttpRequest at 'file:///xxxxx/PQ.BaseInfo.proto' from origin 'null' has been blocked by ...
- go build 编译报错 missing go.sum entry for module providing package
go build 编译报错 missing go.sum entry for module providing package 解决方法 // 移除未使用的依赖 go mod tidy 再次编译,就可 ...
- docker报错:driver failed programming external connectivity on endpoint, iptables:No chain by that name
docker 报错: Error response from daemon: Cannot restart container hello: driver failed programming ext ...
- mybatis查询报错:com.mysql.cj.exceptions.DataConversionException: Cannot determine value type from string
mybatis查询报错: com.mysql.cj.exceptions.DataConversionException: Cannot determine value type from strin ...
最新文章
- 当年轻人开始谈论AI伦理
- Perl默认的内部变量
- 为什么计算机存储的是二进制补码?
- JSP页面如何调用自己写的.class文件
- 自从我上了数据结构课之后……
- [20160704]Block recover using RMAN.txt
- 数据结构与算法之-----队列(Queue)
- ios 时间选择器月份设置中文_iOS 自定义日期选择器 PickerView优化版
- JAVA多线程之男朋友和女朋友之间的故事
- c++循环输入一个字符可录入空格和换行_如何快速录入数据?
- 分享一个超厉害的网站,几乎解决一切command not found问题
- 使用BPEL创建Web服务组合
- IAM之Tivoli Identity Manager(二)
- python在输出中间加空行_python输出空行
- Little_Women2.txt
- odoo中的消息通知
- 《麦肯锡意识》前言 解决问题的战略模型-思维导图
- 计算机专业可以谈手机吗,手机真能取代电脑吗?谈手机/PC系统的大一统
- msra数据集_干货下载 | 中文自然语言处理 语料/数据集
- Vue中实现电子围栏/围栏(高德地图)功能: