inquiry about Loading scATAC-seq matrices into R

[jiangzh-coder]

Hi

i want to analyze some scATAC-seq data. And after unzip, i got 10 folders (1 patient per folder). In folder, there are 2 subfolders in whcih scRNAseq data and scATAC-seq data exist seperately. Within these 2 subfolders, there are files generated by cell ranger ( i attached 2 pic.)
How may i loading scATAC-seq matrices as well as scRNAseq data into R? Could you kindly provide some codes?

[jiangzh-coder]

jzhou@jiang:/mnt/d/##files/#atac/HD2$ tail atac_peaks.bed
GL000218.1 83275 84106
KI270726.1 27131 28058
KI270726.1 41490 42368
KI270711.1 7979 8731
KI270711.1 8887 9380
KI270713.1 15800 16459
KI270713.1 17290 18022
KI270713.1 21445 22340
KI270713.1 32734 33486
KI270713.1 36862 37780

jzhou@jiang:/mnt/d/##files/#atac/HD2$ head atac_peak_annotation.tsv
chrom start end gene distance peak_type
chr1 9778 10667 MIR1302-2HG -18887 distal
chr1 180732 181004 AL627309.5 -6871 distal
chr1 181116 181809 AL627309.5 -7255 distal
chr1 183935 184770 AL627309.5 -10074 distal
chr1 191103 192028 AL627309.5 -17242 distal
chr1 267627 268482 AP006222.2 773 distal
chr1 629498 630379 AC114498.1 41870 distal
chr1 633577 634503 AC114498.1 45949 distal
chr1 778280 779196 LINC01409 0 promoter

jzhou@jiang:/mnt/d/##files/#atac/HD2$ tail -5 atac_fragments.tsv
KI270713.1 39073 39449 TAAGTAGCACAGGATG-1 1
KI270713.1 39075 39208 GCTTAACAGTTCCCGT-1 2
KI270713.1 39697 39805 TGTTATGAGGGCTTTG-1 2
KI270713.1 40366 40677 GACCTAAGTTCCGGCT-1 2
KI270713.1 40639 40708 GGATGGCCACACCAAC-1 2

[jiangzh-coder]

Hello,

i analyzed upstream data in this dataset (GSE199994) . i download SRA file from EBI, then i change then into fastq file, then i transform their names to standard name and run cellranger_atac. i got error as following:

2.7% (< 10%) of read pairs have a valid 10x barcode. This could be a result of poor sequencing quality,
a sample mixup, or running the wrong pipeline, for example, running cellranger-atac on Multiome ATAC + GEX data, or vice versa.

The whole code is as following :

ascp -QT -l 300m -P33001 -i ~/miniconda3/envs/my10x/etc/asperaweb_id_dsa.openssh \era-fasp@fasp.sra.ebi.ac.uk:/vol1/srr/SRR186/006/SRR18613306 .

mv SRR18613306 SRR18613306.sra
parallel-fastq-dump -t 12 -O ./ --split-files --gzip -s SRR18613306.sra

mv *.fastq.gz /home/ubuntu/GSE199994/scATAC/2.raw_fastq

mv SRR18613295_1.fastq.gz SRR18613295-P5_S9_L001_I1_001.fastq.gz
mv SRR18613295_2.fastq.gz SRR18613295-P5_S9_L001_R1_001.fastq.gz
mv SRR18613295_3.fastq.gz SRR18613295-P5_S9_L001_R2_001.fastq.gz
mv SRR18613295_4.fastq.gz SRR18613295-P5_S9_L001_R3_001.fastq.gz

cellranger-atac count --id=SRR18613295-P5
--reference=/home/ubuntu/biosoftware/refdata-cellranger-arc-GRCh38-2020-A-2.0.0
--fastqs=/home/ubuntu/GSE199994/scATAC/2.raw_fastq
--sample=SRR18613295-P5
--localcores=24
--localmem=96

could you help me?