Strand inconsistency between Long read tool and minimap2?
0
0
Entering edit mode
9 weeks ago
shinyjj ▴ 50

I found there is an inconsistency on LR strand between LR gene assignment table by IsoQuant and LR alignment bam file by minimap. Do you know what could cause this inconsistency? I see this NM_006988.5 is - strand according to refseq annotation so IsoQuant gave it correctly but not sure why minimap2 is giving + strand. I gave refseq to IsoQuant when running and used the following command for minimap2. I also gave junctions_bed which gives where splicing junctions are to be more precise. Is this because long read is error-prone?

./minimap2 -a -Y --junc-bed <junctions.bed> -k 14 -x splice -t 20

enter image description here

minimap2 • 227 views
ADD COMMENT

Login before adding your answer.

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