Hi, I have 2 bam files from which the unmapped sequences have to be extracted and aligned against each other to identify unaligned sequences. The purpose is to identify the presence of a contamination in one sample by comparing with other.
Now I have fasta sequences of unmapped reads which am planning to blast against each other considering one as subject and other as query. Both files are of about ~3 gb size. But this is getting killed in my server which has 64 gb ram and 6 tb space.
Can somebody suggest how I can do this effectively in the available ram and server space ? I read about splitting the fasta files but couldn't find a proper tutorial. Hope someone would guide me on this.
Thank you.
0 answers
No answers yet.
Log in to answer this question.
-
GLnexus import issues with 4764 samples, WGS
written by Claire Weber •
-
"sortmerna" quitting early--is this is a memory issue?
written by Andrew •
-
Questions about per base sequence quality and GC content.
written by bioinfo •
-
blastn error : mdb_env_open: No such file or directory
written by Sowmya Pulapet •
-
BLASTn produces empty output + output is 'Killed'
written by Shraddha •
-
Fastest way to perform BLAST search using a multi-FASTA file against a remote database
written by accibio •
-
Blastn, need help to increase speed
written by chiachoong_leong93 •
-
Faster way to get sequences from large fasta files?
written by emilywynn6 •
-
which database to blast against genomes.
written by Jacob •
-
Buying computational resources vs. renting allocation on cluster
written by ebrown1955
