Home > OS >  Problem when using hisat2 ratio
Problem when using hisat2 ratio

Time:11-22

Rookie, I was born a letter when using hisat2 comparing to the genome sequence appeared many problems, the following simple said I tread pit,
1, personal feel, unless you want to build contains exons and splice site index, if just compare to the genome, can be directly downloaded from hisat2 website reference index in the genome, in a folder after download (human is hg19 if download hg19 mm10) in mice, which are genome. 123,,,. Ht2 file
2, when using hisat2 comparison, the best index, x.h t2, fq. Gz file, advance deposit. Sam into the same folder, or may be wrong, and I don't know why because I very vegetables
3, before use hisat2 must hisat2 - h see the command, or you'll, pound-foolish! Hisat2 than general parameters of the following
Hisat2 - p - dta - 8 x hg19/genome - 1 sample_1. Fq. Gz - 2 sample_2. Fq. Gz - S sample. Sam -p 8 general how many nuclear to run, this watch your computer memory, 16 g I choose eight nuclear slightly point card, - dta is report, - x is your identity, hg19 is your deposit index directory, the genome is your prefix index file! Interrupt at this time, if you don't have to see hisat2 requirement, you may be behind the hg19 directly to the folder in - x, or will be one of the genome. X.h t2 using regular expressions genome. * held behind - x, or you could use the cat> the whole into a genome, these will cause hisat2 don't know you of the index file, see hisat2 request, you will understand, the somebody else just want prefix of the index, so you need to put hg19/genome - x can back!
- 1-2 behind fastq file must be written on the path, or like I said to put them into the same folder, folder named
I like to use alignedThis path can't go wrong and hisat2 won't convulsions error!
4, if your fastq file a lot, you can use for the terminal... do... Done to write a simple loop for I in xx yy 'do' seq hisat2-8 - x hg19/p genome (all of your index file prefix) - 1 sample ${I} 1. Fq. Gz - 2 sample ${I} 2. Fq. Gz - S sample ${I}. Sam
If you want to write a sh in vim, pay attention to don't forget to configure the PATH
5, not must ask more and more don't behind closed doors, or a waste of time! Wish everyone born letter as a master!
Finally append a train diagrams (my Ubuntu is so lovely, and the penguin although this is the label but I still want to use centos haha)

CodePudding user response:

Hisat2 I used for the construction of a genome index file, has been running for two days and two nights, but still didn't finish, I was thinking, the need for such a long time?

CodePudding user response:

reference 1st floor weixin_44194636 response:
I use to build the index file of the genome, hisat2 has been running for two days and two nights, but still didn't finish, I was thinking, the need for such a long time?

Whether the card, I don't build so long, 16 gb of memory built more than two hours

CodePudding user response:

Problem has been resolved, because 1) construct the index needed within about 200 g; 2) 8 or more processors,

CodePudding user response:

The
reference 3 floor weixin_44194636 response:
problem has been solved, because 1) construct the index is about 200 g in need; 2) 8 or more processors,

哈哈,还是硬伤,

CodePudding user response:

The
reference 3 floor weixin_44194636 response:
problem has been solved, because 1) construct the index is about 200 g in need; 2) 8 or more processors,


Excuse me, would you please tell me what you said is 200 gb of memory disk space? How can build 16 g, upstairs to 200 gb of memory, you, too scary

CodePudding user response:

Than as shown in figure error when how to do file does exist!

CodePudding user response:

Recently I use Hisat2 processing of mulberry silkworm transcriptome data, but I can't debug mistakes, also please teach a great god, the paper is as follows:
1. The silkworm genome is about 468.3 Mb, 28 chromosomes;
2. The server 288 nuclear, 1 t memory, almost as far as I'm a person with,
3. Hisat2 - build to build the index has not found the problem;
4. Hissat2 alignment, please have a look at an example below, use the top command to check the content usage, with the extension of time, the memory climbed all the way,

30 - dta hisat2 - t - p - x/home/RNAseq_2/source/silkworm/index//data/storage04 silkworm_tran - 1/RNAseq_2 silkworm/majorbio/data4antivirus cleandata/306 d3d1a_r1 - clean. Fastq. Gz - 2/data/storage04/RNAseq_2/silkworm/majorbio/data4antivirus cleandata/306 d3d1a_r2 - clean. Fastq./data/storage04/gz - S RNAseq_2/silkworm/majorbio/data4antivirus alignedFromHisat2Results/306 d3d1a. Sam

5. 22 G, is used to estimate target Sam file currently is 7.8 G, big % mem is 55%, feel the process endlessly in memory,
6. Try to use before 8 cores, 5 similar programs running at the same time, but running mistakes, the error information is as follows:

(ERR) : hisat2 - align died with signal 9 (KILL)

7. I try to Google, may also be a biggest memory burst, before no bugs, Sam file is not in, completely stop,
8. In addition to the above operation failure, is also an unusually long running time,

  • Related