shendurelab / MPRAflow

A portable, flexible, parallelized tool for complete processing of massively parallel reporter assay data
Apache License 2.0
31 stars 16 forks source link

Problems with barcode association analysis #71

Closed clytiaxu closed 2 years ago

clytiaxu commented 2 years ago

Hi there, We prepared the CRS plasmid library followed your nature protocoland. The library was then subjected to MGISeq 2000 to get the CRS-barcode association (CRS length=220bp, complexity=8Million, PE150) . I run the MPRAflow association.nf followed and found that only very few reads () were kept after join the Read1 and Read2 fastq files (fastq-join -R1 -R2 -o xxx.joint). I am confused about why we need to join the R1 and R2 and then subject to bwa alignment. Is this step critical for the correct assignment of barcode to CRS? Many thanks, Clytia