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

Excluding undetermined index reads

parent 699d2afd
No related branches found
No related tags found
No related merge requests found
...@@ -52,7 +52,7 @@ perl -pi -e 's/\r//' $PROJECT_ID.ped ...@@ -52,7 +52,7 @@ perl -pi -e 's/\r//' $PROJECT_ID.ped
mkdir -p $READS_DIR/$PROJECT_ID mkdir -p $READS_DIR/$PROJECT_ID
# generate the family_ids list - makes strong assumption about relative paths! # generate the family_ids list - makes strong assumption about relative paths!
ls ../../data/$PROJECT_ID/*.gz | cut -d '/' -f 5 | cut -f 1,3 -d '_' | sort -u > $PROJECT_ID.family_ids.txt ls ../../data/$PROJECT_ID/*.gz | grep -v Undetermined | cut -d '/' -f 5 | cut -f 1,3 -d '_' | sort -u > $PROJECT_ID.family_ids.txt
PLATE_ID=`cut -f 1 -d '_' $PROJECT_ID.family_ids.txt | sort -u` PLATE_ID=`cut -f 1 -d '_' $PROJECT_ID.family_ids.txt | sort -u`
for FAMILY_ID in `cut -f 2 -d '_' $PROJECT_ID.family_ids.txt` for FAMILY_ID in `cut -f 2 -d '_' $PROJECT_ID.family_ids.txt`
......
...@@ -66,7 +66,7 @@ do ...@@ -66,7 +66,7 @@ do
for FILE in `ls $DOWNLOAD_DIR/$PROJECT_ID/*${SAMPLE}*.gz` for FILE in `ls $DOWNLOAD_DIR/$PROJECT_ID/*${SAMPLE}*.gz`
do do
echo "$FILE,${SAMPLE}_${FAMILY_ID},$FAMILY_ID,$SEX,$PHENOTYPE,$TARGET" >> ${VERSION}_${PROJECT_ID}_${FAMILY_ID}.csv echo "$FILE,${SAMPLE},$FAMILY_ID,$SEX,$PHENOTYPE,$TARGET" >> ${VERSION}_${PROJECT_ID}_${FAMILY_ID}.csv
done done
done 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