New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4498181 =NA (NA)3 (0.080) 3/208 4.0 NA noncoding (1/1331 nt) IS2 repeat region
?NC_000913 = 4499512 NA (NA)intergenic (+12/+81) insD1/yjgX IS2 transposase TnpB/pseudogene fragment
Rejected: Coverage evenness skew score above cutoff.

GCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGT  >  NC_000913/4499375‑4499644
                                                                                                                                         |                                                                                                                                    
gCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAcc                                                                                                                                     <  1:719623/139‑1 (MQ=11)
          tCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGat                                                                                                                           <  1:642021/139‑1 (MQ=39)
                    taACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                                 <  2:600593/139‑1 (MQ=255)
                    taACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCGGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                                 <  1:462758/139‑1 (MQ=255)
                              cATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTtata                                                                                                       >  1:815811/1‑139 (MQ=255)
                                tCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACa                                                                                                     >  1:17437/1‑139 (MQ=255)
                                   gCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGaa                                                                                                  >  2:22424/1‑139 (MQ=255)
                                         tGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACttt                                                                                            >  1:467645/1‑139 (MQ=255)
                                          gcgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACtttt                                                                                           <  1:520773/139‑1 (MQ=255)
                                              tGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAcc                                                                                                                                     <  1:461486/93‑1 (MQ=31)
                                              tGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAcc                                                                                                                                     >  2:461486/1‑93 (MQ=31)
                                                ggTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                                 >  1:605219/1‑111 (MQ=255)
                                                ggTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                                 <  2:605219/111‑1 (MQ=255)
                                                    aTCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGc                                                                                 <  2:405897/139‑1 (MQ=255)
                                                      cgctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                                               <  2:467645/139‑1 (MQ=255)
                                                                gggAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACa                                                                     <  2:815811/139‑1 (MQ=255)
                                                                       tCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                                               <  2:19040/122‑1 (MQ=255)
                                                                       tCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                                               >  1:19040/1‑122 (MQ=255)
                                                                           cGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGat                                                                                                               >  2:196383/1‑86 (MQ=255)
                                                                           cGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGat                                                                                                               <  1:196383/86‑1 (MQ=255)
                                                                                 cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAg                                                                                  >  1:301575/1‑109 (MQ=255)
                                                                                 cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAg                                                                                  <  2:301575/109‑1 (MQ=255)
                                                                                     cTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCt                                                <  1:22424/139‑1 (MQ=255)
                                                                                          aagagacagattGATAACAGATGTCTTGAAATATAGGGGCTAATCCACCTTGTTCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGAt                                           <  1:50613/128‑1 (MQ=255)
                                                                                              ggTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAgg                                       >  1:263126/1‑139 (MQ=255)
                                                                                                ttAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGgaa                                     <  2:17437/139‑1 (MQ=255)
                                                                                                        tAACAGATGTCTGGAAATATAGGGGCAAATCCACCTtgtg                                                                                                                                >  2:908907/1‑40 (MQ=255)
                                                                                                        tAACAGATGTCTGGAAATATAGGGGCAAATCCACCTtgtg                                                                                                                                <  1:908907/40‑1 (MQ=255)
                                                                                                        tAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                                 >  2:810910/1‑55 (MQ=255)
                                                                                                        tAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                                 <  1:810910/55‑1 (MQ=255)
                                                                                                             gATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGc                                                                                 >  2:197676/1‑82 (MQ=255)
                                                                                                             gATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGc                                                                                 <  1:197676/82‑1 (MQ=255)
                                                                                                              aTGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATactgtctctt                       >  2:50613/1‑130 (MQ=255)
                                                                                                                             ggggCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTAcc                                                 >  2:40611/1‑98 (MQ=255)
                                                                                                                             ggggCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTAcc                                                 <  1:40611/98‑1 (MQ=255)
                                                                                                                                gCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGc                                                                                 <  2:822648/63‑1 (MQ=255)
                                                                                                                                gCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGc                                                                                 >  1:822648/1‑63 (MQ=255)
                                                                                                                                gCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGAt                                                     <  2:115504/91‑1 (MQ=255)
                                                                                                                                gCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGAt                                                     >  1:115504/1‑91 (MQ=255)
                                                                                                                                gCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTg       >  2:94074/1‑137 (MQ=255)
                                                                                                                                gCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTg       <  1:94074/137‑1 (MQ=255)
                                                                                                                                   aaTCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGt  <  2:263126/139‑1 (MQ=255)
                                                                                                                                         ccTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTg                        <  2:241057/111‑1 (MQ=255)
                                                                                                                                         ccTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTg                        >  1:241057/1‑111 (MQ=255)
                                                                                                                                         |                                                                                                                                    
GCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGT  >  NC_000913/4499375‑4499644

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