Giter VIP home page Giter VIP logo

bgi_to_illumina's Introduction

Convert BGI to Illumina pipeline

A example sample Info

Sample_ID Name Index 7 Lane
Truseq06 1_31GEX GCAACAAA Lane3
Truseq06 2_31GEX TAGTTGTC Lane3
Truseq06 3_31GEX CGCCATCG Lane3
Truseq06 4_31GEX ATTGGCGT Lane3

FileNames are like:

ID reads
1 V300079632_L03_1_31GEX_1.fq.gz
V300079632_L03_1_31GEX_2.fq.gz
2 V300079632_L03_2_31GEX_1.fq.gz
V300079632_L03_2_31GEX_2.fq.gz
3 V300079632_L03_3_31GEX_1.fq.gz
V300079632_L03_3_31GEX_2.fq.gz
4 V300079632_L03_4_31GEX_1.fq.gz
V300079632_L03_4_31GEX_2.fq.gz

Folders

CRData
├── BC_Attached ## output of 01_batch_attachBarcodes.py
│   └── 31GEX
│       ├── 31GEX_ATTGGCGT_L003_R1.fastq.gz
│       ├── 31GEX_ATTGGCGT_L003_R2.fastq.gz
│       ├── 31GEX_CGCCATCG_L003_R1.fastq.gz
│       ├── 31GEX_CGCCATCG_L003_R2.fastq.gz
│       ├── 31GEX_GCAACAAA_L003_R1.fastq.gz
│       ├── 31GEX_GCAACAAA_L003_R2.fastq.gz
│       ├── 31GEX_TAGTTGTC_L003_R1.fastq.gz
│       └── 31GEX_TAGTTGTC_L003_R2.fastq.gz
├── BGI ## original BGI data
│   └── 31GEX
│       ├── V300079632_L03_1_31GEX_1.fq.gz
│       ├── V300079632_L03_1_31GEX_2.fq.gz
│       ├── V300079632_L03_2_31GEX_1.fq.gz
│       ├── V300079632_L03_2_31GEX_2.fq.gz
│       ├── V300079632_L03_3_31GEX_1.fq.gz
│       ├── V300079632_L03_3_31GEX_2.fq.gz
│       ├── V300079632_L03_4_31GEX_1.fq.gz
│       └── V300079632_L03_4_31GEX_2.fq.gz
├── Converted ## output of 02_run_convertHeaders.py
│   └── 31GEX ## can take this as input of cellranger
│       ├── 31GEX_S1_L003_R1_001.fastq.gz
│       ├── 31GEX_S1_L003_R2_001.fastq.gz
│       ├── 31GEX_S2_L003_R1_001.fastq.gz
│       ├── 31GEX_S2_L003_R2_001.fastq.gz
│       ├── 31GEX_S3_L003_R1_001.fastq.gz
│       ├── 31GEX_S3_L003_R2_001.fastq.gz
│       ├── 31GEX_S4_L003_R1_001.fastq.gz
│       └── 31GEX_S4_L003_R2_001.fastq.gz

Pipeline step

  1. python3 01_batch_attachBarcodes.py to attach barcode to fastq file

    Please set BGI_prefix(V300079632 in our case) & lst before running the script

  2. python 3 02_run_convertHeaders.py to convert to cellranger input format

    Please set dict d before running the script

  3. bash 03_homosapien_cellranger.sh to run cellranger

Reference: https://github.com/powellgenomicslab/BGI_vs_Illumina_Benchmark https://github.com/IMB-Computational-Genomics-Lab/BGIvsIllumina_scRNASeq

bgi_to_illumina's People

Contributors

chengmingbo avatar

Watchers

 avatar

Recommend Projects

  • React photo React

    A declarative, efficient, and flexible JavaScript library for building user interfaces.

  • Vue.js photo Vue.js

    🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.

  • Typescript photo Typescript

    TypeScript is a superset of JavaScript that compiles to clean JavaScript output.

  • TensorFlow photo TensorFlow

    An Open Source Machine Learning Framework for Everyone

  • Django photo Django

    The Web framework for perfectionists with deadlines.

  • D3 photo 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.

  • Game

    Some thing interesting about game, make everyone happy.

Recommend Org

  • Facebook photo Facebook

    We are working to build community through open source technology. NB: members must have two-factor auth.

  • Microsoft photo Microsoft

    Open source projects and samples from Microsoft.

  • Google photo Google

    Google ❤️ Open Source for everyone.

  • D3 photo D3

    Data-Driven Documents codes.