Reanalyze all our datasets including the trimming of the restriction sites as a preprocessing step

[leilaicruz]

No description provided.

[leilaicruz]

• dbem3_a
  Settings for the trimming :

trim to the right of the first restriction site, and discard reads that after trimming are below 50 bp to ensure decent confidence when aligning them.
trimming_settings="ktrim=r k=4 mink=2 hdist=0 qin=33 minlen=50 tossbrokenreads" #-t

• reads on genes with the trimmed data (with the restriction sites)=2640422
• reads on genes with the NON trimmed data (with the restriction sites)=5041605
• insertions on genes with the trimmed data (with the restriction sites)=215235
• insertions on genes with the NON trimmed data (with the restriction sites)=342014
• the number of reads in the ADE2 for the trimmed one is around 7% of the total amount of reads with 187507 and for the non trimmed one it represents 20% of the data with 1003795.