Entering edit mode
2.0 years ago
shaoqian ma
•
0
I mapped the 10x data (v3.1) using STARsolo but the strange thing is that output sam file does not have the tags including CR CY UR UY.... (STAR 2.7.3a)
and here is my setting:
STAR --runThreadN 10 --genomeDir mm10_STARidx_GTF --readFilesIn ../big_R2.fastq.gz ../big_R1.fastq.gz --outSAMattributes NH HI nM AS CR CY UR UY --outFilterMultimapNmax 100 --winAnchorMultimapNmax 100 --outMultimapperOrder Random \ --readFilesCommand zcat --soloType CB_UMI_Simple --soloCBwhitelist 3M-february-2018.txt \ --runRNGseed 777 --outSAMmultNmax 1 --soloUMIlen 12 --soloStrand Forward --soloFeatures Gene GeneFull Velocyto \ --soloOutFileNames solo.out/ genes.tsv barcodes.tsv matrix.mtx --soloBarcodeReadLength 0 --soloCBlen 16 --limitOutSJcollapsed 10000000 --limitIObufferSize=1500000000
I can not find anything wrong with my setting, I plan to use the STARsolo output to quantify TE expression for my 10x single cell data using the scTE tool Identifying transposable element expression dynamics and heterogeneity during development at the single-cell level with a processing pipeline scTE , and it requires bam file with CR tag
I have solved this problem.
I changed this as a comment. Also, could you add how you solved the issue?