New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 301593923 (1.120)4 (0.210) 3/154 1.6 17.1% coding (1265/2871 nt) xdhD molybdopterin‑dependent oxidoreductase Mo/Fe‑S‑binding subunit
?NZ_CP009273 = 3015996 18 (0.970)coding (1322/2871 nt) xdhD molybdopterin‑dependent oxidoreductase Mo/Fe‑S‑binding subunit
Rejected: Frequency below/above cutoff threshold.

ATATCTGCTTTACTCGTATGGACGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAG  >  NZ_CP009273/3015861‑3016021
                                                                              |                                                                                  
atatCTGCTTTACTCGTATGGACGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGcc                                                                         <  2:253818/90‑1 (MQ=255)
     tGCTTTACTCGTATGGACGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTg                                                                    <  2:421489/90‑1 (MQ=255)
      gCTTTACTCGTATGGACGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTgg                                                                   >  1:20418/1‑90 (MQ=255)
         ttACTCGTATGGACGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCa                                                                >  1:124643/1‑90 (MQ=255)
         ttACTCGTATGGACGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCa                                                                >  2:546288/1‑90 (MQ=255)
               gTATGGACGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAAc                                                          <  1:539947/90‑1 (MQ=255)
                tATGGACGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACa                                                         >  2:314495/1‑90 (MQ=255)
                     aCGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATa                                                    >  1:224847/1‑90 (MQ=255)
                      cGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATaa                                                   <  2:144115/90‑1 (MQ=255)
                        gCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAg                                                 >  2:172986/1‑90 (MQ=255)
                        gCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGGGGGCATGAAACAGCATAAAg                                                 >  2:73853/1‑90 (MQ=255)
                           aTCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTc                                              <  1:229151/90‑1 (MQ=255)
                             cGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCAt                                            <  2:219151/90‑1 (MQ=255)
                             cGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCAt                                            <  1:190276/90‑1 (MQ=255)
                               tCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGt                                          <  2:404778/90‑1 (MQ=255)
                                      aTCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTaaa                                   >  1:388106/1‑90 (MQ=255)
                                        ccACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAactcgctctttaatgacatgaactttatgctgtttcatgcccacgaggcgc                                 >  2:323145‑M1/1‑39 (MQ=255)
                                        ccACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAaga                                 >  1:148716/1‑90 (MQ=255)
                                            gCCTCCACCCAGGTACCATGGCACTTACGCCGCCAactcgctctttaatgacatgaactttatgctgtttcatgcccacgaggcgcgcga                             >  2:73723‑M1/1‑35 (MQ=255)
                                            gCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACATCATAAAGTTCATGTCATTAAAGAGCGa                             >  2:152442/1‑90 (MQ=255)
                                                       ggTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGcggcgg                  >  2:324385/1‑90 (MQ=255)
                                                           ccATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTtggcggcgtaagt              >  2:20418‑M2/1‑77 (MQ=255)
                                                           ccATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTtggcggcgtaagt              >  1:314495‑M2/1‑77 (MQ=255)
                                                           ccATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGtttt              <  1:546288/90‑1 (MQ=255)
                                                             aTGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTgg            >  1:287205/1‑90 (MQ=255)
                                                                  aCTTACGCCGCCAactcgctctttaatgacatgaactttatgctgtttcatgcccacgaggcgcgcgacctggcggcgtaagtgccatgg       >  2:20418‑M1/1‑13 (MQ=255)
                                                                  aCTTACGCCGCCAactcgctctttaatgacatgaactttatgctgtttcatgcccacgaggcgcgcgacctggcggcgtaagtgccatgg       >  1:314495‑M1/1‑13 (MQ=255)
                                                                  aCTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCa       >  2:301854/1‑90 (MQ=255)
                                                                       cgccgcCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTCGGCGGCGGTTTTGGTTCCAAACAg  <  1:99497/90‑1 (MQ=255)
                                                                              |                                                                                  
ATATCTGCTTTACTCGTATGGACGGCGATCGTCTGGTTATCCACGCCTCCACCCAGGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAG  >  NZ_CP009273/3015861‑3016021

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