Weird Fragment Length distribution multiQC
1
0
Entering edit mode
7 weeks ago
_deb • 0

I have several samples coming from an RNAseq experiment generated with Illumina sequencing (Pair-End sequencing of 151bp).

I used Salmon for indexing the de novo transcriptome and then I quantified in mapped-based mode against the created index. However, checking the output from MultiQC I noted a very weird thing. The fragment length distribution suddenly drops around 145-150bp and then goes up again... Salmon MultiQC drop image

Is there someone that could have an explanation for this behaviour?

Any putative explanation is welcome!

Salmon RNA-seq paired-end MultiQC • 569 views
ADD COMMENT
2
Entering edit mode
6 weeks ago

Thats weird... could it have something to do with the way Salmon deals with dovetail reads ? For instance, these reads could be excluded when the insert about the size of the read length (150), but not when the insert is shorter (<145). I wonder what would happen if you use option --allowDovetail with Salmon in mapping mode. My guess is that it would solve the issue.

ADD COMMENT
0
Entering edit mode

Thanks a lot for your answer! In your opinion, could this affect the downstream analyses (e.g. DE analysis)?

ADD REPLY
0
Entering edit mode

Can you tell us if the option mentioned above mitigates your observation? Have you tried it? How many cycles were sequenced for your data (300 each way?)

ADD REPLY
0
Entering edit mode

Yes, I tried adding the option --allowDovetail, and it solved the problem. It seems a normal default behaviour of Salmon, so there is nothing to be too worried about apparently. Thank you again!

ADD REPLY
0
Entering edit mode

Please accept the answer by @Carlo to provide closure to this thread.

ADD REPLY

Login before adding your answer.

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