Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	2,737,118	IS1 (+) +8 bp	intergenic (+235/‑29)	bamD → / → raiA	outer membrane protein assembly factor BamD/stationary phase translation inhibitor and ribosome stability factor

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1978503 =	NA (NA)	22 (0.730)	20/160	0.5	100%	noncoding (768/768 nt)	IS1	repeat region
	?	NC_000913	2737118 =	0 (0.000)					intergenic (+235/‑36)	bamD/raiA	outer membrane protein assembly factor BamD/stationary phase translation inhibitor and ribosome stability factor
*	?	NC_000913	= 1979270	NA (NA)	32 (1.060)	27/160	0.2	100%	noncoding (1/768 nt)	IS1	repeat region
	?	NC_000913	= 2737125	0 (0.000)					intergenic (+242/‑29)	bamD/raiA	outer membrane protein assembly factor BamD/stationary phase translation inhibitor and ribosome stability factor

GAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAA  >  NC_000913/1979202‑1979349                                                                     |                                                                                gAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGtat                                                                      <  2:326887/80‑1 (MQ=255)     cAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGt                                                                  >  2:326910/1‑80 (MQ=255)        aaGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATa                                                              <  2:326888/81‑1 (MQ=255)             aCTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACtttt                                                         <  1:326889/81‑1 (MQ=255)             aCTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACtttt                                                         >  1:326917/1‑81 (MQ=255)                     aCCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTggg                                                 >  1:326911/1‑81 (MQ=255)                     aCCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTggg                                                 >  1:326914/1‑81 (MQ=255)                          aaaaaCACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGgtgt                                              <  1:326891/79‑1 (MQ=255)                          aaaaaCACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGgtgt                                              <  2:326890/79‑1 (MQ=255)                           aaaaCACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTAtt                                           >  2:326909/1‑81 (MQ=255)                           aaaaCACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTAtt                                           <  1:326892/81‑1 (MQ=255)                            aaaCACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTAttt                                          <  1:326893/81‑1 (MQ=255)                            aaaCACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTAttt                                          <  2:326894/81‑1 (MQ=255)                               cacCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCaa                                       <  2:326895/81‑1 (MQ=255)                                      atACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTa                                <  2:326896/81‑1 (MQ=255)                                                    tAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgct                  <  1:326897/81‑1 (MQ=255)                                                      aGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctac                <  1:326898/81‑1 (MQ=255)                                                       gTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacg               <  1:326899/81‑1 (MQ=255)                                                          ggCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgc              >  1:326913/1‑79 (MQ=255)                                                             agcaTCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagc          >  2:326916/1‑80 (MQ=255)                                                                    cccATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaa  <  2:326900/81‑1 (MQ=255)                                                                    cccATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaa   >  2:326912/1‑80 (MQ=255)                                                                     |                                                                                GAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAA  >  NC_000913/1979202‑1979349

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