This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Which Kinds Of Pair-End Reads Direction Do We Have Now?

About the pair-end / mate-pair read, because of different technology they use, the normal read direction may be different.

As far as I know, there are directions --> <-- and <-- --> for different tech.

Is there any other possible direction of today's sequencing technology?

• 2,913 views • link
updated by Biostar Biostar • written by Chen Sun Chen Sun

1 answer

SOLiD mate pair reads are oriented in the same direction: R3 ------> F3 ------>

Bowtie provides several options for paired end reads. Below is taken from bowtie2 manual:


--fr/--rf/--ff
The upstream/downstream mate orientations for a valid paired-end alignment against the forward reference strand. E.g., if --fr is specified and there is a candidate paired-end alignment where mate 1 appears upstream of the reverse complement of mate 2 and the fragment length constraints (-I and -X) are met, that alignment is valid. Also, if mate 2 appears upstream of the reverse complement of mate 1 and all other constraints are met, that too is valid. --rf likewise requires that an upstream mate1 be reverse-complemented and a downstream mate2 be forward-oriented. --ff requires both an upstream mate 1 and a downstream mate 2 to be forward-oriented. Default: --fr (appropriate for Illumina's Paired-end Sequencing Assay).

• 102 views • link
written by Ashutosh Pandey Ashutosh Pandey

Log in to answer this question.

More posts like this
  • In Sam File Flag Information, How Could Unmapped Read Have A Direction?
    written by Chen Sun Chen Sun
    <p>In <a href='http://samtools.sourceforge.net/SAM1.pdf'>SAM</a> file I get, there are some flag value that dose not make sense:</p> <pre><code>101 1+4+32+64 means first read, unmapped with direction+ ; …
  • How To Find The Most Frequent Allel Snp In Vcf File
    written by Chen Sun Chen Sun
    <p>When I use <code>samtools mpileup</code>, there are usually several allel SNP in ALT column according to REF column. just as follows:</p> <pre><code>##fileformat=VCFv4.1 ##samtoolsVersion=0.1.19-44428cd ##reference=file:///gpfs/home/cxs1031/backup/ref/chr22.fa ##contig=&lt;ID=chr22,length=51304566&gt; …
  • Samtools Mpileup And Overlapping Paired-End Reads
    written by alpha2zee • alpha2zee
    <p>I am using <a href='http://samtools.sourceforge.net/'>samtools</a> mpileup to generate a pileup for paired-end RNA sequencing data. I am curious about how <a href='http://samtools.sourceforge.net/'>samtools</a> handles pair mates …
  • Reverse Read Direction In A Bam File
    written by lin.barnum • lin.barnum
    <p>I have a <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> alignment file with SOLiD mate-pair data (R3/F3). I want to reverse the second read so that it points towards the …
  • Solid Mate-Pair Alignment With Bwa
    written by lin.barnum • lin.barnum
    <p>I have aligned SOLiD mate-pair reads using <a href='http://bio-bwa.sourceforge.net/'>bwa</a>-0.5.9. Now, it is my understanding that SOLiD mate-pair reads are oriented in the same direction</p> <pre><code>R3 …
  • Paired-End Sam Flag Mismatch During Filtering
    written by Andreas Andreas
    <p>Hi BioStar users,</p> <p>while removing human contamination from paired-end metagenomic Illumina-shotgun data I ran into a strange problem. In a first step I mapped the …
  • Mpileup Output And Quality Scores
    written by NextGenSeek • NextGenSeek
    <p>I am looking at mpileup output from RNA-Seq data, specifically at a location that looks like this.</p> <pre><code>&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&lt;&gt;&gt;&gt;&gt;&gt;&gt;&lt;&lt;&lt;&gt;&gt;&lt;&gt;&gt;&lt;&gt;&gt;&gt;&gt;&lt;&lt;&lt;&lt;&lt; G&gt;BB#GHJCCD#@#5#E;F##ICFEBIHDBDD;BB?IJGGGHJC##?EC </code></pre> <p>My understanding is that the …
  • Tophat Parameters Misunderstanding
    written by Nicolas Rosewick Nicolas Rosewick
    <p>Hi,</p> <p>I have a little understanding problem on several tophat parameters :</p> <p>-r/--mate-inner-dist <int> This is the expected (mean) inner distance between mate pairs. For, …
  • Unmapped Reads And Sequence Name
    written by Stefano Berri Stefano Berri
    <p>Hi. I have two questions that might or might not be related.</p> <p>I want to filter all and only unmapped reads from a pair end …
  • Rescuing Orphaned Reads With Local Alignment Within A Radius Of Mapped Mate
    written by Abhi Abhi
    <p>PS: this message has also been cross posted on seqanswers. I just want to reach out to more bioinfo guys so thought of posting it …