-
Notifications
You must be signed in to change notification settings - Fork 18
Description
I now have a bam file of the space transcriptome. The bar code in it is Cx and the Cy column, as shown in the figure

At the same time, there is also a maternal vcf file and a fetal vcf file. I now need to distinguish the maternal and fetal sources for each spot in the maternal fetal interface.
demuxlet --sam /work/home/luo_funong/Tan_Guang-hui/mutai_qufen/sorted.bam
--vcf mothers.vcf,fetal.vcf
--field GT
--out /work/home/luo_funong/Tan_Guang-hui/mutai_qufen/out
--alpha 0.1
--vcf-list
--tag-UMI null
--tag-group null
--tag-BC1 Cx
--tag-BC2 Cy
--threads 4
--mapq 255
--paired-end false
--per-cell true
--plot true
--plot-prefix /work/home/luo_funong/Tan_Guang-hui/mutai_qufen/out
--plot-format pdf
--ref /work/home/luo_funong/Tan_Guang-hui/mutai_qufen/ARS-UCD1.2_Btau5.0.1Y.fa

I have also tried more methods, but still reported an error. Can you help me solve this problem