You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
I used tmap map4 to realign a aligned bam to another reference fasta, found that tmap will keep original @sq header and just append new reference. But tmap does not consider original @sq header when it assigns the reference id, so when I used samtools to view the new realigned bam, I found many reads's reference are the orignal ones. For example, if original @sq has ['chr1', 'chr2'], and new reference has ['chr3', chr4'], then a read which maps to 'chr3' will be assigned reference id 0 and get 'chr1' as reference.
The text was updated successfully, but these errors were encountered:
I used tmap map4 to realign a aligned bam to another reference fasta, found that tmap will keep original @sq header and just append new reference. But tmap does not consider original @sq header when it assigns the reference id, so when I used samtools to view the new realigned bam, I found many reads's reference are the orignal ones. For example, if original @sq has ['chr1', 'chr2'], and new reference has ['chr3', chr4'], then a read which maps to 'chr3' will be assigned reference id 0 and get 'chr1' as reference.
The text was updated successfully, but these errors were encountered: