Comments (4)
Hi @jwisrael, can you use samtools view to extract the BAM record of "HWI:1:X:1:1216:4127:50166" and "HWI:1:X:1:1216:4180:50180" ?
from rsem.
Hi @bli25wisc here are the records:
HWI:1:X:1:1216:4127:50166 345 ENST00000262493.6 3092 255 77M268435404S * 0 0 AACTTAAAAAAAAAAGGAAAAAAGAAAAAAAAGCCCACGGGTCTTTGTAAG
DA>CDBDDFHJJJIHFGIJJJIHBIJJJJIHBHGHIJIGHHHHFDEFFCCC NH:i:1 HI:i:1
HWI:1:X:1:1216:4180:50180 419 ENST00000467733.1 268 0 50M = 326 109 GAAATTTGTTCTCAAGCACAAAGTAGCACAGAAGAGGGAAGATGCTGTTT B@BDFFFFHHHFHHIJIGIIII?FFDHCHECHG3C>;DDCBG@@gd<?E? NH:i:6 HI:i:1
HWI:1:X:1:1216:4180:50180 339 ENST00000467733.1 326 0 51M = 268 -109 GTGACTCGAAAACTTTCTGAAGCTGATAATAGAAAGATGTCTCGGAAGGAG FFHIHIJIJJJIJJJJJJIIHGIEIIJJJIJJIJJJIGHHFHDDDDFFCCC NH:i:6 HI:i:1
HWI:1:X:1:1216:4180:50180 419 ENST00000496050.1 240 0 50M = 298 109 GAAATTTGTTCTCAAGCACAAAGTAGCACAGAAGAGGGAAGATGCTGTTT B@BDFFFFHHHFHHIJIGIIII?FFDHCHECHG3C>;DDCBG@@gd<?E? NH:i:6 HI:i:2
HWI:1:X:1:1216:4180:50180 339 ENST00000496050.1 298 0 51M = 240 -109 GTGACTCGAAAACTTTCTGAAGCTGATAATAGAAAGATGTCTCGGAAGGAG FFHIHIJIJJJIJJJJJJIIHGIEIIJJJIJJIJJJIGHHFHDDDDFFCCC NH:i:6 HI:i:2
HWI:1:X:1:1216:4180:50180 419 ENST00000463503.1 216 0 50M = 274 109 GAAATTTGTTCTCAAGCACAAAGTAGCACAGAAGAGGGAAGATGCTGTTT B@BDFFFFHHHFHHIJIGIIII?FFDHCHECHG3C>;DDCBG@@gd<?E? NH:i:6 HI:i:3
HWI:1:X:1:1216:4180:50180 339 ENST00000463503.1 274 0 51M = 216 -109 GTGACTCGAAAACTTTCTGAAGCTGATAATAGAAAGATGTCTCGGAAGGAG FFHIHIJIJJJIJJJJJJIIHGIEIIJJJIJJIJJJIGHHFHDDDDFFCCC NH:i:6 HI:i:3
HWI:1:X:1:1216:4180:50180 419 ENST00000476217.1 276 0 50M = 334 109 GAAATTTGTTCTCAAGCACAAAGTAGCACAGAAGAGGGAAGATGCTGTTT B@BDFFFFHHHFHHIJIGIIII?FFDHCHECHG3C>;DDCBG@@gd<?E? NH:i:6 HI:i:4
HWI:1:X:1:1216:4180:50180 339 ENST00000476217.1 334 0 51M = 276 -109 GTGACTCGAAAACTTTCTGAAGCTGATAATAGAAAGATGTCTCGGAAGGAG FFHIHIJIJJJIJJJJJJIIHGIEIIJJJIJJIJJJIGHHFHDDDDFFCCC NH:i:6 HI:i:4
HWI:1:X:1:1216:4180:50180 419 ENST00000467789.1 300 0 50M = 358 109 GAAATTTGTTCTCAAGCACAAAGTAGCACAGAAGAGGGAAGATGCTGTTT B@BDFFFFHHHFHHIJIGIIII?FFDHCHECHG3C>;DDCBG@@gd<?E? NH:i:6 HI:i:5
HWI:1:X:1:1216:4180:50180 339 ENST00000467789.1 358 0 51M = 300 -109 GTGACTCGAAAACTTTCTGAAGCTGATAATAGAAAGATGTCTCGGAAGGAG FFHIHIJIJJJIJJJJJJIIHGIEIIJJJIJJIJJJIGHHFHDDDDFFCCC NH:i:6 HI:i:5
HWI:1:X:1:1216:4180:50180 419 ENST00000265044.2 335 0 50M = 393 109 GAAATTTGTTCTCAAGCACAAAGTAGCACAGAAGAGGGAAGATGCTGTTT B@BDFFFFHHHFHHIJIGIIII?FFDHCHECHG3C>;DDCBG@@gd<?E? NH:i:6 HI:i:6
HWI:1:X:1:1216:4180:50180 339 ENST00000265044.2 393 0 51M = 335 -109 GTGACTCGAAAACTTTCTGAAGCTGATAATAGAAAGATGTCTCGGAAGGAG FFHIHIJIJJJIJJJJJJIIHGIEIIJJJIJJIJJJIGHHFHDDDDFFCCC NH:i:6 HI:i:6
from rsem.
Hi @jwisrael , I think this is a STAR aligner bug. The first alignment has CIGAR string "77M268435404S", which should not appear for transcript-based alignments. If you are not using the latest version of STAR, maybe you can consider update your program.
Hope it helps,
Bo
from rsem.
Hi @bli25wisc yes I think that fixed the problem. Thanks!
from rsem.
Related Issues (20)
- 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
- Not sure if rsem is using --num-threads option...
- fatal error: buffer size for SJ output is too small Solution: increase input parameter --limitOutSJcollapsed HOT 3
- Single-end fragment length mean and SD
- rsem-extract-reference-transcripts fails with "Error Message: Strand is neither '+' nor '-'!"
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.