Skip to content
Snippets Groups Projects
Commit e785f110 authored by ameyner2's avatar ameyner2
Browse files

Using symlinks to get around bcbio issues with file merging

parent a29ba755
No related branches found
No related tags found
No related merge requests found
......@@ -64,7 +64,29 @@ do
SEX=`head -n $i ${PROJECT_ID}_${FAMILY_ID}.ped | tail -n 1 | cut -f 5`
PHENOTYPE=`head -n $i ${PROJECT_ID}_${FAMILY_ID}.ped | tail -n 1 | cut -f 6`
for FILE in `ls $DOWNLOAD_DIR/$PROJECT_ID/raw_data/*/*${SAMPLE}*/*.gz`
# create symlinks for problematic filenames
for FILE in `ls $DOWNLOAD_DIR/$PROJECT_ID/raw_data/*/*${SAMPLE}*/*_1_*_1.fastq.gz`
do
newname=`echo $FILE | sed -e 's/_1_/_one_/'`
ln -s $FILE ${newname%1.fastq.gz}R1.fastq.gz
done
for FILE in `ls $DOWNLOAD_DIR/$PROJECT_ID/raw_data/*/*${SAMPLE}*/*_1_*_2.fastq.gz`
do
newname=`echo $FILE | sed -e 's/_1_/_one_/'`
ln -s $FILE ${newname%2.fastq.gz}R2.fastq.gz
done
for FILE in `ls $DOWNLOAD_DIR/$PROJECT_ID/raw_data/*/*${SAMPLE}*/*_2_*_1.fastq.gz`
do
newname=`echo $FILE | sed -e 's/_2_/_two_/'`
ln -s $FILE ${newname%1.fastq.gz}R1.fastq.gz
done
for FILE in `ls $DOWNLOAD_DIR/$PROJECT_ID/raw_data/*/*${SAMPLE}*/*_2_*_2.fastq.gz`
do
newname=`echo $FILE | sed -e 's/_2_/_two_/'`
ln -s $FILE ${newname%2.fastq.gz}R2.fastq.gz
done
for FILE in `ls $DOWNLOAD_DIR/$PROJECT_ID/raw_data/*/*${SAMPLE}*/*_R[1,2].fastq.gz`
do
echo "$FILE,$SAMPLE,$FAMILY_ID,$SEX,$PHENOTYPE,$TARGET" >> ${PROJECT_ID}_${FAMILY_ID}.csv
done
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment