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

Fastq reader clips long reads at 10,000 bp #1992

Closed
heuermh opened this issue May 24, 2018 · 0 comments
Closed

Fastq reader clips long reads at 10,000 bp #1992

heuermh opened this issue May 24, 2018 · 0 comments
Milestone

Comments

@heuermh
Copy link
Member

heuermh commented May 24, 2018

E.g. ftp://ftp-trace.ncbi.nlm.nih.gov/giab/ftp/data/AshkenazimTrio/HG002_NA24385_son/Ultralong_OxfordNanopore/

$ dsh-bio fastq-sequence-length -i combined_2018-05-18.fastq.gz | sort -n -r | head
322099
298433
283861
262525
247501
scala> val reads = sc.loadAlignments("combined_2018-05-18.fastq.gz")
reads: org.bdgenomics.adam.rdd.read.AlignmentRecordRDD =
ParquetUnboundAlignmentRecordRDD with 0 reference sequences, 0 read groups,
and 0 processing steps

scala> val lengths = reads.rdd.map(_.sequence.length()).collect()
lengths: Array[Int] = Array(641, 3082, 10000, 10000, 10000, ...

scala> lengths.sortBy(-1 * _).take(10)
res3: Array[Int] = Array(10000, 10000, 10000, 10000, 10000, 10000, 10000, 10000, ...
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