Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

How process non-10X RNA-seq library #8

Open
yuanlizhanshi opened this issue Jan 9, 2025 · 4 comments
Open

How process non-10X RNA-seq library #8

yuanlizhanshi opened this issue Jan 9, 2025 · 4 comments

Comments

@yuanlizhanshi
Copy link

I have a batch of self-constructed test data, where UMI and barcode information is located at fixed positions in R1, and R2 contains the sequencing data. I need to append the barcode and UMI information to R2, but in some cases, the library is quite short, and R2 can sequence through into R1. I need to trim the polyA tail in R2 before alignment. Would this type of sequencing data be compatible with preCellar?
The custom library is below.
1736389884180

@yuanlizhanshi
Copy link
Author

I should ask in seqspec

@kaizhang kaizhang reopened this Jan 9, 2025
@kaizhang
Copy link
Member

kaizhang commented Jan 9, 2025

This can be handled by precellar. But more detailed information about the library is needed, such as the "(9) Final library structure" in this example: https://teichlab.github.io/scg_lib_structs/methods_html/SMART-seq_family.html.

@yuanlizhanshi
Copy link
Author

This can be handled by precellar. But more detailed information about the library is needed, such as the "(9) Final library structure" in this example: https://teichlab.github.io/scg_lib_structs/methods_html/SMART-seq_family.html.

I encountered an issue while using the seqspec tool. When attempting to read the local spec.yaml file, an error occurred, preventing the file from being properly accessed. Here are the steps I followed to generate the file:

First, I initialized the sequence specification file with the following command:
seqspec init -o spec_test.yaml -n myassay -m rna -r "rna,R1.fastq.gz,r1_primer,25,pos:rna,R2.fastq.gz,r2_primer,0,neg" "((r1_primer:25,barcode:10,umi:10,cdna:150,r2_primer:0)rna)"
seqspec format spec_test.yaml > spec.yaml

image

@kaizhang
Copy link
Member

I suggest you use the seqspec templates here to start with: https://github.com/regulatory-genomics/precellar/tree/main/seqspec_templates.

Many current seqspec you can find do not conform the seqspec specification/schema: https://github.com/pachterlab/seqspec/blob/main/seqspec/schema/seqspec.schema.json

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants