diff --git a/rnaseq/step5_count/snippet.sh b/rnaseq/step5_count/snippet.sh index 6aff43a..fd3270c 100644 --- a/rnaseq/step5_count/snippet.sh +++ b/rnaseq/step5_count/snippet.sh @@ -22,6 +22,7 @@ BAM="${PROJECT_DIRECTORY}/step3/alignment/sample1_Aligned.sortedByCoord.out.bam" # `--stranded` flag. htseq-count \ --order pos \ + -f bam \ ${BAM} \ ${GTF_FILE} \ - > ${COUNT_OUTPUT}/counts.txt \ No newline at end of file + > ${COUNT_OUTPUT}/counts.txt