New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 4501747 =41 (0.770)3 (0.060) 3/150 4.0 6.7% coding (485/957 nt) fecD Fe(3+) dicitrate ABC transporter permease subunit FecD
?NZ_CP009273 4501885 = 46 (0.960)coding (347/957 nt) fecD Fe(3+) dicitrate ABC transporter permease subunit FecD
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCG  >  NZ_CP009273/4501802‑4501974
                                                                                   |                                                                                         
ccACGCCGGTGAGCGCCAGCTTCATCGGCTGGGGGGGCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCca                                                                                     >  1:1050887/1‑90 (MQ=255)
 cACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCcag                                                                                    >  2:1567647/1‑90 (MQ=255)
  ccGCCGGTGAGCGCCACCTTCATCGGCGGGTGGGTTTTTGCCAGCATCTTAAGCAATACCAACCCCGCCATGCCGCCCGCAAAGGCcagc                                                                                   <  2:1255427/89‑1 (MQ=255)
        gTGAGCGCCAGCTTCACCGGCTGGGGGGTCTTTGCCAGCATCTTCAGCAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGgc                                                                             <  1:1477843/90‑1 (MQ=255)
          gAGCGCCAGCTTCATCGGCTGGTGGGTCTTTACCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGgcag                                                                           >  2:915039/1‑90 (MQ=255)
               ccAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCa                                                                      >  1:858296/1‑90 (MQ=255)
                                 gggTTTTTGCCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCa                                                    <  1:742455/90‑1 (MQ=255)
                                 gggTCTTTGCCAGCATCTTCAGCAATCCCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCa                                                    <  1:53595/90‑1 (MQ=255)
                                 gggTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCa                                                    >  1:1411316/1‑90 (MQ=255)
                                   gtatttGCCAGCACCTTCAGCAATATAACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCGTa                                                  <  1:665723/87‑1 (MQ=255)
                                   gTCTTTGCCGCTATCTTCCGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATa                                                  <  1:39491/90‑1 (MQ=255)
                                   gTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATa                                                  >  2:818554/1‑90 (MQ=255)
                                   gTCTTTGCCAGCATCTTCAGCAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATa                                                  <  1:1383461/90‑1 (MQ=255)
                                     cTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATaag                                                >  2:79414/1‑90 (MQ=255)
                                      tttGCCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATaaga                                               <  1:915039/90‑1 (MQ=255)
                                       ttGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATa                                                  >  1:844050/1‑86 (MQ=255)
                                       ttGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATa                                                  <  2:844050/86‑1 (MQ=255)
                                           cAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTag                                          <  1:1560838/90‑1 (MQ=255)
                                            aGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTaga                                         <  2:1050887/90‑1 (MQ=255)
                                             gCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGAc                                                        <  2:6134/74‑1 (MQ=255)
                                             gCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGAc                                                        >  1:6134/1‑74 (MQ=255)
                                             gCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTagag                                        >  2:1190700/1‑90 (MQ=255)
                                             gCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTagag                                        >  2:686749/1‑90 (MQ=255)
                                             gCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTagag                                        >  1:493006/1‑90 (MQ=255)
                                                 cTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGcccc                                    <  2:483576/90‑1 (MQ=255)
                                                         atatCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACag                               <  1:353731/87‑1 (MQ=255)
                                                         atatCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACag                               >  2:353731/1‑87 (MQ=255)
                                                           atCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGcc                          >  2:1194805/1‑90 (MQ=255)
                                                              aaCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACag                               <  2:526664/82‑1 (MQ=255)
                                                              aaCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACag                               >  1:526664/1‑82 (MQ=255)
                                                                ccccGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCcaca                                <  1:1255377/79‑1 (MQ=255)
                                                                ccccGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCcaca                                >  2:1255377/1‑79 (MQ=255)
                                                                ccccGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCt                     >  1:486301/1‑90 (MQ=255)
                                                                   cGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTagag                                        >  2:936750/1‑68 (MQ=255)
                                                                   cGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTagag                                        <  1:936750/68‑1 (MQ=255)
                                                                    gCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCAc                                                                        >  1:925892/1‑35 (MQ=255)
                                                                    gCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCAc                                                                        <  2:925892/35‑1 (MQ=255)
                                                                    gCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCa                                                                      <  1:229254/37‑1 (MQ=255)
                                                                    gCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCa                                                                      >  2:229254/1‑37 (MQ=255)
                                                                     ccATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTggcgg                >  2:1158866/1‑90 (MQ=255)
                                                                       aTGGCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCa              <  1:632316/90‑1 (MQ=255)
                                                                       aTGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCa              <  1:1266619/90‑1 (MQ=255)
                                                                       aTGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCa              <  1:1580200/90‑1 (MQ=255)
                                                                       aTGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCa              <  2:214854/90‑1 (MQ=255)
                                                                       aTGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCa              <  2:1545680/90‑1 (MQ=255)
                                                                       aTGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCa              <  1:1246486/90‑1 (MQ=255)
                                                                        tGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCAt             <  1:829719/90‑1 (MQ=255)
                                                                            gcccgcAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGtt         <  2:808909/90‑1 (MQ=255)
                                                                             cccgcAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTaga                                         <  1:714077/57‑1 (MQ=255)
                                                                             cccgcAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTaga                                         >  2:714077/1‑57 (MQ=255)
                                                                             cccgcAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTa        <  1:1118227/90‑1 (MQ=255)
                                                                              ccgcAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTaa       <  1:1032184/90‑1 (MQ=255)
                                                                                 cAAAGGCCAGCAGCGGCAGCCCCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAacac    <  1:916557/90‑1 (MQ=255)
                                                                                 cAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAacac    <  1:818554/90‑1 (MQ=255)
                                                                                   aaGGCCAGCAGCGGCAGCCCCATCCCGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCg  <  1:1085386/90‑1 (MQ=255)
                                                                                   |                                                                                         
CCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCG  >  NZ_CP009273/4501802‑4501974

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