A script to reformat FASTQ files from ASAP-seq for downstream processing.
Note that it was originally designed for kite (kallisto | bustools), but can be used for something else as well.
docker run -it --rm \
-v $(pwd)/test:/data \
-w /data \
asap2kite:0.0.1 \
--fastqs /tmp/asap_to_kite-0.0.1/test/data1/ \
--sample test1 \
--id asapseq-test \
--cores 2
./build.sh
Either you can use the docker push
command or run push.sh
(requires SCING):
./push.sh
Illumina demultiplexed:
R1
: UMI / UBI (10 nts) at position 0R2
: CB (16 nts)R3
: HTO (15 nts)I
: Illumina indexing
After reformatted:
R1
: CB (16nts) + UMI (10 nts)R2
: HTO (15 nts) at position 0
Conversion code: https://github.com/hisplan/asap_to_kite/blob/master/asap_to_kite_v2.py#L122
R2
from sequencing must have exactly 16 nts. It appears that >16 could cause issues during the reformatting.