Comments (2)
Yes, RSEM does require the HI:i: tag in order to correctly process multimappers. The HI tag (HI:i:) is a field in SAM/BAM files that is used to indicate the index of a read within the input file. When processing multimapped reads, RSEM uses the HI tag to keep track of which reads are duplicates, so that they are not counted multiple times.
If the BAM files generated by Hisat2 do not contain the HI tag, you can add it using the samtools command-line tool. The samtools addreplacerg command allows you to add or modify the values of any field in the SAM/BAM header, including the HI tag.
Here's an example of how to add the HI tag to a BAM file using samtools:
samtools addreplacerg -r 'ID:1' -r 'SM:sample' -r 'PL:illumina' -r 'PU:unit' -r 'LB:library' -r 'DT:2023-05-03T12:34:56-0700' -R 'HiSeq2500' -i 'HI' -o output.bam input.bam
In this example, the ID, SM, PL, PU, LB, and DT tags are arbitrary and can be changed to match your specific sample information. The -R option specifies the sequencing platform used to generate the data, and the -i option specifies the name of the index tag to be added (in this case, "HI"). The input BAM file is "input.bam", and the output BAM file with the added HI tag is "output.bam".
from rsem.
Hi Aryan,
Thank you for your answer! I am confused, however. RSEM appears to support hisat2 as a mapping tool - yet the options used by RSEM do not allow for HI tag to be outputted.
The issue I was having was due to sorting of the BAM file; apparently, transcriptomic files are always read name-sorted by default, and both RSEM and Salmon expect this to be true in order to work correctly. Adding HI tag does not seem to influence the result.
from rsem.
Related Issues (20)
- single-cell-prior is not giving me any values HOT 1
- rsem-calculate-expression
- gene id
- import using tximport for the rsem count table to get integers?
- how to simulate single-end reads?
- The SAM/BAM file declares less than one reference sequence! HOT 2
- genomeSAsparseD option for RSEM Star?
- Problem implementing a reverse-strand specific analysis HOT 2
- about gffcompare
- RSEM deal with multimapping? HOT 1
- Pearl error
- error when using convert-sam-for-rsem to a RSEM compatible BAM file.
- rsem-prepare-reference arguments
- Using STAR BAMs from genomic alignment HOT 2
- rsem-parse-alignments and rsem-build-read-index
- Gtf file
- cleaned.fq.gz: command not found HOT 3
- Two mates are aligned to two different transcripts HOT 1
- rsem-generate-ngvector command not found HOT 1
- bam only containing aligned reads produced different expression profile
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from rsem.