Tag: Rsubread

Possible bugs in Rsubread/stad-alone featureCounts options fracOverlap and largestOverlap with fractional counts

Hi, running Rsubread 2.8.2/2.12.0 or featureCounts 2.0.3/2.0.1, I stumbled over two issues when allowing ambiguous read assignment (-O/allowMultiOverlap) 1) regarding assignment via minimum fractional overlap (–fracOverlap) using featureCounts stand-alone binary. 2) when combined with –/largestOverlap and –/fraction using Rsubread featureCounts function or the stand-alone binary. to 1) Assume a read…

Continue Reading Possible bugs in Rsubread/stad-alone featureCounts options fracOverlap and largestOverlap with fractional counts

Rsubread featureCounts outputs dozens of temp files, no counts

Rsubread featureCounts outputs dozens of temp files, no counts 1 @83165de1 Last seen 16 hours ago United States Hello, I am having trouble getting an output file in Rsubreads using featureCounts. I want to set up my data to run analysis of differential expresssion in EdgeR. I’m running about 40…

Continue Reading Rsubread featureCounts outputs dozens of temp files, no counts

Index of /Rsubread/annot

Name Last modified Size Description Parent Directory   –   chm13v2.0_RefSeqStri..> 2022-08-22 01:51 2.0M   chm13v2.0_RefSeq_12A..> 2022-08-22 01:51 2.5M   code/ 2022-08-22 01:53 –   hg38_RefSeqStrict_22..> 2022-08-22 01:51 2.3M   hg38_RefSeq_22Apr202..> 2022-08-22 01:51 3.0M   mm10_RefSeqStrict_22..> 2022-08-22 01:51 2.0M   mm10_RefSeq_22Jun202..> 2022-08-22 01:51 2.5M   mm39_RefSeqStrict_23..> 2022-08-22 01:51 1.9M  …

Continue Reading Index of /Rsubread/annot

Rsubread featurecounts

Rsubread featurecounts 1 Hi there, I seem to be getting this error when reading in a BAM file which was generated by PBMM2 align on pacbio data. I have tried to google the error message but there are no results. I wonder if anyone has ideas on what the error…

Continue Reading Rsubread featurecounts

Using featureCounts and downloading Rsubread

Using featureCounts and downloading Rsubread 1 @4769e097 Last seen 23 hours ago United Kingdom I am trying to perform a count per gene analysis using featureCounts in R. I have downloaded the gtf file and edited it within R to only contain the gene ID, chr, start, end, and strand,…

Continue Reading Using featureCounts and downloading Rsubread

Accepted r-bioc-rsubread 2.10.0-1 (source) into unstable

—–BEGIN PGP SIGNED MESSAGE—– Hash: SHA512 Format: 1.8 Date: Thu, 12 May 2022 14:00:29 +0200 Source: r-bioc-rsubread Architecture: source Version: 2.10.0-1 Distribution: unstable Urgency: medium Maintainer: Debian R Packages Maintainers <r-pkg-t…@alioth-lists.debian.net> Changed-By: Andreas Tille <ti…@debian.org> Changes: r-bioc-rsubread (2.10.0-1) unstable; urgency=medium . * Team upload. * New upstream version * Standards-Version:…

Continue Reading Accepted r-bioc-rsubread 2.10.0-1 (source) into unstable

Annotated file with gene ID (instead of gene symbol)

Annotated file with gene ID (instead of gene symbol) 0 @9cb59de3 Last seen 14 hours ago United States Hello, I am using “featureCounts” in Rsubread package for analyzing bulk RNA-seq of drosophila. Since there is no inbuilt annotations of drosophila, I am using a gtf file in the homepage of…

Continue Reading Annotated file with gene ID (instead of gene symbol)

Using Rsubread buildindex with GRCh37.p13.genome.fa.gz gives me an error

Using Rsubread buildindex with GRCh37.p13.genome.fa.gz gives me an error 0 @efernandez-22025 Last seen 1 day ago Argentina Hi I am triying to build the human index using ftp.ebi.ac.uk/pub/databases/gencode/Gencode_human/release_19/GRCh37.p13.genome.fa.gz I am using Rsubread 2.4.3 an it gives me the following error //================================= Running ==================================\ || || || Check the integrity of…

Continue Reading Using Rsubread buildindex with GRCh37.p13.genome.fa.gz gives me an error

Bioconductor – Rsubread

    This package is for version 2.13 of Bioconductor; for the stable, up-to-date release version, see Rsubread. Rsubread: high-performance read alignment, quantification and mutation discovery Bioconductor version: 2.13 This R package provides easy-to-use tools for analyzing next-gen sequencing read data. Functions of these tools include quality assessment, read alignment,…

Continue Reading Bioconductor – Rsubread

Accepted r-bioc-rsubread 2.8.2-1 (source) into unstable

—–BEGIN PGP SIGNED MESSAGE—– Hash: SHA512 Format: 1.8 Date: Mon, 21 Mar 2022 21:39:43 +0100 Source: r-bioc-rsubread Architecture: source Version: 2.8.2-1 Distribution: unstable Urgency: medium Maintainer: Debian R Packages Maintainers <r-pkg-t…@alioth-lists.debian.net> Changed-By: Andreas Tille <ti…@debian.org> Changes: r-bioc-rsubread (2.8.2-1) unstable; urgency=medium . * Team upload. * New upstream version * Maintainer:…

Continue Reading Accepted r-bioc-rsubread 2.8.2-1 (source) into unstable

Error in Rsubread featureCounts

Hi there, Excellent package! I am using it to do RNA-seq. But I encountered a small problem when using featureCounts(). The code is as follows: featureCounts( “A1.raw_1.fastq.gz.subjunc.BAM”, annot.inbuilt = NULL, annot.ext = “GCF_015227675.2_mRatBN7.2_genomic.gtf”, isGTFAnnotationFile=TRUE, isPairedEnd=TRUE, nthreads = 8 ) And it returns this: ========== _____ _ _ ____ _____ ______…

Continue Reading Error in Rsubread featureCounts

Aligning multiple single and paired-end reads from multiple files (lanes)

Rsubread: Aligning multiple single and paired-end reads from multiple files (lanes) 0 Hello, I am new to bioinformatics and looking for some help. I have 27 files from an Illumina output. There are 4 paired end and 23 single read files. I am trying to align them using Rsubread in…

Continue Reading Aligning multiple single and paired-end reads from multiple files (lanes)

“Paired-end reads were detected in single-end read library”

“Paired-end reads were detected in single-end read library” 0 @9cb59de3 Last seen 12 hours ago United States Hello, I am using “featureCounts” in Rsubread package for analyzing bulk RNA-seq of drosophila. Since there is no inbuilt annotations of drosophila, I am trying to use a gtf file in the homepage…

Continue Reading “Paired-end reads were detected in single-end read library”

Rsubread align maximum nthreads

Hi Experts, I am using Rsubread align using following comand- align (index=”my_index”, readfile1 = “SRR123456_1.fastq” ,readfile2= “SRR123456_2.fastq”, type=”rna”,input_format = “FASTQ”, minFragLength=35,maxFragLength=151,useAnnotation=”TRUE”, nthreads=64, annot.ext = “my_annotation.gtf.gz”, isGTF = “TRUE”, sortReadsByCoordinates = “TRUE”, output_format = “BAM”) here i have asigned 64 threads but in console, i see only 40 threads, I dont…

Continue Reading Rsubread align maximum nthreads