User Tools

Site Tools


custom_scripts_for_16s_rrna_sequence_analysis

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Next revision
Previous revision
custom_scripts_for_16s_rrna_sequence_analysis [2019/09/09 21:51] – created hyjeongcustom_scripts_for_16s_rrna_sequence_analysis [2021/03/17 13:09] (current) – external edit 127.0.0.1
Line 7: Line 7:
   # https://github.com/edamame-course/Amplicon_Analysis/blob/master/resources/Merged_Reads_Script.sh   # https://github.com/edamame-course/Amplicon_Analysis/blob/master/resources/Merged_Reads_Script.sh
   #   #
 +  
   out=merged_dir   out=merged_dir
 +  
   if [ ! -d merged_reads ]; then   if [ ! -d merged_reads ]; then
       echo Creating ${out} directory...       echo Creating ${out} directory...
       mkdir ${out}       mkdir ${out}
   fi   fi
 +  
   for f in *F_sub.*   for f in *F_sub.*
   do   do
 +  
       if [ ! -e ${f} ]; then       if [ ! -e ${f} ]; then
           echo Can\'t find 'F_sub' reads!           echo Can\'t find 'F_sub' reads!
Line 23: Line 23:
           exit           exit
       fi       fi
 +  
       r=$(sed -e "s/F_sub/R_sub/" <<< "$f")       r=$(sed -e "s/F_sub/R_sub/" <<< "$f")
       file=$(cut -d_ -f1 <<< "$f")       file=$(cut -d_ -f1 <<< "$f")
custom_scripts_for_16s_rrna_sequence_analysis.1568033481.txt.gz · Last modified: 2021/03/17 13:09 (external edit)