SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Segmentation fault, cuffdiff labrat73 Bioinformatics 12 03-22-2018 11:00 AM
soapsnp segmentation fault - help! jtjli Bioinformatics 6 02-29-2016 08:54 AM
Segmentation fault (core dumped) at contig step during SOAP denovo assembly tangzhonghui Bioinformatics 1 10-09-2012 05:32 PM
soap segmentation fault scami Bioinformatics 6 04-17-2012 06:08 AM
Newbler segmentation fault flobpf Bioinformatics 4 04-18-2011 11:45 AM

Reply
 
Thread Tools
Old 09-09-2010, 01:40 PM   #1
Mansequencer
Member
 
Location: Canada

Join Date: Jul 2010
Posts: 16
Default SOAP:segmentation fault

Hi All,
I am trying to map solexa PE reads (fastq) from different strains onto reference scaffolds file (fasta) using soap2. The questions I have are:
1. The alignment works but it seems that soap is not able to read seq data from all the 8 strains I have in the command as the number of reads that it processes is too low to encompass all the 8 strains.
2. When I try to msort (-k 8,n9) the PE output file from the alignment, it comes back with error (segmentation fault) however, the SE output gets sorted with the same command.
3. The snp calling through soapsnp (-d <ref> -i <SEoutput.sort> -r 0.00005 e 0.0001 -t -u -L <100>) generates a huge file many times larger than original file. Is that usual?
Thanks for the patient reading.
Mansequencer is offline   Reply With Quote
Old 11-19-2010, 10:47 AM   #2
Awesome
Junior Member
 
Location: california

Join Date: Aug 2009
Posts: 7
Thumbs up

Hi Mansequencer,

For some reason, MSORT doesn't work with files over a certain size. It works just fine for small files using the command:
Code:
$ msort -k 8,n9 mapped.out > mapped.out.sort
To get around this, you can use the SORT command that is native to most Linux and Unix distributions:
Code:
sort -t $'\t' -k 8f,8 -k 9n,9 mapped.out > mapped.out.sort
To check to see if your files are indeed sorted, you can always isolate just the two columns in question, the chromosome and position:
Code:
cut -f 8,9 mapped.out.sort > col_89.txt
SOAPSNP generates a huge file on purpose. Either of the other formats GLFv2 and GPFv2 may give a smaller footprint output. Try those.

Best of luck,
-Awesome
Awesome is offline   Reply With Quote
Reply

Tags
msort, segmentation fault, soap, soapsnp

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:19 PM.


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