Biostar Beta. Not for public use.
Using parallel with Mem
0
Entering edit mode
2.9 years ago
@always_learning8364

Hi All.

We are trying to us MEM with GNU parallel with paired fastQ files, but its saying that can't accept more then one file. What could be the possible reasons? What if we concatante both file as single one and then use that with mem? Does Mem command support multi threading too ?

Thanks

Syed

mem BWA Parallel • 3.2k views
ADD COMMENT • link
0
Entering edit mode

show us the cmd-line please.

ADD REPLY • link
0
Entering edit mode

this exact commant we are using:

${bwa} mem -t 8 -R \\"@RG\\tID:${ursId}_pairend\\tSM:${sampleName}\\tLB:${libraryId}\\tPL:${platform}\\" ${refGenomeFastaFile} \${seqReadsDir}/${ursId}_pairend1.fastq.gz ${seqReadsDir}/${ursId}_pairend2.fastq.gz

ADD REPLY • link
0
Entering edit mode

but how do you invoke parallel ?

ADD REPLY • link
0
Entering edit mode

Sorry as I didn't able to post complete command. Here it goes as we are trying to convert SAM files to BAM files at same time only

${bwa} mem -t 8 -R "@RG\\tID:${ursId}_pairend\\tSM:${sampleName}\\tLB:${libraryId}\\tPL:${platform}" ${refGenomeFastaFile} \
${seqReadsDir}/${ursId}_pairend1.fastq.gz ${seqReadsDir}/${ursId}_pairend2.fastq.gz | ${samtools} view -bhS - > ${bamDir}/${ursId}.bam" > \
${seqReadsDir}/${ursId}_pairend.txt; parallel -a ${seqReadsDir}/${ursId}_pairend.txt&#39
ADD REPLY • link
4
Entering edit mode
22 months ago
@Matt Shirley1681

This command is telling bwa mem to use 8 threads (-t) so you will have some level of parallelism. Your command above does not invoke GNU parallel, and to do so you would use something like:

parallel bwa mem -t 8 -R "@RG\t@ID:blahblahblah" reference.fasta {1} {2} > {1.}.sam ::: R1.fastq.gz ::: R2.fastq.gz See here: https://www.gnu.org/software/parallel/man.html#example__use_multiple_inputs_in_one_command and here: https://www.biostars.org/p/63816/

ADD COMMENT • link
0
Entering edit mode

Hi Matt,

I am using gnu parallel as explained by you for cufflinks.Here is my post.

https://www.biostars.org/p/122961/

I want to add the prefix to my cufflink output files using filenames from bam samples.

ADD REPLY • link
0
Entering edit mode
2.9 years ago
Korsocius • 110
@Korsocius14475

There is snytax for parallel with bwa mem, but I can't load ID and SM from INPUT. Do you know, where could be a problem? 

ls *R1_001.fastq | parallel 'bwa mem -k 19 -A 1 -B 4 -O 6 -L 5 -R '@RG\tID:'{.}'\tSM:'{.}'\tLB:'Trusight_custom_amplicon_CARR'\tPL:'ILLUMINA'\tPI:150' $REFERENCE {} {= s/_R1_001/_R2_001/ =} > {= s/_R1_001.fastq/.sam/ =}'
ADD COMMENT • link
1
Entering edit mode

ask this as a new question please.

ADD REPLY • link

Login before adding your answer.

Similar Posts
Loading Similar Posts
Content Search
Users
Tags
Badges
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.3