Comments (7)
https://github.com/singleron-RD/CeleScope/blob/master/doc/chemistry.md
singleron kit v1 corresponds to chemistry scopeV2.2.1, which requires at least 69bp R1. If you have 68bp R1, the last 1bp UMI will be thrown away.
change -soloUMIposition 0_57_0_68 --soloUMIlen 12
to -soloUMIposition 0_57_0_67 --soloUMIlen 11
should work.
from celescope.
https://github.com/singleron-RD/CeleScope/blob/master/doc/chemistry.md singleron kit v1 corresponds to chemistry scopeV2.2.1, which requires at least 69bp R1. If you have 68bp R1, the last 1bp UMI will be thrown away. change
-soloUMIposition 0_57_0_68 --soloUMIlen 12
to-soloUMIposition 0_57_0_67 --soloUMIlen 11
should work.
In multi_rna, use '--chemistry customized --pattern C8L16C8L16C8L1U11T18 --whitelist "./scopeV2.2.1/bclist ./scopeV2.2.1/bclist ./scopeV2.2.1/bclist"
from celescope.
Thanks for your advice! I found the linker in scopeV2.2.1 was "TCGGTGACAGCCATATCGTAGTCAGAAGCTGAC", and the last base "C" was not sequenced in BGI data. So, the "L1" in '--pattern C8L16C8L16C8L1U11T18' was not exact I think.
from celescope.
https://github.com/singleron-RD/CeleScope/blob/master/doc/chemistry.md singleron kit v1 corresponds to chemistry scopeV2.2.1, which requires at least 69bp R1. If you have 68bp R1, the last 1bp UMI will be thrown away. change
-soloUMIposition 0_57_0_68 --soloUMIlen 12
to-soloUMIposition 0_57_0_67 --soloUMIlen 11
should work.In multi_rna, use '--chemistry customized --pattern C8L16C8L16C8L1U11T18 --whitelist "./scopeV2.2.1/bclist ./scopeV2.2.1/bclist ./scopeV2.2.1/bclist"
Thanks for your advice! I found the linker in scopeV2.2.1 was "TCGGTGACAGCCATATCGTAGTCAGAAGCTGAC", and the last base "C" was not sequenced in BGI data. So, the "L1" in '--pattern C8L16C8L16C8L1U11T18' was not exact I think.
from celescope.
from celescope.
if so,just remove the L1 and use U12
…
---------- 该邮件从移动设备发送
--------------原始邮件-------------- 发件人:"biozhao77 @.>; 发送时间:2024年3月5日(星期二) 晚上9:19 收件人:"singleron-RD/CeleScope" @.>; 抄送:"周义其 @.>;"Comment @.>; 主题:Re: [singleron-RD/CeleScope] BGI data was not suitable for version 2.0.7 (Issue #276) ----------------------------------- https://github.com/singleron-RD/CeleScope/blob/master/doc/chemistry.md singleron kit v1 corresponds to chemistry scopeV2.2.1, which requires at least 69bp R1. If you have 68bp R1, the last 1bp UMI will be thrown away. change -soloUMIposition 0_57_0_68 --soloUMIlen 12 to -soloUMIposition 0_57_0_67 --soloUMIlen 11 should work. In multi_rna, use '--chemistry customized --pattern C8L16C8L16C8L1U11T18 --whitelist "./scopeV2.2.1/bclist ./scopeV2.2.1/bclist ./scopeV2.2.1/bclist" Thanks for your advice! I found the linker in scopeV2.2.1 was "TCGGTGACAGCCATATCGTAGTCAGAAGCTGAC", and the last base "C" was not sequenced in BGI data. So, the "L1" in '--pattern C8L16C8L16C8L1U11T18' was not exact I think. — Reply to this email directly, view it on GitHub, or unsubscribe. You are receiving this because you commented.Message ID: @.***>
I'm trying. But I have another question, if the last base "C" will be thrown when L1 was removed? Because I noticed that the linker was 33bp including "C".
from celescope.
The lastest pipeline actually does NOT use the sequence of the linker. It just use the --pattern
to locate the barcode segments and check if the barcode segments match the whitelist or not.
from celescope.
Related Issues (20)
- Subread and cutadapt working as Single-end reads BUT I have Paired end reads HOT 1
- Some ISSUEs concerning about constructing REFERENCE on vdj sequencing and full length immune receptor sequencing HOT 6
- No bclist found for chemistry scopeV1 under HOT 2
- aggr function HOT 1
- Problem in mkref command HOT 2
- Problem in running Celescope for RNA HOT 4
- Splitting VDJ fails HOT 3
- Problem installing celescope, cannot build wheel for pysam HOT 2
- Auto Chemistry Detection Failing and FIFO File Generation Problem HOT 7
- Celescope needs an improved documentation/manual HOT 3
- Individual steps are not working. HOT 4
- celescope rna mkref but no folder output HOT 2
- 1.12.0 can be installed, but failed to be run. HOT 3
- Problem with mkref
- IndexError: sequence length is not enough in R1 read HOT 7
- featureCounts subprocess.CalledProcessError HOT 1
- Multiple errors trying to analyze flv_trust4 data HOT 5
- FileNotFoundError: Did not find file SMA/outs/filtered/matrix.mtx.gz. HOT 1
- numpy.dtype size changed HOT 4
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from celescope.