Search
Question: FastqStreamer and qualitytype
1
gravatar for David
14 months ago by
David860
David860 wrote:

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")

 

ADD COMMENTlink modified 14 months ago • written 14 months ago by David860
0
gravatar for Martin Morgan
14 months ago by
Martin Morgan ♦♦ 20k
United States
Martin Morgan ♦♦ 20k wrote:

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 COMMENTlink modified 14 months ago • written 14 months ago by Martin Morgan ♦♦ 20k
0
gravatar for David
14 months ago by
David860
David860 wrote:

Ok great,

Thanks so much.

david

 

ADD COMMENTlink written 14 months ago by David860
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 299 users visited in the last hour