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!