MorrellLAB / sequence_handling

A series of scripts to automate sequence workflows
19 stars 8 forks source link

Quality_Assessment outputs incorrect read depth when read length is variable #33

Closed Aerin13 closed 7 years ago

Aerin13 commented 7 years ago

If the sequence reads have a variable length according to FastQC (ex: ReadLength = 1-100), then the computation of the read depth fails. The PROJECT_quality_summary.txt output from Quality_Assessment has the number of reads listed again under the Read Depth column instead of a depth. This should be fixed to give depth estimates using the longest read length.