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

A few questions about using nanopolish to estimate poly(A) tail length #1125

Open
lopp1Y opened this issue Jan 15, 2024 · 1 comment
Open

Comments

@lopp1Y
Copy link

lopp1Y commented Jan 15, 2024

Hi,I'm running the module how to estimate poly(A) tail lengths from nanopore native RNA reads, and I've got the result file, but I don't understand the result file right now, I think this command line: nanopolish polya --threads=8 --reads=data/30xpolyA.fastq -- bam=data/30xpolyA.sorted.bam --genome=data/enolase_reference.fas > polya_results.tsv

The data in this command line is fastq, bam, fas, the explanation of this module in the literature is to use the hidden Markov model for segmentation of the raw sequencing signal , I would like to ask which of these has the information that provides signals for the hidden Markov model to analyze? What I understand is that fastq is the base sequence and BAM is the sequence alignment result. Or does this provide additional signal information? Can you help me with that?
Thank you .

@jts
Copy link
Owner

jts commented Jan 16, 2024

Hi @lopp1Y,

Before running nanopolish polya you should have run nanopolish index, which builds a map from the basecalled reads (fastq) to the signal reads (fast5/pod5). This allows nanopolish to access the signal data.

Jared

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

No branches or pull requests

2 participants