This is a test version of Biostars. For the public version, visit https://www.biostars.org.
problem with bedtools

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' \
    | bed_to_juncs | flankBed -g '/home/mra/Downloads/genome_size.txt' -b 1 \
    | slopBed -g '/home/mra/Downloads/genome_size.txt' -b 5 \
    > '/media/mra/4A3B44B64CA65A93/SNP/DATA/drone-larvae/SRR1571716/Trimmomatic/tophat_out/window.bed'

error is : 

Warning: malformed line 1, missing columns
     track name=junctions description="TopHat junctions"
ERROR: chrom "1" not found in genome file. Exiting.
Traceback (most recent call last):
  File "/home/mra/miniconda3/bin/bed_to_juncs", line 69, in <module>
    sys.exit(main())
  File "/home/mra/miniconda3/bin/bed_to_juncs", line 59, in main
    print("%s\t%d\t%d\t%s\t%s" % (chromosome, left_pos, right_pos, orientation, counts))
BrokenPipeError: [Errno 32] Broken pipe
snp rna-seq software error
• 1,984 views • link
updated by Ram Ram • written by mra8187 • mra8187

1 answer

Are you sure your input files are built for the same reference genome? This looks like it could be a case of GRCh37 naming vs ucsc hg19 naming conventions (>1,>2,>3,... vs >chr1, >chr2,>chr3,...)

• 0 views • link
written by Ram Ram

thanks for your help .... very thanks ..........

• 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 …
  • GBS pipeline POD Error: does not exist or is unreadable
    written by csjlee08 • csjlee08
    Hello, I am trying to use GBS-POD pipeline, having trouble with the system to find my files although I did specify the file paths in …
  • 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 for snp filtering
    written by mra8187 • mra8187
    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 …
  • 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 …