Current work on simplifying the inventory tracking page.

Its still a work in progress and I need to implement the compontents to
browse inventory items by category.

In addition, somehow I need to re-integrate the component
that prints labels for archived runfolders on disk.

Add utility to fix the off-by-3 error in the amplified from field.

The bug appeared because the library name and library primary key
drifted apart because of some accidental deletions.

I eventually moved us back to using raw primary keys as that
made constructing django urls simplier & easier. However
I'd apparently forgotten to adjust the amplified from field.

Make the inventory login page look a bit nicer.
This involved adjusting some of the base django templates.
Whenever given a choice I try to make our pages look more like the
admin site.

Adjust some whitespace

This mildly adjusts the css to be a bit prettier

Add a new ini-file option to point to the flowcell repository.
It's in the [frontend] section called "results_dir".

Add config option notification_bcc for always attaching to notify emails.
Its called notification_bcc in the htsworkflow.ini file and
NOTIFICATION_BCC in the django settings.py file.

It ignores the "send bcc" option on the email_started form.

Cache the attributes for each type of fastq file.
Previously it was just using whatever was last as the
set of attributes. Which led to everything being flagged as
FastqRd2

Map junctions.bed to the Junctions view.
Hopefully it'll pass validation and not be treated as a bed file,
even though it ends in the .bed extension.

(There were lots of problems with a previous subission because the
cufflinks bed file was almost but not quite a match to ucsc's bed
file validator. Or maybe they caved on the definition of a score.)

Autodetect location of *2fastq scripts
Uses the .__file__ variable of the python modules.

Don't count failed flowcells when guessing if a library is paired or single ended
In addition the NameToView map was extended to cache the is_paired flag
so I don't have to keep searching through the dictionary of library information.

Merge branch 'master' of mus.cacr.caltech.edu:htsworkflow

Very basic non-ExtJS version of the inventory page
It needs to group by object type. Possible with an index of types
followed by index of objects before getting to the detail page.

The component where it forces you to login before viewing the
inventory page is broken. It shows a super-plain form, but
fails on submit with some kind of cross site forgery prevention
cookie fails

Remove another ExtJS panel reference.
This version has a working two column independantly scrollable
library browse with no ExtJS dependency.

Update to the new UCSC DAF file.
We are submitting 3 different types of expression data,
one for Gencode v3c, one for Gencode v4, and one for the
de novo cufflinks assembly.

As a result I needed to update the file name to ucsc view map
to describe these new, more specific datasets.

Also I modified the name of the archival and upload condor scripts
to include a username.

Create a class to handle mapping extension to ucsc view attributes.

This will attempt to guess if a library is paired end by seeing
if there are more paired end lanes than single read lanes.

The file exention matching was changed to use fnmatch instead of
ends with.

I still haven't added the ability to define extensions to ucsc DAF
view maps in a config file, as I don't know how to handle the case
of the .bam file that goes to a different view depending on if its
a paired end vs single end.

Also the ucsc_gather script is too long and parts of it need
to migrate into the rest of the htsworkflow tree.

Two lines between functions (trivial change)

Only complane about missing genomes if we _wanted_ to align the lane.
If it was flagged for sequencing, don't complain.

Check for s_${lane}_02_matrix.txt as well as s_${lane}_1_matrix.txt
It still could use more test coverage

Remove dependency on ExtJS for library & lanes_for pages.
The inventory page still needs extjs.

In the process I cleaned up some of the CSS as well.

The downside is you'll need to link
htsworkflow/static/js/jquery.min.js
to a valid jquery sourcefile.

The WITH_SEQUENCE flag needs to be lower-case "true".

Report version number derived from git tag.
This patch includes the necessary infrastructure to support that
feature and its been added to qseq2fastq and srf2fastq.

Additionally to improve testability of qseq2fastq and srf2fastq, the
original standalone module was moved into htsworkflow.pipelines and a
small stub module was placed in scripts.

Set WITH_SEQUENCE as both a per-lane AND global parameter
the per-lane version set $(lane)_WITH_SEQUENCE := TRUE in the makefiles
the gerald WITH_SEQUENCE TRUE changed the setting in the config.txt file

