Read length mismatch between bam and query reads
0
0
Entering edit mode
6 months ago
蓝寒 • 0

Hi all,

I met a problem when mapping Pacbio Iso-seq reads.

I got ccs reads and corrected by NGS data, but when I use GMAP located in SMRTlink, I found that:

The length of reads I extracted from bam file 'SEQ' list are longer than query sequences I submitted.

Which means that, for example, I use 1k read 'AA' to map reference but the 'SEQ' length of AA in bam is 1.2k.

I just wonder why, I have checked previous questions but didn't get accurate result, is that due to full-length? Some answer speak of Insertion or trim, but that should be query longer than bam ?

Thanks all!

GMAP Pacbio • 267 views
ADD COMMENT

Login before adding your answer.

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