Skip to content

Issue when running mm10 using custom bed file #89

Open
@saimukund20

Description

@saimukund20

I uploaded my inputs to google drive. https://drive.google.com/drive/folders/1jJ2QvBTdKuUFfnQS6769IHGLSqXelcVB?usp=sharing

I tried running the ngsplot.r on my downsampled bam as a test using my custom bed file with the following command:

ngs.plot.r -G mm10 -R bed -E Analysis7/WGBS_Analysis3_dmrs.bed -C downsampled.bam -O out -T 45417_h3K27_tri -P 4 -L 10000

I dont get any errors but warnings. Here is the output i am getting:

Configuring variables...Done
Loading R libraries.....Done
Analyze bam files and calculate coverageWarning message:
In headerIndexBam(bam.list) :
Aligner for: downsampled.bam cannot be determined. Style of
standard SAM mapping score will be used. Would you mind submitting an issue
report to us on Github? This will benefit people using the same aligner.
..................................................................................Done
Plotting figures...Done
Saving results...Warning message:
In data.row.names(row.names, rowsi, i) :
some row.names duplicated: 2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100,101,102,103,104,105,106,107,108,109,110,111,112,113,114,115,116,117,118,119,120,121,122,123,124,125,126,127,128,129,130,131,132,133,134,135,136,137,138,139,140,141,142,143,144,145,146,147,148,149,150,151,152,153,154,155,156,157,158,159,160,161,162,163,164,165,166,167,168,169,170,171,172,173,174,175,176,177,178,179,180,181,182,183,184,185,186,187,188,189,190,191,192,193,194,195,196,197,198,199,200,201,202,203,204,205,206,207,208,209,210,211,212,213,214,215,216,217,218,219,220,221,222,223,224,225,226,227,228,229,230,231,232,233,234,235,236,237,238,239,240,241,242,243,244,245,246,247,248,249,250,251,252,253,254,255,256,257,258,259,260,261,262,263,264,265,266,267,268,269,270,271 [... truncated]
Done
Wrapping results up...Done
All done. Cheers!

The outputs seem to be empty. I have placed them also in the google drive folder.

I would appreciate your help in troubleshooting this error.

Attached are all the inputs and outputs for your verification
Archive.zip

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions