Advertisement
Guest User

Untitled

a guest
Jun 18th, 2019
80
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.74 KB | None | 0 0
  1. # creating index :
  2. bowtie2-build reference_file_name.fasta prefix_of_output_index_file
  3.  
  4. # mapping reads :
  5. # -x path to index, only prefix needed
  6. # -1 first file with paired end reads
  7. # -2 second file with paired end reads
  8. # -U files with unpaired reads
  9. # -S output file name
  10.  
  11. ./bowtie2 -x prefix_of_output_index_file \
  12. -1 file_1P.fastq \
  13. -2 file_2P.fastq \
  14. -U file_1U.fastq \
  15. -U file_2U.fastq \
  16. -S mapping.sam
  17.  
  18. # NOTE:
  19. # if number of reads in paired-end files is uneven, bowtiew will throw and exception
  20. # you can get uneven number of files, when you additionally recrop one of them in SE mode of Trimmomatic
  21. # in this case use this approach :
  22.  
  23. bowtie2 -x prefix_of_output_index_file \
  24. -U file_1P.fastq \
  25. -U file_2P.fastq \
  26. -U file_1U.fastq \
  27. -S mapping.sam
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement