FastqStreamer and qualitytype
2
1
Entering edit mode
David ▴ 860
@david-3335
Last seen 6.1 years ago

Hello,

I have a very large fastq file  in old illumina format (ascii 64 illumina 1.5). I cannot use FastqStreamer with the qualityType. Is there any workaround ? 

FastqStreamer(file, 1e+06, qualityType = "SFastqQuality") : 
  unused argument (qualityType = "SFastqQuality")

 

shortread • 1.2k views
ADD COMMENT
0
Entering edit mode
@martin-morgan-1513
Last seen 3 days ago
United States

I think the quality is being inferred, via the ShortReadQ() function; if the inference is not to your liking, you can pass an argument qualityType= to yield(), for instance after example(FastqStreamer)

> class(quality(yield(FastqStreamer(fl), qualityType="SFastqQuality")))
[1] "SFastqQuality"
attr(,"package")
[1] "ShortRead"
> class(quality(yield(FastqStreamer(fl), qualityType="FastqQuality")))
[1] "FastqQuality"
attr(,"package")
[1] "ShortRead"

 

ADD COMMENT
0
Entering edit mode
David ▴ 860
@david-3335
Last seen 6.1 years ago

Ok great,

Thanks so much.

david

 

ADD COMMENT

Login before adding your answer.

Traffic: 517 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6