This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Make Paired File Task-Scripture

Hello, I am trying to prepare the input file for Scripture and have stumbled upon the fact that I need to have the TopHat alignment output split into right and left reads. Can somebody tell me how to do this? I obviously had a bam file,containing both right and left reads. Thanks, G.

• 2,324 views • link
updated by Biostar Biostar • written by GPR GPR

1 answer

I think you could get that information from the SAM file. the spec says that column 8 PNEXT indicates the position of the mate/next segment.

So for each line if the position is larger than the current position then it is a left read otherwise a right read.

• 0 views • link
written by Istvan Albert Istvan Albert

Thanks a lot, I will look into this.

• 0 views • link
written by GPR GPR

Log in to answer this question.

More posts like this
  • Tophat, low alignement left read
    written by claire.morandin • claire.morandin
    Hello, I downloaded a bunch of RNA-seq dataset from Genbank and used tophat (with Bowtie2) to map back to the genome (to create a gtf …
  • Strand Specific Tophat Output
    written by GPR GPR
    <p>Hello, I am having a hard time coming about the possible bad quality of a data-set I just got. The library is paired-end and strand-specific …
  • Tophat2, Strand-Specific And Genes.Gtf
    written by GPR GPR
    <p>Hello, I have been running Tophat2 for a while and haven't had issues with it until now that my input is a paired-end, strand-specific data …
  • Varscan Mpileup2Snp Output Missing Varfreq Values
    written by GPR GPR
    <p>Hello, I am running an SNP analysis and have two questions about the output.</p> <p>My command-line is the following: samtools mpileup -f genome.fa left-reads.bam right …
  • How To Generate A Paired End File Suitable For Bwa
    written by GPR GPR
    <p>A question about BWA and paired-end reads. In TopHat, one inputs the paired-end fastq file twice and the tool understands you have paired-end reads. I …
  • Running Tophat Over Multiple Files
    written by narges • narges
    <p>Hello all</p> <p>I have a list containing left mate sample reads and a list for right mate sample reads like: list1 = sample1<em>1.fastq , ... …
  • Write Diffgenes/Diffgeneids In A Tab File
    written by GPR GPR
    <p>Hello, I am relatively new to CummeRbund. Can somebody tell me how to write a tab file containing diffGenes/diffGeneIDs? Thanks G.</p>
  • Paire-End Alignments And Scripture
    written by GPR GPR
    <p>Hello, I am slowly getting my hands on Scripture and have stumbled on yet another unclear (at least to me )detail. I have paired-end data …
  • Scripture: One Chromosome At A Time?
    written by GPR GPR
    <p>Hello, I am new to Scripture and am trying to understand the command-line arguments. Am I right in that Scripture only reconstructs transcripts one chromosome …
  • Scripture And Transcript Abundance
    written by GPR GPR
    <p>Hello, After gaped read alignment with TopHat, I plan to compare the transcript reconstruction outputs of Cufflinks and Scripture, and potentially combine them. My question …