problem running featurecounts using hg38 in-built annotation
0
0
Entering edit mode
raf4 ▴ 30
@raf4-8249
Last seen 23 months ago
United States

Dear List,

I am having trouble running feature counts as implemented in Rsubread on my mac

My script is:

library(Rsubread)

b1.res<-featureCounts(files="b1.res.bam",annot.inbuilt="hg38",isPairedEnd=TRUE,requireBothEndsMapped=FALSE)

I get an output file with numbers instead of gene names and no counts.

I have 2 warning messages.

WARNING: reads from the same pair were found not adjacent to each       ||
||             other in the input (due to read sorting by location or         ||
||             reporting of multi-mapping read pairs).                        ||
||                                                                            ||
||    Read re-ordering is performed.           

and

warning message:
package ‘Rsubread’ was built under R version 3.3.3

(I am using 3.3.2)

I am not reproducing the whole output file because I get a

"language "af" is not one of the supported languages ['en']!"

from the site.

Thanks and best wishes,

Rich

Richard Friedman

 

 

 

featurecounts hg38 • 1.7k views
ADD COMMENT
0
Entering edit mode

I believe I found the problem. I should  not have used the patched sequence.

I am currently downloading


ftp://ftp.ncbi.nlm.nih.gov/genomes/genbank/vertebrate_mammalian/Homo_sapiens/all_assembly_versions/GCA_000001405.15_GRCh38/seqs_for_alignment_pipelines.ucsc_ids/ GCA_000001405.15_GRCh38_no_alt_analysis_set.fna.gz

If this is not the right version of the sequence, please let me know.

Thanks and best wishes,

Rich

 

 

ADD REPLY
0
Entering edit mode

Hi Rich, could you please provide the complete output from your featureCounts run? Also you should update your copy of Rsubread and Bioc to their latest version as well.

ADD REPLY

Login before adding your answer.

Traffic: 907 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