Clustering RBH hits into files by reference gene

mfort001

Junior Member

Join Date: Oct 2011
Posts: 2

Clustering RBH hits into files by reference gene

06-20-2016, 07:36 AM

I have performed RBH for each of five transcriptomes against a reference genome, producing five separate output files. I figured, for ease, it was best to combine these output files into one using 'cat' command. The output looks like so, e.g.,:

Code:

    Acur_1001484	1751	OFAS000562-RA-EXON01	467	1	1039	1240	1	467	264	4.00E-62	237	87.75	179	204	1	264	467	1	1240	1039	1.00E-62	237	87.75	179	204
    Apil_comp17655_c0_seq1	2446	OFAS000562-RA-EXON01	467	1	288	453	1	289	454	2.00E-66	252	93.98	156	166	1	289	454	1	288	453	1.00E-66	252	93.98	156	166
    Atri_comp16713_c0_seq1	2069	OFAS000562-RA-EXON01	467	1	1576	1754	1	467	289	5.00E-67	254	92.18	165	179	1	289	467	1	1754	1576	4.00E-67	254	92.18	165	179
    Btri_comp12490_c0_seq1	2547	OFAS000562-RA-EXON01	467	1	2069	2267	1	467	268	2.00E-66	252	89.5	179	200	1	268	467	1	2267	2069	1.00E-66	252	89.5	179	200
    Atri_comp10477_c0_seq1	2608	OFAS000604-RA-EXON01	1887	1	4	1159	1	298	1537	0	701	78.5	997	1270	1	298	1537	1	4	1159	0	701	78.54	999	1272
    Ctom_1008131	415	OFAS000604-RA-EXON01	1887	1	101	399	1	1260	959	1.00E-81	300	85.06	262	308	1	959	1260	1	399	101	4.00E-81	300	85.06	262	308
    Acur_1001075	997	OFAS000604-RA-EXON01	1887	1	469	954	1	756	1258	1.00E-125	448	83.69	431	515	1	247	650	1	16	417	1.00E-121	435	86.54	360	416
    Btri_comp18753_c0_seq1	2786	OFAS000604-RA-EXON01	1887	1	420	1181	1	761	1593	1.00E-108	392	77.03	654	849	1	297	650	1	4	336	1.00E-87	324	84.49	305	361
    Apil_comp16297_c0_seq1	1235	OFAS000628-RA-EXON01	453	1	182	553	1	372	1	4.00E-66	250	78.93	296	375	1	1	372	1	553	182	4.00E-66	250	78.93	296	375
    Atri_comp15731_c0_seq1	657	OFAS000628-RA-EXON01	453	1	131	517	1	387	1	9.00E-55	211	76.73	300	391	1	1	387	1	517	131	2.00E-54	211	76.73	300	391
    Apil_comp14961_c1_seq1	1230	OFAS000635-RA-EXON02	212	1	254	452	1	1	198	2.00E-54	211	85.93	171	199	1	1	198	1	254	452	8.00E-55	211	85.93	171	199

For a given reference gene (e.g., OFAS000562-RA-EXON01), I would now like to retrieve the query IDs and the corresponding reference gene ID and copy this data into a new text file named as the reference gene ID (e.g., OFAS000562-RA-EXON01_IDs.txt). Below are four examples of what I'm wanting in terms of output from this approach:

File 1:

Code:

   OFAS000562-RA-EXON01
    Acur_1001484
    Apil_comp17655_c0_seq1
    Atri_comp16713_c0_seq1
    Btri_comp12490_c0_seq1

File 2:

Code:

   OFAS000604-RA-EXON01
    Acur_1001075
    Btri_comp18753_c0_seq1
    Atri_comp10477_c0_seq1
    Ctom_1008131

File 3:

Code:

   OFAS000628-RA-EXON01
    Apil_comp16297_c0_seq1
    Atri_comp15731_c0_seq1

File 4:

Code:

   OFAS000635-RA-EXON02
    Apil_comp14961_c1_seq1

I don't have strong bioinformatic/scripting skills, so I appreciate any feedback on how to do this.

Last edited by mfort001; 06-20-2016, 10:54 AM.

Tags: genome, reciprocal blast hits, transcriptomes

GenoMax

Senior Member

Join Date: Feb 2008

Posts: 7142
- Share
- Tweet
#2

06-20-2016, 07:42 AM

For reference cross-posted: https://www.biostars.org/p/196945/
Comment

Previous template Next

From Collection to Sequencing: Why Sample Preparation and Preservation Define Sequencing Data

by SEQadmin2

Data variability is still an issue in sequencing technologies despite the advances in reproducibility and accuracy of these platforms. But the problem does not originate in the sequencing itself, but in the previous steps, before the sample reaches the sequencer.

The first step is collection, followed by preservation and sample preparation for analysis. Most scientists overlook those steps, but not being careful might just be skewing the experiment’s results.
...
- Channel: Articles
Yesterday, 10:05 AM
Single-Cell Sequencing at an Inflection Point: Early Impacts of New Platforms and Emerging Trends

by SEQadmin2

With the launch of new single-cell sequencing platforms in 2026, the field stands at an exciting inflection point. This article surveys the most impactful advances in the field and discusses how they’re reshaping research in cancer, immunology, and beyond.

Introduction

Single-cell sequencing technologies have undergone remarkable advances over the past decade, transitioning from low-throughput experimental approaches to highly scalable platforms capable of...
- Channel: Articles
05-22-2026, 06:42 AM
Environmental Genomics in the Age of NGS: From Microbes to Conservation Strategies

by SEQadmin2

Studying ecosystems means dealing with complex, multi-species communities that are hard to observe at scale. This complexity, however, hides many important questions to be answered, from how biogeochemical cycles work and how climate change can affect species distribution to how conservation strategies can work best.

Genomics, particularly since the expansion of NGS, has transformed ecosystem ecology. By sequencing environmental DNA, we can now assess biodiversity without direct...
- Channel: Articles
05-06-2026, 09:04 AM

Topics	Statistics	Last Post
Long-Read RNA Sequencing Uncovers a Hidden Layer of Immune Cell Regulation by SEQadmin2 Started by SEQadmin2, Yesterday, 12:03 PM	0 responses 19 views 0 reactions	Last Post by SEQadmin2 Yesterday, 12:03 PM
DNA Methylation Study Reveals How Epigenetic Changes Pass Between Generations by SEQadmin2 Started by SEQadmin2, Yesterday, 11:40 AM	0 responses 14 views 0 reactions	Last Post by SEQadmin2 Yesterday, 11:40 AM
MetaBeeAI Helps Scientists Process Research Literature Faster by SEQadmin2 Started by SEQadmin2, 05-28-2026, 11:40 AM	0 responses 29 views 0 reactions	Last Post by SEQadmin2 05-28-2026, 11:40 AM
Scientists Solve a 25-Year Mystery in RNA Interference by SEQadmin2 Started by SEQadmin2, 05-26-2026, 10:12 AM	0 responses 31 views 0 reactions	Last Post by SEQadmin2 05-26-2026, 10:12 AM

Unconfigured Ad

Clustering RBH hits into files by reference gene

Comment

Latest Articles

ad_right_rmr

News