New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	= 322269	38 (2.480)	3 (0.210)	3/238	1.3	8.9%	coding (690/1671 nt)	betA	choline dehydrogenase
	?	NZ_CP009273	= 322352	26 (1.840)					coding (607/1671 nt)	betA	choline dehydrogenase
Rejected: Frequency below/above cutoff threshold.

CGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAA  >  NZ_CP009273/322219‑322487                                                                                                                                      |                                                                                                                                        cgGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCacgcg                                                                                                                                    >  2:85192‑M2/1‑134 (MQ=255) cgGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCacgcg                                                                                                                                    >  2:138837‑M2/1‑134 (MQ=255) cgGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGAGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCacgcg                                                                                                                                    >  2:84589‑M2/1‑134 (MQ=255)  gggTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCacgcgc                                                                                                                                   >  1:145134‑M2/1‑133 (MQ=255)              tGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTCCGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGgcgc                                                                                                                        >  2:70920/1‑138 (MQ=255)              tGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCAAGCGCGGTGCTGGcgcg                                                                                                                       >  2:133816/1‑139 (MQ=255)                ctttCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGa                                                                                                                      <  2:180259/135‑1 (MQ=255)                cTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGa                                                                                                                      <  1:133816/138‑1 (MQ=255)                cTGGCGCCTTCCAGCCAGTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGa                                                                                                                      <  1:85793/138‑1 (MQ=255)                     gCCTTCCAGCCATTCGACGCCCACCGGGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGccc                                                                                                                <  1:247782/139‑1 (MQ=255)                        ttCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCacgcg                                                                                                                                    <  1:141233‑M2/116‑6 (MQ=255)                        ttCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCacgcg                                                                                                                                    >  2:141233‑M2/1‑111 (MQ=255)                              ccAGTCGACGCCCGCCGCGCGTTTGACGTCAAAGATGATGGGATCGGTCATAGCGCGAGTGCGGAGGGTCAGGTTAGGACGCGATGTGGCCTGATCGAGATAGCCACGCGCGGTGCCGGCGCGACGGCCCTGCGGCGTg                                                                                                       <  2:167344/139‑1 (MQ=255)                               cATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGa                                                                                                      <  1:58188/139‑1 (MQ=255)                                 ttCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACg                                                                                                    >  1:193100/1‑139 (MQ=255)                                               cgcgTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCg                                                                                       >  1:183566/1‑138 (MQ=255)                                                       ccGTCAAAAATGATGTGATCGGTCAGAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCacgcgcggtgggcgtcgaatggctggaaggcgacagcaccatcccaacccgcgcaacg                                                                               >  1:258001‑M2/1‑80 (MQ=255)                                                              aaaTGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCacgcgcggtgggcgtcgaatggctggaaggcgacagca                                                                                                   <  1:152029‑M2/111‑39 (MQ=255)                                                              aaaTGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCacgcgcggtgggcgtcgaatggctggaaggcgacagca                                                                                                   >  2:152029‑M2/1‑73 (MQ=255)                                                                   atgTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGAt                                                                  <  2:183566/139‑1 (MQ=255)                                                                        aTCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCa                                                                                                                                        <  2:240162‑M2/64‑2 (MQ=255)                                                                        aTCGGGCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCa                                                                                                                                        >  1:240162‑M2/1‑63 (MQ=255)                                                                                tAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGtt                                                           >  2:143630/1‑133 (MQ=255)                                                                                   cGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCacgcg                                                                                                                                    >  1:66043‑M2/1‑52 (MQ=255)                                                                                   cGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCacgcg                                                                                                                                    <  2:66043‑M2/57‑6 (MQ=255)                                                                                       aGTGCGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTg                                              <  2:11391/139‑1 (MQ=255)                                                                                       aGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATGGTCCATg                                              >  2:183567/1‑139 (MQ=255)                                                                                       aGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGAGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTg                                              >  2:185130/1‑139 (MQ=255)                                                                                             aaTGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGAt                                                                  >  1:133956/1‑113 (MQ=255)                                                                                             aaTGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGAt                                                                  <  2:133956/113‑1 (MQ=255)                                                                                             aaTGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCggg                                        >  1:196417/1‑139 (MQ=255)                                                                                                    aGGTTAGGACGCGATTTGGCCGGATCGAGATGGCCACGCGCGGTGCTGGCGCGACCGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGccc                                  <  2:232263/138‑1 (MQ=255)                                                                                                       ttAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCt                              <  2:194458/139‑1 (MQ=255)                                                                                                           gACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCAc                          >  1:192228/1‑139 (MQ=255)                                                                                                            aCGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGCGATCGTCCGTGCGCGGGTAGCCCGCCTGCACg                         >  2:61709/1‑139 (MQ=255)                                                                                                                gATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCg                     >  1:78797/1‑139 (MQ=255)                                                                                                                 aTTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAACCCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGc                    >  1:229882/1‑139 (MQ=255)                                                                                                                           aTCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAatcatc          >  1:25889/1‑139 (MQ=255)                                                                                                                                     ccACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCaaa  >  2:199323/1‑137 (MQ=255)                                                                                                                                      |                                                                                                                                        CGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTTTGCCGTCAAAAATGATGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAA  >  NZ_CP009273/322219‑322487

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