New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	= 729706	NA (NA)	5 (0.100)	5/150	3.2	NA	pseudogene (103/1506 nt)	BW25113_RS03635	RHS repeat protein
	?	NZ_CP009273	= 729729	NA (NA)					pseudogene (126/1506 nt)	BW25113_RS03635	RHS repeat protein
Rejected: Coverage evenness skew score above cutoff.

GCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATG  >  NZ_CP009273/729665‑729790                                          |                                                                                     gCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGaaaa                                      >  1:216751/1‑90 (MQ=33)       tCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGAcc                                <  1:1576483/90‑1 (MQ=35)       tCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGAcc                                <  2:58554/90‑1 (MQ=33)            aTGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATc                           <  1:1343171/90‑1 (MQ=35)              gtgGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATccc                         >  2:1385956/1‑90 (MQ=33)              gtgGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATccc                         >  2:313480/1‑90 (MQ=33)              gtgGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATccc                         >  1:117546/1‑90 (MQ=33)               tgGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCg                        <  1:1230001/90‑1 (MQ=35)                ggCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCgg                       >  2:831011/1‑90 (MQ=33)                ggCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCgg                       >  2:1104382/1‑90 (MQ=33)                  ccGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGaa                     <  1:1104382/90‑1 (MQ=35)                            gTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATcc           <  1:1389342/90‑1 (MQ=33)                              aTCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGa                      >  1:1288587/1‑77 (MQ=34)                              aTCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGa                      <  2:1288587/77‑1 (MQ=35)                              aTCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGc         <  1:1048218/90‑1 (MQ=33)                               tCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGAcc                                <  2:405464/66‑1 (MQ=35)                               tCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGAcc                                >  1:405464/1‑66 (MQ=35)                                cGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCAc       <  1:981081/90‑1 (MQ=33)                                cGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCAc       <  1:1503964/90‑1 (MQ=33)                                cGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCAc       <  2:768260/90‑1 (MQ=33)                                cGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCAc       <  2:95361/90‑1 (MQ=33)                                   ccGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGga    >  1:669584/1‑90 (MQ=33)                                     gTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCACCCCGGAAGGGGTTATCCGCACGgatg  >  1:1360792/1‑90 (MQ=18)                                          |                                                                                     GCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATG  >  NZ_CP009273/729665‑729790

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