--- title: "Alissa Doucet- Amsonia tharpii Project" breaks: false tags: RADSeq --- # Buxbaumia ## Login ### Away from CBG: Login to Forticlient- VPN to access Fabronmia and Buxbaumia (See email) ### At CBG: Need to be at CBG to download files onto my computer: Open Windows Powershell Login to fabromia (Jeremie's login info), then buxbaumia (Dylan's login info) ### Get to where data is: **Change directory three times:** cd .. cd .. cd .. **List what's in folder:** ls **To go to specific folder:** (cd folder/) cd data cd Alissa (my data folder) Stacks === ### Overview **1. Radtags:** Clean-up and demultiplex (do this with each sublibrary idividually) **2. Denovomap:** Combine all sublibraries into one folder; assemble loci, call SNPs, based on parameter choices; output vcf files and basic pop stats; can refine parameters with python script **3. Clean and use this output for analyses** ### 1. Process Radtags 1. nano barsub1.txt 2. save - control x 3. Remove files 4. rm barcode.txt 5. Remove folder 6. rm -rf test3 7. Make folder 8. mkdir (make directory) 9. Create new folder 10. mkdir sublib1 11. Note: 12. Folders (blue) 13. gunzipped files (red) 14.text files, scripts (white) 15.nano process_radtag (copy code) ### Run Script sh process_radtags.sh enter ### Check to see if it's working cat process_radtags1.out cp = copy/paste cp process_radtags.sh (duplicate file under new name) cp process_radtags.sh (file name I want to copy) process_radtags2.sh (new file name) top = see who is using Buxbaumia can process 2-3 radtags at a time q = get out ### How to know whena job is done If run at noon, check around 4-5pm. If done, will shoe in output file when its odne running; will summarize info from run at the bottom ### Notes When copy/pasting barcodes into file, must use up arrow rather than scroll to see all 48 samples 01/26/2023 Processed Radtags for sublibaries #1-#10 ### DeNovo Assembly denovo_map.[d]1-h (denovo help file) Navigate to my folder cd sublib1 [process_radtags.log] (can check reads here) 1. Remove files ending in rem.1.fq.gz & rem.2.gq.gz **zless** file name (view in zip file) or raw sequence data **less**