New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 419593421 (0.660)3 (0.110) 3/146 2.4 12.2% coding (1527/1590 nt) purH bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase
?NZ_CP009273 4196004 = 25 (0.900)coding (1457/1590 nt) purH bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase
Rejected: Frequency below/above cutoff threshold.

GTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGCCGCTTTAATACCGG  >  NZ_CP009273/4195916‑4196091
                                                                                        |                                                                                       
gTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATc                                                                                        >  2:761631/1‑90 (MQ=255)
gTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATc                                                                                        >  1:546416/1‑90 (MQ=255)
gTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATc                                                                                        >  1:706553/1‑90 (MQ=255)
gTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGCCAGGCTGGATTACGCAAGTCACGCCCGCAGCGGCGGCGGCATc                                                                                        >  2:900344/1‑90 (MQ=255)
  gCTCGTCGGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCaa                                                                                      >  1:808648/1‑90 (MQ=255)
     cgtcgGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATAc                                                                                   >  1:578922/1‑90 (MQ=255)
              cggcAATCCCTCCGTCATCACGGATAGACCCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGa                                                                          <  2:294365/90‑1 (MQ=255)
                       cTTCGTCATCACGGATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGgaaga                                                                 >  1:639842/1‑90 (MQ=255)
                                     aTAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGcc                                                   >  2:317509/1‑90 (MQ=255)
                                     aTAGAACCGCCAGGCTGGATTACGCAGGCCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGcc                                                   <  2:706553/90‑1 (MQ=255)
                                            cGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAgg                                            >  1:291616/1‑90 (MQ=255)
                                            cGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAgg                                            >  2:744783/1‑90 (MQ=255)
                                            cGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAgg                                            <  1:744783/90‑1 (MQ=255)
                                            cGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAgg                                            >  1:655771/1‑90 (MQ=255)
                                             gCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGa                                           >  1:823579/1‑90 (MQ=255)
                                             gCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGa                                           <  2:171724/90‑1 (MQ=255)
                                              ccAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGaa                                          >  2:8894/1‑90 (MQ=255)
                                                     ggATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCa                                   <  1:317509/90‑1 (MQ=255)
                                                      gagcatcgcaataccgtgctcggcggcggcggcaTCAATACCGTCGCGGAAAGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCAc                                  >  2:487673‑M2/35‑90 (MQ=255)
                                                        gcatcgcaataccgtgctcgtcggcggcggcaTCAATACCGTCGCGGAACGGGAAGAAc                                                               >  1:307790‑M2/33‑59 (MQ=255)
                                                        gcatcgcaataccgtgctcgtcggcggcggcaTCAATACCGTCGCGGAACGGGAAGAAc                                                               <  2:307790‑M2/27‑1 (MQ=255)
                                                          aCGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTcc                              <  2:209836/90‑1 (MQ=255)
                                                              aGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGc                          <  2:663138/90‑1 (MQ=255)
                                                               ggTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGcc                         >  1:884232/1‑90 (MQ=255)
                                                                         gCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGcc               <  1:516712/90‑1 (MQ=255)
                                                                         gCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGcc               <  1:761631/90‑1 (MQ=255)
                                                                         gCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGcc               <  2:28/90‑1 (MQ=255)
                                                                         gCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGcc               >  2:269682/1‑90 (MQ=255)
                                                                              ggcggcggcaTCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCAc                                  >  2:423142‑M2/11‑66 (MQ=255)
                                                                              ggcggcggcaTCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCAc                                  <  1:423142‑M2/56‑1 (MQ=255)
                                                                                      caTCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGCCGCTTTAATACCgg  >  2:381586‑M2/3‑90 (MQ=255)
                                                                                      caTCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGCCGCTTTAATACCgg  >  1:187010‑M2/3‑90 (MQ=255)
                                                                                        |                                                                                       
GTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGCCGCTTTAATACCGG  >  NZ_CP009273/4195916‑4196091

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