text: user989: I recently had some Iso-Seq sequencing done on my organism catfish on the new Sequel platform and got weird alignments for a size selected 4 Kilobase and up fraction after running the isoseq3 pipeline.  When I investigated the exon intron junctions using data from (Pac Bio Iso-Seq, Illumina RNA-seq) against our private reference genome, I get PacBio Iso-Seq reads where the sequence length is greater than the seq length on the reference for my bam alignments:
  Seq length: 6751
  Seq length on ref: 5267
  When blasting the sequence in NCBI it shows 5% similarity to the publicly-available genome.  I'll attach a JBrowse image with the same RNA sequenced with Illumina paired RNA-seq.  I also included a image and circled the pacbio contig and Illumina is above.
  
  Question 1:
  Is this a problem with JBrowse?
  Question 2:
  Does this mean that some DNA was sequenced with the catfish cDNA?
  Question 3:
  Why is this aligned sequence's length longer than the length of the reference genome region to which it is aligned?
  Thanks
  
  user2085: Q1
  Probably not - this is probably just displaying the contents of the input (I presume a bam file) exactly as they are encoded.
  Q2
  The splicing looks a little weird, but the exon-mapping pattern in the 5' part of the gene makes it seem like RNA rather than DNA. I think that the large intron is what makes you think that it is DNA, correct?
  
  Explanations for weird splicing
  
  The weird read is a rare splice variant but a true mRNA.
  The weird read is not completely matured and has an intron that has not yet been spliced. The IsoSeq library prep protocol picked this up somehow and it ended up being sequenced.
  
  
  
  Q3
  
  Explanations for why Seq length: 6751 and Seq length on ref: 5267
  
  There is probably a large amount of clipping on the 5' or 3' end of the read. This could be polyA or some other transcribed-but-not-yet-completely-matured seq. Or, the extra sequence could be sequencing error, PacBio adapters, or a chimeric molecule.
  
  
  
  things to try
  How to check what is going on:
  I'm guessing that you had to upload a bam file to JBrowse to make the IsoSeq track. Use the sequence name that you got from Jbrowse and run the following command to extract the alignment info from the bam file:
  samtools view /path/to/your/isoseq/file.bam | grep 'weird_read_id'
  
  Look at the cigar string to see what it says about 5' or 3' clipping and corroborate with the sequence data.
  
original_text: 989: <p>I recently had some Iso-Seq sequencing done on my organism catfish on the new Sequel platform and got weird alignments for a size selected 4 Kilobase and up fraction after running the isoseq3 pipeline.  When I investigated the exon intron junctions using data from (Pac Bio Iso-Seq, Illumina RNA-seq) against our private reference genome, I get PacBio Iso-Seq reads where the sequence length is greater than the seq length on the reference for my bam alignments:</p>
  
  <p><strong>Seq length:</strong> 6751</p>
  
  <p><strong>Seq length on ref:</strong> 5267</p>
  
  <p>When blasting the sequence in NCBI it shows 5% similarity to the publicly-available genome.  I'll attach a JBrowse image with the same RNA sequenced with Illumina paired RNA-seq.  I also included a image and circled the pacbio contig and Illumina is above.</p>
  
  <p><a href="https://i.stack.imgur.com/vtf1j.jpg" rel="noreferrer"><img src="https://i.stack.imgur.com/vtf1j.jpg" alt="Pac Bio Iso-seq and Illumina RNA-seq"></a></p>
  
  <h2>Question 1:</h2>
  
  <p>Is this a problem with JBrowse?</p>
  
  <h2>Question 2:</h2>
  
  <p>Does this mean that some DNA was sequenced with the catfish cDNA?</p>
  
  <h2>Question 3:</h2>
  
  <p>Why is this aligned sequence's length longer than the length of the reference genome region to which it is aligned?</p>
  
  <p>Thanks</p>
   2085: <h3>Q1</h3>
  <p>Probably not - this is probably just displaying the contents of the input (I presume a <code>bam</code> file) exactly as they are encoded.</p>
  <h3>Q2</h3>
  <p>The splicing looks a little weird, but <strong>the exon-mapping pattern in the 5' part of the gene makes it seem like RNA rather than DNA</strong>. I think that the large intron is what makes you think that it is DNA, correct?</p>
  <ul>
  <li>Explanations for weird splicing
  <ul>
  <li>The weird read is a rare splice variant but a true mRNA.</li>
  <li>The weird read is not completely matured and has an intron that has not yet been spliced. The IsoSeq library prep protocol picked this up somehow and it ended up being sequenced.</li>
  </ul>
  </li>
  </ul>
  <h2>Q3</h2>
  <ul>
  <li>Explanations for why <strong>Seq length: 6751</strong> and <strong>Seq length on ref: 5267</strong>
  <ul>
  <li>There is probably a large amount of clipping on the 5' or 3' end of the read. This could be polyA or some other transcribed-but-not-yet-completely-matured seq. <strong>Or</strong>, the extra sequence could be sequencing error, PacBio adapters, or a chimeric molecule.</li>
  </ul>
  </li>
  </ul>
  <h3>things to try</h3>
  <p>How to check what is going on:</p>
  <p>I'm guessing that you had to upload a <code>bam</code> file to JBrowse to make the IsoSeq track. Use the sequence name that you got from Jbrowse and run the following command to extract the alignment info from the bam file:</p>
  <pre><code>samtools view /path/to/your/isoseq/file.bam | grep 'weird_read_id'
  </code></pre>
  <p>Look at the <a href="http://bioinformatics.cvr.ac.uk/blog/tag/cigar-string/" rel="noreferrer">cigar string</a> to see what it says about 5' or 3' clipping and corroborate with the sequence data.</p>
  