Skip to content
Snippets Groups Projects
BamView.java 89.5 KiB
Newer Older
  • Learn to ignore specific revisions
  • tjc's avatar
    tjc committed
              if(snpPanel != null)
                snpPanel.repaint();
    
    tjc's avatar
    tjc committed
            }
          });
          bottomPanel.add(scrollBar, BorderLayout.SOUTH);
        }
    
        mainPanel.add(bottomPanel, BorderLayout.SOUTH);
    
        coveragePanel.setPreferredSize(new Dimension(900, 100));
        coveragePanel.setVisible(false);
    
    tjc's avatar
    tjc committed
        snpPanel.setPreferredSize(new Dimension(900, 100));
        snpPanel.setVisible(false);
    
    
    tjc's avatar
    tjc committed
        jspView.getVerticalScrollBar().setValue(
            jspView.getVerticalScrollBar().getMaximum());
    
    tjc's avatar
    tjc committed
        jspView.getVerticalScrollBar().setUnitIncrement(maxUnitIncrement);
        
    
    tjc's avatar
    tjc committed
        if(feature_display == null)
    
    tjc's avatar
    tjc committed
        {
    
          addKeyListener(new KeyAdapter()
    
    tjc's avatar
    tjc committed
          {
    
            public void keyPressed(final KeyEvent event)
    
    tjc's avatar
    tjc committed
            {
    
              switch (event.getKeyCode())
              {
    
    tjc's avatar
    tjc committed
              case KeyEvent.VK_UP:
    
                setZoomLevel((int) (BamView.this.nbasesInView * 1.1));
    
    tjc's avatar
    tjc committed
                break;
              case KeyEvent.VK_DOWN:
    
                if (showBaseAlignment)
                  break;
    
                setZoomLevel((int) (BamView.this.nbasesInView * .9));
    
    tjc's avatar
    tjc committed
                break;
              default:
                break;
    
    tjc's avatar
    tjc committed
            }
    
        addMouseListener(new PopupListener());
    
    tjc's avatar
    tjc committed
        setFocusable(true);
        requestFocusInWindow();
    
      
      private void addToViewMenu(final int thisBamIndex)
      {
        File f = new File(bamList.get(thisBamIndex));
        final JCheckBoxMenuItem cbBam = new JCheckBoxMenuItem(
                                         f.getName(), true);
        bamFilesMenu.add(cbBam);
        cbBam.addActionListener(new ActionListener()
        {
          public void actionPerformed(ActionEvent e)
          {
            if(cbBam.isSelected())
              hideBamList.remove(new Integer(thisBamIndex));
            else
              hideBamList.add(new Integer(thisBamIndex));
            laststart = -1;
            lastend = -1;
            repaint();
          } 
        });
      }
    
    tjc's avatar
    tjc committed
    
    
      private void createMenus(JComponent menu)
    
      {
        final JMenuItem addBam = new JMenuItem("Add Bam ...");
        menu.add(addBam);
        addBam.addActionListener(new ActionListener()
        {
          public void actionPerformed(ActionEvent e)
          {
    
            FileSelectionDialog bamFileSelection = new FileSelectionDialog(
                null, false, "BamView", "BAM");
    
            List<String> bamFiles = bamFileSelection.getFiles(".*\\.bam$");
    
            int count = bamList.size();
           
    
            bamList.addAll(bamFileSelection.getFiles(".*\\.bam$"));
    
            
            for(int i=0; i<bamFiles.size(); i++)
              addToViewMenu(i+count);
            laststart = -1;
            lastend = -1;
            
            repaint();
          } 
        });
        
        bamFilesMenu.setFont(addBam.getFont());
        menu.add(bamFilesMenu);
        for(int i=0; i<bamList.size(); i++)
          addToViewMenu(i);
        
        menu.add(new JSeparator());
        
    
        JMenu viewMenu = new JMenu("Views");
    
    tjc's avatar
    tjc committed
        ButtonGroup group = new ButtonGroup();
    
        final JCheckBoxMenuItem checkBoxPairedStackView = new JCheckBoxMenuItem("Paired Stack View");
    
    tjc's avatar
    tjc committed
        final JCheckBoxMenuItem checkBoxStrandStackView = new JCheckBoxMenuItem("Strand Stack View");
    
    tjc's avatar
    tjc committed
        final JCheckBoxMenuItem checkIsizeStackView = new JCheckBoxMenuItem("Inferred Size View", false);
    
        checkBoxStackView.setFont(checkIsizeStackView.getFont());
    
        baseQualityColour.setFont(checkIsizeStackView.getFont());
    
        colourByCoverageColour.setFont(checkIsizeStackView.getFont());
    
        markInsertions.setFont(checkIsizeStackView.getFont());
        
    
    tjc's avatar
    tjc committed
        group.add(checkBoxStackView);
        group.add(checkBoxPairedStackView);
        group.add(checkBoxStrandStackView);
        group.add(checkIsizeStackView);
        
        checkIsizeStackView.addActionListener(new ActionListener()
        {
          public void actionPerformed(ActionEvent e)
          {
            laststart = -1;
            lastend = -1;
            
            if(checkIsizeStackView.isSelected())
            {
              isStackView = false;
              isPairedStackView = false;
              isStrandStackView = false;
    
              logMenuItem.setEnabled(true);
    
    tjc's avatar
    tjc committed
            }
            repaint();
          }
        });
    
        viewMenu.add(checkIsizeStackView);
    
    tjc's avatar
    tjc committed
        checkBoxStackView.addActionListener(new ActionListener()
        {
          public void actionPerformed(ActionEvent e)
          {
            laststart = -1;
            lastend = -1;
            isStackView = !isStackView;
    
            
            if(isStackView)
            {
              isPairedStackView = !isStackView;
    
    tjc's avatar
    tjc committed
              isStrandStackView = !isStackView;
    
              checkBoxPairedStackView.setSelected(!isStackView);
    
    tjc's avatar
    tjc committed
              checkBoxStrandStackView.setSelected(!isStackView);
    
              logMenuItem.setEnabled(false);
    
    tjc's avatar
    tjc committed
            repaint();
          }
        });
    
        viewMenu.add(checkBoxStackView);
    
    
        checkBoxPairedStackView.addActionListener(new ActionListener()
        {
          public void actionPerformed(ActionEvent e)
          {
            laststart = -1;
            lastend = -1;
            isPairedStackView = !isPairedStackView;
            
            if(isPairedStackView)
            {
              isStackView = !isPairedStackView;
    
    tjc's avatar
    tjc committed
              isStrandStackView = !isPairedStackView;
    
              checkBoxStackView.setSelected(!isPairedStackView);
    
    tjc's avatar
    tjc committed
              checkBoxStrandStackView.setSelected(!isPairedStackView);
    
              logMenuItem.setEnabled(false);
    
        viewMenu.add(checkBoxPairedStackView);
    
    tjc's avatar
    tjc committed
        checkBoxStrandStackView.addActionListener(new ActionListener()
        {
          public void actionPerformed(ActionEvent e)
          {
            laststart = -1;
            lastend = -1;
            isStrandStackView = !isStrandStackView;
            
            if(isStrandStackView)
            {
              isStackView = !isStrandStackView;
              isPairedStackView = !isStrandStackView;
              checkBoxStackView.setSelected(!isStrandStackView);
              checkBoxPairedStackView.setSelected(!isStrandStackView);
              setViewportMidPoint();
    
              logMenuItem.setEnabled(false);
    
    tjc's avatar
    tjc committed
            }
            repaint();
          }
        });
    
        viewMenu.add(checkBoxStrandStackView);  
        menu.add(viewMenu);
    
    tjc's avatar
    tjc committed
        final JCheckBoxMenuItem checkBoxSNPs = new JCheckBoxMenuItem("SNP marks");
    
        // 
        JMenu colourMenu = new JMenu("Colour By");
        colourMenu.add(colourByCoverageColour);
        
        baseQualityColour.addActionListener(new ActionListener()
        {
          public void actionPerformed(ActionEvent e)
          {
            if(baseQualityColour.isSelected())
            {
              checkBoxSNPs.setSelected(false);
              isSNPs = false;
            }
            repaint();
          }
        });
        colourMenu.add(baseQualityColour);
        menu.add(colourMenu);
        
        //
    
        JMenu showMenu = new JMenu("Show");
        JCheckBoxMenuItem checkBoxSingle = new JCheckBoxMenuItem("Single Reads");
        checkBoxSingle.addActionListener(new ActionListener()
        {
          public void actionPerformed(ActionEvent e)
          {
            repaint();
            isSingle = !isSingle;
          }
        });
        showMenu.add(checkBoxSingle);
    
        checkBoxSNPs.addActionListener(new ActionListener()
        {
          public void actionPerformed(ActionEvent e)
          {
            if (isSNPs && bases == null)
            {
              JOptionPane.showMessageDialog(null,
                  "No reference sequence supplied to identify SNPs.", "SNPs",
                  JOptionPane.INFORMATION_MESSAGE);
            }
            isSNPs = !isSNPs;
            
            if(isSNPs)
              baseQualityColour.setSelected(false);
            repaint();
          }
        });
        showMenu.add(checkBoxSNPs);
    
        
        markInsertions.addActionListener(new ActionListener()
        {
          public void actionPerformed(ActionEvent e)
          {
            repaint();
          }
        });
        showMenu.add(markInsertions);
    
    tjc's avatar
    tjc committed
        menu.add(showMenu);
    
    tjc's avatar
    tjc committed
        //
        JMenu graphMenu = new JMenu("Graph");
    
        JCheckBoxMenuItem checkBoxCoverage = new JCheckBoxMenuItem("Coverage");
        checkBoxCoverage.addActionListener(new ActionListener()
        {
          public void actionPerformed(ActionEvent e)
          {
            isCoverage = !isCoverage;
            coveragePanel.setVisible(isCoverage);
            repaint();
          }
        });
    
    tjc's avatar
    tjc committed
        graphMenu.add(checkBoxCoverage);
        
        JCheckBoxMenuItem checkBoxSNP = new JCheckBoxMenuItem("SNP");
        checkBoxSNP.addActionListener(new ActionListener()
        {
          public void actionPerformed(ActionEvent e)
          {
            isSNPplot = !isSNPplot;
            snpPanel.setVisible(isSNPplot);
            repaint();
          }
        });
        graphMenu.add(checkBoxSNP);
        menu.add(graphMenu);
        
    
        if(feature_display != null)
    
          final JCheckBoxMenuItem checkBoxSync =
            new JCheckBoxMenuItem("Asynchronous", asynchronous);
          checkBoxSync.addActionListener(new ActionListener()
    
            public void actionPerformed(ActionEvent e)
            {
              asynchronous = checkBoxSync.isSelected();
            }
          });
    
          menu.add(checkBoxSync);
    
        menu.add(new JSeparator());
    
        JMenu maxHeightMenu = new JMenu("BamView Height");
    
        menu.add(maxHeightMenu);
    
        final String hgts[] =
           {"500", "800", "1000", "1500", "2500", "5000", "50000"};
    
        
        ButtonGroup bgroup = new ButtonGroup();
        for(int i=0; i<hgts.length; i++)
        {
          final String hgt = hgts[i];
          final JCheckBoxMenuItem maxHeightMenuItem = new JCheckBoxMenuItem(hgt);
          bgroup.add(maxHeightMenuItem);
          maxHeightMenuItem.setSelected(hgts[i].equals(Integer.toString(maxHeight)));
          maxHeightMenu.add(maxHeightMenuItem);
          maxHeightMenuItem.addActionListener(new ActionListener()
          {
            public void actionPerformed(ActionEvent e)
            {
              if(maxHeightMenuItem.isSelected())
                maxHeight = Integer.parseInt(hgt);
              int start = getBaseAtStartOfView();
              setDisplay(start, nbasesInView+start, null);
            }
          });
        }
    
        menu.add(new JSeparator());
        logMenuItem.setFont(checkIsizeStackView.getFont());
        menu.add(logMenuItem);
        logMenuItem.addActionListener(new ActionListener()
        {
          public void actionPerformed(ActionEvent e)
          {
            logScale = logMenuItem.isSelected();
            repaint();
          }
        });
        
    
        JMenuItem filter = new JMenuItem("Filter by Flag...");
        menu.add(filter);
        filter.addActionListener(new ActionListener()
        {
          public void actionPerformed(ActionEvent e)
          {
            new SAMRecordFilter(BamView.this);
          } 
        });
    
    tjc's avatar
    tjc committed
        menu.add(new JSeparator());
    
      }
      
      public void setVisible(boolean visible)
      {
        super.setVisible(visible);
        mainPanel.setVisible(visible);
    
    tjc's avatar
    tjc committed
      }
      
    
    tjc's avatar
    tjc committed
      private void setViewportMidPoint()
      {
        Point p = jspView.getViewport().getLocation();
        p.y = (getHeight() - jspView.getViewport().getViewRect().height)/2;
        jspView.getViewport().setViewPosition(p);
      }
      
    
    tjc's avatar
    tjc committed
      protected int getBaseAtStartOfView()
    
    tjc's avatar
    tjc committed
      {
    
    tjc's avatar
    tjc committed
        if(feature_display != null)
          return feature_display.getForwardBaseAtLeftEdge();
        else
          return scrollBar.getValue();
    
    tjc's avatar
    tjc committed
      }
      
      /**
       * Set the panel size based on the number of bases visible
       * and repaint.
       * @param nbasesInView
       */
    
    tjc's avatar
    tjc committed
      private void setZoomLevel(final int nbasesInView)
      {
    
    tjc's avatar
    tjc committed
        int startValue = getBaseAtStartOfView();
    
    tjc's avatar
    tjc committed
        this.nbasesInView = nbasesInView;
    
    tjc's avatar
    tjc committed
        float pixPerBase = getPixPerBaseByWidth(); 
    
    tjc's avatar
    tjc committed
    
    
        if(pixPerBase*1.08f >= ALIGNMENT_PIX_PER_BASE)
    
    tjc's avatar
    tjc committed
        {
          pixPerBase = ALIGNMENT_PIX_PER_BASE;
    
          this.nbasesInView = (int)(mainPanel.getWidth()/pixPerBase);
    
    tjc's avatar
    tjc committed
          jspView.getVerticalScrollBar().setValue(0);
    
          
          if(ruler == null)
            ruler = new Ruler();
    
    tjc's avatar
    tjc committed
          jspView.setColumnHeaderView(ruler);
          showBaseAlignment = true;
    
          baseQualityColour.setEnabled(true);
    
          markInsertions.setEnabled(true);
    
    tjc's avatar
    tjc committed
        }
        else if(jspView != null)
        {
          jspView.setColumnHeaderView(null);
    
    tjc's avatar
    tjc committed
          
          if(!isStrandStackView)
            jspView.getVerticalScrollBar().setValue(
                jspView.getVerticalScrollBar().getMaximum());
    
    tjc's avatar
    tjc committed
          showBaseAlignment = false;
    
          baseQualityColour.setEnabled(false);
    
          markInsertions.setEnabled(false);
    
    tjc's avatar
    tjc committed
        }
        
    
    tjc's avatar
    tjc committed
        if(scrollBar != null)
    
    tjc's avatar
    tjc committed
          scrollBar.setValues(startValue, nbasesInView, 1, 
                 getMaxBasesInPanel(getSequenceLength()));
    
          scrollBar.setUnitIncrement(nbasesInView/20);
    
          scrollBar.setBlockIncrement(nbasesInView);
    
    tjc's avatar
    tjc committed
      }
    
    tjc's avatar
    tjc committed
    
    
    tjc's avatar
    tjc committed
      /**
    
       * Set the start and end base positions to display.
       * @param start
       * @param end
       * @param event
    
    tjc's avatar
    tjc committed
       */
    
      public void setDisplay(int start,
                             int end,
                             DisplayAdjustmentEvent event)
    
    tjc's avatar
    tjc committed
      {
    
        this.startBase = start;
        this.endBase   = end;
    
    tjc's avatar
    tjc committed
        this.nbasesInView = end-start+1;
    
        lastMousePoint = null;
    
    tjc's avatar
    tjc committed
        float pixPerBase;
        if(jspView.getViewport().getViewRect().width > 0)
    
    tjc's avatar
    tjc committed
          pixPerBase = getPixPerBaseByWidth();
    
    tjc's avatar
    tjc committed
        else
        {
          if(feature_display == null)
            pixPerBase = 1000.f/(float)(end-start+1);
          else
            pixPerBase = feature_display.getWidth()/(float)(end-start+1);
        }
    
    tjc's avatar
    tjc committed
        Dimension d = new Dimension();
    
    tjc's avatar
    tjc committed
        d.setSize(nbasesInView*pixPerBase, maxHeight);
    
    tjc's avatar
    tjc committed
        setPreferredSize(d);
    
    tjc's avatar
    tjc committed
        if(event == null)
        {
          this.startBase = -1;
          this.endBase   = -1;
        }
    
    tjc's avatar
    tjc committed
      /**
       * Return an Artemis entry from a file 
       * @param entryFileName
       * @param entryGroup
       * @return
       * @throws NoSequenceException
       */
      private Entry getEntry(final String entryFileName, final EntryGroup entryGroup) 
                       throws NoSequenceException
      {
        final Document entry_document = DocumentFactory.makeDocument(entryFileName);
        final EntryInformation artemis_entry_information =
          Options.getArtemisEntryInformation();
        
    
        //System.out.println(entryFileName);
    
    tjc's avatar
    tjc committed
        final uk.ac.sanger.artemis.io.Entry new_embl_entry =
          EntryFileDialog.getEntryFromFile(null, entry_document,
                                           artemis_entry_information,
                                           false);
    
        if(new_embl_entry == null)  // the read failed
          return null;
    
        Entry entry = null;
        try
        {
          if(entryGroup.getSequenceEntry() != null)
            bases = entryGroup.getSequenceEntry().getBases();
    
    tjc's avatar
    tjc committed
          
    
    tjc's avatar
    tjc committed
          if(bases == null)
    
    tjc's avatar
    tjc committed
          {
    
    tjc's avatar
    tjc committed
            entry = new Entry(new_embl_entry);
    
    tjc's avatar
    tjc committed
            bases = entry.getBases();
          }
    
    tjc's avatar
    tjc committed
          else
            entry = new Entry(bases,new_embl_entry);
          
          entryGroup.add(entry);
        } 
        catch(OutOfRangeException e) 
        {
          new MessageDialog(null, "read failed: one of the features in " +
              entryFileName + " has an out of range " +
                            "location: " + e.getMessage());
        }
        return entry;
      }
    
      
      private boolean isShowScale()
      {
        return showScale;
      }
    
      public void setShowScale(boolean showScale)
      {
        this.showScale = showScale;
      }
      
      public JScrollPane getJspView()
      {
        return jspView;
      }
      
      public void setBases(Bases bases)
      {
        this.bases = bases;
      }
      
      /**
       * Remove JScrollPane border
       */
      public void removeBorder()
      {
        Border empty = new EmptyBorder(0,0,0,0);
        jspView.setBorder(empty);
      }
      
      /**
       *  Handle a mouse drag event on the drawing canvas.
       **/
    
    tjc's avatar
    tjc committed
      private void handleCanvasMouseDrag(final MouseEvent event)
    
        if(event.getButton() == MouseEvent.BUTTON3 || bases == null) 
    
        if(event.getClickCount() > 1)
        {
          getSelection().clear();
          repaint();
          return;  
        }
        
    
        highlightRange(event, 
            MouseEvent.BUTTON1_DOWN_MASK & MouseEvent.BUTTON2_DOWN_MASK);
      }
      
      /**
       * 
       * @param event
       * @param onmask
       */
      private void highlightRange(MouseEvent event, int onmask)
      {
    
    tjc's avatar
    tjc committed
        int seqLength = getSequenceLength();
    
    tjc's avatar
    tjc committed
        float pixPerBase = getPixPerBaseByWidth();
    
    tjc's avatar
    tjc committed
        int start = (int) ( ( (event.getPoint().getX())/pixPerBase) + getBaseAtStartOfView() );
    
        if(start < 1)
          start = 1;
        if(start > seqLength)
          start = seqLength;
    
        if (dragStart < 0 && (event.getModifiersEx() & onmask) == onmask)
          dragStart = start;
        else if((event.getModifiersEx() & onmask) != onmask)
          dragStart = -1;
        
    
        MarkerRange drag_range;
        try
    
          if(dragStart < 0)
            drag_range = new MarkerRange (bases.getForwardStrand(), start, start);
          else
            drag_range = new MarkerRange (bases.getForwardStrand(), dragStart, start);
    
          getSelection().setMarkerRange(drag_range);
          repaint();
    
        catch (OutOfRangeException e)
    
          e.printStackTrace();
    
      /**
       * Get the colour for the given read given to it by the coverage plot.
       * @param samRecord
       * @return
       */
      private Color getColourByCoverageColour(SAMRecord samRecord)
      {
        LineAttributes lines[] = CoveragePanel.getLineAttributes(bamList.size());
        int fileIndex = 0;
        if(bamList.size()>1)
          fileIndex = (Integer) samRecord.getAttribute("FL");
        return lines[fileIndex].getLineColour(); 
      }
      
    
      private Selection getSelection()
      {
        return selection;
      }
      
    
    tjc's avatar
    tjc committed
      protected List<SAMRecord> getReadsInView()
    
    tjc's avatar
    tjc committed
        return readsInView;
      }
      
    
    tjc's avatar
    tjc committed
      protected int getBasesInView()
      {
        return nbasesInView;
      }
      
    
    tjc's avatar
    tjc committed
      private String getVersion()
      {
        final ClassLoader cl = this.getClass().getClassLoader();
        try
        {
          String line;
          InputStream in = cl.getResourceAsStream("etc/versions");
          BufferedReader reader = new BufferedReader(new InputStreamReader(in));
          while((line = reader.readLine()) != null)
          {
            if(line.startsWith("BamView"))
              return line.substring( "BamView".length() ).trim();
          }
          reader.close();
          in.close();
        }
        catch (Exception ex)
        {
        }
        return null;
      }
      
    
    tjc's avatar
    tjc committed
      /**
       * Artemis event notification
       */
    
      public void displayAdjustmentValueChanged(final DisplayAdjustmentEvent event)
    
        if(event.getType() == DisplayAdjustmentEvent.REV_COMP_EVENT &&
           event.isRevCompDisplay())
          JOptionPane.showMessageDialog(this, 
              "Flipping the display is not supported by BamView.", "Warning", 
              JOptionPane.WARNING_MESSAGE);
    
    
        if(!asynchronous)
        {
          // if not asynchronous
          displayAdjustmentWork(event);
          return;
        }
        
    
        SwingWorker worker = new SwingWorker()
    
    tjc's avatar
    tjc committed
        {
    
          public Object construct()
          {
            try
            {
              Thread.sleep(500);
            }
            catch (InterruptedException e)
            {
              e.printStackTrace();
            }
            
            if(event.getStart() != ((FeatureDisplay)event.getSource()).getForwardBaseAtLeftEdge())
            {
    
              waitingFrame.showWaiting("waiting...", mainPanel);
    
            displayAdjustmentWork(event);
    
            waitingFrame.setVisible(false);
    
            return null;
          }
        };
        worker.start();
    
      /**
       * Carry out the display agjustment event action.
       * @param event
       */
      private void displayAdjustmentWork(final DisplayAdjustmentEvent event)
      {
        if(event.getType() == DisplayAdjustmentEvent.SCALE_ADJUST_EVENT)
        {
          laststart = -1;
          lastend = -1;
          BamView.this.startBase = event.getStart();
          BamView.this.endBase   = event.getEnd();
    
          int width = feature_display.getMaxVisibleBases();
          setZoomLevel(width);
          repaint();
    
          if(coveragePanel != null && coveragePanel.isVisible()) 
          {
            coveragePanel.setStartAndEnd(BamView.this.startBase, BamView.this.endBase);
            coveragePanel.setPixPerBase(getPixPerBaseByWidth());
    
    tjc's avatar
    tjc committed
            coveragePanel.repaint();
    
    tjc's avatar
    tjc committed
          if(snpPanel != null && snpPanel.isVisible())
    
          {
            snpPanel.setStartAndEnd(BamView.this.startBase, BamView.this.endBase);
            snpPanel.setPixPerBase(getPixPerBaseByWidth());
    
    tjc's avatar
    tjc committed
            snpPanel.repaint();
    
        }
        else
        {
          setDisplay(event.getStart(), 
            event.getStart()+feature_display.getMaxVisibleBases(), event);
          repaint();
    
    tjc's avatar
    tjc committed
          if(coveragePanel != null && coveragePanel.isVisible())
            coveragePanel.repaint();
          if(snpPanel != null && snpPanel.isVisible())
            snpPanel.repaint();
    
    tjc's avatar
    tjc committed
      public void selectionChanged(SelectionChangeEvent event)
      {
        repaint();
    
    tjc's avatar
    tjc committed
      private class Ruler extends JPanel
    
    tjc's avatar
    tjc committed
      {
    
    tjc's avatar
    tjc committed
        private static final long serialVersionUID = 1L;
    
    tjc's avatar
    tjc committed
        int start;
        int end;
    
    tjc's avatar
    tjc committed
    
    
    tjc's avatar
    tjc committed
        public Ruler()
        {
          super();
    
          setPreferredSize(new Dimension(mainPanel.getWidth(), 15));
    
    tjc's avatar
    tjc committed
          setBackground(Color.white);
          setFont(getFont().deriveFont(11.f));
        }
    
    tjc's avatar
    tjc committed
    
    
    tjc's avatar
    tjc committed
        public void paintComponent(Graphics g)
        {
          super.paintComponent(g);
          Graphics2D g2 = (Graphics2D)g;
          drawBaseScale(g2, start, end, 12);
        }
    
    tjc's avatar
    tjc committed
    
        private void drawBaseScale(Graphics2D g2, int start, int end, int ypos)
        {
          int startMark = (((int)(start/10))*10)+1;
    
    
          if(end > getSequenceLength())
            end = getSequenceLength();
          
    
    tjc's avatar
    tjc committed
          for(int i=startMark; i<end; i+=10)
          {
    
    tjc's avatar
    tjc committed
            int xpos = (i-start)*ALIGNMENT_PIX_PER_BASE;
    
    tjc's avatar
    tjc committed
            g2.drawString(Integer.toString(i), xpos, ypos);
            
            xpos+=(ALIGNMENT_PIX_PER_BASE/2);
            g2.drawLine(xpos, ypos+1, xpos, ypos+5);
          }
        }
    
    tjc's avatar
    tjc committed
      }
      
    
      /**
      * Popup menu listener
      */
      class PopupListener extends MouseAdapter
      {
    
        public void mouseClicked(MouseEvent e)
        {
    
          if(e.isPopupTrigger() ||
             e.getButton() == MouseEvent.BUTTON3)
            return;
          
    
          
          if(e.getClickCount() > 1)
            getSelection().clear(); 
          else if(e.getButton() == MouseEvent.BUTTON1)
            highlightSAMRecord = mouseOverSAMRecord;
    
          else
            highlightRange(e, MouseEvent.BUTTON2_DOWN_MASK);
    
        }
        
        public void mousePressed(MouseEvent e)
        {
          maybeShowPopup(e);
        }
    
        public void mouseReleased(MouseEvent e)
        {
    
          maybeShowPopup(e);
        }
    
        private void maybeShowPopup(MouseEvent e)
        {
          if(e.isPopupTrigger())
    
    tjc's avatar
    tjc committed
          {
            if(popup == null)
            {
              popup = new JPopupMenu();
    
              createMenus(popup);
    
    tjc's avatar
    tjc committed
            }
    
            
            if(gotoMateMenuItem != null)
              popup.remove(gotoMateMenuItem);
    
    
            if(showDetails != null)
              popup.remove(showDetails);
    
            
            if( mouseOverSAMRecord != null && 
                mouseOverSAMRecord.getReadPairedFlag() &&
               !mouseOverSAMRecord.getMateUnmappedFlag() )
            {
              final SAMRecord thisSAMRecord = mouseOverSAMRecord;
              gotoMateMenuItem = new JMenuItem("Go to mate of : "+
                  thisSAMRecord.getReadName());
              gotoMateMenuItem.addActionListener(new ActionListener()
              {
    			public void actionPerformed(ActionEvent e) 
    			{
    			  String name = thisSAMRecord.getMateReferenceName();
    			  if(name.equals("="))
    			    name = thisSAMRecord.getReferenceName();
    			  int offset = getSequenceOffset(name);
    			  if(feature_display != null)
    			    feature_display.makeBaseVisible(
    			        thisSAMRecord.getMateAlignmentStart()+offset);
    			  else
    			    scrollBar.setValue(
    			        thisSAMRecord.getMateAlignmentStart()+offset-
    			        (nbasesInView/2));
    			  
    			  highlightSAMRecord = thisSAMRecord; 
    			}  
              });
              popup.add(gotoMateMenuItem);
    
            }  
              
            if( mouseOverSAMRecord != null)
            {
              final SAMRecord thisSAMRecord = mouseOverSAMRecord;
              showDetails = new JMenuItem("Show details of : "+
                  thisSAMRecord.getReadName());
              showDetails.addActionListener(new ActionListener()
              {
                public void actionPerformed(ActionEvent e) 
                {
                  FileViewer viewDetail = new FileViewer(thisSAMRecord.getReadName(), true, false);
                  appendToDetailView(thisSAMRecord, viewDetail);
                }  
              });
              popup.add(showDetails);
    
            popup.show(e.getComponent(),
                    e.getX(), e.getY());
    
    tjc's avatar
    tjc committed
          }
    
      private void appendToDetailView(SAMRecord thisSAMRecord, FileViewer viewDetail)
      {
        viewDetail.appendString("Read Name             "+thisSAMRecord.getReadName()+"\n", Level.INFO);
        viewDetail.appendString("Coordinates           "+thisSAMRecord.getAlignmentStart()+".."+
                                                         thisSAMRecord.getAlignmentEnd()+"\n", Level.DEBUG);
        viewDetail.appendString("Length                "+thisSAMRecord.getReadLength()+"\n", Level.DEBUG);
        viewDetail.appendString("Reference Name        "+thisSAMRecord.getReferenceName()+"\n", Level.DEBUG);
        viewDetail.appendString("Inferred Size         "+thisSAMRecord.getInferredInsertSize()+"\n", Level.DEBUG);
    
    tjc's avatar
    tjc committed
        viewDetail.appendString("Mapping Quality       "+thisSAMRecord.getMappingQuality()+"\n", Level.DEBUG);
    
    tjc's avatar
    tjc committed
        if(thisSAMRecord.getReadPairedFlag() && thisSAMRecord.getProperPairFlag() && !thisSAMRecord.getMateUnmappedFlag())
    
        {
          viewDetail.appendString("Mate Reference Name   "+thisSAMRecord.getMateReferenceName()+"\n", Level.DEBUG);
          viewDetail.appendString("Mate Start Coordinate "+thisSAMRecord.getMateAlignmentStart()+"\n", Level.DEBUG);
          viewDetail.appendString("Strand (read/mate)    "+
            (thisSAMRecord.getReadNegativeStrandFlag() ? "-" : "+")+" / "+
            (thisSAMRecord.getMateNegativeStrandFlag() ? "-" : "+"), Level.DEBUG);
        }
        else
        {
          viewDetail.appendString("Strand (read)         "+
              (thisSAMRecord.getReadNegativeStrandFlag() ? "-" : "+"), Level.DEBUG);
        }
        
        viewDetail.appendString("\n\nCigar String          "+thisSAMRecord.getCigarString(), Level.DEBUG);
        
        viewDetail.appendString("\n\nFlags:", Level.INFO);
        viewDetail.appendString("\nDuplicate Read    "+
            (thisSAMRecord.getDuplicateReadFlag() ? "yes" : "no"), Level.DEBUG);
    
    tjc's avatar
    tjc committed
        
        viewDetail.appendString("\nRead Paired       "+
            (thisSAMRecord.getReadPairedFlag() ? "yes" : "no"), Level.DEBUG);
        if(thisSAMRecord.getReadPairedFlag())
        {
          viewDetail.appendString("\nFirst of Pair     "+
    
            (thisSAMRecord.getFirstOfPairFlag() ? "yes" : "no"), Level.DEBUG);
    
    tjc's avatar
    tjc committed
          viewDetail.appendString("\nMate Unmapped     "+
    
            (thisSAMRecord.getMateUnmappedFlag() ? "yes" : "no"), Level.DEBUG);  
    
    tjc's avatar
    tjc committed
          viewDetail.appendString("\nProper Pair       "+
    
            (thisSAMRecord.getProperPairFlag() ? "yes" : "no"), Level.DEBUG);
    
        viewDetail.appendString("\nRead Fails Vendor\nQuality Check     "+
            (thisSAMRecord.getReadFailsVendorQualityCheckFlag() ? "yes" : "no"), Level.DEBUG);
        viewDetail.appendString("\nRead Unmapped     "+
            (thisSAMRecord.getReadUnmappedFlag() ? "yes" : "no"), Level.DEBUG);
    
    tjc's avatar
    tjc committed
        
        if(thisSAMRecord.getReadPairedFlag())
          viewDetail.appendString("\nSecond Of Pair    "+
    
            (thisSAMRecord.getSecondOfPairFlag() ? "yes" : "no"), Level.DEBUG);
    
        
        viewDetail.appendString("\n\nRead Bases:\n", Level.INFO);
        viewDetail.appendString(new String(thisSAMRecord.getReadBases()), Level.DEBUG);
    
    
      protected SAMRecordFlagPredicate getSamRecordFlagPredicate()
      {
        return samRecordFlagPredicate;
      }
    
      protected void setSamRecordFlagPredicate(
          SAMRecordFlagPredicate samRecordFlagPredicate)
      {
        laststart = -1;
        lastend = -1;
        this.samRecordFlagPredicate = samRecordFlagPredicate;
      }
    
    tjc's avatar
    tjc committed
      protected SAMRecordMapQPredicate getSamRecordMapQPredicate()
      {
        return samRecordMapQPredicate;
      }
    
      protected void setSamRecordMapQPredicate(
          SAMRecordMapQPredicate samRecordMapQPredicate)
      {
        laststart = -1;
        lastend = -1;
        this.samRecordMapQPredicate = samRecordMapQPredicate;
      }
      
    
      class PairedRead
      {
        SAMRecord sam1;
        SAMRecord sam2;
    
    tjc's avatar
    tjc committed
      public static void main(String[] args)
      {
    
    tjc's avatar
    tjc committed
        BamFrame frame = new BamFrame();
    
    tjc's avatar
    tjc committed
        
        if(args.length == 0 && BamFrame.isMac())
        {
          try
          {
            Thread.sleep(1000);
          }
          catch (InterruptedException e1) {}
    
    tjc's avatar
    tjc committed
          if(frame.getBamFile() != null)
    
    tjc's avatar
    tjc committed
            args = new String[]{ frame.getBamFile() };
        }
    
    tjc's avatar
    tjc committed
          
        List<String> bam = new Vector<String>();
    
    tjc's avatar
    tjc committed
        String reference = null;
    
        if(args.length == 0)
        {
    
    tjc's avatar
    tjc committed
          System.setProperty("default_directory", System.getProperty("user.dir"));
    
          FileSelectionDialog fileSelection = new FileSelectionDialog(
              null, true, "BamView", "BAM");
    
          bam = fileSelection.getFiles(".*\\.bam$");
    
          reference = fileSelection.getReferenceFile();
          if(reference == null || reference.equals(""))
            reference = null;
    
    tjc's avatar
    tjc committed
          
          if(bam == null || bam.size() < 1)
    
          {
            System.err.println("No files found.");
    
    tjc's avatar
    tjc committed
            System.exit(0);
    
        else if(!args[0].startsWith("-"))
        {
          for(int i=0; i< args.length; i++)
            bam.add(args[i]);
        }