New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 931034	9 (0.100)	4 (0.050)	4/184	NT	31.6%	intergenic (‑72/+51)	cydD/trxB	glutathione/cysteine ABC transporter export permease/ATPase/thioredoxin reductase, FAD/NAD(P)‑binding
	?	NC_000913	= 931049	9 (0.100)					intergenic (‑87/+36)	cydD/trxB	glutathione/cysteine ABC transporter export permease/ATPase/thioredoxin reductase, FAD/NAD(P)‑binding

AATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTTGCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAA  >  NC_000913/930989‑931176                                                             |                                                                                                                                aaTTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTg                                                                                                        >  1:49627/1‑86 (MQ=255) aaTTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTg                                                                                                        <  2:49627/86‑1 (MQ=255)            gagCAATGTTACAACGGGGCAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACGTTTGTTACTg                                                                                                        <  4:100442/75‑1 (MQ=255)            gagCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTgtcgcct                                                                                                                          <  1:953931‑M2/57‑8 (MQ=255)            gagCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTgtcgcct                                                                                                                          >  2:953931‑M2/1‑50 (MQ=255)            gagCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTg                                                                                                        >  3:100442/1‑75 (MQ=255)                                         ggCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTTGCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATaaa  >  2:775771/1‑148 (MQ=255)                                          gCGACCCATAGTCGCATGGTgtcgcctttattttttccccgttgtaacattgctctgcaaataattctgataactcacctgctaag                                                               <  3:193264‑M2/86‑67 (MQ=255)                                          gCGACCCATAGTCGCATGGTgtcgcctttattttttccccgttgtaacattgctctgcaaataattctgataactcacctgctaag                                                               >  4:193264‑M2/1‑20 (MQ=255)                                               ccATAGTCGCATGGTgtcgcctttattttttccccgttgtaacattgctctgcaa                                                                                         <  1:1042404‑M2/55‑41 (MQ=255)                                               ccATAGTCGCATGGTgtcgcctttattttttccccgttgtaacattgctctgcaa                                                                                         >  2:1042404‑M2/1‑15 (MQ=255)                                               ccATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTa                                                                                                           >  3:277438/1‑37 (MQ=255)                                               ccATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTa                                                                                                           <  4:277438/37‑1 (MQ=255)                                                     tCGCATGGTGTCGCCTTCTCTACTTTTGTTACTGATTTGTAAAATTATTTTGCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGtgc                                              <  1:590974/92‑1 (MQ=255)                                                     tCGCATGGTGTCGCCTTCTCTACTTTTGTTACTGATTTGTAAAATTATTTTGCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGtgc                                              >  2:590974/1‑92 (MQ=255)                                                             |                                                                                                                                AATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTTGCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAA  >  NC_000913/930989‑931176

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ 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.