Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

SquigglePlot.py with -n 20 option produces more than one plot #53

Closed
waltergallegog opened this issue Jan 7, 2022 · 3 comments
Closed
Assignees
Labels
question Further information is requested

Comments

@waltergallegog
Copy link

Hello,
I'm following your tutorial and wanted to plot a .fast5 I got from here, given that the one under the examples folder was giving me an error (same as issue 51 )

I used the -n NUMoption, which if I understood correctly should plot only the first NUM data points.
However what I observed is that multiple plots of size NUM are generated.

Same behavior observed if using -n START,FINISH or the save options.
Is this the intended behavior?
Any way to limit to just one plot?

Some of the commands I used:

python SquigglePlot.py -i ../datasets/Chip137_IVT_NA12878_Data_reads_0.fast5 -n 20
Generates a plot window with 20 datapoints. As soon as I close it, a new one opens, with different data ( I suppose the next 20 points).
I assume it would continue to generate new plots until reaching the end of the file ( I closed maybe 15 windows, then just killed the process)

python SquigglePlot.py -i ../datasets/Chip137_IVT_NA12878_Data_reads_0.fast5 -n 40,50
Generates a plot window with 10 datapoints. As soon as I close it, a new one opens, with different data ( I suppose points 50,60)

python SquigglePlot.py -i ../datasets/Chip137_IVT_NA12878_Data_reads_0.fast5 -n 20 --dpi 300 --no_show --save test.pdf --save_path example/pngs20
Generates a bunch of pdf files, each with a graph of 20 datapoints.

Thanks
Walter

@Psy-Fer
Copy link
Owner

Psy-Fer commented Jan 7, 2022

Hello,

It's not the NEXT 20 datapoints, it's the 20 datapoints of the next read in the multifast5 file.

If you wish to limit which read is being analysed, you can use the --readID flag to say which read.

I hope that helps

@Psy-Fer Psy-Fer self-assigned this Jan 7, 2022
@Psy-Fer Psy-Fer added the question Further information is requested label Jan 7, 2022
@waltergallegog
Copy link
Author

Hello,
I didn't realized it was a multfast5 file. Thanks for the quick response and sorry for the trouble.
Walter

@Psy-Fer
Copy link
Owner

Psy-Fer commented Jan 7, 2022

All good :) let me know if you have any other questions or something isn't working how you want it to.

James

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
question Further information is requested
Projects
None yet
Development

No branches or pull requests

2 participants