Question about ISSAC seq computational pipeline

[teichlab]

Question about ISSAC seq computational pipeline

[dKlee99]

Also, was there any reasons why you --clip3pNbases 116 option and make reads compaing that 10x r1 uses around 90 bp ?

[dbrg77]

Hi @dKlee99

Thanks for your question. Back in the time, we were not very familiar with all the STARsolo options, so we kept to the default most of the time. There was not any particular reasons. We probably overlooked that option.

Indeed, GeneFull is the preferred option considering the method is mainly profiling the poly-A RNA from the nucleus. We are now using GeneFull as well.

In terms of --clip3pNbases 116, this is due to the nature of the SHERRY library, which is generally short. We are doing 151bp PE sequencing. The majority of the bases at the 3' end are just adapters, so we only used the first 36 bp for the mapping. It seems the first 50 bp also work well.

I hope this helps.

Xi

[dKlee99]

Dear Dr. Chen,

Thank you so much for your reply!