New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	= 729706	NA (NA)	12 (0.110)	10/154	3.5	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.

CCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAG  >  NZ_CP009273/729658‑729795                                                 |                                                                                          ccGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACg                                                  <  2:1898884/90‑1 (MQ=31)  cGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACgg                                                 >  1:123501/1‑90 (MQ=32)  cGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACgg                                                 >  2:2799434/1‑90 (MQ=32)  cGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACgg                                                 >  2:1944394/1‑90 (MQ=32)      cAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGaa                                             <  1:2779538/90‑1 (MQ=35)       aGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGaaa                                            >  2:1271086/1‑90 (MQ=33)       aGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGaaa                                            >  1:1913237/1‑90 (MQ=33)          aCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTTAGGGAGaaaaac                                         >  2:2444836/1‑88 (MQ=255)              tCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGAcc                                     >  1:2407065/1‑90 (MQ=33)              tCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGAcc                                     <  2:1695663/90‑1 (MQ=35)              tCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGAcc                                     <  1:164586/90‑1 (MQ=35)                aGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTc                                   <  1:1626294/90‑1 (MQ=35)                  cATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCAt                                 >  2:531917/1‑90 (MQ=33)                  cATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCAt                                 >  1:2957603/1‑90 (MQ=33)                  cATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCAt                                 >  1:3291646/1‑90 (MQ=33)                   aTGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATc                                <  2:3112950/90‑1 (MQ=35)                    tgtgGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGTAGGCTGACGGAGAAAACCGACCTCATcc                               <  2:65815/90‑1 (MQ=25)                    tgtgGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATcc                               <  2:2829458/90‑1 (MQ=35)                     gtgGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATccc                              >  2:3049523/1‑90 (MQ=33)                     gtgGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATccc                              >  2:3110434/1‑90 (MQ=33)                     gtgGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATcc                               >  2:688032/1‑89 (MQ=33)                     gtgGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATcc                               <  1:688032/89‑1 (MQ=35)                        gCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGa                           <  1:2649147/90‑1 (MQ=35)                         ccGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGaa                          <  2:1364704/90‑1 (MQ=35)                         ccGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGaa                          <  2:391178/90‑1 (MQ=35)                         ccGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGaa                          >  1:2684829/1‑90 (MQ=33)                          cGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAg                         <  1:2937041/90‑1 (MQ=35)                           ggATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAG‑AAACCGACCTCATCCCGGAAggg                       >  1:743520/1‑90 (MQ=16)                           ggATAACCGTATCGACCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAgg                        >  1:511965/1‑90 (MQ=18)                               aaCCGTATCGCCCGTGACGGGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGtt                    <  1:1942814/90‑1 (MQ=25)                                   gTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATcc                <  1:3316686/90‑1 (MQ=33)                                   gTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATcc                <  1:2300812/90‑1 (MQ=33)                                     aTCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGc              >  2:940561/1‑90 (MQ=33)                                      tCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGGTGACGGAGAAAACCCACCGCTAACCGGTAGGGTTTATTCGCa             >  2:1578830/1‑90 (MQ=11)                                       cGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCAc            <  2:2913645/90‑1 (MQ=33)                                       cGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCAc            <  2:3124899/90‑1 (MQ=33)                                       cGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCAc            <  2:3193559/90‑1 (MQ=33)                                       cGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCAc            <  1:1815014/90‑1 (MQ=33)                                       cGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCAc            <  1:1205962/90‑1 (MQ=33)                                        gCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACg           <  2:2968755/90‑1 (MQ=33)                                        gCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACg           <  1:1201941/90‑1 (MQ=33)                                          ccGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGGGGAAACCGACCTCATCCCGGAAGGGGTTATCCACACggg         >  2:193398/1‑89 (MQ=11)                                          ccGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGga         >  1:2322676/1‑90 (MQ=33)                                           cGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGgat        >  1:3002712/1‑90 (MQ=33)                                                 gcgcACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGGGAAAACCGACCTCATCCCGGAAGGGGTTATCCCCCCGGATgaggag  >  2:451786/1‑90 (MQ=11)                                                 gcgcACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAg  >  2:1607857/1‑90 (MQ=33)                                                 |                                                                                          CCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAG  >  NZ_CP009273/729658‑729795

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