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

GTTTTAACCAGCGGGTTAACTCTTTTTGACGAGATTTATTCATTGCACGCTTAGCAGGTGAGTTATCAGAATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTTGCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAA  >  NC_000913/930920‑931176                                                                                                                   |                                                                                                                                               gTTTTAACCAGCGGGTTAACTCTTTTTGACGAGATTTATTCATTGCACGCTTAGCAGGTGAGTTATCAGAATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGc                                                                                                                                      <  1:392751/125‑1 (MQ=255) gTTTTAACCAGCGGGTTAACTCTTTTTGACGAGATTTATTCATTGCACGCTTAGCAGGTGAGTTATCAGAATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGc                                                                                                                                      >  2:392751/1‑125 (MQ=255)                                                  cTTAGCAGGTGAGTTATCAGAATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACaccatgcgactatgggtcgc                                                                                                                            >  4:193264‑M1/1‑66 (MQ=255)                                                  cTTAGCAGGTGAGTTATCAGAATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACaccatgcgactatgggtcgc                                                                                                                            <  3:193264‑M1/86‑21 (MQ=255)                                                                   cAGAATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCa                                                                                                                                     <  1:163098/60‑1 (MQ=255)                                                                   cAGAATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCa                                                                                                                                     >  2:163098/1‑60 (MQ=255)                                                                      aaTTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTg                                                                                                        <  2:49627/86‑1 (MQ=255)                                                                      aaTTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTg                                                                                                        >  1:49627/1‑86 (MQ=255)                                                                            ttGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACaccatgcgactatgg                                                                                                                                 >  2:1042404‑M1/1‑40 (MQ=255)                                                                            ttGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACaccatgcgactatgg                                                                                                                                 <  1:1042404‑M1/55‑16 (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)                                                                                                             aGGCGACaccatgcgactatgggtcgcctttattttttccccgttgtaacattgctc                                                                                              >  2:953931‑M1/1‑7 (MQ=255)                                                                                                             aGGCGACaccatgcgactatgggtcgcctttattttttccccgttgtaacattgctc                                                                                              <  1:953931‑M1/57‑51 (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)                                                                                                                   |                                                                                                                                               GTTTTAACCAGCGGGTTAACTCTTTTTGACGAGATTTATTCATTGCACGCTTAGCAGGTGAGTTATCAGAATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTTGCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAA  >  NC_000913/930920‑931176

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