So lets just set both versions and hopefully we'll get sequence files
out the other end.

Let me show you my hammer...

WITH_SEQUENCE is actually a global GERALD option, remove the lane specifier.

Update gerald config file generation.
Include ELAND_SET_SIZE, WITH_SEQUENCE options
Stop commenting SEQUENCE_FORMAT out

Minor changes to display

Removed some models from the admin index view

Changed Frequently Used links to be properly relative

Removed from admin index page

Removed unused DataRun model from admin index page

Removed all analysis models from main admin index

Limited which models show up on admin index page
Removed some fields from view in the Library Change/Add page
Rearranged the remaining Library fields

Ignore some files generated by python
.coverage
*.egg-info

Adjust spoolwatchers use of inotify to correspond to pyinotify 0.8.9

Handle the case when a sequencing lane lacks any yield information.
(For instance if we only run one lane of analysis on a flowcell)

Remind myself that dry-run doesn't work with the extract results code

Add support for CASAVA 1.7

They stopped generating eland_extended files so I needed to parse
the export files instead.

Also more carefully test how I'm computing the U0-2 and R0-2
scores.

Better detect which program ids generated CNF4 srf files.
apparently illumina2srf v1.11.5Illumina.1.3 also genrated CNF4s.
It's only illumina2srf v1.11.6... that generated CNF1s.

Split parts of build_fastqs into a seperate function.
Use the database to flag if a lane is bad instead of hard coding it.
Support comments and spaces in the library_id to target directory file.

Add a "All Lanes" option to the control-lane field for the cases
where we didn't set a control lane, but instead let the pipeline
estimate from the entire flowcell.

Shrink flowcell notes field.
See [ticket:167]

Don't offer to track flowcell config files in svn

Remove debug print statement

Reduce the size of the bioanalyzer summary text box
to 3 lines.

