This is a test version of Biostars. For the public version, visit https://www.biostars.org.
error for snp filtering

Dear all i am using this script to SNP filtering and get error ... what 's wrong .. best regards 

**bowtie-inspect -s /media/mra/4A3B44B64CA65A93/SNP/DATA/Genom/refrense_genome/Bowtie2/indexed-refrence.fa  | tail -n +4 | cut -f 2- > '/home/mra/Desktop/genome_size.txt' cat junctions.bed | bed_to_juncs | flankBed -g '/home/mra/Desktop/genome_size.txt' -b 1 | s lopBed -g '/home/mra/Desktop/genome_size.txt' -b 2 > '/home/mra/Desktop/windows.bed'** # my scipt

Error

cut: cat: No such file or directory
cut: junctions.bed: No such file or directory
Could not locate a Bowtie index corresponding to basename "/media/mra/4A3B44B64CA65A93/SNP/DATA/Genom/refrense_genome/Bowtie2/indexed-refrence.fa"
Command: bowtie-inspect --wrapper basic-0 -s /media/mra/4A3B44B64CA65A93/SNP/DATA/Genom/refrense_genome/Bowtie2/indexed-refrence.fa 
s: command not found
rna-seq snp software error
• 1,827 views • link
updated by Pierre Lindenbaum Pierre Lindenbaum • written by mra8187 • mra8187

1 answer

missing pipe or semicolon before cat:

not

... nome_size.txt' cat junctio ...

but something like

nome_size.txt' |  cat junctio

or 

nome_size.txt' ;  cat junctio

furthermore:

not:

| s lopBed

but

| slopBed
• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

thanks a lot. now .. my problem is :

my script : bowtie-inspect -s indexed-refrence.fa | tail -n +4 | cut -f 2- > '/home /mra/Desktop/genome_size.txt' cat junctions.bed | bed_to_juncs | flankBed -g '/home/mra/Desktop/genome_size.txt' -b 1 | slopBed -g '/home/mra/Desktop/genome_size.txt' -b 2 > '/home/mra/Desktop/windows.bed'

error is :

cut: cat: No such file or directory cut: junctions.bed: No such file or directory Could not locate a Bowtie index corresponding to basename "indexed-refrence.fa" Command: bowtie-inspect --wrapper basic-0 -s indexed-refrence.fa

• 0 views • link
written by mra8187 • mra8187

Log in to answer this question.

More posts like this
  • SNP Discovery error
    written by a.hasankhani74 • a.hasankhani74
    hi I want to create a bed file from bowtie-inspect summary and the TopHat output file with the bed extension, so a new bed-format file …
  • problem with bedtools
    written by mra8187 • mra8187
    Dear all i am trying to use bedtools to snp filtering by get the error . what's wrong ? my script is : cat '/media/mra/4A3B44B64CA65A93/SNP/DATA/drone-larvae/SRR1571716/Trimmomatic/tophat_out/junctions.bed' …
  • how to use bowtie-inspect
    written by mra8187 • mra8187
    Dear all how can i use bowtie-inspect .. my script is : bowtie-inspect -s indexed-refrence.fa | tail -n +4 | cut -f 2- > '/home …
  • error with samtools
    written by mra8187 • mra8187
    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‬‬ …
  • problem with trinity
    written by mra8187 • mra8187
    dear all when i am using Trinity v 2.3.2 by this command .`/Trinity --seqType fq --max_memory 5G --single '/media/mra/4A3B44B64CA65A93/SNP/DATA/Queen-Larvae/SRR1571722/Trimmomatic/SE/SRR1571722-1.trimmed.fastq' --CPU 4` i see error ... …
  • Question: Unable to find Bowtie2 index
    written by mra8187 • mra8187
    hello i indexed genome by bowtie2 and then use them to tophat2 but still give error ..... tophat2 -p 4 -r 200 --mate-std-dev 50 -G …
  • .gtf file error in Tophat2 - Error at parsing .tlst line (invalid strand):
    written by 4galaxy77 4galaxy77
    I'm trying to align my reads to a genome + transcriptome, which is in the form of a .gtf file. It gets half way through …