ATACseqQC function fragSizeDist returns error
0
0
Entering edit mode
@xuyao15927402897-20826
Last seen 3.2 years ago
Germany

Species: wheat

There is no error message when using bam files of other species as input, The following is the error message.

Error in value[[3L]](cond) : 'end' must be <= 536870912
  file: ../mapping/CT_2_1.uni.dedup.bam
  index: ../mapping/CT_2_1.uni.dedup.bam
Calls: fragSizeDist ... tryCatch -> tryCatchList -> tryCatchOne -> <Anonymous>


# include your problematic code here with any corresponding output 
# please also include the results of running the following in an R session 

sessionInfo( )
R version 3.5.1 (2018-07-02)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: CentOS Linux 7 (Core)

Matrix products: default
BLAS: /PUBLIC/software/RNA/R/R-3.5.1/lib64/R/lib/libRblas.so
LAPACK: /PUBLIC/software/RNA/R/R-3.5.1/lib64/R/lib/libRlapack.so

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
 [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8    
 [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8   
 [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C            
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       

attached base packages:
[1] stats4    parallel  stats     graphics  grDevices utils     datasets 
[8] methods   base     

other attached packages:
[1] ATACseqQC_1.6.4     S4Vectors_0.20.1    BiocGenerics_0.28.0

loaded via a namespace (and not attached):
 [1] ProtGenerics_1.14.0           bitops_1.0-6                 
 [3] matrixStats_0.55.0            bit64_0.9-7                  
 [5] progress_1.2.2                httr_1.4.1                   
 [7] GenomeInfoDb_1.18.2           tools_3.5.1                  
 [9] R6_2.4.1                      rGADEM_2.30.0                
[11] KernSmooth_2.23-16            seqLogo_1.48.0               
[13] DBI_1.1.0                     lazyeval_0.2.2               
[15] colorspace_1.4-0              ade4_1.7-15                  
[17] motifStack_1.26.0             prettyunits_1.1.1            
[19] bit_1.1-15.2                  curl_4.3                     
[21] compiler_3.5.1                VennDiagram_1.6.20           
[23] graph_1.60.0                  Biobase_2.42.0               
[25] formatR_1.7                   DelayedArray_0.8.0           
[27] grImport_0.9-3                rtracklayer_1.42.2           
[29] scales_1.0.0                  randomForest_4.6-14          
[31] RBGL_1.58.2                   stringr_1.4.0                
[33] digest_0.6.18                 Rsamtools_1.34.1             
[35] XVector_0.22.0                pkgconfig_2.0.3              
[37] htmltools_0.4.0               fastmap_1.0.1                
[39] ensembldb_2.6.8               limma_3.38.3                 
[41] BSgenome_1.50.0               regioneR_1.14.0              
[43] htmlwidgets_1.5.1             rlang_0.4.8                  
[45] RSQLite_2.2.0                 shiny_1.4.0                  
[47] BiocParallel_1.16.6           RCurl_1.98-1.1               
[49] magrittr_1.5                  polynom_1.4-0                
[51] GO.db_3.7.0                   GenomeInfoDbData_1.2.0       
[53] futile.logger_1.4.3           Matrix_1.2-15                
[55] Rcpp_1.0.3                    munsell_0.5.0                
[57] stringi_1.4.6                 yaml_2.2.1                   
[59] MASS_7.3-51.1                 SummarizedExperiment_1.12.0  
[61] zlibbioc_1.28.0               AnnotationHub_2.14.5         
[63] grid_3.5.1                    blob_1.2.1                   
[65] promises_1.1.0                crayon_1.3.4                 
[67] lattice_0.20-38               Biostrings_2.50.2            
[69] splines_3.5.1                 multtest_2.38.0              
[71] GenomicFeatures_1.34.8        hms_0.5.3                    
[73] GenomicScores_1.6.0           MotIV_1.38.0                 
[75] GenomicRanges_1.34.0          seqinr_3.6-1                 
[77] biomaRt_2.38.0                futile.options_1.0.1         
[79] XML_3.98-1.16                 lambda.r_1.2.4               
[81] BiocManager_1.30.10           idr_1.2                      
[83] vctrs_0.3.4                   httpuv_1.5.2                 
[85] mime_0.9                      preseqR_4.0.0                
[87] xtable_1.8-4                  AnnotationFilter_1.6.0       
[89] later_1.0.0                   survival_3.1-8               
[91] ChIPpeakAnno_3.16.1           GenomicAlignments_1.18.1     
[93] AnnotationDbi_1.44.0          memoise_1.1.0                
[95] IRanges_2.16.0                interactiveDisplayBase_1.20.0

Thanks,

Yao

2021-8-28

fragSizeDist ATACseqQC • 555 views
ADD COMMENT

Login before adding your answer.

Traffic: 577 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6