Search
Question: FastqStreamer and qualitytype
1
gravatar for David
22 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 22 months ago • written 22 months ago by David860
0
gravatar for Martin Morgan
22 months ago by
Martin Morgan ♦♦ 21k
United States
Martin Morgan ♦♦ 21k 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 22 months ago • written 22 months ago by Martin Morgan ♦♦ 21k
0
gravatar for David
22 months ago by
David860
David860 wrote:

Ok great,

Thanks so much.

david

 

ADD COMMENTlink written 22 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: 158 users visited in the last hour