sorting bowtie2 results using samtools sort
0
0
Entering edit mode
6.8 years ago
ggman ▴ 80

Greetings friends,

I have successfully ran bowtie2 on my file and I am having trouble sorting my data based on a genome file that I have genome.fa.fai

It appears samtools does not use import anymore so I do not know how to go about using this file as a reference. This is the command I have been trying to use.

samtools import ~/genome_data/ref_genome.fa.fai - - | samtools sort - bowtie_run.pair

I am going to eventually merge the files, which I do not have a problem with but the command above is giving me trouble when running it through a shell script, bowtie2 is running successfully so I am almost certain that is not the problem.

bowtie2 samtools sort • 2.4k views
ADD COMMENT
1
Entering edit mode
 import    samtools import <in.ref_list> <in.sam> <out.bam>

                 Since 0.1.4, this command is an alias of:

                 samtools view -bt <in.ref_list> -o <out.bam> <in.sam>
ADD REPLY
0
Entering edit mode

what is samtools import ?

ADD REPLY
0
Entering edit mode

it ran on an older version of samtools at some point or something like that. Actually... I think I found a solution....

 import    samtools import <in.ref_list> <in.sam> <out.bam>

                 Since 0.1.4, this command is an alias of:

                 samtools view -bt <in.ref_list> -o <out.bam> <in.sam>

can be found here

ADD REPLY

Login before adding your answer.

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