New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	3209524 =	19 (0.650)	3 (0.120)	3/150	2.0	12.8%	coding (327/765 nt)	nfeF	NADPH‑dependent ferric chelate reductase
	?	NZ_CP009273	3209606 =	24 (0.930)					coding (245/765 nt)	nfeF	NADPH‑dependent ferric chelate reductase
Rejected: Frequency below/above cutoff threshold.

AATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGCGCCTGCATCGCCCAGCCGCTGGCGACCCCACCGTCGTGAATAAAGAAATCAATCGCCAGTTCATGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGT  >  NZ_CP009273/3209436‑3209611                                                                                         |                                                                                        aaTCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGCGCCTGCATCGCCCAGCCGCTGGCGAccc                                                                                        <  1:109043/90‑1 (MQ=255)    cTTCCGGCCCCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGCGCCTGCATCGCCCAGCCGCTGGCGACCCCAc                                                                                     <  1:647172/90‑1 (MQ=255)    cTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGCGCCTGCATCGCCCAGCCGCTGGCGACCCCAc                                                                                     <  2:790563/90‑1 (MQ=255)       ccGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGCGCCTGCATCGCCCAGCCGCTGGCGACCCCACcgt                                                                                  >  2:752863/1‑90 (MQ=255)         ctcttctttcctcaacctgacgctcactttgtgccgccaacggtaacggaagagggcatcgtctggccggaaggaccacgCCCACcgtcg                                                                                <  2:307154‑M1/10‑1 (MQ=255)         ctcttctttcctcaacctgacgctcactttgtgccgccaacggtaacggaagagggcatcgtctggccggaaggaccacgCCCACcgtcg                                                                                >  2:264315‑M1/81‑90 (MQ=255)                    cGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGCGCCTGCATCGCCCAGCCGCTGGCGACCCCACCGTCGTGAATAAAGaa                                                                     <  1:206572/90‑1 (MQ=255)                          gcgcgGACCTGCCACCGTAAGTTTATCGCCAGGTTGCCCCTGCATCGCCCCGCCGCCGGCGACCCCACCGTCGTGAATAAAGAaatcaat                                                               >  1:637895/1‑90 (MQ=255)                              ggACCTGCCACCGTAAGTTTATCGCCCGGTTGCGCCTGCATCGCCCAGCCGCTGGCGACCCCACCGTCGTGAATAAAGAAATCAATCGcc                                                           <  1:700132/90‑1 (MQ=255)                              ggACCTGCCACCGTAAGTTTATCGCCCGGTTGCGCCTGCATCGCCCAGCCGCTGGCGACCCCACCGTCGTGAATAAAGAAATCAATCGcc                                                           <  1:752863/90‑1 (MQ=255)                                          gTAAGTTTATCGCCCGGTTGCGCCTGCATCGCCCAGCCGCTGGCGACCCCCCCGTCGTGAATAAAGAAATCAATCGCCAGTTCATGGCGt                                               <  2:326063/90‑1 (MQ=255)                                           cgccaacggtaacggaagagggcatcgtctggccggaaggaccacgCCCACCGCCGTGAATAAAGAAATCAATCGCCAGTTCATGgcgca                                              <  2:219899‑M1/44‑3 (MQ=255)                                            gccaacggtaacggaagagggcatcgtctggccggaaggaccacgCCCACCGTCGGGAATAAAGAAATCAATCGCCAGTTCATGGCGTAg                                             <  2:433599‑M1/45‑1 (MQ=255)                                            aaGTTTATCGCCCGGTTGCGCCTGCATCGCCCAGCCGCTGGCGACCCCACCGTCGTGAATAAAGAAATCAATCGCCAGTTCATGGCGTAg                                             <  2:255367/90‑1 (MQ=255)                                            aaGTTTATCGCCCGGTTGCGCCGGCATCGCCCAGCCGCTGGCGACCCCACCGTCGTGAATAAAGAAATCAATCGCCAGTTCATGGCGTAg                                             <  1:81985/90‑1 (MQ=255)                                                      cccGGTTGCGCCTGCATCGCCCAGCCGCTGGCGACCCCACCGTCGTGAATAAAGAAATCAATCGCCAGTTCATGGCGTAGTTCGTCATAc                                   <  2:200920/90‑1 (MQ=255)                                                        cGGTTGCGCCTGCATCGCCCAGCCGCTGGCGACCCCACCGTCGTGAATAAAGAAATCAATCGCCAGTTCATGGCGTAGTTCGTCATACAg                                 >  1:321452/1‑90 (MQ=255)                                                         ggTTGCGCCTGCATCGCCCAGCCGCTGGCGACCCCACCGTCGTGAATAAAGAAATCAATCGCCAGTTCATGGCGTAGTTCGTCATACAgc                                <  2:199264/90‑1 (MQ=255)                                                           agagggcatcgtctggccggaaggaccaccCCCACCGTCGTGAAAAAAAAAAACAATCGCCAGGTCATGGCGTAGTTCGTCATACAgcgg                              <  1:554671‑M1/60‑1 (MQ=255)                                                                   gCATCGCCCAGCCGCTGGCGACCCCACCGTCGTGAATAAAGAAATCAATCGCCAGTTCATGGCGTAGTTCGTCATACAGCGGCGTATAGt                      >  1:530382/1‑90 (MQ=255)                                                                        gCCCAGCCGCTGGCGACCCCACCGTCGTGAATAAAGAAATCAATCGCCAGTTCATGGCGTAGTTCGTCATACAgcgg                              <  1:718703/77‑1 (MQ=255)                                                                        gCCCAGCCGCTGGCGACCCCACCGTCGTGAATAAAGAAATCAATCGCCAGTTCATGGCGTAGTTCGTCATACAgcgg                              >  2:718703/1‑77 (MQ=255)                                                                                        ccccACCGTCGTGAATAAAGAAATCAATCGCCAGTTCATGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGt  <  1:566783/89‑1 (MQ=255)                                                                                        ccccACCGTCGTGAATAAAGAAATCAATCGCCAGTTCATGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGt  >  2:566783/1‑89 (MQ=255)                                                                                         |                                                                                        AATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGCGCCTGCATCGCCCAGCCGCTGGCGACCCCACCGTCGTGAATAAAGAAATCAATCGCCAGTTCATGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGT  >  NZ_CP009273/3209436‑3209611

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