New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 168811	26 (0.830)	3 (0.100)	3/480	3.4	10.3%	coding (1328/2244 nt)	fhuA	ferrichrome outer membrane transporter/phage receptor
	?	NC_000913	= 258675	NA (NA)					noncoding (1/768 nt)	IS1	repeat region
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCAAACTCCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTAC  >  NC_000913/168568‑169034                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                gtgCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGAtttggtgat                                                                                                                                                                                                                           <  2:138803/250‑7 (MQ=255)           aGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCaa                                                                                                                                                                                                                  >  1:235372/1‑249 (MQ=255)                         aaGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCAAACTcc                                                                                                                                                                                                  >  2:109800/1‑251 (MQ=255)                         aaGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTACTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCAAACTcc                                                                                                                                                                                                  <  2:64973/251‑1 (MQ=255)                              cATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCAAACTCCGGccc                                                                                                                                                                                             >  1:296891/1‑251 (MQ=255)                                  aTCTGGCTCTTAAATACGTTGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATGCCCAGTTGCAGAGCAATGTGGACACTGGCGAGAGCGGCCACACCCTGCTGACCGGTGTCGTCTTTATGCTTATGCGTAATGACATCATCGCCTGGTTTGGTTACGACGTCTCTGTGCCACTGCTCAATCTGTACAACCCGTTGAATACCGATTTCGACTTCAAGGCCAAAGATCCGGCAAACTCCGGCCCTTa                                                                                                                                                                                          <  2:204764/250‑1 (MQ=255)                                        cACGTAAATACGTCGTTGATGATGAGAAGCTGCACAACTTCTCCGTTTATACCCAGTGGCATAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCAAACTCCGGCCCTTACCGCAtt                                                                                                                                                                                   <  1:340368/251‑1 (MQ=255)                                           gTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATTACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAGTATGTACAATCCGGTGAATACCGATTTCGACTTCAATACCAAAGATCCGGCAACCTCCGGCCCTTACCGCAttatg                                                                                                                                                                                >  2:92534/1‑248 (MQ=255)                                                     cgtTGATGATGAGAAGCTTCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACGGGCGATATCGACCACACCCTGCTGACCGGTTGCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCAAACTCCGTCCCTTACCGCATTCTGAATAAACAg                                                                                                                                                                       <  2:177150/250‑1 (MQ=255)                                                            gatgaGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCAAACTCCGGCCCTTACCGCATTCTGAATAAACAGaaacaaac                                                                                                                                                               <  1:147468/251‑1 (MQ=255)                                                               gagaAGCTGCAAAACTTCTCCGTTGATACACAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCGGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTTCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCAAACTCCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACggg                                                                                                                                                            <  1:113353/251‑1 (MQ=255)                                                                               tctcCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGccaa                                                                                                                                                                                                      >  2:9908/1‑190 (MQ=255)                                                                               tctcCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCaaa                                                                                                                                                                                                      <  1:9908/193‑1 (MQ=255)                                                                                  cctgtGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCAGCCCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCAAACTCCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGTATc                                                                                                                                           <  2:116314/245‑1 (MQ=255)                                                                                                    agagCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCAAACTCCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTg                                                                                                                       >  1:78834/1‑251 (MQ=255)                                                                                                                           atCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCAAACTCCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTa                                                                                                <  1:109800/251‑1 (MQ=255)                                                                                                                             cGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCAAACTCCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATg                                                                                              <  1:257256/251‑1 (MQ=255)                                                                                                                                       ctgctgACCGGGGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGTCAAACTCCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCCGGCGCAGTGGGATAAAGTGCTGGTCCCCCACTGCGGTCGTTATGACTGGGc                                                                                       >  2:271995/1‑248 (MQ=255)                                                                                                                                         gctgCCCGGTGTCGTCTTTATGAGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACCACTCTGTGCCACTGCTCAATCTGTACAAGCCGGTGAATACCGATTTCGACTTCAATGCCAAATATCCGGCAAACTCCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGAt                                                                                   <  1:62825/250‑1 (MQ=255)                                                                                                                                                            atgcgtatgcgtAATGACATCAACGCCTTGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCAAACTCCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACcgcg                                                                <  1:26102/250‑1 (MQ=255)                                                                                                                                                                      gtaATGACATCAACGCCTGGTTTGGTTACGACGACTCTGGGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCAAACTCCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGgacg                                                     <  2:193065/251‑1 (MQ=255)                                                                                                                                                                             cATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCAAACTCCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATa                                              <  2:78834/251‑1 (MQ=255)                                                                                                                                                                                                 acgacgACTTGGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCAAACTCCGGCCCGTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGGTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCTGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGttt                          <  1:336303/251‑1 (MQ=255)                                                                                                                                                                                                    acgacTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCAAACTCCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTAcc                       >  2:135808/1‑251 (MQ=255)                                                                                                                                                                                                           ggtgCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCAAACTCCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCgtg                <  2:296891/250‑1 (MQ=255)                                                                                                                                                                                                              gCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTGCGACTTCAATGCCAAAGATCCGGCAAACTCCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCgtggtg             <  1:350801/251‑1 (MQ=255)                                                                                                                                                                                                                             gTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCAAACTCCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTAc  <  1:135808/247‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                GTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCAAACTCCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGTGATGACAAACAGTTTACCTGGCGTGGTGGTGTTAACTAC  >  NC_000913/168568‑169034

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