arshajii / ema

Fast & accurate alignment of barcoded short-reads
http://ema.csail.mit.edu
MIT License
32 stars 7 forks source link

Incorrect template length output #54

Open pdimens opened 4 months ago

pdimens commented 4 months ago

Looking over the output of some test [haplotag] data, column 9 (TLEN) is showing 0 for many alignments, even if the alignment positions are present.

Example record:

A00470:481:HNYFWDRX2:1:2239:18150:28354_A61C72B93D96_GAACGACTACCACAG    113 2R  1555    54  150M    =   1555    0   AAAGAATAAATTGACTAAGTTATGTCATTTTAAGCGGTCAAAATGGGTGAATTTCCGATTTCAAGTACCAGGCGAACAGAAGACACCTTCTAGAGATTCTGTTCAACTGGTAAGCAAAAACAGTAAATTGCCTAAGTTTTACATTATAAG  FFFFFFFFFFFFFFFFFF:FFFFFFFFFFFF,FFFFFFFFFF,FFF:F:FF:FF:FFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF  NM:i:2  BX:Z:A61C72B93D96   XG:f:1  MI:i:66006  XF:i:1  RG:Z:sample1

TLEN = 0 alignment start = 1554 alignemnt end = 1704