RNA-Seqのカウントツールを色々試しているのですが、同じgeneIDなのにツールごとに全然違う数値が出てきます。
手順としては、trimgaloreでトリミング&QCし、その後hisat2でmapping, samtoolsでソートBAMファイルに変換、そしてstringtieまたはfeatureCountsまたはHTSeqでカウントしました。あとhisat2とRSEMでもカウントしてみました。
大体のgeneIDでは傾向は似ているのですが、たまにツールによって全然違う(HTseqでは数万のカウントなのにrsemやfeaturecountでは0に近いカウント)を示すgeneがあります。
この原因はなんなのでしょうか?またこうなってしまった場合はどう対応したらよいのでしょうか? |
|