# Fq2bam rel3-nanopore-wgs-288418386-FAB39088.fastq.gz

**URL:** <https://forums.developer.nvidia.com/t/fq2bam-rel3-nanopore-wgs-288418386-fab39088-fastq-gz/166666>\
**Category:** Parabricks\
**Created:** [January 24, 2021, 8:00pm UTC](https://forums.developer.nvidia.com/t/fq2bam-rel3-nanopore-wgs-288418386-fab39088-fastq-gz/166666 "2021-01-24T20:00:49Z")\
**Posts on this page:** 5\
**Page:** 1

<div class="post-metadata">

**Author:** ![alex-dna-tech](https://developer.download.nvidia.com/images/forums/profile-default-devtalk-84.png) [@alex-dna-tech](https://forums.developer.nvidia.com/u/alex-dna-tech)\
**Post date:** [January 24, 2021, 8:00pm UTC](https://forums.developer.nvidia.com/t/fq2bam-rel3-nanopore-wgs-288418386-fab39088-fastq-gz/166666/1 "2021-01-24T20:00:49Z")

</div>

Hi. Please help to debug.  
Sample downloaded from  
`aws s3 cp s3://nanopore-human-wgs/rel3-nanopore-wgs-288418386-FAB39088.fastq.gz .`

```
sudo pbrun fq2bam --ref Ref/Homo_sapiens_assembly38.fasta \
	--in-se-fq rel3-nanopore-wgs-288418386-FAB39088.fastq.gz \
	--knownSites Ref/Homo_sapiens_assembly38.known_indels.vcf.gz \
	--out-recal-file recal_gpu.txt \
	--out-bam mark_dups_gpu.bam \
	--tmp-dir /raid
Please visit https://docs.nvidia.com/clara/#parabricks for detailed documentation

[Parabricks Options Mesg]: Checking argument compatibility
[Parabricks Options Mesg]: Automatically generating ID prefix
[Parabricks Options Mesg]: Read group created for /raid/rel3-nanopore-wgs-288418386-FAB39088.fastq.gz
[Parabricks Options Mesg]:
@RG\tID:0cf3beca-6611-4690-b756-32f304.1\tLB:lib1\tPL:bar\tSM:sample\tPU:0cf3beca-6611-4690-b756-32f304.1
------------------------------------------------------------------------------
|| Parabricks accelerated Genomics Pipeline ||
|| Version v3.2.0_ampere ||
|| GPU-BWA mem, Sorting Phase-I ||
|| Contact: Parabricks-Support@nvidia.com ||
------------------------------------------------------------------------------
[M::bwa_idx_load_from_disk] read 0 ALT contigs

GPU-BWA mem
ProgressMeter	Reads Base Pairs Aligned
PARABRICKS: ParaBricks/src/CReadWrite.cpp:197: bool FastQReader::getSingleGZRead(bseq1_t*, int, char*, int&, int&): Assertion `(sizeOfSeq < MAX_READ_SIZE) || (!printf("Size of seq is %d\n", sizeOfSeq))' failed.
Received signal: 6
Please make sure you have enough disk space in the output and temp directory
Contact Parabricks-Support@nvidia.com for support
2493: a6701b1e-67ed-4a56-8225-ad4234af3eef_Basecall_Alignment_template
Size of seq is 2493
Please contact Parabricks-Support@nvidia.com for any questions
There is a forum for Q&A as well at https://forums.developer.nvidia.com/c/healthcare/Parabricks/290
Exiting...

Could not run fq2bam
Exiting pbrun ...

```

---

<div class="post-metadata">

**Author:** ![mdemouth](https://sea2.discourse-cdn.com/nvidia/user_avatar/forums.developer.nvidia.com/mdemouth/32/16591_2.png) [@mdemouth](https://forums.developer.nvidia.com/u/mdemouth)\
**Post date:** [January 25, 2021, 3:17pm UTC](https://forums.developer.nvidia.com/t/fq2bam-rel3-nanopore-wgs-288418386-fab39088-fastq-gz/166666/2 "2021-01-25T15:17:00Z")

</div>

Hi Alex,

We are sorry that you are facing this error.  
Unfortunately, fq2bam is only intended for short reads, and here your you use a long reads fastq file.

Myrieme

---

<div class="post-metadata">

**Author:** ![alex-dna-tech](https://developer.download.nvidia.com/images/forums/profile-default-devtalk-84.png) [@alex-dna-tech](https://forums.developer.nvidia.com/u/alex-dna-tech)\
**Post date:** [January 28, 2021, 7:01am UTC](https://forums.developer.nvidia.com/t/fq2bam-rel3-nanopore-wgs-288418386-fab39088-fastq-gz/166666/3 "2021-01-28T07:01:15Z")

</div>

Hi, mdemouth.

Thank you for your reply.  
Another question. Originally [BWA MEM works with](https://github.com/lh3/bwa#1-what-types-of-data-does-bwa-work-with) Oxford Nanopore reads. Is this option `-x` implemented in Parabriks? Or Parabricks does work exclusively with Illumina?

---

<div class="post-metadata">

**Author:** ![jamesdalg](https://developer.download.nvidia.com/images/forums/profile-default-devtalk-84.png) [@jamesdalg](https://forums.developer.nvidia.com/u/jamesdalg)\
**Post date:** [January 29, 2024, 3:20pm UTC](https://forums.developer.nvidia.com/t/fq2bam-rel3-nanopore-wgs-288418386-fab39088-fastq-gz/166666/4 "2024-01-29T15:20:45Z")

</div>

I wish to report the exact same thing happens with PacBio reads when using the minimap2 from version 4.2.1 parabricks. There is a hardcoded limit of 100,000bp. If a bam contains a read that long or longer, then the entire application collapses like a house of cards and ends the alignment. It would be better if it would just skip that read or, better yet, increase the maximum read size to be larger than this. There are some reads that are better than 30kb and sometimes it happens. There should be a way to deal with it gracefully.

---

<div class="post-metadata">

**Author:** ![gburnett](https://sea2.discourse-cdn.com/nvidia/user_avatar/forums.developer.nvidia.com/gburnett/32/14047_2.png) [@gburnett](https://forums.developer.nvidia.com/u/gburnett)\
**Post date:** [January 29, 2024, 5:57pm UTC](https://forums.developer.nvidia.com/t/fq2bam-rel3-nanopore-wgs-288418386-fab39088-fastq-gz/166666/5 "2024-01-29T17:57:09Z")

</div>

Hi @jamesdalg,

This issue will be fixed in the next version of Parabricks 4.3.

Thank you for your patience.