HT-Seq Batch
0
0
Entering edit mode
3 months ago
Nai • 0

Will this shell script for batch ..Related HTSeq-count. will it create separate file for each sample or merge count file. Kindly interpret the scrpt also. Because I am new in batch scripting.

for i in *.sam
do
echo $i
./local/bin/htseq-count  $i ./hg19.gtf > $i.count
done
Batch script HT-Seq • 429 views
ADD COMMENT
2
Entering edit mode

why don't you run it and figure it out?

ADD REPLY
1
Entering edit mode

you can try this, it will generate a count matrix for all .sam files in the current folder

./local/bin/htseq-count *.sam ./hg19.gtf >count.mat.tsv
ADD REPLY
1
Entering edit mode

the script you show here will create one count file per input sam file (so no merged output, see Nitin Narwade comment to achieve that merged kind of output)

ADD REPLY
0
Entering edit mode

Thank you...

  1. Here I should not call .sam in a variable $i.
  2. I have one more query regarding to convert count.mat.tsv to TPM matrix. How can I do that? I have to find eQTL . ....
ADD REPLY
0
Entering edit mode

Thank you Lieven and Nitin! Can you guide me how to convert read count matrix to TPM matrix.

ADD REPLY
0
Entering edit mode

This post discussed TPM and other normalized count matrices and also has a link for a code to convert raw count to TPM.

ADD REPLY

Login before adding your answer.

Traffic: 1734 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6