New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 392902 =37 (0.170)46 (0.330)
+TTAGCCGCGTGTCCACTATTGCTGGGTAAGATCA
13/118 17.1 40.6% noncoding (65/1258 nt) IS3 repeat region
?NC_000913 = 572471 175 (0.800)coding (6/456 nt) ybcN DLP12 prophage; uncharacterized protein
Rejected: Coverage evenness skew score above cutoff.

TTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAAC  >  NC_000913/392874‑393001
                            |                                                                                                   
ttttACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACa                              <  1:1021591/100‑1 (MQ=31)
 tgtACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACaa                             <  1:4706405/98‑1 (MQ=17)
  ttACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAAc                            <  1:5424620/100‑1 (MQ=33)
   tACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAAcc                           >  1:2040995/1‑100 (MQ=34)
     cTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTg                         <  1:3487359/100‑1 (MQ=37)
       ggTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTa                       >  1:1034687/1‑100 (MQ=39)
       ggTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTa                       >  1:1844231/1‑100 (MQ=39)
        gTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTAc                      <  1:3382970/100‑1 (MQ=40)
         ttGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACg                     <  1:765604/100‑1 (MQ=255)
         ttGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACg                     <  1:472716/100‑1 (MQ=255)
           taTACTGTTTTTGTCATGTGAGTCACTTCTGTCTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGcc                   <  1:702144/99‑1 (MQ=255)
           gATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACa                              >  2:3196629/1‑89 (MQ=31)
             tACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAAc                            >  2:199118/1‑89 (MQ=33)
             tACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCaa                 >  1:3257832/1‑100 (MQ=255)
              aCTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAAcc                           <  2:2677939/89‑1 (MQ=35)
               cTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCt                          >  2:2843713/1‑89 (MQ=37)
               cTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATa               <  1:341614/100‑1 (MQ=255)
               cTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATa               >  1:5095275/1‑100 (MQ=255)
                  tttttGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTa            >  1:474489/1‑100 (MQ=255)
                   ttttGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTAt           <  1:3776056/100‑1 (MQ=255)
                    tttGTCATGTGTGTCACCTCTGACTGAGTGTTTACACACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACg                     <  2:2469098/89‑1 (MQ=14)
                    tttGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATc          >  1:4164959/1‑100 (MQ=255)
                      tGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGcc                   <  2:22958/89‑1 (MQ=255)
                      tGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGcc                   <  2:37681/89‑1 (MQ=255)
                      tGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCAc        >  1:5212161/1‑100 (MQ=255)
                      tGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCAc        <  1:473145/100‑1 (MQ=255)
                      tGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCAc        >  1:2631155/1‑100 (MQ=255)
                      tGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCAc        >  1:1263985/1‑100 (MQ=255)
                       gTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCa                  <  2:1464389/89‑1 (MQ=255)
                       gTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCAcc       <  1:5351809/100‑1 (MQ=255)
                        tCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCg      >  1:696825/1‑100 (MQ=255)
                          aTGTGAGTCACCTCTGACTGAGTGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACCACCTCTACGCCaaca               <  2:1295336/89‑3 (MQ=14)
                          aTGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATa               >  2:268004/1‑89 (MQ=255)
                           tgtgAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATAc              <  2:2135322/89‑1 (MQ=255)
                           tgtgAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATAc              >  2:835522/1‑89 (MQ=255)
                           tgtgAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGtaa   >  1:927822/1‑100 (MQ=255)
                           tgtgAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGtaa   >  1:4082541/1‑100 (MQ=255)
                            gtgAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGtaac  >  1:4735512/1‑100 (MQ=255)
                            |                                                                                                   
TTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAAC  >  NC_000913/392874‑393001

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