error with samtools
0
0
Entering edit mode
7.9 years ago
mra8187 ▴ 20

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 • 2.6k views
ADD COMMENT
3
Entering edit mode

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.

ADD REPLY
2
Entering edit mode

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.

ADD REPLY
1
Entering edit mode

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.

ADD REPLY
0
Entering edit mode

so how i can edit my script ?

ADD REPLY
4
Entering edit mode

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

ADD REPLY

Login before adding your answer.

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