New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 3950787 =44 (0.740)4 (0.080) 3/144 3.6 9.5% coding (394/894 nt) ilvY HTH‑type transcriptional activator IlvY
?NZ_CP009273 3950869 = 38 (0.730)coding (312/894 nt) ilvY HTH‑type transcriptional activator IlvY
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CATCGCATCTGCCGCATCACCAGTAGTAAGTTTAATCTCCACCGACGGGTGTTCCGCGCGGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTCACCGAGCAGAAGATATGTAATTCGCCAGAGAGCGACGGCCCTTGCTGATCGATGGTGTGGCGCAACTGCTGAT  >  NZ_CP009273/3950782‑3950957
                                                                                       |                                                                                        
cATCGCATCTGCCGCATCACCAGTAGTAAGTTTAATCTCCACCGACGGGTGTTCCGCGCGGAAGCGATCCAGAATCGGCGGCAGATGGCt                                                                                        >  1:1348520/1‑90 (MQ=255)
 aTCGCATCTGCCGCATCACCAGTAGTAAGTTTAATCTCCACCGACGGGTGTTCCGCGCGGAAGCGATCCAGAATCGGCGGCAGATGGCTg                                                                                       <  2:1308582/90‑1 (MQ=255)
    gCATCTGCCGCATCACCAGTAGTAAGTTTAATCTCCACCTACGGGTGTTCCGCGCGGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAg                                                                                    >  1:1477644/1‑90 (MQ=255)
     cATCTGCCGCATCACCAGTAGTAAGTTTAATCTCCACCGACGGGTGTTCCGCGCGGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAgg                                                                                   <  2:74417/90‑1 (MQ=255)
                 cACCAGTAGTAAGTTTAATCTCCACCGACGGGTGTTCCGCGCGGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTCACCg                                                                       <  2:524304/90‑1 (MQ=255)
                    cAGTAGTAAGTTTAATCTCCACCGACGGGTGTTCCGCGCGGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTCACCGAGc                                                                    <  2:559043/90‑1 (MQ=255)
                       tagtaAGTTTAATCTCCACCGACGGGTGTTCCGCGCGGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTCACCGAGCaga                                                                 <  2:998211/90‑1 (MQ=255)
                       tagtaAGTTTAATCTCCACCGACGGGTGTTCCGCGCGGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTCACCGAGCaga                                                                 <  2:750076/90‑1 (MQ=255)
                       tagtaAGTTTAATCTCCACCGACGGGTGTTCCGCGCGGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTCACCGAGCaga                                                                 <  2:1116148/90‑1 (MQ=255)
                        agtaAGTTTAATCTCCACCGACGGGTGTTCCGCGCGGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTCACCGAGCagaa                                                                >  2:1739655/1‑90 (MQ=255)
                           aaGTTTAATCTCCACCGACGGGTGTTCCGCGCGGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTCACCGAGCaga                                                                 >  1:1085961/1‑86 (MQ=255)
                           aaGTTTAATCTCCACCGACGGGTGTTCCGCGCGGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTCACCGAGCaga                                                                 <  2:1085961/86‑1 (MQ=255)
                           aaGTTTAATCTCCACCGACGGGTGTTCCGCGCGAAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTc                                                                           >  2:1614717/1‑76 (MQ=255)
                           aaGTTTAATCTCCACCGACGGGTGTTCCGCGCGAAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTc                                                                           <  1:1614717/76‑1 (MQ=255)
                                tAATCTCCACCGACGGGTGTTCCGCGCGGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTCACCGAGc                                                                    >  2:336177/1‑78 (MQ=255)
                                tAATCTCCACCGACGGGTGTTCCGCGCGGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTCACCGAGc                                                                    <  1:336177/78‑1 (MQ=255)
                                tAATCTCCACCGACGGGTGTTCCGCGCGGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTCACCGAGCAGAAGATATGTa                                                        >  2:1705180/1‑90 (MQ=255)
                                tAATCTCCACCGACGGGTGTTCCGCGCGGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTCACCGAGCAGAAGATATGTa                                                        >  2:1613964/1‑90 (MQ=255)
                                tAATCTCCACCGACGGGTGTTCCGCGCGGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGGCACCGAGCAGAAGATATGTa                                                        >  2:50139/1‑90 (MQ=255)
                                    ctcCACCGACGGGTGTTCCGCGCGGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTCACCGAGCAGAAGATATGTAATTc                                                    <  1:1705180/90‑1 (MQ=255)
                                       cACCGACGGGTGTTCCGCGCGGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTCACCGAGCAGAAGATATGTAATTCGcc                                                 >  2:1409442/1‑90 (MQ=255)
                                       cACCGACGGGTGTTCCGCGCGGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTCACCGAGCAGAAGATATGTAATTCGcc                                                 >  1:772544/1‑90 (MQ=255)
                                           gACGGGTGTTCCGCGCGGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTCACCGAGc                                                                    >  1:463892/1‑67 (MQ=255)
                                           gACGGGTGTTCCGCGCGGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTCACCGAGc                                                                    <  2:463892/67‑1 (MQ=255)
                                                gtgtTCCGCGCGGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTCACCGAGCAGAAGATATGTAATTCGCCAGAGAGCGa                                        <  2:234807/90‑1 (MQ=255)
                                                 tgtTCCGCGCGGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTCACCGAGCAGAAGATATGTAATTCGCCAGAGAGCGAc                                       <  2:1254930/90‑1 (MQ=255)
                                                  gtTCCGCGCGGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTCACCGAGCAGAAGATATGTAATTCGCCAGAGAGCGACg                                      <  2:286427/90‑1 (MQ=255)
                                                   ttCCGCGCGGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTCACCGAGCAGAAGATATGTAATTCGCCAGAGAGCGACgg                                     >  2:1413696/1‑90 (MQ=255)
                                                     ccGCGCGGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTCACCGAGCAGAAGATATGTAATTCGCCAGAGAGCGACGGcc                                   >  1:1205622/1‑90 (MQ=255)
                                                       gcgcgGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTCACCGAGCAGAAGATATGTAATTCGCCAGAGAGCGACGGCCCt                                 <  2:1477644/90‑1 (MQ=255)
                                                       gcgcgGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTCACCGAGCAGAAGATATGTAATTCGCCAGAGAGCGACGGCCCt                                 <  2:899749/90‑1 (MQ=255)
                                                        cgcgGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTCACCGATCAGAAGATATGTAATTCGCCAGAGAGCGACGGCCCtt                                <  1:285342/90‑1 (MQ=255)
                                                          cgGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTCACCGAGCAGAAGATATGTAATTCGCCAGAGAGCGACGGCCCTTGc                              >  1:1158476/1‑90 (MQ=255)
                                                            gAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTCACCGGGCAGAAGATATGTAATTCGCCAGAGAGCGGCGGCCCctgctg                            >  2:209307/1‑90 (MQ=255)
                                                                 gATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTCACCGAGCAGAAGATATGTAATTCGCCAGAGAGCGACGGCCCTTGCTGATCGa                       >  1:20431/1‑90 (MQ=255)
                                                                            ggcggcAGATGGCTGTAGGCAGCGGTCACCGAGCaga                                                                 >  1:362601/1‑37 (MQ=255)
                                                                            ggcggcAGATGGCTGTAGGCAGCGGTCACCGAGCaga                                                                 <  2:362601/37‑1 (MQ=255)
                                                                            ggcggcAGATGGCTGTAGGCAGCGGTCACCGAGCAGAAGATATGTAATTCGCCAGAGAGCGACGGCCCTTGCTGATCGATGGTGTGgcgc            <  2:416339/90‑1 (MQ=255)
                                                                                 cAGATGGCTGTAGGCAGCGGTCACCGAGCAGAAGATATGTAATTCGCCAGAGAGCGACg                                      >  1:1377691/1‑59 (MQ=255)
                                                                                 cAGATGGCTGTAGGCAGCGGTCACCGAGCAGAAGATATGTAATTCGCCAGAGAGCGACg                                      <  2:1377691/59‑1 (MQ=255)
                                                                                      ggCTGTAGGCAGCGGTCACCGAGCAGAAGATATGTAATTCGCCAGAGAGCGACGGCCCTTGCTGATCGATGGTGTGGCGCAACTGCTGAt  <  1:693313/90‑1 (MQ=255)
                                                                                       |                                                                                        
CATCGCATCTGCCGCATCACCAGTAGTAAGTTTAATCTCCACCGACGGGTGTTCCGCGCGGAAGCGATCCAGAATCGGCGGCAGATGGCTGTAGGCAGCGGTCACCGAGCAGAAGATATGTAATTCGCCAGAGAGCGACGGCCCTTGCTGATCGATGGTGTGGCGCAACTGCTGAT  >  NZ_CP009273/3950782‑3950957

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