This is a test version of Biostars. For the public version, visit https://www.biostars.org.
mapping transcriptome to genome to find TSS/TES/(+/-)

I have a transcriptome fasta file that looks like this: 

>loc.00001.1 gene=loc.00001
TCATTTTCCTGGCATGGATAATTTCGCTGCAAAAGGTGTCGCACCTGAAACAGGAGAAATGCATTAGCAG
CTGTCAAATTACCACGAATTTCGGAGATTATGACGACATTTTAGAATGCGACCAATCGAAACTCGAATCG
AGAAATCTCTATGAGATATTTTCATTTCTCACATTGATAACGGGAATCGAGGAATTCCATGAGCAGAACA
TTTGGATTAAATGCTTTCACGATTTGGTGACCGATTTCATCATCTACATGCCGAAAAGGCTCAAGGATTT

I also have a genome fasta file that is around -500 contigs long.

  > contig_1  
attgacaatttaactcgagcattgttctcgtttacataacattatctcatatgtcctttaccacatctttgtagtttttctggctactcctgggataat

I want to map the transcriptome to the genome such that I can produce a resultant bed file that has the transcript ID, contig ID, TSS, TES, and +/- orientation?

Would I have to use a program like gmap? Would that give me the information I need?

mapping genome alignment
• 1,792 views • link
written by a.rex a.rex

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • Adding gene names to fasta file
    written by mm2568 • mm2568
    I have file 1 (a FASTA file): >dmel_X type=golden_path_region; loc=X:2270008..2271068; ID=X; dbxref=GB:AE014298,GB:AE014298,REFSEQ:NC_004354; release=r6.37; species=Dmel; CTCTTTGTTAGCCTACGCTTTCTGCTGAGTTTGTTATTTTTGTCTGCTCCCCACAAGGATATTGTTACAGAGAAAAAGCT CGAATTGAAGGGAAAATGGAGACAAATAAGAAAACCCATGACAAAGAGGAAAGTTTCAAATATGGGCAATCGAAAAAATC GAGAAGTGAGCCAATTTTTTTTTCGCCGAGGCTCCACTGTTCCCAGCTGCATAACTGTTTTCCCTCGGCACCTCTCTTTT >dmel_3L type=golden_path_region; loc=3L:20341634..20342694; ID=3L; dbxref=GB:AE014296,GB:AE014296,REFSEQ:NT_037436; release=r6.37; species=Dmel; ATTAGTATATAGGCATATGCTTAAGTCTTAGGGTCTTATGGATATGTCACTATATATATATATAATTGCATAAATAGAGA …
  • WGS upload prep - table2asn
    written by ian.will • ian.will
    Hi, I've wrapped up an assembly and will soon be uploading a new genome and annotations to NCBI - but am having a little trouble …
  • Perform GO analysis on genome
    written by a.rex a.rex
    I have a protein fasta file, for each longest ORF for an annotated locus in my genome. Could anyone advise how best to perform a …
  • gmap_build database problem
    written by Kenny • Kenny
    Hi all, I have a scaffold sequence and transcriptome. I am currently running GMAP to map and align the transcriptome to the scaffold genome. First …
  • select -500 from the TSS to +500 after the TES strand +/-
    written by Lila M Lila M
    Hi everybody! Hi have a terrible doubt, I want to select the coverage over -500 from the TSS to +500 after the TES, from my …
  • ncbi blast results in multi-sequence fasta file
    written by a.rex a.rex
    I have a fasta file that has 1000s of sequences, i.e.: >loc.1 ATGGTTTTTCCGTATACTTCACTGACTGCTGCTTGTAATTTTTCAGCATCCTCAAATTTTCAATTGTCAGAAACATGTCTAAACAATGAGGTCCCTCCAGTATTTATAAATAGCAACGGCCAAAATTGTTCCCGGCCCAATATGTCTTATAATGCTGACTTTCAAGTTTCAAATAGTCATTGTAATGAAACCACTGACAGTATTGGTAGTGGTCAAAATACTACTACTGATATGAATTACGATCCAAATAATAGTCAGAACTTTTCATTTTCCTCAAATGTTCTTGGTAATCTACAAAACTGGAATGGTAAAAGATCTAATTATTTCAGTTACAAACTCAATGACATGAAACAATTTTATAATCAAGAAATACCGTTAGTGGACAATTCCGTACCGATTTACACAAATGG >loc.2 CAGGTAAAAATCTGGTTTCAAAATCGTCGGTCCAAGTATAAAAAGCTTATTAAGCAAGGTCAGGATCCAAGCATCCTGATGAATGGAGAATTTAATGACAGCATGGATGAAATGACGGAAGATCAAATCGACGAAGATAATTGTATCAAACCAAAAGCAGAAATGCTATTAACAAGCGATCCAAATAATCCTAGAGGCGATAGTTCTGATATTCCAACTGAAA I am using ncbi blastn in terminal against this file …
  • append to column in bedfile
    written by a.rex a.rex
    I have a bed file that looks like this: scaffold1 206 262 loc.00001 + . scaffold1 749 1391 loc.00001 + . scaffold1 1549 1595 loc.00001 …
  • annotating transposons in transcriptome
    written by a.rex a.rex
    I have recently run a relatively obscure genome through RepeatMasker and RepeatModeler, and produced a bedfile whereby I have all the instances of a TE …
  • parsing bed file
    written by a.rex a.rex
    Apologies if a similar question has been posed before - I was unable to find anything helpful. I have a bed file that looks like …
  • parsing through gtf file
    written by a.rex a.rex
    I have the following gtf file layout, with the 'features' transcript (i.e. full length of the transcript) and the exons within that transcript. For example: …