New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	2837257 =	27 (1.240)	3 (0.150)	3/152	1.6	12.5%	coding (313/768 nt)	hycG	formate hydrogenlyase subunit HycG
	?	NZ_CP009273	2837317 =	18 (0.920)					coding (253/768 nt)	hycG	formate hydrogenlyase subunit HycG
Rejected: Frequency below/above cutoff threshold.

TTTATCCGTACCGCCCCACACGCAGTAGAGATCGTGGAAGATCCCGCCACTGTTACCGCAGGCACCGTAGGAGATACAAATTTTCGGGTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCG  >  NZ_CP009273/2837168‑2837343                                                                                          |                                                                                       tttACCCGTACCGCCCCACACGCAGTAGAGACCGTGGAAGATCCCGCCACTGTTACCGCAGGCACCGTAGGAGATACAAATTTTCGGGTc                                                                                        <  2:194800/90‑1 (MQ=255)    aTCCGTACCGCCCCACACGCAGTAGAGATCGTGGAAGATCCCGCCACTGTTACCGCAGGCACCGTAGGAGATACAAATTTTCGGGTCCgg                                                                                     >  1:350692/1‑90 (MQ=255)         tACCGCCCCACACGCAGTAGAGATCGTGGAAGATCCCGCCACTGTTACCGCAGGCACCGTAGGAGATACAAATTTTCGGGTCCGGCGCgg                                                                                >  1:365579/1‑90 (MQ=255)                       cAGTAGAGATCGTGGAAGATCCCGCCACTGTTACCGCAGGCACCGTAGGAGATACAAATTTTCGGGTCCGGCGCGGACTGCCacgcacgc                                                                  <  2:570562/90‑1 (MQ=255)                       cAGTAGAGATCGTGGAAGATCCCGCCACTGTTACCGCAGGCACCGTAGGAGATACAAATTTTCGGGTCCGGCGCGGACTGCCacgcacgc                                                                  >  1:570562/1‑90 (MQ=255)                         gTAGAGATCGTGGAAGATCCCGCCACTGTTACCGCAGGCACCGTAGGAGATACAAATTTTCGGGTCCGGCGCGGACTGCCACGCAcgcag                                                                >  2:386070/1‑90 (MQ=255)                            gagaTCGTGGAAGATCCCGCCACTGTTACCGCAGGCACCGTAGGAGATACAAATTTTCGGGTCCGGCGCGGACTGCCACGCAcgcagcgc                                                             >  2:172687/1‑90 (MQ=255)                              gaTCGTGGAAGATCCCGCCACTGTTACCGCAGGCACCGTAGGAGATACAAATTTTCGGGTCCGGCGCGGACTGCCACGCAcgcagcgcag                                                           >  2:556668/1‑90 (MQ=255)                                 cGTGGAAGATCCCGCCACTGTTACCGCAGGCACCGTAGGAGATACAAATTTTCGGGTCCGGCGCGGACTGCCacgcacgc                                                                  >  1:19115/1‑80 (MQ=255)                                 cGTGGAAGATCCCGCCACTGTTACCGCAGGCACCGTAGGAGATACAAATTTTCGGGTCCGGCGCGGACTGCCacgcacgc                                                                  <  2:19115/80‑1 (MQ=255)                                   tGGAAGATCCCGCCACTGTTACCGCAGGCACCGTAGGAGATACAAATTTTCGGGTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGAt                                                      >  2:259355/1‑90 (MQ=255)                                     attcggcattaaagtcgttccttcaccgcgtcatgcggatattttactgtttaCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCg                                                    <  1:94389‑M1/37‑1 (MQ=255)                                      ttcggcattaaagtcgttccttcaccgcgtcatgcggatattttactgtttaCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGc                                                   >  2:94389‑M1/53‑90 (MQ=255)                                      aaGACCCCGCCACTGTTACCGCAGGCACCGTAGGAGATACAAATTTTCGGGTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGc                                                   <  2:319736/90‑1 (MQ=255)                                        gATCCCGCCACTGTTACCGCAGGCACCGTAGGAGATACAAATTTTCGGGTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCAt                                                 <  1:556668/90‑1 (MQ=255)                                            ccGCCACTGTTACCGCAGGCACCGTAGGAGATACAAATTTTCGGGTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCa                                             >  1:515258/1‑90 (MQ=255)                                                    gTTACCGCAGGCACCGTAGGAGATACAAATTTTCGGGTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGAc                                     >  2:358822/1‑90 (MQ=255)                                                    gTTACCGCAGGCACCGTAGGAGATACAAATTTTCGGGTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGAc                                     >  1:516953/1‑90 (MQ=255)                                                       aCCGCAGGCACCGTAGGAGATACAAATTTTCGGGTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACcgc                                  <  1:386070/90‑1 (MQ=255)                                                       aCCGCAGGCACCGTAGGAGATACAAATTTTCGGGTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACcgc                                  <  2:318872/90‑1 (MQ=255)                                                       aCCGCAGGCACCGTAGGAGATACAAATTTTCGGGTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACcgc                                  <  2:515258/90‑1 (MQ=255)                                                       aCCGCAGGCACCGTAGGAGAGCCAAATTTTCGGGTCCGGCGCGGACTGCCAGGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACcgc                                  <  1:313077/90‑1 (MQ=255)                                                              gCACCGTAGGAGATACAAATTTTCGGGTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTa                           <  1:388799/90‑1 (MQ=255)                                                               cACCGTAGGAGATACAAATTTTCGGGTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTaa                          <  2:366687/90‑1 (MQ=255)                                                                      ggAGATACAAATTTTCGGGTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTaa                   <  1:274495/90‑1 (MQ=255)                                                                           tACAAATTTTCGGGTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAatat              >  2:300526/1‑90 (MQ=255)                                                                              aaaTTTTCGGGTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCg           >  1:176042/1‑90 (MQ=255)                                                                                     gtttaCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCCGTAAACAGTAAAATATCCGCATGAcg    <  1:395859‑M1/85‑1 (MQ=255)                                                                                      gggTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGAcgc   <  1:98979/90‑1 (MQ=255)                                                                                      gggTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGAcgc   <  1:259355/90‑1 (MQ=255)                                                                                       ttaccGGCGCGGTCTGCCAAGCACTCCGCGACGGGGGTCGCGTTGCACGGGTGACCGCGCCGGTAAACAGGAAAATGTCCGCATGGcgcg  >  1:48447/4‑90 (MQ=255)                                                                                          |                                                                                       TTTATCCGTACCGCCCCACACGCAGTAGAGATCGTGGAAGATCCCGCCACTGTTACCGCAGGCACCGTAGGAGATACAAATTTTCGGGTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCG  >  NZ_CP009273/2837168‑2837343

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