error with samtools

dear all

i have some problem with samtools ....

samtools mpileup -u -q 10 -f '/media/mra/4A3B44B64CA65A93/SNP/DATA/Genom/refrense_genome/Apis_mellifera.GCA_000002195.1.dna.chromosome_group.1.fa' '/media/mra/4A3B44B64CA65A93/SNP/DATA/drone-larvae/SRR1571717/segemehl/f-sort.bam' ‫|‬ ‫‪bcftools‬‬ ‫‪view‬‬ ‫‪-g‬‬ ‫‪-‬‬ ‫|‬ ‫‪vcfutils.pl‬‬ ‫‪varFilter‬‬ ‫|‬ ‫‪awk‬‬ ‫‪'($6‬‬ ‫>= 50)' > '/media/mra/4A3B44B64CA65A93/SNP/DATA/drone-larvae/SRR1571717/segemehl/‫‪file.raw.pileup‬‬'  


[E::hts_open_format] fail to open file '‫'
[mpileup] failed to open ‫: No such file or directory
: command not found
: command not found
: command not found

what's my script problem ..?

rna-seq

• 3,083 views • link

updated 9.5 years ago by Biostar

• written 9.6 years ago by mra8187 •

Unless you know how each component works, avoid pipes because those will just confuse you. Based on a quick look (while preparing food) it looks like quite a number of things in your command are wrong.

• 0 views • link

written 9.6 years ago by WouterDeCoster

I am trying to format the content OP posted but the content is confusing both the website and me - what I see above is not what is being shown when I hit edit.

• 0 views • link

written 9.6 years ago by Ram

Then let's wait for OP to clarify. But looking at the error messages I get the idea OP is piping things into "empty" commands.

• 0 views • link

written 9.6 years ago by WouterDeCoster

so how i can edit my script ?

• 0 views • link

written 9.6 years ago by mra8187 •

First doing things step by step, if all parts work you can chain the commands using pipes.

• 0 views • link

written 9.6 years ago by WouterDeCoster

0 answers

No answers yet.