RNA-seq : Quantification problem after mapped with BBMap
2
0
Entering edit mode
5.6 years ago
k.kathirvel93 ▴ 300

Hi EveryOne,

I have used BBMap for RNA-seq reads mapping. After mapping, i tried HTSeq, FeatureCounts and STRINGtie for quantification. But, i didn't get any results from featureCounts and HTSeq (readcount for all genes 0), STRINGtie showing that " your gtf file doesn't matched with your reference genome ". The strange is the same reference genome and annotation gtf files were tried with STAR aligner, but i didn't get any errors , working fine. Any suggestions, Thanks .

Here is the command used for BBMap :

./bbmap.sh threads=16 in1=/home/1.fastq in2=/home/2.fastq out=/home/output.sam ref=/home/Homo_sapiens.GRCh38.77.fa nodisk
RNA-Seq next-gen software error • 2.0k views
ADD COMMENT
0
Entering edit mode

please list some header of sequences and screenshot of gtf file to give us better idea.

ADD REPLY
1
Entering edit mode
5.6 years ago
michael.ante ★ 3.8k

Hi k.kathirvel93,

Next to trausch's suggestion that bbmap's cigar string is per default in sam 1.4 notation (which can be switched by setting sam=1.3), BBmap uses the full fasta header for the chromosome name. Instead of having just e.g. "21" in the sam file, you have "21 dna:chromosome chromosome:GRCh38:21:1:46709983:1 REF".

To switch this behaviour of, use the parameter trimreaddescriptions=t.

Cheers,

Michael

ADD COMMENT
0
Entering edit mode

where should I use this trimreaddescriptions=t option? Thanking you

ADD REPLY
1
Entering edit mode

AFAIK, you can use it either during index generation or during mapping.

ADD REPLY
0
Entering edit mode
5.6 years ago
trausch ★ 1.8k

I think, BBMap uses the new CIGAR characters = (sequence match) and X (sequence mismatch) and that might not be supported yet in your quantification methods. The latest version of Alfred v0.1.8 can handle these and also does RNA-Seq quantification:

alfred count_rna -g Homo_sapiens.GRCh38.77.gtf.gz <input.bam>
ADD COMMENT

Login before adding your answer.

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