New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4417460	132 (1.070)	3 (0.030)	3/264	NT	2.3%	coding (509/750 nt)	yjfP	acyl CoA esterase
	?	NC_000913	= 4417479	128 (1.090)					coding (528/750 nt)	yjfP	acyl CoA esterase
Rejected: Frequency below/above cutoff threshold.

GATGACGGCACTGGGGATTACCGCTCGCCACCCCACGGTGAGATGTACCGCCAGCATGATGGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4417305‑4417460 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttgtcgcTTCCCACTCTGCCAGTGGCGCGACAATGTTATTGAATTCATTCTGCTGTGCTGCCGTTTCAGGTATCAGCGGTGGAAACAGTGAACGGGCGAGTGATGTAAAATAGCCCGATCCCATCATGCTGGCGGTACATCTCACCGTGGGGTG  <  NC_000913/4417479‑4417333                                                                                                                                                                                                                                                                                                                 GATGACGGCACTGGGGATTACCGCTCGCCACCCCACGGTGAGATGTACCGCCAGCATGATGGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACAT                                                                                                                                                           >  1:163555/1‑150  ATGACGGCACTGGGGATTAGCGCTCGCCACCCCACGGTGAGATGTACCGCCAGCATGATGGGATCGGGCTATTTTACATCTCTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATT                                                                                                                                                          <  2:991551/150‑1   TGACGGCACTGGGGATTACCGCTCGCCACCCCACGGTGAGATGTACCGCCAGCATGATGGGATCGGGCTATTTTACCTCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTG                                                                                                                                                         >  1:1465568/1‑150     cCGGCACTGGGGATTACCGCTCGCCACCCCACGGTGAGATGTACCGCCAGCATGATGGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCCCCGCTGAGAGCTGAAACGGCCGCACAGCAGAAGGAATTCAATAACATTGTC                                                                                                                                                       >  2:201569/2‑150     ACGGCACTGGGGATTACCGCTCGCCACCCCACGGTGAGATGTACCGCCAGCATGATGGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTC                                                                                                                                                       >  1:1104671/1‑150      CGGCACTGGGGATTACCGCTCGCCACCCCACGGTGAGATGTACCGCCAGCATGATGGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCG                                                                                                                                                      >  1:51202/1‑150       GGCACTGGGGATTACCGCTCGCCACCCCACGGTGAGATGTACCGCCTGCATGATGGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGC                                                                                                                                                     <  2:1728395/150‑1       GGCACTGGGGATTACCGCTCGCCACCCCACGGTGAGATGTACCGCCAGCATGATGGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGC                                                                                                                                                     <  2:284003/150‑1                                    CGGTGAGATGTACCGCCAGCATGATGGGCTCGGGCTATTTAACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCTTCCCACTCTGCCAGTGGCGCGACAATGT                                                                                                                        <  1:617415/150‑1                                              TACCGCCAGCATGATGGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCTTCCCACTCTGCCAGTGGCGCGACAATGTTATTGAATTC                                                                                                              <  1:1560457/150‑1                                                CCGCCAGCATGATGGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCTTCCCACTCTGCCAGTGGCGCGACAATGTTATTGAATTCAT                                                                                                            <  2:417629/150‑1                                                                                                                                                      TTGTCGCTTCCCACTCTGCCAGTGGCGCGACAATGTTATTGAATTCATTCTGCTGTGCTGCCGTTTCAGGTATCAGCGGTGGAAACAGTGAACGGGCGAGTGATGTAAAATAGCCCGATCCCATCATGCTGGCGGTACATCTCACCGTGG      <  2:1950049/150‑1                                                                                                                                                          ctCTTCCCACTCTGCCAGTGGCGCGACAATGTTATTGAATTCATTCTGCTGTGCTTCCGTTTCAGGTATCAGCGGTGGAAACAGTGAACGGGCGAGGGATGTAAAATAGCCCGATCCCATCATGCTGGCGGTACATCTCACCGTGGGGTG  <  2:72378/148‑1                                                                                                                                                          CGCTTCCCACTCTGCCAGTCGCGCGACAATGTTATTGAATTCATTCTGCTGTGCTGCCGTTTCAGGTATCAGCGGTGGAAACAGTGAACGGGCGAGTGATGTAAAATAGCCCTATCCCATCATGCTGGCGGTACATCTCACCGTGGGGTG  <  1:973887/150‑1                                                                                                                                                                                                                                                                                                                 GATGACGGCACTGGGGATTACCGCTCGCCACCCCACGGTGAGATGTACCGCCAGCATGATGGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4417305‑4417460 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttgtcgcTTCCCACTCTGCCAGTGGCGCGACAATGTTATTGAATTCATTCTGCTGTGCTGCCGTTTCAGGTATCAGCGGTGGAAACAGTGAACGGGCGAGTGATGTAAAATAGCCCGATCCCATCATGCTGGCGGTACATCTCACCGTGGGGTG  <  NC_000913/4417479‑4417333

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