Header Leaderboard Ad

Collapse

samtools faidx Segmentation fault

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • totalnew
    replied
    I had this problem too, my solution was to use unix command to trim and fold the fasta file. You would have to cut the header first, and catenate it with your sorted fasta file. It perfectly solves the problem.

    Writing a shell script may be a good idea to make things easlier.

    Leave a comment:


  • webbrewer
    replied
    This bioperl snippet fixes the fasta:

    Code:
    use Bio::SeqIO;
    $in  = Bio::SeqIO->new(-file => "inputfilename",
                           -format => 'Fasta');
    $out = Bio::SeqIO->new(-file => ">outputfilename",
                           -format => 'Fasta');
    while ( my $seq = $in->next_seq() ) {$out->write_seq($seq); }

    Leave a comment:


  • webbrewer
    replied
    I am having this same issue. While I have verified that there are lines of different length in the sequence file, however, why should this matter?

    Leave a comment:


  • nilshomer
    replied
    Originally posted by baohua100 View Post
    [email protected]:~/samtools-0.1.5c_x86_64-linux$ ./samtools faidx /media/Poplar/baohua/genome/poplar_genome.fa
    [fai_build_core] different line length in sequence 'scaffold_28'.
    Segmentation fault

    What's the meaning of defferent line length ?

    It's the standard fasta file.

    I download it from:
    ftp://ftp.jgi-psf.org/pub/JGI_data/P...asked.fasta.gz
    Can you verify that every sequence line in the FASTA has the same length?

    Leave a comment:


  • baohua100
    started a topic samtools faidx Segmentation fault

    samtools faidx Segmentation fault

    [email protected]:~/samtools-0.1.5c_x86_64-linux$ ./samtools faidx /media/Poplar/baohua/genome/poplar_genome.fa
    [fai_build_core] different line length in sequence 'scaffold_28'.
    Segmentation fault

    What's the meaning of defferent line length ?

    It's the standard fasta file.

    I download it from:
    ftp://ftp.jgi-psf.org/pub/JGI_data/P...asked.fasta.gz
Working...
X