This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Circos input format

Does anyone know how to convert *.fasta format contig data into a circos input data format? What tools do I need to use?

How could I find the start position and end position of the sequence?

The input karyotype file of circos usually takes inputs as name, label, start and end position and a color in order.

Any help or suggestion would be appreciated!

data-format-conversion alignment
• 6,455 views • link
updated by Ram Ram • written by cdwilliam524 • cdwilliam524

Your first link does not work; can you please post a working link to your contig data?

• 0 views • link
written by Alex Reynolds Alex Reynolds

1 answer

Perhaps use command-line tools like faToTwoBit to build an indexed reference genome of interest, and then command-line BLAT and your 2bit and FASTA files to query that reference genome. Where it can find matches, BLAT will yield hits that include chromosome name and start/stop positions, which you can parse into input to feed into Circos. These tools are part of the Kent Tools source code package.

• 0 views • link
updated by Ram Ram • written by Alex Reynolds Alex Reynolds

I am currently trying tools like bowtie and bwa to index the reference genome, Candidatus Kuenenia stuttgartiensis; however, the output file .sam is not in its supposed result. I suppose to see actual genome information but I see mass code like NNNNNNNNNNNNCNNNNNNNNNNNNNNNNNN...

I don't think BLAT and BLAST have Candidatus Kuenenia stuttgartiensis (bacteria).

Thanks!

• 0 views • link
written by cdwilliam524 • cdwilliam524

The web version of BLAT does not, very probably, but you can definitely build your own index files and query against them with the command-line tools, if you have the reference genome somewhere.

• 0 views • link
written by Alex Reynolds Alex Reynolds

Yes, I have the reference genome. How could I build my own index and query against it? Could you suggest me some tools and how to do it? What are the command-line tools? It is part of Xcode on Mac?

• 0 views • link
written by cdwilliam524 • cdwilliam524

See: http://genome.ucsc.edu/FAQ/FAQblat.html#blat3

You will need a compiler installed to build these tools. If you are using OS X, you will need to install Xcode and then install the command-line tools via that app. Then you can download the Kent Tools source code and compile it to get blat and faToTwoBit and other tools.

You would take your FASTA-formatted reference genome and convert it to 2bit format.

$ faToTwoBit myAssembly.fa myAssembly.2bit

Then, you might do something like:

$ blat myAssembly.2bit -oneOff=0 -noHead myQueryContigSeqs.fa myHits.psl

The output file myHits.psl is in PSL format. You can convert to BED with psl2bed and cut -f1-3 to grab the first three columns, or just read the specs for PSL and use cut or awk etc. to grab those columns.

• 0 views • link
updated by Ram Ram • written by Alex Reynolds Alex Reynolds

Thanks!

I am working on an ANAMMOX project. I don't need to compare the Kuenenia stuttgartiensis genome against human genome sequence.

Have you had experiences in bwa and Samtools? I generated the *.sam file but it has so many gaps. I need to remove the gaps and find the start and end positions of each alignment.

• 0 views • link
updated by Ram Ram • written by cdwilliam524 • cdwilliam524

Log in to answer this question.

More posts like this