SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Convert csfasta to fastq yksikaksi Bioinformatics 2 10-30-2011 08:36 PM
Error Bowtie with fastq files (from Solid .csfasta & .qual) pepperoni Bioinformatics 0 10-12-2011 08:20 PM
Problem converting solid2fastq with Barcode m_elena_bioinfo Bioinformatics 7 09-15-2011 07:17 AM
Bowtie cant read my csfasta file ! lakshmaa Bioinformatics 4 06-07-2010 07:22 AM
fastq to csfasta and .qual samt SOLiD 15 10-29-2009 09:11 AM

Reply
 
Thread Tools
Old 11-06-2011, 05:53 AM   #1
yksikaksi
Member
 
Location: Scandinavia

Join Date: Dec 2009
Posts: 20
Question Unmatched read name error when converting csfasta to fastq with solid2fastq.pl

I'm encountered with "unmatched read name" error when I converting paired SOLiD csfasta files to fastq by using solid2fastq.pl from bwa-0.5.9

The error message:
================================
** unmatched read name: '>4_26_52_F3
' != '>4_26_52_F3
'
================================

I read there was suggestion to modify the script as mentioned by

1. drio's post in PE SOLiD reads alignment by bwa thread.

2. javijevi's post in bwa samse segmentation fault thread.

I tired the modification and the "unmatched read name" error remained.

Could anyone shares your solution for this?

Thank you very much.
yksikaksi is offline   Reply With Quote
Old 11-06-2011, 07:48 AM   #2
Richard Finney
Senior Member
 
Location: bethesda

Join Date: Feb 2009
Posts: 700
Default

It's dieing in this routine:
Code:
sub read1 {
  my $i = shift(@_);
  my $j = ($i-1)<<1;
  my ($key, $seq);
  my ($fhs, $fhq) = ($fhr[$j], $fhr[$j|1]);
  while (<$fhs>) {
        my $t = <$fhq>;
        if (/^>(\d+)_(\d+)_(\d+)_[FR]3/) {
          $key = sprintf("%.4d_%.4d_%.4d", $1, $2, $3); # this line could be improved on 64-bit machines
          die(qq/** unmatched read name: '$_' != '$_'\n/) unless ($_ eq $t);
          my $name = "$pre:$1_$2_$3/$i";
          $_ = substr(<$fhs>, 2);
          tr/0123./ACGTN/;
          my $s = $_;
          $_ = <$fhq>;
          s/-1\b/0/eg;
          s/^(\d+)\s*//;
          s/(\d+)\s*/chr($1+33)/eg;
          $seq = qq/\@$name\n$s+\n$_\n/;
          last;
        }
  }
  return defined($seq)? ($key, $seq) : ();
}
Things you might try are modifying this to print out $t and the line number. Please do it on a copy of solid2fastq.pl so you keep the original. You can look at that line number in both files to see what's up. I'm not really sure what's up with that extra carriage return(s) in the error message.

Last edited by Richard Finney; 11-06-2011 at 07:51 AM. Reason: format
Richard Finney is offline   Reply With Quote
Old 11-06-2011, 08:18 PM   #3
yksikaksi
Member
 
Location: Scandinavia

Join Date: Dec 2009
Posts: 20
Default

Thanks Richard Finney.

I came across to this piece of modified solid2fastq.pl script by bgulko in bwa color-space index thread.

I did minor changed as suggested by drio's post in PE SOLiD reads alignment by bwa thread for pairend reads on this modified solid2fasta.pl script.

Without that minor changed on that script, it only converted the F3 csfasta to single end reads fastq.

So far, the script worked with no complaint. I got read1.fastq and read2.fastq after the converting task done. However, I'm not yet carry out the mapping task. I will do it later.
yksikaksi is offline   Reply With Quote
Reply

Tags
bwa, csfasta, fastq, solid, solid2fastq

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 09:10 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO