Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC JC	REL606	2,405,131	IS1 (–) +9 bp	100%	coding (168‑176/1347 nt)	fadL →	long‑chain fatty acid outer membrane transporter

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 2263507	NA (NA)	35 (1.310)	26/278	NT	100%	noncoding (1/768 nt)	IS1	repeat region
	?	REL606	2405131 =	0 (0.000)					coding (168/1347 nt)	fadL	long‑chain fatty acid outer membrane transporter
*	?	REL606	= 2405139	0 (0.000)	32 (1.200)	23/278	NT	100%	coding (176/1347 nt)	fadL	long‑chain fatty acid outer membrane transporter
	?	REL606	2407442 =	NA (NA)					noncoding (768/768 nt)	IS1	repeat region

CTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGTGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTCCGATCATCAGGTTGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGC  >  REL606/2407294‑2407580                                                                                                                                                     |                                                                                                                                           caaaaacaTCCAAAAAGAGATTAGTATTAATGGTGATAGTCTCACCCAATATTGTAATGAAAACGCCAGGAGAGCTTCATCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTCCGATCATCAGGTTGTAgg                                                                                                                                           <  1:52284/145‑1 (MQ=255) cTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTATCACCCAATATTGTAGTGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATAATGGCAGAATCATTTGCAACTGGTTCCGATCATCAGGTTGTggg                                                                                                                                           <  1:353391/150‑1 (MQ=255)  tACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGTGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTCCGATCATCAGGTTGTAGGt                                                                                                                                          <  1:373381/150‑1 (MQ=255)               gagGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGTGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTCCGATCATCAGGTTGTAGGTAATGACTCCAACt                                                                                                                             >  2:182500/1‑150 (MQ=255)                      tCGTATTAATTGTGATAGTATCACCCAATATTGTAGTGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTCCGATCATCAGGTTGTAGGTAATGACTCCAACTTATTGAt                                                                                                                      <  1:65985/150‑1 (MQ=255)                      tCGTATTAATGGTGATAGTCTCACCCAATATTGTAGTGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTCCGATCATCAGGTTGTAGGTAATGACTCCAACTTATTGAt                                                                                                                      <  1:35011/150‑1 (MQ=255)                               tGGTGATAGTCTCACCCAATATTGTAGTGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTCCGATCATCAGGTTGTAGGTAATGACTCCAACTTATTGATAGTGTTTTa                                                                                                             <  2:578316/150‑1 (MQ=255)                               tGGTGATAGTCTCACCCAATATTGTAGTGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTCCGATCATCAGGTTGTAGGTAATGACTCCAACTTATTGATAGTGTTTTa                                                                                                             >  2:6629/1‑150 (MQ=255)                                       gTCTCACCCAATATTGTAGTGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCACGGCAGAATCATTTTCAACTGGTTCCGATCATCAGGTTGTAGGTACTGACTCCAACTTATTGATAGTGTTTTCAGTCCAGa                                                                                                     >  1:673367/1‑150 (MQ=255)                                                 atatTGTAGTGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATGATTTGCAACTGGTTCCGATCATCAGGTTGTAGGTAATGACGCCAAATTGTTGATAGTGTTTTATGTTCAGATGATGCCCgc                                                                                           >  1:680071/1‑149 (MQ=39)                                                      gTAGTGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTCCGATCATCAGGTTGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACt                                                                                      >  2:474248/1‑150 (MQ=255)                                                                cGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTCCGATCATCAGGTTGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCa                                                                            >  2:71299/1‑150 (MQ=255)                                                                cGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTCCGATCATCAGGTTGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCa                                                                            >  2:327282/1‑150 (MQ=255)                                                                          gCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTCCGATCATCAGGTTGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGa                                                                  >  2:7925/1‑150 (MQ=255)                                                                                                                     cATTTGCAACTGGTTCCGATCATCAGGTTGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGtgc                       <  1:213282/150‑1 (MQ=255)                                                                                                                                ggTTCCGATCATCAGGTTGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGAtt            >  1:552620/1‑150 (MQ=255)                                                                                                                                          atcaGGTTATATATAAAGACTACAACTTCTTGATTGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGACTTTGAGAACGACAGCGACTTCCGTCACAGCCGTGCCAGGTGCTGTCTCAGATTTAGGTTATGc  <  1:360989/150‑1 (MQ=11)                                                                                                                                                     |                                                                                                                                           CTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGTGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTCCGATCATCAGGTTGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGC  >  REL606/2407294‑2407580

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