For one thing, you're missing a space in your script, as you pasted it.
I=$FILESO=/raid/development/anusha/DGFDATA_Tissue_DATA/brain/brain.merge.bam
---
Additional comments.
Picard-tools 1.53 dates back to 2011.
The current version is Picard-tools 1.114.
I find Python scripts easier to read, and maintain, but perhaps bash is appropriate for such a short script.
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Picard tools merge sam isssues ;;;;;;
Hi All,
I am trying to use merge sam with script so that I can automate this for several tissues . Command is working when I am giving this in command line
################################################
!/bin/bash
FILES=''
for file in $(find /raid/development/anusha/DGFDATA_Tissue_DATA/brain -name \*.sorted.bam)
do
FILES="$FILES $file"
done
#echo $FILES
java -Xmx2g -jar /raid/software/src/picard-tools-1.53/MergeSamFiles.jar I=$FILESO=/raid/development/anusha/DGFDATA_Tissue_DATA/brain/brain.merge.bam AS=FALSE USE_THREADING=TRUE VALIDATION_STRINGENCY=SILENT&
#############################################
anusha@cn1:/raid/development/anusha/DGFDATA_Tissue_DATA> echo $FILES
/raid/development/anusha/DGFDATA_Tissue_DATA/brain/SRR203405.sorted.bam /raid/development/anusha/DGFDATA_Tissue_DATA/brain/SRR203400.sorted.bam /raid/development/anusha/DGFDATA_Tissue_DATA/brain/SRR203406.sorted.bam /raid/development/anusha/DGFDATA_Tissue_DATA/brain/SRR203407.sorted.bam /raid/development/anusha/DGFDATA_Tissue_DATA/brain/SRR203408.sorted.bam /raid/development/anusha/DGFDATA_Tissue_DATA/brain/SRR203404.sorted.bam /raid/development/anusha/DGFDATA_Tissue_DATA/brain/SRR203401.sorted.bam /raid/development/anusha/DGFDATA_Tissue_DATA/brain/SRR203402.sorted.bam /raid/development/anusha/DGFDATA_Tissue_DATA/brain/SRR203403.sorted.bam /raid/development/anusha/DGFDATA_Tissue_DATA/brain/SRR203409.sorted.bam
###########################################
ERROR: Invalid argument '/raid/development/anusha/DGFDATA_Tissue_DATA/brain/SRR203400.sorted.bam'
Merges multiple SAM/BAM files into one file.
######################################
I am getting error like this
Hope to help on this . Thanks in advance .Tags: None
Latest Articles
Collapse
-
by seqadmin
Innovations in next-generation sequencing technologies and techniques are driving more precise and comprehensive exploration of complex biological systems. Current advancements include improved accessibility for long-read sequencing and significant progress in single-cell and 3D genomics. This article explores some of the most impactful developments in the field over the past year.
Long-Read Sequencing
Long-read sequencing has seen remarkable advancements,...-
Channel: Articles
12-02-2024, 01:49 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Today, 08:24 AM
|
0 responses
10 views
0 likes
|
Last Post
by seqadmin
Today, 08:24 AM
|
||
Started by seqadmin, Yesterday, 07:41 AM
|
0 responses
8 views
0 likes
|
Last Post
by seqadmin
Yesterday, 07:41 AM
|
||
Started by seqadmin, 12-11-2024, 07:45 AM
|
0 responses
15 views
0 likes
|
Last Post
by seqadmin
12-11-2024, 07:45 AM
|
||
Started by seqadmin, 12-10-2024, 07:59 AM
|
0 responses
14 views
0 likes
|
Last Post
by seqadmin
12-10-2024, 07:59 AM
|
Leave a comment: