Home > Community > Running htseq-count over BAM files
Upvote

26

Downvote
+ Bioinformatics
+ Rna-seq
Posted by
Monterino Overson

Running htseq-count over BAM files

Dave Aheto  Follow

If your BAM file is using UCSC chromosome names (chr1, chr2, etc.) and your GTF file is using Ensembl chromosome names (1, 2, etc.) then htseq-count won't know how to match your alignments to any of the genes. The simplest solution is to not use htseq-count, but instead to use featureCounts, which will handle the chromosome name difference transparently (it's one of the few programs that does) and is much faster anyway. In the future, download your GTF and fasta files from the same source, don't mix and match them.

More

Upvote

VOTE

Downvote
Jessica Allen  Follow
and made sure the chromosome names were corrected; they're now chr1, chr2, etc. However, when I ran the original command with the new file I still got the same problem. Is it maybe an issue with my More
Upvote

VOTE

Downvote
Hobby Logie  Follow
sed /^[[:digit:]]/s/^/chr/ file.gtf > newfile.gtfMore
Upvote

VOTE

Downvote
Douglas Klein  Follow
file.fixed.gtf` for your analysis. However, More
Upvote

VOTE

Downvote
Gail Fisher  Follow
sed 's/^/chr/' file.gtf > file.fixed.gtf`` and then use the More
Upvote

VOTE

Downvote
Gulabo Swat  Follow
Yes, you can edit the GTF, though if you want real consistency you might see if you were using the exact same fasta and GTF files before.More
Upvote

VOTE

Downvote
more replies