Tell git to ignore *.py[co]~ files.
(Useful if you're using git-svn)

This patch adds in fields to support storing bioanalyzer results. [ticket:166]

Georgi used slightly different extensions for the second batch
of data

Add search field to the lanes admin page

Use the htsworkflow API to determine if a flowcell is paired end or not.

Under some conditions urlerror doesn't have a code,
so just report the raw error message in that case

Report if a lane was on a paired end flowcell from the library
json structure.

Since srf2named_fastq detects what type of srf file I can remove
the -c option that was forcing it to read CNF1 formatted srf files

Update srf2named_fastq to try to detect if the srf file is CNF1 or CNF4
and figure out the correct option to pass to srf2fastq.

Add option to force overwriting old fastqs.

It will ignore the current existence of a fastq file when generating
the condor submit script.

In addition it will tell srf2named_fastq to --force as well.

If a quality score started with an @ sign it was treated as a header
which created an invalid fastq file.

This patch fixes that, and introduces some test cases for srf2named_fastq.py

This still isn't ideal as the RESULTS_HOME_DIR is still out of my
home directory, but at least the use of expanduser means it'll work on
both linux and os x.

Save changes needed to submit to UCSC in Jun 2010.

Some of those changes include:
  * modifying the list of variables to include in the ddf
  * making it easier to set the MapAlgorithm.
  * Return information about the condor scripts so I can make a
    condor dagman script to run all the compression jobs

Perhaps some parts of this should be moved into the main
htsworkflow. I can see wanting the code to batch convert
srf/qseqs to fastq being more generally useful.

Fix srf2named_fastq to output the proper /2 in paired end reads

Update the 'flowcell started' email message to include the sequencer
in the url.

Make it easier to run srf2named_fastq.py

Wrapper script that helps convert srf files to fastq files.
It can do the following:
  split the fastq into two files (for paired end reads)
  add in a flowcell id to the header (for either type of read).

Rename avg_lib_size to gel_cut_size, and add insert_size
to clear up the confusion about what was supposed to be
being recorded in avg_lib_size.

In addition, this patch adds in a per lane status field.

Update internal copy of the django admin/templates/admin/index.html
They changed how the css was being imported in the django css
so when running with django 1.1 the admin index page was messed up.

This gets all the borders to show up correctly.

Collect fastqs by read and add them to the configuration ini file as a
single line. (As desired by UCSC).

Also the library to result map file supports a basic comment character.
if # is the _first_ character it will skip that line.

Next I should fix the avg library size / insert length variables.

Put partial support back in for srf files.

Since I don't know if the srf file is supposed to be single or
paired end, this version assumes paired end unless you provide the
--single option.

Currently it'll give up if you try to convert a paired srf file
to a fastq file.

Also I made the code formatting in the make_parser function
look cleaner, and changed it to allow setting the logging verbosity
via command line options --verbose/--debug.

Add script to try and build submissions to the UCSC encode project.

This version supports generating qseq2fastq entries using the
htsworkflow api and scanning the flowcell repository directory.

There was code to generate the ddf files (from ini files).

I need to update the ini generation code to scan the submisison directory
for fastq files and group them by read.

Extend htsworkflow.pipelines.sequences to also try to figure out the cycle count.
In addition there is experimental code to shove the found sequences into a
sql database.

I also needed to bug fix the sequence patterns to catch the fake flowcell
ilmn200901 which wasn't matching my regexp for detecting flowcell ids.

Move the code to scan the sequence file archive to its own module so
I can use it in scripts other than make-library-tree

new api module actually needed logging.

Move the knowledge of the urls for the REST API to one new file
htsworkflow.util.api and then update some of the scripts that were
using the api to import from the new module.

Yes this increases the dependencies, but it does mean it'll be
easier to update the urls if we need to change them

apparently commas are important

make-library-tree should be an installed script too

Add in extra fields lorian asked for to library detail page

Do not insert anything into the header if there is no flowcell info.

(Previously there was a spurious _)

qseq2fastq should also be installed as a script

A patch that allows printing more than 11 labels at a time.

Default to pointing to the Ubuntu python location for django admin templates.

Matches can have trailing AGCT in addition to a number

Always return a count from carefully_make_hardlinks
Be more flexible about which json parser to use

Disabling 'delete selected'.

Added a 'Print Labels' action to the Library Admin Page.
 * Django 1.1 feature.
 * FIXME: Requires a Printer Template (Inventory) with type Library to work...
   * returns a useful error if template does not exist.

Update summary script to read from the GERALD Summary.xml file
instead of depending on randomly changign html code.

Add support for generating fasta files in addition to fastq files
Add an option to add a flowcell ID to the header

WARNING: Django 1.0.2 to Django 1.1.1 compatibility patch... There's not going back now!

Extend qseq2fastq to write to two fastq files,
one for files that pass filter and one that doesn.

Report hidden field in the library API

Use the HTS workflow API to figure out the library tree.
This also needed to search flow flowcell id by the starting name
because we still have the status of a flowcell being part of the
name in a few places.

there is no such thing as sequence_extended. I was using the wrong
suffix generator for paired end sequencing

Return affiliation, library name, and comment in the lanes_for json
api call

Actually implement the code to loop over a list of runfolders
on the command line.

Update the inventory tracker code for the split from lanes being
in the flowcell table to their own stand-alone model.

Also I made the mark_archived_data script take a list of
runfolder archives so I can archive a whole hard disk in one go.

Adds a json api 'lanes_for' feature

Force auth_backend error messages to sys stderr, as
wsgi hates stdout

Remove debugging code that breaks mod_wsgi

Added 'lanes_for' which will show recent flowcell lanes ordered by date,
and allows filtering by username.

In addition I modified the library index to bin runs into
small (<40), medium (<100), and large (>=100) runs seperated by single
and paired end reads.

Refine user handling.
The sysadmins need username to match up with the unix accounts,
The site manager needs a meaningful name to attach users to samples.
So the HTSUser string representation is first/last name and then username
in the corner.
In addition I modified the add user popup form to allow setting the
first/last name during the user creation.

Update test code to deal with the switch to storing archive
qseq files instead of srf files

Don't throw an error if library.cell_line is None.

The API was having problems where if the cell_line wasn't set it was
trying to do None.cellline_name, which didn't work so well.

In addition there were a few other type conversion issues, such as
unicode(None) != None.

So I added unicode_or_none

Modify qseq2fastq to also read from compressed tar files containing qseq files

Update the usage string for qseq2fastq