How to solve the problem related to unmatched chromosome identifiers?
While using Featurecounts tool, I am getting zero counts for RNA seq reads, It can be due to unmatched chromosome identifiers in BAM and GTF file. Can someone please help me out in making the identifiers same for both the file?
Thank You
• 1,457 views
•
link
updated
by
Biostar
•
written
by
priya.saxena •
•
written
by
priya.saxena •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
How to create a command that matches 13 characters in each column to be extracted into a new file, …
written by MetagenomicsPlease •
-
Help with ERROR: Paired-end reads were detected in single-end read library
written by Chris •
-
featureCounts strandedness option
written by Noah •
-
Transcriptome assembly using stringtie and stringtie-merge
written by priya.saxena •
-
Formatting BAM file and GTF file for feature counts
written by abhilashtripathi10 •
-
featureCounts on antisense genes
written by till.schwaemmle •
-
counting reads for genes
written by kemuntobilia •
-
Error in piRNA gtf file while featurecounts step in STAR aligner?
written by geethapriyanka7 •
-
Very low percentage of assigned reads using featureCounts
written by Meghdad •
-
About the ERCC spike-in sequence reads quantification
written by younglin113
How can this be? Usually one uses a reference genome that matches the GTF. What species is this and where are the files from?
Yes, I used the reference genome and then I did the mapping with Bowtie2 to generate a BAM file. This BAM file and the gtf file from NCBI RefSeq was used as input for Featurecounts. The species used is Desulfovibrio alaskensis G20 from NCBI RefSeq. The FTP link is from where I uploaded GTF and genome sequence to galaxy is the following
https://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/000/012/665/GCF_000012665.1_ASM1266v1/
Which files did you exactly use?
Try using this for genome: GCF_000012665.1_ASM1266v1_genomic.fna.gz
And this for gtf: GCF_000012665.1_ASM1266v1_genomic.gtf.gz
Both are coming from: https://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/000/012/665/GCF_000012665.1_ASM1266v1/
This should work, good luck!
Thank You so much I will try this and let you know.