Use of Hisat2 with single end data
Respected all Greetings can anyone let me know what is the command line for single end data mapping on refernce genome in hisat2 2.what is the difference between single and double end refernce genome index
• 3,417 views
•
link
updated
by
Ram
•
written
by
usman.mehmood •
•
written
by
usman.mehmood •
1 answer
If you're asking about how to align with single end data, then use the -U flag instead of flagging your first pair as -1 and second as -2.
hisat2 [Flags you want] -U /path/to/reads.fq -x /path/to/genomeIndex -S out.sam
If you're asking about genome indexing, it doesn't technically matter, just use hisat2-build on your reference genome fastA.
hisat2-build [Flags you want] reference.fastA genomeIndex
Log in to answer this question.
More posts like this
-
FastQC for paired end data
written by Dk •
-
HISAT2 index difference GRCh38 vs hg38
written by Atsunop •
-
Differences in HISAT2 indexes for download
written by vitor •
-
Hisat2 error: --read-lengths arg must be at least 20
written by aindap
-
HISAT2 reference genomes and starting code
written by Morris_Chair
-
What is the difference between HISAT2 index "genome.tran" and "genome"
written by Chevy Xu •
-
Tophat alignment against BALB_cJ_R
written by anu014
-
Python scripts to map multiple pair end files to the genome
written by Bioinfonext
-
mapping read is low for miRNA data after trimming
written by cyn.liu •
-
Tophat unable to produce accepted_hits.bam file for single end read data
written by nikhilvgbt •
I will check it. Thankyou for the quick reply
Use
ADD COMMENT/ADD REPLYwhen responding to existing posts to keep threads logically organized.SUBMIT ANSWERis for new answers to original question.