[mira_talk] Same singlets,

  • From: wakamoto5959@xxxxxxxxxxx
  • To: "mira_talk@xxxxxxxxxxxxx" <mira_talk@xxxxxxxxxxxxx>
  • Date: Wed, 31 Aug 2016 09:23:06 +0900 (JST)

Dear Bastien and All, 

I really appreciate Mira. Thank you very much! 

Today, I have a question about singlets. 

After an assembly, if I check singlets in out.unpadded.fasta file, many of 
those singlet reads are exactly same sequence. 
So I am wondering why this is happening. 
I have pasted my manifest information below. 
I also uploaded input fastq files. 

I would very much appreciate your help. 

Best regards, 

Lyu  


Manifest file below 

project = 10K
job = est,denovo,accurate
parameters = COMMON_SETTINGS  

readgroup = SomePairedEndIlluminaReadsIGotFromTheLab
data = /mira_4.0.2_linux-gnu_x86_64_static/bin/10KR1.fastq 
/mira_4.0.2_linux-gnu_x86_64_static/bin/10KR2.fastq

technology = solexa
autopairing  
segment_naming = solexa 

parameters = -SK:pr=90   

parameters = -SK:mmhr=10 
parameters = -SK:bph=32 
parameters = -HS:mnr=no  
parameters = -HS:ldn=no 
parameters = -NW:cmrnl=no 
parameters = -CO:mr=yes  
parameters = -CL:ascdc=yes 

parameters = SOLEXA_SETTINGS 
parameters = -AL:mrs=90 
parameters = -AS:mrpc=1  
parameters = -AL:mo=15 
parameters = -OUT:sssip=yes 

Attachment: 10KR1.fastq.tar.gz
Description: application/gzip

Attachment: 10KR2.fastq.tar.gz
Description: application/gzip

Other related posts: