Question: Error in processAmplicons "hairpin sequence length not with specified length"
gravatar for lucia.caceres
20 months ago by
lucia.caceres0 wrote:


When I type the following code:

> x = processAmplicons("Index2.Plate10.fastq", barcodefile="Samples1.txt",
+     hairpinfile="Hairpins sense.txt", hairpinStart=29, hairpinEnd=50, verbose=TRUE)
I get:

Error in processAmplicons("Index2.Plate10.fastq", barcodefile = "Samples1.txt",  :
  Hairpin sequence length is set to 22, there are hairpin sequences not with specified length.

I checked the length of my hairpin sequences, in my hairpin file/reference library and they are all 22nt long.

I don't understand what is going wrong. If you have any ideas of what I can do/check please let me know.



ADD COMMENTlink modified 20 months ago by Matthew Ritchie690 • written 20 months ago by lucia.caceres0
gravatar for Matthew Ritchie
20 months ago by
Matthew Ritchie690 wrote:

This error message suggests you have hairpins of different lengths in your library.

To check if this is the case, try running something like: 

hp = read.delim("Hairpins sense.txt", stringsAsFactors=FALSE)
summary(sapply(hp$Seq, FUN="nchar"))
which(sapply(hp$Seq, FUN="nchar")!=22)
hp[sapply(hp$Seq, FUN="nchar")!=22,]
ADD COMMENTlink written 20 months ago by Matthew Ritchie690

Thank you for the code. It helped a lot.

ADD REPLYlink written 20 months ago by lucia.caceres0
Please log in to add an answer.


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