New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	3642840 =	39 (1.260)	3 (0.110)	3/218	2.5	8.4%	coding (546/1290 nt)	arsB	arsenite/antimonite:H(+) antiporter ArsB
	?	NZ_CP009273	3644036 =	30 (1.080)					intergenic (+14/‑615)	arsC/BW25113_RS18155	glutaredoxin‑dependent arsenate reductase/hypothetical protein
Rejected: Frequency below/above cutoff threshold.

CTCGCCTGTGCCGCCCTTCAGAAGTGGTGCTGGAAATTCTGCCAGATGCGCAAAAAGGCGCATTCTCCAAGGAAGATGGCGAGAAAGTGGTTGATGAAGCGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACAAAGTTATGAAAATTCGATGAATAGTATGTTATTTGTAGGCCTGACAGGCGTAGCGCATCAGG  >  NZ_CP009273/3643904‑3644170                                                                                                                                     |                                                                                                                                       cTCGCCTGTGCCGCCCTTCAGAAGTGGTGCTGGAAATTCTGCCAGATGCGCAAAAAGGCGCATTCTCCAAGGAAGATGGCGAGAAAGTGGTTGATGAAGCGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCcac                                                                                                                                  <  2:118632/139‑1 (MQ=255)     ccTGTGCCGCCCTTCAGAAGTGGTGCTGGAAATTCTGCCAGATGCGCAAAAAGGCGCATTCTCCAAGGAAGATGGCGAGAAAGTGGTTGATGAAGCGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGgt                                                                                                                              <  1:581959/139‑1 (MQ=255)                           gtgCTGGAAATTCTGCCAGATGCGCAAAAAGGCGCATTCTCCAAGGAAGATGGCGAGAAAGTGGTTGATGAAGCGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCtttatt                                                                                                        >  1:414954/1‑139 (MQ=255)                              cTGGAAATTCTGCCAGATGCGCAAAAAGGCGCATTCTCCAAGGAAGATGGCGAGAAAGTGGTTGATGAAGCGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCtttattta                                                                                                      <  1:581216/138‑1 (MQ=255)                                      tCTGCCAGATGCGCAAAAAGGCGCATTCTCCAAGGAAGATGGCGAGAAAGTGGTTGATGAAGCGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATg                                                                                             >  2:212968/1‑139 (MQ=255)                                                 cgcAAAAAGGCGCATTCTCCAAGGAAGATGGCGAGAAAGTGGTTGATGAAGCGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGgc                                                                                   <  1:20105/138‑1 (MQ=255)                                                           cgcATTCTCCAAGGAAGATGGCGAGAAAGTGGTTGATGAAGCGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGt                                                                               >  2:59441/1‑132 (MQ=255)                                                           cgcATTCTCCAAGGAAGATGGCGAGAAAGTGGTTGATGAAGCGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGt                                                                               <  1:59441/132‑1 (MQ=255)                                                                  ctgcggaatatctttgcgaaaatagagatgtaacatcaccagcgtggcaacaatcgcggcgatatccCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTa                                                                  >  2:212601‑M2/68‑138 (MQ=255)                                                                            agggCGAGAAAGTGGTTGATGAAGCGGGTAAGCGCCTGAAATAAAGCGGCAATATCCCCCCCACAGGTTGTTAGAACAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACAAGGTTATg                                                        >  2:586940/3‑138 (MQ=255)                                                                            aTGGCGAGAAAGTGGTTGATGAAGCGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACg                                                                                        <  2:431341/106‑1 (MQ=255)                                                                            aTGGCGAGAAAGTGGTTGATGAAGCGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACAAAGTTATg                                                        >  2:119454/1‑138 (MQ=255)                                                                            aTGGCGAGAAAGTGGTTGATGAAGCGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACAAAGTTATg                                                        >  2:528210/1‑138 (MQ=255)                                                                            aTGGCGAGAAAGTGGTTGATGAAGCGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACAAAGTTATg                                                        >  2:479166/1‑138 (MQ=255)                                                                            aTGGCGAGAAAGTGGTTGATGAAGCGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACAAAGTTATg                                                        >  2:425794/1‑138 (MQ=255)                                                                            aTGGCGAGAAAGTGGTTGATGAAGCGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACAAAGTTATg                                                        >  2:338836/1‑138 (MQ=255)                                                                            aTGGCGAGAAAGTGGTTGATGAAGCGGGGAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACg                                                                                        >  1:431341/1‑106 (MQ=255)                                                                               gCGAGAAAGTGGTTGATGAAGCGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACAAAGTTATGaaa                                                     <  1:212968/138‑1 (MQ=255)                                                                                cGAGAAAGTGGTTGATGAAGCGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACAAAGTTATGAAAAt                                                   <  1:87788/139‑1 (MQ=255)                                                                                   gaAAGCGGTTGCTGAAGCGGGCAAGCGCCTGAAATAACGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCGCCCTCCAAAGTTATGAAAATTCg                                                <  2:129741/139‑1 (MQ=39)                                                                                        tagagatgtaacatcaccagcgtggcaacaatcgcggcgatatccCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACAAAGTTATGAAAATTCGATGaa                                           >  2:402790‑M2/46‑139 (MQ=255)                                                                                        tGGTTGATGAAGCGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACAAAGTTATGAAAATTCGATGaa                                           <  2:591519/139‑1 (MQ=255)                                                                                        tGGTTGATGAAGCGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACAAAGTTATGAAAATTCGATGaa                                           <  1:449041/139‑1 (MQ=255)                                                                                        tGGTTGATGAAGCGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACAAAGTTATGAAAATTCGATGaa                                           <  1:338836/139‑1 (MQ=255)                                                                                         ggTTGATGAAGCGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGGTTATAGAAAAGTGCGCTTTATTTATGGCGTATGCGACGCTGGCGCGTCTTATCCGACCTACAAAGTTATGAAAATTCGATGaa                                           >  1:129741/1‑138 (MQ=255)                                                                                           ttGATGAAGCGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACAAAGTTATGAAAATTCGATGAATAg                                        <  2:61580/139‑1 (MQ=255)                                                                                           ttGATGAAGCGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACAAAGTTATGAAAATTCGATGAATAg                                        <  1:363732/139‑1 (MQ=255)                                                                                             gatgaAGCGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACAAAGTTATGAAAATTCGATGAATAGt                                       <  1:425794/138‑1 (MQ=255)                                                                                             gatgaAGCGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACAAAGTTATGAAAATTCGATGAATAGt                                       <  1:479166/138‑1 (MQ=255)                                                                                                  acatcaccagcgtggcaacaatcgcggcgatatccCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACAAAGTTATGAAAATTCGATGAATAGTATGtt                                  <  1:212601‑M2/103‑1 (MQ=255)                                                                                                      ggTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACAAAGTTATGAAAATTCGATGAATAGTATGTTATTTg                             >  2:47628/1‑139 (MQ=255)                                                                                                        tAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACAAAGTTATGAAAATTCGATGAATAGTATGTTATTTGTa                           <  2:456341/139‑1 (MQ=255)                                                                                                               cTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACAAAGTTATGAAAATTCGATGAATAGTATGTTATTTGTAGGCCTGa                    <  2:462665/139‑1 (MQ=255)                                                                                                                                  tccCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACAAAGTTATGAAAATTCGATGAATAGTATGTTATTTGTAGGCCTGACAGGCGTAGCGCATCAgg  <  1:586940‑M2/135‑1 (MQ=255)                                                                                                                                     |                                                                                                                                       CTCGCCTGTGCCGCCCTTCAGAAGTGGTGCTGGAAATTCTGCCAGATGCGCAAAAAGGCGCATTCTCCAAGGAAGATGGCGAGAAAGTGGTTGATGAAGCGGGTAAGCGCCTGAAATAAAGCGGCGATATCCCCCCCACAGGTTGTTAGAAAAGTGCGCTTTATTTATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACAAAGTTATGAAAATTCGATGAATAGTATGTTATTTGTAGGCCTGACAGGCGTAGCGCATCAGG  >  NZ_CP009273/3643904‑3644170

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.