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.

GGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTC  >  NZ_CP009273/3015916‑3016084
                                                                                |                                                                                        
ggTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGcggcgg                                                                                 >  2:324385/1‑90 (MQ=255)
    ccATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTtggcggcgtaagt                                                                             >  1:314495‑M2/1‑77 (MQ=255)
    ccATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTtggcggcgtaagt                                                                             >  2:20418‑M2/1‑77 (MQ=255)
    ccATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGtttt                                                                             <  1:546288/90‑1 (MQ=255)
      aTGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTgg                                                                           >  1:287205/1‑90 (MQ=255)
           aCTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCa                                                                      >  2:301854/1‑90 (MQ=255)
                cgccgcCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTCGGCGGCGGTTTTGGTTCCAAACAg                                                                 <  1:99497/90‑1 (MQ=255)
                          tCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTtggcggcgtaagtgccatggtacctgggtggaggc                                                       >  2:73723‑M2/1‑55 (MQ=255)
                              gcgcCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTtggcggcgtaagtgccatggtacctgggtggaggcgtgg                                                   >  2:323145‑M2/1‑51 (MQ=255)
                                           aTGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGcc                                      >  1:79930/1‑90 (MQ=255)
                                              aaaCAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTgg                                   <  2:124643/90‑1 (MQ=255)
                                                  aGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCaa                               >  1:545988/1‑90 (MQ=255)
                                                  aGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCaa                               <  1:172986/90‑1 (MQ=255)
                                                             cATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGAcc                    >  2:234802/1‑90 (MQ=255)
                                                             cATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGAcc                    >  2:416292/1‑90 (MQ=255)
                                                             cATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGAcc                    >  2:526535/1‑90 (MQ=255)
                                                                gTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCggg                 <  1:320458/90‑1 (MQ=255)
                                                                     tAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTcc            <  2:388106/90‑1 (MQ=255)
                                                                          agCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTAc       <  1:429081/90‑1 (MQ=255)
                                                                           gCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACt      <  2:375667/90‑1 (MQ=255)
                                                                            cGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTg     >  2:57007/1‑90 (MQ=255)
                                                                               gTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTc  <  1:258256/90‑1 (MQ=255)
                                                                                |                                                                                        
GGTACCATGGCACTTACGCCGCCAGGTCGCGCGCCTCGTGGGCATGAAACAGCATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGCGGTTTTGGTTCCAAACAGGACATCCTGCTGGAAGAAGTGTGCGCCTGGGCAACCTGCGTGACCGGGCGTCCGGTACTGTTC  >  NZ_CP009273/3015916‑3016084

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