search

Teichlab / tracer

TraCeR - reconstruction of T cell receptor sequences from single-cell RNAseq data
Other
124 stars 48 forks source link

Adding assembled_file option assembly process #77

Open JohnMCMa opened 6 years ago

JohnMCMa commented 6 years ago

As we discussed in Teichlab/bracer#21 , when using 10X Genomics outputs, reads should be split by cell and then called their VDJ types using bracer assembly's --assembled_file option.

I suppose the main idea also applies to TraCeR, but tracer assembly doesn't have a similar function. This means I'll have to split the fastq, calculate fragment length and SD and run the entire VDJ assembly pipeline despite the cellranger vdj output files are already assembled. My rough calculation shows this increases the running time by 3-4X.

mstubb commented 6 years ago

Yes, this is a direct consequence of TraCeR being written prior to release of the 10x Genomics VDJ product. We’ll put this functionality on the to-do list but I’d be happy to accept a pull request that does this if you’d like to write one.

All the best,

Mike

On 6 Sep 2018, at 15:23, John Ma notifications@github.com wrote:

As we discussed in issue 21 in BraCeR, when using 10X Genomics outputs, reads should be split by cell and then called their VDJ types using bracer assembly's --assembled_file option.

I suppose the main idea also applies to TraCeR, but tracer assembly doesn't have a similar function. This means I'll have to split the fastq, calculate fragment length and SD and run the entire VDJ assembly pipeline despite the cellranger vdj output files are already assembled. My rough calculation shows this increases the running time by 3-4X.

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/Teichlab/tracer/issues/77, or mute the thread https://github.com/notifications/unsubscribe-auth/ABFwhk3vFO-dRy-6-zq4wKri9egJHaQzks5uYS_JgaJpZM4WdCt4.