The overall quality of the sequencing reads was assessed using FastQC (v0.10.1) [45]. Raw reads were processed to remove adapter sequence using FASTA/Q Clipper from the FASTX_Toolkit [46] and sequences of poor quality using Sickle [47]. Reference indexes of the Ae. aegypti genome assembly AaegL2 [48] (obtained from VectorBase) and the C. capitata genome assembly Ccap_1.0 [49] (obtained from NCBI) were constructed using Bowtie2 [50]. Trimmed reads were aligned to these indexes using TopHat2 (v2.0.9) [51]. Transcript assemblies were created from the alignments using the reference annotation based transcript assembly method [52] with Cufflinks (v2.1.1) [53] followed by Cuffmerge (v1.0.0) [54]. Transcript expression in each sample was quantified using Cuffdiff2 (v2.1.1) [55].
Do you have any questions about this protocol?
Post your question to gather feedback from the community. We will also invite the authors of this article to respond.