This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Concatenating all unmapped or alternative contigs in a reference genome

Hi, I want to use ConTExt to analyse TEs in Drosophila. To use the tool, I need to concatenate all the unmapped and alternative contigs in the reference genome to reduce the number of output files. How should I do this? Thank you

contigs fasta reference
• 1,892 views • link
updated by Semir • Semir • written by manikin_python9f • manikin_python9f

cat contig1.fa contig2.fa .. contigN.fa > one_file.fa will concatenate the files into one.

• 0 views • link
written by GenoMax GenoMax

Or if you can find any sort of common identifier (here assuming its *.fa) you can automate find . -maxdepth 1 -name '*fa' | xargs cat > one_file.fa.

• 0 views • link
written by ATpoint ATpoint

Thank you. I don’t want to cat some files to one file. I want to put all the unmapped and alternative contigs together, under one header, and add it to the reference genome.

• 0 views • link
written by manikin_python9f • manikin_python9f

2 answers

Then use one of the solutions here: HOw to merge multifasta sequence into a single sequence having only one header?
https://stackoverflow.com/questions/69471751/how-to-concatenate-sequences-in-the-same-multifasta-files-and-then-print-result

Once you have this file you can cat it at the end of original reference genome.

• 0 views • link
written by GenoMax GenoMax

I'd also like to share an alternative answer that might be less cryptic than the sed command:

# First concatenate all fa files in the directory
seqkit concat -j ${task.cpus} --full ${fasta_dir}/*.fa > temp_concat.fa

# Replace all headers with a single header and remove line breaks in sequence
seqkit replace -j ${task.cpus} -p "^>.+" -r ">combined_sequence" temp_concat.fa | \
seqkit replace -s -p "\n" -r "" > combined.fa

# Combine with reference genome
cat ${ref_genome} combined.fa > final_combined.fa

# Clean up temporary files
rm temp_concat.fa combined.fa

I generated this answer using amplicon.ai, a tool I've been building to iteratively write and execute pipelines easier. Feel free to try it out

enter image description here

• 0 views • link
written by Semir • Semir

Log in to answer this question.

More posts like this