metatran_LBE/20200811_count_htseq.sh
b4a73f19d755master

20200811_count_htseq.sh
View Options

			conda activate count_env

			export WORKING_DIR=//media/bioinfoteam/SCRATCH_850GB/temp_metatranscriptomic/bowtie_output/SAMTools
			export DATABASE=/media/bioinfoteam/SCRATCH_850GB/databases/annotation_hybridAB/d427AB_annotation_removeFASTA.gff

			cd $WORKING_DIR

			#Count number of reads mapping
			for filename in $(ls *'.sam')

			do

			#echo $filename

			fileID=`expr match "$filename" '$[0-9]*C_map_on_d427_hybridAB_sort$'`

			#echo $fileID

			export OUTPUT=/media/bioinfoteam/SCRATCH_850GB/temp_metatranscriptomic/counting_htseq/${fileID}

			htseq-count -o $OUTPUT.sam --nonunique all $filename $DATABASE > $OUTPUT.txt

			done