Question: Hisat2 splicesite file vs Tophat
0
gravatar for shrutidabral
19 months ago by
shrutidabral10
shrutidabral10 wrote:

I am using hisat2 for RNA seq data i have some queries --

1) There is option --phred33 ----Input qualities are ASCII chars equal to the Phred quality plus 33. This is also called the "Phred+33" encoding, which is used by the very latest Illumina pipelines.

Is this means Ididn't have to run trimmomatic and using this score for alignment the selection quality of reads will be phred 33 ? Please explain.

2) Command I am using --

hisat2 -p 24 --dta --known-splicesite-infile splicesites.txt  --phred33  -x hg38  -1 31_1.fastq.gz -2 /media31_2.fastq.gz -S H031.sam --summary-file H031_summary.txt --novel-splicesite-outfile H031_outfile --novel-splicesite-infile H031_infile

In H031_outfile file my getting result like

> 1 15037   15795   -

Anyone can explain me how to interpret this result and this result is totally different form tophat junction file or its gives us appropriate information.How can I get the information of reads skipped ?

rna-seq tool assembly • 536 views
ADD COMMENTlink modified 19 months ago by Bastien Hervé4.9k • written 19 months ago by shrutidabral10
0
gravatar for Bastien Hervé
19 months ago by
Bastien Hervé4.9k
Karolinska Institutet, Sweden
Bastien Hervé4.9k wrote:

1) No the phred33 option is use by default and will not filter out any reads based on their quality. It is a table cotation for score, on old sequencing machine phred64 table where used to report scores. Nowadays, sequencing machine use phred33 score table. I suggest you to not use this option at all if you get reads from recent machine.

See also : http://drive5.com/usearch/manual/quality_score.html

2) The result of --novel-splicesite-outfile option is explained in the HISAT2 manual

chromosome name <tab> genomic position of the flanking base on the left side of an intron <tab> genomic position of the flanking base on the right <tab> strand (+, -, and .) '.' indicates an unknown strand for non-canonical splice sites.

ADD COMMENTlink modified 19 months ago • written 19 months ago by Bastien Hervé4.9k

Thankyou for explaination. For Second question I know the defination of result ouput of command --novel-splicesite-outfile but i did not give the number of reads that contain the junction as tophat give in junction file thats why i ask that How to interpret these result. Further more I used regtools get the results as in tophat format but i am not satisfy because by comparing the position of chromosome from output file of hisat2 and regtools is totally different .

Any suggestions?

ADD REPLYlink written 19 months ago by shrutidabral10

The information you're looking for isn't present in the file you're looking in. If the novel sites are that disparate then it's likely that neither tool should be relied upon. Have a look at STAR, which generally performs better in comparisons for finding novel splice sites.

ADD REPLYlink written 19 months ago by Devon Ryan97k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2237 users visited in the last hour