Question

How Do You Convert A Sam/Bam File To A Psl Or Axt Format?

2

Entering edit mode

10.8 years ago

lp364 ▴ 20

I am trying to make a chain file to use in liftOver and these can be made using PSL or AXT files. My current data is in SAM/BAM. What's the best way to convert this to a psl or axt format?

samtools conversion • 3.8k views

ADD COMMENT • link updated 10.8 years ago by Alex Reynolds 35k • written 10.8 years ago by lp364 ▴ 20

0

Entering edit mode

cross posted on seqanswers: http://seqanswers.com/forums/showthread.php?p=112080#post112080

ADD REPLY • link 10.8 years ago by Pierre Lindenbaum 161k

score 0 · Answer 1 · 2013-07-30

0

Entering edit mode

10.8 years ago

Alex Reynolds 35k

Perhaps you could try passing the output from the sam2bed script to the bed2psl.py script. I'd take care to make sure coordinate indexing is correct at each step of the process.

ADD COMMENT • link 10.8 years ago by Alex Reynolds 35k

2

Entering edit mode

where can you download the bed2psl.py? It seems that the sam2bed doesn't process the cigar string that is needed to generate the psl.

ADD REPLY • link 10.8 years ago by lenore.pipes ▴ 20