New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	= 3209533	20 (0.740)	3 (0.130)	3/150	2.1	14.1%	coding (318/765 nt)	nfeF	NADPH‑dependent ferric chelate reductase
	?	NZ_CP009273	= 3209595	19 (0.790)					coding (256/765 nt)	nfeF	NADPH‑dependent ferric chelate reductase
Rejected: Frequency below/above cutoff threshold.

GCCCAGCCGCTGGCGACCCCACCGTCGTGAATAAAGAAATCAATCGCCAGTTCATGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAG  >  NZ_CP009273/3209507‑3209675                                                                                         |                                                                                 gCCCAGCCGCTGGCGACCCCACCGTCGTGAATAAAGAAATCAATCGCCAGTTCATGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGc                                                                                 >  2:348234‑M2/1‑89 (MQ=255)      gCCGCTGGCGACCCCACCGTCGTGAATAAAGAAATCAATCGCCAGTTCATGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGcgacgg                                                                            <  1:118396‑M2/90‑7 (MQ=255)            ggCGACCCCCCCGTCGTGAATAAAGAAATCAATCGCCAGTCCATGGCGTAGTCCGCCATCCAGCGGCGTATAGCCACGcgacggtggggt                                                                      >  1:654260‑M2/1‑78 (MQ=255)                    cACCGTCGTGAATAAAGAAATCAATCGCCAGTTCATGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCtt                                                              >  2:695752/1‑90 (MQ=255)                         tcgtGAATAAAGAAATCAATCGCCAGTTCATGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGc                                                         <  2:383540/90‑1 (MQ=255)                             gAATAAAGAAATCAATCGCCAGTTCATGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTc                                                             >  1:210554/1‑82 (MQ=255)                             gAATAAAGAAATCAATCGCCAGTTCATGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTc                                                             <  2:210554/82‑1 (MQ=255)                                tAAAGAAATCAATCGCCAGTTCATGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGa                                                  <  1:21230/90‑1 (MQ=255)                                tAAAGAAATCAATCGCCAGTTCATGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGa                                                  <  1:428248/90‑1 (MQ=255)                                tAAAGAAATCAATCGCCAGTTCAGGGCGTGGTTCGTCATACAGCGGCGTATAGTCACGcgacggtggggtcgccagcggctgggcgatgc                                                  >  1:424548‑M2/1‑58 (MQ=255)                                      aatcaatcGCCAGTTCATGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCct                                            <  1:688199/90‑1 (MQ=255)                                        tcaatcGCCAGTTCATGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCctct                                          <  2:118396/90‑1 (MQ=255)                                              gctaGTTCATGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGtt                                    >  2:743643/4‑90 (MQ=255)                                                 aGTTCATGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGcgacggtggggtcgccagcggctgggcgatgcaggcgcaaccgggcgat                                 >  1:573671‑M2/1‑41 (MQ=255)                                                       tGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGAAGATGCCCTCTTCCGTTACCGTTggc                           >  1:604411/1‑90 (MQ=255)                                                            tAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTTCCGTTGGCGGcac                      <  2:770886/90‑1 (MQ=255)                                                            tAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGcac                      <  1:739069/90‑1 (MQ=255)                                                                  gTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTg                <  1:465544/90‑1 (MQ=255)                                                                   tCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGa               >  2:251041/1‑90 (MQ=255)                                                                          gcggcgTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAg        >  2:61175/1‑90 (MQ=255)                                                                          gcggcgTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAg        >  1:383835/1‑90 (MQ=255)                                                                            ggcgTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGGGAGCGTCAGGt      >  2:567108/1‑90 (MQ=255)                                                                             gcgTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGtt     >  2:469524/1‑90 (MQ=255)                                                                                tataGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAg  <  1:220486/90‑1 (MQ=255)                                                                                         |                                                                                 GCCCAGCCGCTGGCGACCCCACCGTCGTGAATAAAGAAATCAATCGCCAGTTCATGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAG  >  NZ_CP009273/3209507‑3209675

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.