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.

CATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGG  >  NZ_CP009273/4501662‑4501836
                                                                                     |                                                                                         
catcaGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCgg                                                                                       >  2:742455/1‑90 (MQ=255)
catcaGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCgg                                                                                       >  2:721728/1‑90 (MQ=255)
       gggATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCgag                                                                                <  1:1128029/90‑1 (MQ=255)
       gggATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCgag                                                                                <  1:159757/90‑1 (MQ=255)
       gggATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCgag                                                                                <  2:1274000/90‑1 (MQ=255)
       gggATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCgag                                                                                <  2:931642/90‑1 (MQ=255)
                      acaAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATcc                                                                 >  1:97639/1‑90 (MQ=255)
                       caAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCg                                                                >  1:176760/1‑90 (MQ=255)
                           gCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAg                                                            >  2:1351481/1‑90 (MQ=255)
                            cTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAgg                                                           >  2:165151/1‑90 (MQ=255)
                             tCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGc                                                          <  2:944534/90‑1 (MQ=255)
                              ccAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCa                                                                            <  2:507090/71‑1 (MQ=255)
                              ccAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCa                                                                            >  1:507090/1‑71 (MQ=255)
                              ccAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCt                                                         >  1:475461/1‑90 (MQ=255)
                              ccAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCt                                                         >  2:1583776/1‑90 (MQ=255)
                                        ccccATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCAt                                               >  1:230132/1‑90 (MQ=255)
                                           cATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCgag                                                                                >  1:850306/1‑54 (MQ=255)
                                           cATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCgag                                                                                <  2:850306/54‑1 (MQ=255)
                                           cATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCa                                            >  2:162064/1‑90 (MQ=255)
                                           cATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCa                                            >  2:1118227/1‑90 (MQ=255)
                                               aGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGaaa                                        >  1:490922/1‑90 (MQ=255)
                                               aGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGaaa                                        >  1:394341/1‑90 (MQ=255)
                                               aGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGaaa                                        >  1:1300008/1‑90 (MQ=255)
                                               aGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGaaa                                        <  2:866572/90‑1 (MQ=255)
                                                 cTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAgc                                      <  2:230132/90‑1 (MQ=255)
                                                   gCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGagag                                                                              >  1:4491/1‑48 (MQ=255)
                                                   gCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGagag                                                                              <  2:4491/48‑1 (MQ=255)
                                                   gCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAgcgc                                    >  1:54136/1‑90 (MQ=255)
                                                   gCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAgcgc                                    >  1:771755/1‑90 (MQ=255)
                                                   gCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAgcgc                                    >  1:379273/1‑90 (MQ=255)
                                                      ggTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCAc                                 >  2:1580200/1‑90 (MQ=255)
                                                          aGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCg                             >  2:227280/1‑90 (MQ=255)
                                                                     ggCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAg                  >  2:260398/1‑90 (MQ=255)
                                                                     ggCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAg                  >  2:1046702/1‑90 (MQ=255)
                                                                     ggCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAg                  <  2:177708/90‑1 (MQ=255)
                                                                     ggCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAg                  >  1:368408/1‑90 (MQ=255)
                                                                        gttgttCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCtt               <  1:1046702/90‑1 (MQ=255)
                                                                         ttgttCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTc              >  2:1128398/1‑90 (MQ=255)
                                                                               acaTCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGc        <  2:1466943/90‑1 (MQ=255)
                                                                                caTCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCt       >  2:1070311/1‑90 (MQ=255)
                                                                                    cTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCtggtg   <  1:162064/90‑1 (MQ=255)
                                                                                     tGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCtggtgg  >  2:945158/1‑90 (MQ=255)
                                                                                     tGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCtggggg  >  1:420441/1‑90 (MQ=255)
                                                                                     tGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCtggggg  >  1:1157436/1‑90 (MQ=255)
                                                                                     tGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCtggggg  >  1:1121580/1‑90 (MQ=255)
                                                                                     |                                                                                         
CATCAGCGGGATGGCAATCTTCACAAAGCTCCAGTCACGGCCCCATAAGCTGCCGGTCAGCCACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGG  >  NZ_CP009273/4501662‑4501836

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