New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4200145 =	28 (0.910)	4 (0.160)	4/110	2.1	18.9%	intergenic (+51/‑136)	yjaG/hupA	DUF416 domain protein/HU, DNA‑binding transcriptional regulator, alpha subunit
	?	NC_000913	4200182 =	11 (0.430)					intergenic (+88/‑99)	yjaG/hupA	DUF416 domain protein/HU, DNA‑binding transcriptional regulator, alpha subunit
Rejected: Frequency below/above cutoff threshold.

GATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGAT  >  NC_000913/4200120‑4200244                                                               |                                                               catacacgcaccgatagccactttttttcgccccccaaatgtagaccagacggggggcgaaaAAAAGt                                                           <  2:792517‑M2/6‑1 (MQ=255)       cgcaccgatagccacttttTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTAt                                                     <  2:913129‑M1/49‑1 (MQ=255)        cGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATc                                                    >  1:1094645/1‑68 (MQ=255)          accgatagccacttttTTTCGCAACCCGTCAGGTCTACAGTTGGGGGGCGAAAAAAAGTGGCTATCgg                                                  <  2:49498‑M1/52‑1 (MQ=255)           ccgatagccactttttttcgccccccaactgtagacctgacgggttgcgaaaAAAAGTGGCTATCGGt                                                 <  2:49498‑M2/16‑1 (MQ=255)             cTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGc                                               >  2:15055/1‑68 (MQ=255)              cggAGTCTTCCCTTTCGCCCCCCGTCTTGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGcgcg                                              <  2:670808/67‑1 (MQ=255)              atagccactttttttcgccccccaaatgtagaccagacggggggcgaaaAAAAGTGGCTATCGGTGCg                                              <  2:913129‑M2/19‑1 (MQ=255)                 aGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCgtgt                                           >  2:721509/1‑68 (MQ=255)                 aGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGCGCgcgt                                           <  2:30093/68‑1 (MQ=255)                   tCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTAt                                         <  1:1062266/68‑1 (MQ=255)                   tCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTAt                                         <  2:717562/68‑1 (MQ=255)                   tCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTAt                                         >  1:1090300/1‑68 (MQ=255)                   tCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGCGCGTGTAt                                         <  2:1099555/68‑1 (MQ=255)                   tCTTCCCTTTCGCCCCCCGGCTGTTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTAt                                         <  2:294464/68‑1 (MQ=255)                    acttttTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATg                                        <  2:792517‑M1/62‑1 (MQ=255)                          TTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGgag                                  >  1:1059639‑M1/1‑68 (MQ=255)                           TTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAggagg                                 <  1:159577‑M1/68‑2 (MQ=255)                             CGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGgagagc                               >  2:1094640‑M1/1‑67 (MQ=255)                             CGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGt                               >  2:468329‑M1/1‑68 (MQ=255)                             CGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGt                               >  1:367176‑M1/1‑68 (MQ=255)                                CCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCt                            <  1:242615‑M1/68‑1 (MQ=255)                                  CCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCTGTTCGTGTATGCAGGAGAGTGCTAt                          >  1:802258‑M1/1‑68 (MQ=255)                                           agaccagacggggggcgaaaAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCAtt                 >  1:1128879‑M2/21‑68 (MQ=255)                                                aTTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGt            >  2:1145754/1‑68 (MQ=255)                                                     gggggcgaaaAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCAc       >  2:900710‑M2/11‑68 (MQ=255)                                                     gggggcgaaaAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCAc       >  1:315605‑M2/11‑68 (MQ=255)                                                          cgaaaAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGAt  <  1:15055‑M2/63‑1 (MQ=255)                                                          cgaaaAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGAt  <  2:499941‑M2/63‑1 (MQ=255)                                                          cgaaaAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGAt  <  2:315605‑M2/63‑1 (MQ=255)                                                               |                                                               GATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGAT  >  NC_000913/4200120‑4200244

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