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

What is readsize? #5

Open
yangyxt opened this issue Dec 2, 2020 · 0 comments
Open

What is readsize? #5

yangyxt opened this issue Dec 2, 2020 · 0 comments

Comments

@yangyxt
Copy link

yangyxt commented Dec 2, 2020

Dear Firtina,

THanks for the update. I still have some confusion about the definition of the input argument READSIZE.

Is this metric reflecting the number of pairs of reads in my prepared fastq file? or the total number of reads (e.g. 2 times number of pairs, in my case I don't have singleton reads, every read is paired. ) or the read length?

In extract reads.sh, you typed:
Uploading image.png…

Can I ask why comparing the PNEXT in SAM with readsize? array[8] is supposed to be the 8th field delimited by tab in a line of SAM format, which is PNEXT instead of CIGAR string.

From this command, I suppose readsize is read length. While not every read has the exact same length. So could you pls help explain a bit of this argument's usage? Much appreciated.

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

1 participant