Skip to content

generat data #7

@atefeh-f

Description

@atefeh-f

hello
I faced a problem in my implementation phase. Would you please guide me to solve it? Please find the details in bellow:
I downloaded dataset from REMC and I did run the readme file instructions on this dataset. I run this instruction:

bedtools bedtobam -i E128-H3K4me3.tagAlign -g hg19chrom.sizes > ft2.bam

but I faced with this error message :

Error: The requested genome file (hg19chrom.sizes) could not be opened. Exiting!

Thus I downloaded the file “hg19.chrom.sizes” and I run this instruction:

bedtools bedtobam -i E128-H3K4me3.tagAlign -g hg19.chrom.sizes > ft2.bam

then the error message omited. Then the above instruction produced a bam file but when I run “bedtools multicov” I faced with this error message:

bedtools multicov -bams ft2.bam -bed E128-H3K4me3.tagAlign
Could not open input BAM files.

Thus i install "samtools" and i run this instruction:

samtools sort ft2.bam > ft2.sort
samtools index ft2.sort
bedtools multicov -bams ft2.sort -bed E128-H3K4me3.tagAlign

then the error message omited.
I get RPKM from : http://egg2.wustl.edu/roadmap/data/byDataType/rna/expression/
How could I run RPKM on my files? Plus, where the instructions which I used above correct?

Best Regards

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type
    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions