-
Notifications
You must be signed in to change notification settings - Fork 8
Description
********** The command line looks like this in the new syntax:
********** AddOrReplaceReadGroups -I G1_FRAS210279619-2r_dir/G1_tmp/3-0-0Combine_bam/G1_combine.bam -O G1_FRAS210279619-2r_dir/G1_tmp/3-0-0Combine_bam/G1_combine_rgadd.bam -SO coordinate -RGID 1 -RGLB lib1 -RGPL illumina -RGPU unit1 -RGSM 20
23:24:23.901 INFO NativeLibraryLoader - Loading libgkl_compression.so from jar:file:/training/nong/pp/RADAR/tools/picard.jar!/com/intel/gkl/native/libgkl_compression.so
[Tue Jan 25 23:24:23 CST 2022] AddOrReplaceReadGroups INPUT=G1_FRAS210279619-2r_dir/G1_tmp/3-0-0Combine_bam/G1_combine.bam OUTPUT=G1_FRAS210279619-2r_dir/G1_tmp/3-0-0Combine_bam/G1_combine_rgadd.bam SORT_ORDER=coordinate RGID=1 RGLB=lib1 RGPL=illumina RGPU=unit1 RGSM=20 VERBOSITY=INFO QUIET=false VALIDATION_STRINGENCY=STRICT 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
[Tue Jan 25 23:24:23 CST 2022] Executing as nong@amax on Linux 5.11.0-46-generic amd64; OpenJDK 64-Bit Server VM 1.8.0_312-b07; Deflater: Intel; Inflater: Intel; Provider GCS is not available; Picard version: 2.18.27-SNAPSHOT
[Tue Jan 25 23:24:23 CST 2022] picard.sam.AddOrReplaceReadGroups done. Elapsed time: 0.00 minutes.
Runtime.totalMemory()=2058354688
To get help, see http://broadinstitute.github.io/picard/index.html#GettingHelp
Exception in thread "main" htsjdk.samtools.SAMFormatException: Error parsing SAM header. Problem parsing @pg key:value pair. Line:
@pg ID:hisat2 PN:hisat2 VN: CL:"/home/nong/anaconda3/envs/RADAR/bin/hisat2-align-s --wrapper basic-0 --rna-strandness RF --no-mixed --secondary --no-temp-splicesite --known-splicesite-infile /dat1/dat/ref/hg38/RADAR/reference/Human/GRCh38/annotation_splice_sites.bed --no-softclip --score-min L,-16,0 --mp 7,7 --rfg 0,7 --rdg 0,7 --max-seeds 20 -k 10 --dta -t -p 16 -x /dat1/dat/ref/hg38/RADAR/reference/Human/GRCh38/GRCh38_all.fa --passthrough --read-lengths 150 -1 /tmp/2813493.inpipe1 -2 /tmp/2813493.inpipe2"; File /dat1/dat/ref/hg38/RADAR/reference/Human/GRCh38/G1_FRAS210279619-2r_dir/G1_tmp/3-0-0Combine_bam/G1_combine.bam; Line number 197
at htsjdk.samtools.SAMTextHeaderCodec.reportErrorParsingLine(SAMTextHeaderCodec.java:268)
at htsjdk.samtools.SAMTextHeaderCodec.access$200(SAMTextHeaderCodec.java:45)
at htsjdk.samtools.SAMTextHeaderCodec$ParsedHeaderLine.(SAMTextHeaderCodec.java:317)
at htsjdk.samtools.SAMTextHeaderCodec.decode(SAMTextHeaderCodec.java:98)
at htsjdk.samtools.BAMFileReader.readHeader(BAMFileReader.java:676)
at htsjdk.samtools.BAMFileReader.(BAMFileReader.java:298)
at htsjdk.samtools.BAMFileReader.(BAMFileReader.java:176)
at htsjdk.samtools.SamReaderFactory$SamReaderFactoryImpl.open(SamReaderFactory.java:396)
at picard.sam.AddOrReplaceReadGroups.doWork(AddOrReplaceReadGroups.java:152)
at picard.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:295)
at picard.cmdline.PicardCommandLine.instanceMain(PicardCommandLine.java:103)
at picard.cmdline.PicardCommandLine.main(PicardCommandLine.java:113)
NOTICE: RNA-editing calling is done.