How i can get unalignment read in hisat2 ?
0
0
Entering edit mode
17 months ago
Ahmed ▴ 10

is there any option in hisat2 for viewing unalignment read? first: i make index for my refrence genome with this command:

 hisat2-build -p 3 Danio_rerio.GRCz11.cdna.all.fa /home/ngs/zebra_fish 

second: i make alignment my fastq file with reference genome, using this command:

hisat2 -p 4 -x /home/ngs/zebra_fish/zebra_fish -U SRR1048063_Zebrafish_pineal_gland_dark_replica_3_1.fastq -S SRR1048063_Zebrafish_pineal_gland_dark_replica_3_1_chrl.sam

i get my output in .sam format, then i transfer it to bam file, using this command:

samtools view -bS SRR1048063_Zebrafish_pineal_gland_dark_replica_3_1_chrl.sam > eg2.bam
after that i get my bam file 

i read the manual of hisat2 and i cannot find any option to view the unalignment reads the question i how i can get my unalignment read and view it with readable way? Note: the data is RNA seq

hisat2 • 848 views
ADD COMMENT
0
Entering edit mode

I extract it from bam file but, how i can view it on terminal?

ADD REPLY
0
Entering edit mode

To view the file contents, you can use any unix file command head/tail/cat/more etc if the file is uncompressed. Use zcat if the file is compressed.

ADD REPLY

Login before adding your answer.

Traffic: 1633 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