New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	3022788 =	10 (0.460)	3 (0.150)	3/156	1.6	23.8%	coding (1503/1920 nt)	ygfT	formate‑dependent uric acid utilization protein YgfT
	?	NZ_CP009273	3022819 =	10 (0.500)					coding (1472/1920 nt)	ygfT	formate‑dependent uric acid utilization protein YgfT

CATCTTCGTCACAAGCGATATATTGCGGTTGAACATTGAACTGAAACTCGACACCTTCCTCGCGCGCATTAACCACCTCTTTGCGCGAGCCCGGCATACTGACTTCATCACGACGATACGCGCAGGTCACGCTGGCGGCATTGAGGCGGATGGAAGTCCGCAAACAATCCATTGTCGT  >  NZ_CP009273/3022730‑3022907                                                                                          |                                                                                         cATCTTCGTCACAAGCGATATATTGCGGTTGAACATTGAACTGAAACTCGACACCTTCCTCGCGCGCATTAACCACCTCTTTGCGCGAGc                                                                                          >  2:17143/1‑90 (MQ=255)               gCGATATATTGCGGTTGAACATTGAACTGAAACTCGACACCTTCCTCGCGCGCATTAACCACCTCTTTGCGCGAGCCCGGCATACTGACt                                                                            <  2:73345/90‑1 (MQ=255)               gCGATATATTGCGGTTGAACATTGAACTGAAACTCGACACCTTCCTCGCGCGCATTAACCACCTCTTTGCGCGAGCCCGGCATACTGACt                                                                            <  2:194166/90‑1 (MQ=255)                          cGGTTGAACATTGAACTGAAACTCGACACCTTCCTCGCGCGCATTAACCACCTCTTTGCGCGAGCCCGGCATACTGACTTCATCacgacg                                                                 >  1:210683/1‑90 (MQ=255)                           gcgtatcgtcgtgatgaagtcagtatgccgggctcgcgcaaagaggtggttaatgcgcgcgagCCCGGCATACTGACTTCATCAcgacga                                                                <  1:241563‑M2/27‑1 (MQ=255)                             gtatcgtcgtgatgaagtcagtatgccgggctctcgcaaagaggtggttaatgcgcgcgagCCCGGCATACTGACTTCATCACGACGATa                                                              <  2:332508‑M2/29‑1 (MQ=255)                              tatcgtcgtgatgaagtcagtatgccgggCTCGCGCGCATTAACCACCTCTTTGCGAGAGCCCGGCATACTGACTTCATCACGACGATAc                                                             <  2:332508‑M1/61‑1 (MQ=255)                              tGAACATTGAACTGAAACTCGACACCTTCCTCGCGCGCATTAACCACCTCTTTGCGCGAGCCCGGCATACTGACTTCATCACGACGATAc                                                             >  2:435839/1‑90 (MQ=255)                                tcgtcgtgatgaagtcagtatgccgggCTCGCGCGCATTAACCACCTCTTTGCGCGAGCCCGGCATACTGACTTCATCACGACGATAcgc                                                           <  1:241563‑M1/63‑1 (MQ=255)                                  cATTGACCTGAAACTCGACACCTTCCTCGCGCGCATTAACCTCCTCTTTGCGCGAGCCCGGCATACTGACTTCATCACGACGATACgcgc                                                         <  1:600485/90‑1 (MQ=255)                                           gAAACTCGACACCTTCCTCGCGCGCATTAACCACCTCTTTGCGCGAGCCCGGCATACTGACTTCATCACGACGATACGCGCAGGTCACGc                                                >  2:407648/1‑90 (MQ=255)                                                       cTTCCTCGCGCGCATTAACCACCTCTTTGCGCGAGCCCGGCATACTGACTTCATCACGACGATACGCGCAGGTCACGCTGGCGGCATTGa                                    <  1:573659/90‑1 (MQ=255)                                                           ctcgcgcaaagaggtggttaatgcgcgcgggCCCGGCATACTGACTTCATCACGACGATACGCGCAGGTCACGCTGGCGGCATTGAGGCg                                <  1:375413‑M2/59‑1 (MQ=255)                                                                  gcATTAACCGCCTCTTTGCGCGAGCCCGGCATACTGACTTCATCACGACGATACGCGCAGGTCACGCTGGCGGCATTGAGGCGGATGGaa                         <  1:17143/90‑1 (MQ=255)                                                                             ctctTTGCGCGAGCCCGGCATACTGACTTCATCACGACGATACGCGCAGGTCACGCTGGCGGCATTGAGGCGGATGGAAGTCCGCAAACa              >  1:555165/1‑90 (MQ=255)                                                                                      cgagCCCGGCATACTGACTCCATCACGACGATACGCGCAGGTCACGCTGGCGGCATTGAGGCGGATGGAAGTCCGCAAACAATCCATTGt     >  1:231613‑M2/5‑90 (MQ=255)                                                                                         gCCCGGCATACTGACTTCATCACGACGATACGCGCAGGTCACGCTGGCGGCATTGAGGCGGATGGAAGTCCGCAAACAATCCATTGTCGt  >  1:296887‑M2/2‑90 (MQ=255)                                                                                          |                                                                                         CATCTTCGTCACAAGCGATATATTGCGGTTGAACATTGAACTGAAACTCGACACCTTCCTCGCGCGCATTAACCACCTCTTTGCGCGAGCCCGGCATACTGACTTCATCACGACGATACGCGCAGGTCACGCTGGCGGCATTGAGGCGGATGGAAGTCCGCAAACAATCCATTGTCGT  >  NZ_CP009273/3022730‑3022907

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