Using unaligned reads as input for trinity
0
0
Entering edit mode
5.5 years ago

'I have done HiSAT mapping for my RNAseq sample and I have got unaligned FASTQ sanger file. Can i use this unaligned files as input for running trinity or Should i use the original SRA reads file.

RNA-Seq Assembly trinity • 1.7k views
ADD COMMENT
2
Entering edit mode

What you should do is to take some of those unaligned reads and check them using blast. At a minimum there may be contamination in data you are using. If so there would be no point in wasting further time on this dataset.

If the reads are from the right genome then you should spend some time figuring out why HISAT2 is unable to align them.

ADD REPLY
0
Entering edit mode

Thank you. I'll check that.

ADD REPLY
1
Entering edit mode

In addition to what genomax proposes, I would also just support your approach, or better simply running a full Trinity assembly. It doesn't cost you much except some CPU time and after that you can simply blast all the assembled transcripts. More generally speaking, whenever problems in an analysis arise, one should try out a bunch of different things and make some hypothesis about what is going wrong. Btw, what is your definition of a low mapping rate?

  • wrong genome
  • draft genome with lots of gaps, unmapped regions
  • contamination
  • ribosomal RNA
  • failed sequencing or protocol
ADD REPLY
0
Entering edit mode

Well technically you can anything you wish. Feed the software with something and it will spit out something more but why do want to do that ?

Share your objective here.

ADD REPLY
0
Entering edit mode

Because the result of HISAT mapping was low.

ADD REPLY
1
Entering edit mode

You should always mention that in your post along with mapping statistics. Anyways, look at this old post, there is a long discussion on similar problem I had faced couple of months back!

Too low mapping percentage using HISAT2 on human reference genome.

ADD REPLY
0
Entering edit mode

Thanks. I'll have a look at it.

ADD REPLY

Login before adding your answer.

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