New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	503443 =	40 (1.220)	4 (0.130)	4/256	3.8	11.6%	intergenic (‑205/+33)	ybaL/fsr	inner membrane putative NAD(P)‑binding transporter/putative fosmidomycin efflux system protein
	?	NC_000913	503470 =	24 (0.790)					intergenic (‑232/+6)	ybaL/fsr	inner membrane putative NAD(P)‑binding transporter/putative fosmidomycin efflux system protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TGTTCATGAATAATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATG  >  NC_000913/503327‑503614                                                                                                                                                |                                                                                                                                                 tgTTCATGAATAATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGaaa                                                                                                                                             >  3:76300/1‑149 (MQ=255)   ttCATGAATAATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAAt                                                                                                                                            <  4:48037/148‑1 (MQ=255)     cATGAATAATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCa                                                                                                                                          <  5:103266/148‑1 (MQ=255)       tGAATAATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTc                                                                                                                                       <  3:56297/149‑1 (MQ=255)        gAATAATTATTGAATTTGGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGAGCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTATAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCt                                                                                                                                      <  8:122165/149‑1 (MQ=255)         aataatTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGTAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCtt                                                                                                                                     >  4:112711/1‑149 (MQ=255)              ttattGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGc                                                                                                                                >  7:28579/1‑149 (MQ=255)                 ttGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGTTCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCGTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCgg                                                                                                                             >  2:108215/1‑149 (MQ=255)                  tGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGcctga                                                                                                                            >  3:100391/1‑146 (MQ=255)                      ttttGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATc                                                                                                                        <  4:76300/149‑1 (MQ=255)                              aaaaaTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATcaggcagg                                                                                                                 >  7:65799/1‑148 (MQ=255)                                 aaTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGACATCAGTCTTTATGCCTGATAAccggcataaaga                                                                                                             <  4:100391/149‑7 (MQ=255)                                        ttAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCaa                                                                                                      <  8:65799/149‑1 (MQ=255)                                               cAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATccc                                                                                                <  1:75293/148‑1 (MQ=255)                                                         aCGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCa                                                                                     <  8:96822/149‑1 (MQ=255)                                                             aTCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGa                                                                                  >  7:67319/1‑148 (MQ=255)                                                                       aGAAGGCGGTTTAACAGAGGAAAGGCTTATGGGGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGGTATCAGGCAGGAATATGGTCAACGTCCCCAATAGTGGCAGGAAAGCGCGGAtg                                                                       >  2:114619/1‑148 (MQ=255)                                                                           ggCGGGTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTtat                                                                   <  6:42964/149‑1 (MQ=255)                                                                                       gagGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGa                                                       >  2:107330/1‑149 (MQ=255)                                                                                         ggAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATg                                                     <  1:107330/149‑1 (MQ=255)                                                                                          gAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTTTGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGc                                                    <  1:72605/149‑1 (MQ=255)                                                                                            aaGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTg                                                  >  7:99649/1‑149 (MQ=255)                                                                                                             tAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAg                                 <  3:71882/149‑1 (MQ=255)                                                                                                                              aGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTc                <  8:28579/149‑1 (MQ=255)                                                                                                                                            tGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATg  >  4:76355/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                 TGTTCATGAATAATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATG  >  NC_000913/503327‑503614

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