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
The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...-
Channel: Articles
07-08-2024, 03:19 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 07-25-2024, 06:46 AM
|
0 responses
9 views
0 likes
|
Last Post
by seqadmin
07-25-2024, 06:46 AM
|
||
Started by seqadmin, 07-24-2024, 11:09 AM
|
0 responses
26 views
0 likes
|
Last Post
by seqadmin
07-24-2024, 11:09 AM
|
||
Started by seqadmin, 07-19-2024, 07:20 AM
|
0 responses
160 views
0 likes
|
Last Post
by seqadmin
07-19-2024, 07:20 AM
|
||
Started by seqadmin, 07-16-2024, 05:49 AM
|
0 responses
127 views
0 likes
|
Last Post
by seqadmin
07-16-2024, 05:49 AM
|
Leave a comment: