Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

There are some errors running bam file #51

Open
SevenTea7 opened this issue Nov 13, 2024 · 14 comments
Open

There are some errors running bam file #51

SevenTea7 opened this issue Nov 13, 2024 · 14 comments

Comments

@SevenTea7
Copy link

Hi,

I conducted a test on a dataset and generated my BAM and BAI files using 'samtools'. However, I encountered some errors during the process. I would appreciate any solutions you might have. Additionally, please let me know if there is anything I may have overlooked while executing the workflow."
Best,
Yang

`

gcap.workflow(
tumourseqfile = "/data/liuhui/ysr/ascat/SRR8615228_sorted.sorted_PE.bam",
normalseqfile = "/data/liuhui/ysr/ascat/SRR8615244_sorted.sorted_PE.bam",
tumourname = "wes_395LC",
normalname = "wes_395N",
jobname = "wes_395",
outdir = "/data/liuhui/ysr/ascat",
allelecounter_exe = "/data/liuhui/anaconda3/envs/ysr_cancerit/bin/alleleCounter",
g1000allelesprefix = file.path(
"/data/liuhui/ysr/app/gcap-master/data/1000G_loci_hg38/",
"1kg.phase3.v5a_GRCh38nounref_allele_index_chr"
),
g1000lociprefix = file.path("/data/liuhui/ysr/app/gcap-master/data/1000G_loci_hg38/",
"1kg.phase3.v5a_GRCh38nounref_loci_chrstring_chr"
),
GCcontentfile = "/data/liuhui/ysr/app/gcap-master/data/GC_correction_hg38.txt",
replictimingfile = "/data/liuhui/ysr/app/gcap-master/data/RT_correction_hg38.txt",
skip_finished_ASCAT = TRUE,
skip_ascat_call = FALSE,
result_file_prefix = "wes_395",
genome_build = "hg38",
model = "XGB11"
)
`

`
<gcap> 2024-11-13 14:30:36 info [gcap.workflow]: =====================
<gcap> 2024-11-13 14:30:36 info [gcap.workflow]:    GCAP WORKFLOW
<gcap> 2024-11-13 14:30:36 info [gcap.workflow]: =====================
<gcap> 2024-11-13 14:30:36 info [gcap.workflow]: 
<gcap> 2024-11-13 14:30:36 info [gcap.workflow]: =====================
<gcap> 2024-11-13 14:30:36 info [gcap.workflow]: Step 1: Run ASCAT 3.0
<gcap> 2024-11-13 14:30:36 info [gcap.workflow]: =====================
<gcap> 2024-11-13 14:30:36 info [gcap.runASCAT]: > Run ASCAT on WES data <
<gcap> 2024-11-13 14:30:36 info [gcap.runASCAT]: 
<gcap> 2024-11-13 14:30:36 info [gcap.runASCAT]: Configs:
<gcap> 2024-11-13 14:30:36 info [gcap.runASCAT]:   result path set to /data/liuhui/ysr/ascat
<gcap> 2024-11-13 14:30:36 info [gcap.runASCAT]:   allelecounter_exe set to /data/liuhui/anaconda3/envs/ysr_cancerit/bin/alleleCounter
<gcap> 2024-11-13 14:30:36 info [gcap.runASCAT]:   g1000allelesprefix set to /data/liuhui/ysr/app/gcap-master/data/1000G_loci_hg38//1kg.phase3.v5a_GRCh38nounref_allele_index_chr
<gcap> 2024-11-13 14:30:36 info [gcap.runASCAT]:   g1000lociprefix set to /data/liuhui/ysr/app/gcap-master/data/1000G_loci_hg38//1kg.phase3.v5a_GRCh38nounref_loci_chrstring_chr
<gcap> 2024-11-13 14:30:36 info [gcap.runASCAT]:   GCcontentfile set to /data/liuhui/ysr/app/gcap-master/data/GC_correction_hg38.txt
<gcap> 2024-11-13 14:30:36 info [gcap.runASCAT]:   replictimingfile set to /data/liuhui/ysr/app/gcap-master/data/RT_correction_hg38.txt
<gcap> 2024-11-13 14:30:36 info [gcap.runASCAT]:   nthreads set to 22
<gcap> 2024-11-13 14:30:36 info [gcap.runASCAT]:   minCounts set to 10
<gcap> 2024-11-13 14:30:36 info [gcap.runASCAT]:   BED_file set to NA
<gcap> 2024-11-13 14:30:36 info [gcap.runASCAT]:   probloci_file set to NA
<gcap> 2024-11-13 14:30:36 info [gcap.runASCAT]:   chrom_names set to <1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22>
<gcap> 2024-11-13 14:30:36 info [gcap.runASCAT]:   gender set to <XX>
<gcap> 2024-11-13 14:30:36 info [gcap.runASCAT]:   min_base_qual set to 20
<gcap> 2024-11-13 14:30:36 info [gcap.runASCAT]:   min_map_qual set to 35
<gcap> 2024-11-13 14:30:36 info [gcap.runASCAT]:   penalty set to 70
<gcap> 2024-11-13 14:30:36 info [gcap.runASCAT]:   skip_finished_ASCAT set to TRUE
<gcap> 2024-11-13 14:30:36 info [gcap.runASCAT]: 1 jobs detected
<gcap> 2024-11-13 14:30:36 info [gcap.runASCAT]: No ASCAT job to skip.
<gcap> 2024-11-13 14:30:36 info [FUN]: start submitting job wes_395
<gcap> 2024-11-13 14:30:36 info [FUN]:      tumor data file: /data/liuhui/ysr/ascat/SRR8615228_sorted.sorted_PE.bam
<gcap> 2024-11-13 14:30:36 info [FUN]:     normal data file: /data/liuhui/ysr/ascat/SRR8615244_sorted.sorted_PE.bam
<gcap> 2024-11-13 14:30:36 info [FUN]:    tumor sample name: wes_395LC
<gcap> 2024-11-13 14:30:36 info [FUN]:   normal sample name: wes_395N
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615228_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615228_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615228_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615228_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615228_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615228_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615228_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615228_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615228_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615228_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615228_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615228_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615228_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615228_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615228_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615228_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615228_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615228_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615228_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615228_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615228_sorted.sorted_PE.bam.bai
Reading locis
Done reading locis
Multi pos start:
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615228_sorted.sorted_PE.bam.bai
Reading locis
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615244_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615244_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615244_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615244_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615244_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615244_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615244_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615244_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615244_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615244_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615244_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615244_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615244_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615244_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615244_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615244_sorted.sorted_PE.bam.bai
Reading locis
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615244_sorted.sorted_PE.bam.bai
Done reading locis
Multi pos start:
Reading locis
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615244_sorted.sorted_PE.bam.bai
Reading locis
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615244_sorted.sorted_PE.bam.bai
Reading locis
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615244_sorted.sorted_PE.bam.bai
Done reading locis
Multi pos start:
Reading locis
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615244_sorted.sorted_PE.bam.bai
Reading locis
Done reading locis
Multi pos start:
[W::hts_idx_load3] The index file is older than the data file: /data/liuhui/ysr/ascat/SRR8615244_sorted.sorted_PE.bam.bai
Reading locis
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
Done reading locis
Multi pos start:
[1] Reading Tumor LogR data...
[1] Reading Tumor BAF data...
[1] Reading Germline LogR data...
[1] Reading Germline BAF data...
[1] Registering SNP locations...
[1] Splitting genome in distinct chunks...
                                                                                                                                                                                                                                                                                                                                                                                    [1] Sample wes_395LC (1/1)
GC correlation:  25bp 0.03909 ; 50bp 0.04406 ; 100bp 0.05095 ; 200bp 0.06052 ; 500bp 0.07683 ; 1kb 0.08195 ; 2kb 0.08014 ; 5kb 0.07102 ; 10kb 0.06518 ; 20kb 0.06307 ; 50kb 0.06000 ; 100kb 0.05410 ; 200kb 0.04864 ; 500kb 0.04306 ; 1Mb 0.04496 ; 2Mb 0.04315 ; 5Mb 0.00028 ; 10Mb 0.02476 ; 
Short window size:  1kb 
Long window size:  5kb 
Replication timing correlation:  Bg02es 0.036 ; Bj 0.071 ; Gm06990 0.049 ; Gm12801 0.050 ; Gm12812 0.048 ; Gm12813 0.049 ; Gm12878 0.052 ; Helas3 0.069 ; Hepg2 0.025 ; Huvec 0.059 ; Imr90 0.062 ; K562 0.044 ; Mcf7 0.011 ; Nhek 0.063 ; Sknsh 0.016 ; 
Replication dataset:  Bj 
[1] Plotting tumor data
[1] Plotting germline data
[1] Sample wes_395LC (1/1)
[1] Sample wes_395LC (1/1)
<gcap> 2024-11-13 14:38:14 info [doTryCatch]: job wes_395 done
<gcap> 2024-11-13 14:38:14 info [gcap.runASCAT]: ASCAT analysis done, check /data/liuhui/ysr/ascat for results
<gcap> 2024-11-13 14:38:14 info [gcap.workflow]: checking ASCAT result files
<gcap> 2024-11-13 14:38:14 warn [FUN]: /data/liuhui/ysr/ascat/wes_395.ASCAT.rds contains a failed ASCAT job, will discard it before next step
<gcap> 2024-11-13 14:38:14 fatal [gcap.workflow]: no sucessful ASCAT result file to proceed!
<gcap> 2024-11-13 14:38:14 fatal [gcap.workflow]: check your ASCAT setting before make sure this case could not be used!
Error in gcap.workflow(tumourseqfile = "/data/liuhui/ysr/ascat/SRR8615228_sorted.sorted_PE.bam",  : 
  
In addition: Warning message:
In runASCAT(lrr, baf, lrrsegm, bafsegm, ASCATobj$gender[arraynr],  :
  ASCAT could not find an optimal ploidy and cellularity value for sample wes_395LC.
`

Session info


- [ ] > `
- [ ] > > devtools::session_info()
- [ ] > ─ Session info ────────────────────────────────────────────────────────────────────────────
- [ ] >  setting  value
- [ ] >  version  R version 4.1.3 (2022-03-10)
- [ ] >  os       Ubuntu 22.04.4 LTS
- [ ] >  system   x86_64, linux-gnu
- [ ] >  ui       X11
- [ ] >  language (EN)
- [ ] >  collate  en_US.UTF-8
- [ ] >  ctype    en_US.UTF-8
- [ ] >  tz       Asia/Shanghai
- [ ] >  date     2024-11-13
- [ ] >  pandoc   NA
- [ ] > ─ Packages ────────────────────────────────────────────────────────────────────────────────
- [ ] >  package          * version   date (UTC) lib source
- [ ] >  ASCAT            * 3.0.0     2024-11-11 [1] Github (ShixiangWang/ascat@51fd695)
- [ ] >  bigmemory        * 4.6.4     2024-01-09 [1] CRAN (R 4.1.3)
- [ ] >  bigmemory.sri      0.1.8     2024-01-10 [1] CRAN (R 4.1.3)
- [ ] >  Biobase          * 2.54.0    2021-10-26 [1] Bioconductor
- [ ] >  BiocGenerics     * 0.40.0    2021-10-26 [1] Bioconductor
- [ ] >  BiocManager        1.30.25   2024-08-28 [1] CRAN (R 4.1.3)
- [ ] >  bit                4.5.0     2024-09-20 [1] CRAN (R 4.1.3)
- [ ] >  bit64              4.5.2     2024-09-22 [1] CRAN (R 4.1.3)
- [ ] >  bitops             1.0-9     2024-10-03 [1] CRAN (R 4.1.3)
- [ ] >  brio               1.1.5     2024-04-24 [1] CRAN (R 4.1.3)
- [ ] >  cachem             1.1.0     2024-05-16 [1] CRAN (R 4.1.3)
- [ ] >  Cairo              1.6-2     2023-11-28 [1] CRAN (R 4.1.3)
- [ ] >  cli                3.6.3     2024-06-21 [1] CRAN (R 4.1.3)
- [ ] >  cluster            2.1.6     2023-12-01 [1] CRAN (R 4.1.3)
- [ ] >  codetools          0.2-20    2024-03-31 [1] CRAN (R 4.1.3)
- [ ] >  colorspace         2.1-1     2024-07-26 [1] CRAN (R 4.1.3)
- [ ] >  crayon             1.5.3     2024-06-20 [1] CRAN (R 4.1.3)
- [ ] >  curl               6.0.0     2024-11-05 [1] CRAN (R 4.1.3)
- [ ] >  data.table         1.16.2    2024-10-10 [1] CRAN (R 4.1.3)
- [ ] >  desc               1.4.3     2023-12-10 [1] CRAN (R 4.1.3)
- [ ] >  devtools           2.4.5     2022-10-11 [1] CRAN (R 4.1.3)
- [ ] >  digest             0.6.37    2024-08-19 [1] CRAN (R 4.1.3)
- [ ] >  doParallel       * 1.0.17    2022-02-07 [1] CRAN (R 4.1.3)
- [ ] >  dplyr              1.1.4     2023-11-17 [1] CRAN (R 4.1.3)
- [ ] >  ellipsis           0.3.2     2021-04-29 [1] CRAN (R 4.1.3)
- [ ] >  fansi              1.0.6     2023-12-08 [1] CRAN (R 4.1.3)
- [ ] >  fastmap            1.2.0     2024-05-15 [1] CRAN (R 4.1.3)
- [ ] >  foreach          * 1.5.2     2022-02-02 [1] CRAN (R 4.1.3)
- [ ] >  fs                 1.6.5     2024-10-30 [1] CRAN (R 4.1.3)
- [ ] >  furrr              0.3.1     2022-08-15 [1] CRAN (R 4.1.3)
- [ ] >  future             1.34.0    2024-07-29 [1] CRAN (R 4.1.3)
- [ ] >  gcap             * 1.2.0     2024-11-11 [1] https://shixiangwang.r-universe.dev (R 4.1.3)
- [ ] >  generics           0.1.3     2022-07-05 [1] CRAN (R 4.1.3)
- [ ] >  GenomeInfoDb     * 1.30.1    2022-01-30 [1] Bioconductor
- [ ] >  GenomeInfoDbData   1.2.7     2024-11-08 [1] Bioconductor
- [ ] >  GenomicRanges    * 1.46.1    2021-11-18 [1] Bioconductor
- [ ] >  GetoptLong         1.0.5     2020-12-15 [1] CRAN (R 4.1.3)
- [ ] >  ggplot2            3.5.1     2024-04-23 [1] CRAN (R 4.1.3)
- [ ] >  GlobalOptions      0.1.2     2020-06-10 [1] CRAN (R 4.1.3)
- [ ] >  globals            0.16.3    2024-03-08 [1] CRAN (R 4.1.3)
- [ ] >  glue               1.8.0     2024-09-30 [1] CRAN (R 4.1.3)
- [ ] >  gridBase           0.4-7     2014-02-24 [1] CRAN (R 4.1.3)
- [ ] >  gtable             0.3.6     2024-10-25 [1] CRAN (R 4.1.3)
- [ ] >  hms                1.1.3     2023-03-21 [1] CRAN (R 4.1.3)
- [ ] >  htmltools          0.5.8.1   2024-04-04 [1] CRAN (R 4.1.3)
- [ ] >  htmlwidgets        1.6.4     2023-12-06 [1] CRAN (R 4.1.3)
- [ ] >  httpuv             1.6.15    2024-03-26 [1] CRAN (R 4.1.3)
- [ ] >  IRanges          * 2.28.0    2021-10-26 [1] Bioconductor
- [ ] >  iterators        * 1.0.14    2022-02-05 [1] CRAN (R 4.1.3)
- [ ] >  jsonlite           1.8.9     2024-09-20 [1] CRAN (R 4.1.3)
- [ ] >  later              1.3.2     2023-12-06 [1] CRAN (R 4.1.3)
- [ ] >  lattice            0.22-6    2024-03-20 [1] CRAN (R 4.1.3)
- [ ] >  lgr                0.4.4     2022-09-05 [1] CRAN (R 4.1.3)
- [ ] >  lifecycle          1.0.4     2023-11-07 [1] CRAN (R 4.1.3)
- [ ] >  listenv            0.9.1     2024-01-29 [1] CRAN (R 4.1.3)
- [ ] >  magrittr           2.0.3     2022-03-30 [1] CRAN (R 4.1.3)
- [ ] >  Matrix             1.6-1.1   2023-09-18 [1] CRAN (R 4.1.3)
- [ ] >  memoise            2.0.1     2021-11-26 [1] CRAN (R 4.1.3)
- [ ] >  mime               0.12      2021-09-28 [1] CRAN (R 4.1.3)
- [ ] >  miniUI             0.1.1.1   2018-05-18 [1] CRAN (R 4.1.3)
- [ ] >  munsell            0.5.1     2024-04-01 [1] CRAN (R 4.1.3)
- [ ] >  NMF                0.28      2024-08-22 [1] CRAN (R 4.1.3)
- [ ] >  parallelly         1.39.0    2024-11-07 [1] CRAN (R 4.1.3)
- [ ] >  pillar             1.9.0     2023-03-22 [1] CRAN (R 4.1.3)
- [ ] >  pkgbuild           1.4.5     2024-10-28 [1] CRAN (R 4.1.3)
- [ ] >  pkgconfig          2.0.3     2019-09-22 [1] CRAN (R 4.1.3)
- [ ] >  pkgload            1.4.0     2024-06-28 [1] CRAN (R 4.1.3)
- [ ] >  plyr               1.8.9     2023-10-02 [1] CRAN (R 4.1.3)
- [ ] >  profvis            0.4.0     2024-09-20 [1] CRAN (R 4.1.3)
- [ ] >  promises           1.3.0     2024-04-05 [1] CRAN (R 4.1.3)
- [ ] >  purrr              1.0.2     2023-08-10 [1] CRAN (R 4.1.3)
- [ ] >  quadprog           1.5-8     2019-11-20 [1] CRAN (R 4.1.3)
- [ ] >  R6                 2.5.1     2021-08-19 [1] CRAN (R 4.1.3)
- [ ] >  rappdirs           0.3.3     2021-01-31 [1] CRAN (R 4.1.3)
- [ ] >  RColorBrewer     * 1.1-3     2022-04-03 [1] CRAN (R 4.1.3)
- [ ] >  Rcpp               1.0.13-1  2024-11-02 [1] CRAN (R 4.1.3)
- [ ] >  RCurl              1.98-1.16 2024-07-11 [1] CRAN (R 4.1.3)
- [ ] >  readr            * 2.1.5     2024-01-10 [1] CRAN (R 4.1.3)
- [ ] >  registry           0.5-1     2019-03-05 [1] CRAN (R 4.1.3)
- [ ] >  remotes            2.5.0     2024-03-17 [1] CRAN (R 4.1.3)
- [ ] >  reshape2           1.4.4     2020-04-09 [1] CRAN (R 4.1.3)
- [ ] >  rjson              0.2.23    2024-09-16 [1] CRAN (R 4.1.3)
- [ ] >  rlang              1.1.4     2024-06-04 [1] CRAN (R 4.1.3)
- [ ] >  rngtools           1.5.2     2021-09-20 [1] CRAN (R 4.1.3)
- [ ] >  S4Vectors        * 0.32.4    2022-03-24 [1] Bioconductor
- [ ] >  scales             1.3.0     2023-11-28 [1] CRAN (R 4.1.3)
- [ ] >  sessioninfo        1.2.2     2021-12-06 [1] CRAN (R 4.1.3)
- [ ] >  shiny              1.9.1     2024-08-01 [1] CRAN (R 4.1.3)
- [ ] >  sigminer         * 2.3.2     2024-11-11 [1] https://shixiangwang.r-universe.dev (R 4.1.3)
- [ ] >  stringi            1.8.4     2024-05-06 [1] CRAN (R 4.1.3)
- [ ] >  stringr            1.5.1     2023-11-14 [1] CRAN (R 4.1.3)
- [ ] >  testthat         * 3.2.1.1   2024-04-14 [1] CRAN (R 4.1.3)
- [ ] >  tibble             3.2.1     2023-03-20 [1] CRAN (R 4.1.3)
- [ ] >  tidyr              1.3.1     2024-01-24 [1] CRAN (R 4.1.3)
- [ ] >  tidyselect         1.2.1     2024-03-11 [1] CRAN (R 4.1.3)
- [ ] >  tzdb               0.4.0     2023-05-12 [1] CRAN (R 4.1.3)
- [ ] >  urlchecker         1.0.1     2021-11-30 [1] CRAN (R 4.1.3)
- [ ] >  usethis            3.0.0     2024-07-29 [1] CRAN (R 4.1.3)
- [ ] >  utf8               1.2.4     2023-10-22 [1] CRAN (R 4.1.3)
- [ ] >  uuid               1.2-1     2024-07-29 [1] CRAN (R 4.1.3)
- [ ] >  vctrs              0.6.5     2023-12-01 [1] CRAN (R 4.1.3)
- [ ] >  vroom              1.6.5     2023-12-05 [1] CRAN (R 4.1.3)
- [ ] >  withr              3.0.2     2024-10-28 [1] CRAN (R 4.1.3)
- [ ] >  xgboost            1.5.2.1   2022-02-21 [1] CRAN (R 4.1.3)
- [ ] >  xtable             1.8-4     2019-04-21 [1] CRAN (R 4.1.3)
- [ ] >  XVector            0.34.0    2021-10-26 [1] Bioconductor
- [ ] >  zlibbioc           1.40.0    2021-10-26 [1] Bioconductor
- [ ] >  [1] /data/liuhui/anaconda3/envs/ysr_cancerit/lib/R/library
- [ ] > ─────────────────────────────────────────────────────────────
- [ ] > `

@ShixiangWang
Copy link
Owner

ShixiangWang commented Nov 13, 2024

It seems the ASCAT failed to get results, you can run https://github.com/VanLoo-lab/ascat analysis (it has been maintained and updated) by your own and check this.

@SevenTea7
Copy link
Author

It seems the ASCAT failed to get results, you can run https://github.com/VanLoo-lab/ascat analysis (it has been maintained and updated) by your own and check this.

Hello, could you kindly provide me with two test BAM files? I would like to determine whether the issue lies with my data or if there is a problem with the installation of my software.

@ShixiangWang
Copy link
Owner

I generated a debug path for testing the pipeline, please see https://github.com/ShixiangWang/gcap/tree/master/test-workflow/debug

@SevenTea7
Copy link
Author

I generated a debug path for testing the pipeline, please see https://github.com/ShixiangWang/gcap/tree/master/test-workflow/debug
Hello,
I sincerely appreciate the folder you provided. However, when I attempted to test and execute this program, I encountered the data ERR5243012. Could you kindly share the download link for the required data?

@ShixiangWang
Copy link
Owner

You can refer to goole to see how to download data from SRA, it's quite easy.

@ShixiangWang
Copy link
Owner

Also, if you have multiple samples, you can use other sample to test the workflow.

@SevenTea7
Copy link
Author

Also, if you have multiple samples, you can use other sample to test the workflow.

Hello, thank you very much for your help, I have solved the problem of the program. But now I have a guess, whether I can use inferCNV and other software to analyze the scRNA-seq, obtain the corresponding copy number data, and build a table similar to "ascn" for eccDNA prediction.Thank you very much for your help.

@ShixiangWang
Copy link
Owner

ShixiangWang commented Nov 17, 2024

In logic, it makes sense, however, the single cell data cannot provide accurate absolute copy number profiles, so it will largely fails using gcap, as the copy number is a dominant predictive feature. I beleive there some ongoing developments for sc data, I once wanted to do a job similar to this, but have not been able to do so yet. You have to take a risk for using it, and you'd better add the necessary filtering and verification to the results.

@SevenTea7
Copy link
Author

Hello, I would like to ask, I want to use the data of cancer cell lines to run, and I need to pair normal samples, how should I choose this normal sample, such as the para-cancer data corresponding to cancer, the sequencing data corresponding to tissue, or the data of embryonic stem cells? Do you have any recommendations.Thanks.

@ShixiangWang
Copy link
Owner

For cancer cell line data, it would be better to use WGS. However, if you only got WES, just run the cnvkit to get absolute copy number values and use the integer copy number to infer the results (the accuracy should be lower). On the other hand, if you have your own samples, sequencing with Circle-seq is a good option.

@SevenTea7
Copy link
Author

For cancer cell line data, it would be better to use WGS. However, if you only got WES, just run the cnvkit to get absolute copy number values and use the integer copy number to infer the results (the accuracy should be lower). On the other hand, if you have your own samples, sequencing with Circle-seq is a good option.

Thank you very much for your advice. My data is collected as WGS data, but I'm not sure which normal samples of cancer cell lines to use for comparison. For instance, now that I have obtained the data of the gastric cancer cell line, can I use the data of normal gastric tissue, normal gastric cell culture, or embryonic stem cells as its normal group data? Which data is more in line with practical theory?

@ShixiangWang
Copy link
Owner

ShixiangWang commented Nov 28, 2024

For WGS data, you can use the AmpliconSuite, it uses tumor WGS as input, no normal control is required.

@SevenTea7
Copy link
Author

SevenTea7 commented Dec 23, 2024 via email

@ShixiangWang
Copy link
Owner

Hello, when I used AA for ecc recognition of WGS data, I found that the results obtained were very few, even 0, which was very difficult to analyze in bioinformatics. Do you have any good suggestions for using the software? Thank you very much.

------------------ 原始邮件 ------------------ 发件人: "ShixiangWang/gcap" @.>; 发送时间: 2024年11月28日(星期四) 上午10:35 @.>; @.@.>; 主题: Re: [ShixiangWang/gcap] There are some errors running bam file (Issue #51) For WGS data, you can use the AmpliconSuite, it uses tumor WGS as input, no normal control is required. — Reply to this email directly, view it on GitHub, or unsubscribe. You are receiving this because you commented.Message ID: @.***>

You can tune the parameters of AA.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants