New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 359329830 (0.410)28 (0.390) 23/98 1.1 49.6% coding (140/714 nt) livF branched chain amino acid/phenylalanine ABC transporter ATP binding subunit LivF
?NC_000913 3593309 = 27 (0.370)coding (129/714 nt) livF branched chain amino acid/phenylalanine ABC transporter ATP binding subunit LivF

CACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACAATCTCGCC  >  NC_000913/3593249‑3593347
                                                 |                                                 
cacaATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCa                                                   <  1:2723995/50‑1 (MQ=255)
cacaATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCa                                                   >  1:2863467/1‑50 (MQ=255)
cacaATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCa                                                   <  1:4385612/50‑1 (MQ=255)
 acaATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAc                                                  <  1:2217418‑M1/50‑2 (MQ=255)
    aTTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAggt                                               >  1:3411640/1‑50 (MQ=255)
       cGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAggtgg                                             <  1:1518908/49‑1 (MQ=255)
        gCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAccccgccc                                           >  1:2537747‑M1/1‑42 (MQ=255)
        gCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAccccgccc                                           >  1:2447823‑M1/1‑42 (MQ=255)
        gCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAccccgccc                                           >  1:5404427‑M1/1‑42 (MQ=255)
         cccGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGttt                                          >  1:685842/1‑50 (MQ=255)
          ccGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAccccgccccg                                         >  1:2702752‑M1/1‑40 (MQ=255)
           cGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTc                                        >  1:1358872/1‑50 (MQ=255)
           cGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTc                                        >  1:304850/1‑50 (MQ=255)
            gCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTcc                                       <  1:402088/50‑1 (MQ=255)
             cTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTccc                                      >  1:433028/1‑50 (MQ=255)
             cTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTccc                                      <  1:5141282/50‑1 (MQ=255)
              tggtggCACGCGGATCGCCGCATAACGTGCCGAGCAccccgccccgttcg                                     >  1:3177903‑M1/1‑36 (MQ=255)
              tggtggCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTcccc                                     >  1:5187769/1‑50 (MQ=255)
               ggtggCACGCGGATCGCCGCATAACGTGCCGAGCAccccgccccgttcgc                                    >  1:1740940‑M1/1‑35 (MQ=255)
               ggtggCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTccccg                                    >  1:2768977/1‑50 (MQ=255)
                gtggCACGCGGATCGCCGCATAACGTGCCGAGCAccccgccccgttcg                                     <  1:1809663‑M1/48‑15 (MQ=255)
                gtggCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTcccc                                     <  1:3831095/48‑1 (MQ=255)
                   gCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTccccgcccc                                >  1:1164307/1‑50 (MQ=255)
                    cACGCGGATCGCCGCATAACGTGCCGAGCAccccgccccgttcgcgccaa                               >  1:3725427‑M1/1‑30 (MQ=255)
                     aCGCGGATCGCCGCATAACGTGCCGAGCAccccgccccgttcgcgccaat                              >  1:1684056‑M1/1‑29 (MQ=255)
                     aCGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGt                              >  1:1290533/1‑50 (MQ=255)
                       gcgGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTc                            >  1:4913412/1‑50 (MQ=255)
                       gcgGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTc                            <  1:2396154/50‑1 (MQ=255)
                       gcgGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTc                            >  1:3002905/1‑50 (MQ=255)
                        cgGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTcg                           <  1:2354910/50‑1 (MQ=255)
                         ggATCGCCGCATAACGTGCCGAGCAccccgccccgttcgcgccaatcagc                          <  1:1395923‑M1/50‑26 (MQ=255)
                         ggATCGCCGCATAACGTGCCGAGCAccccgccccgttcgcgccaatcagc                          >  1:3413245‑M1/1‑25 (MQ=255)
                         ggATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTcgc                          >  1:2533881/1‑50 (MQ=255)
                         ggATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTcgc                          >  1:2584030/1‑50 (MQ=255)
                           aTCGCCGCATAACGTGCCGAGCAccccgccccgttcgcgccaatcagcgt                        <  1:2462804‑M1/50‑28 (MQ=255)
                           aTCGCCGCATAACGTGCCGAGCAccccgccccgttcgcgccaatcagcgt                        >  1:1583548‑M1/1‑23 (MQ=255)
                            tCGCCGCATAACGTGCCGAGCAccccgccccgttcgcgccaatcagcgtg                       <  1:3219265‑M1/50‑29 (MQ=255)
                             cgccgcATAACGTGCCGAGCAccccgccccgttcgcgccaatcagcgtga                      >  1:2090133‑M1/1‑21 (MQ=255)
                             cgccgcATAACGTGCCGAGCAccccgccccgttcgcgccaatcagcgtga                      >  1:4355630‑M1/1‑21 (MQ=255)
                              gccgcATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCaa                     >  1:1648837/1‑50 (MQ=255)
                               ccgcATAACGTGCCGAGCAccccgccccgttcgcgccaatcagcgtgaca                    <  1:4440703‑M1/50‑32 (MQ=255)
                                cgcATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAAt                    >  1:1191620/1‑49 (MQ=255)
                                cgcATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATc                   <  1:610121/50‑1 (MQ=255)
                                 gcATAACGTGCCGAGCAccccgccccgttcgcgccaatcagcgtgacaat                  >  1:682507‑M1/1‑17 (MQ=255)
                                 gcATAACGTGCCGAGCAccccgccccgttcgcgccaatcagcgtgacaat                  <  1:3952354‑M1/50‑34 (MQ=255)
                                   aTAACGTGCCGAGCAccccgccccgttcgcgccaatcagcgtgacaatct                >  1:72730‑M1/1‑15 (MQ=255)
                                    tAACGTGCCGAGCAccccgccccgttcgcgccaatcagcgtgacaatctc               >  1:3952148‑M1/1‑14 (MQ=255)
                                    tAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCg               <  1:2532827/50‑1 (MQ=255)
                                       cGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGa            >  1:4756509/1‑50 (MQ=255)
                                        gTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGa            >  1:3382556/1‑49 (MQ=255)
                                          gCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACaa         >  1:1623614/1‑50 (MQ=255)
                                          gCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACaa         >  1:113584/1‑50 (MQ=255)
                                            cGAGCAccccgccccgttcgcgccaatcagcgtgacaatctcgccctgat       >  1:3915632‑M1/1‑6 (MQ=255)
                                            cGAGCAccccgccccgttcgcgccaatcagcgtgacaatctcgccctgat       >  1:129540‑M1/1‑6 (MQ=255)
                                             gAGCAccccgccccgttcgcgccaatcagcgtgacaatctcgccctgatt      <  1:4261602‑M1/50‑46 (MQ=255)
                                             gAGCAccccgccccgttcgcgccaatcagcgtgacaatctcgccctgatt      <  1:3707557‑M1/50‑46 (MQ=255)
                                             gAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACAAtct      >  1:5130974/1‑50 (MQ=255)
                                                cAccccgccccgttcgcgccaatcagcgtgacaatctcgccctgattgat   <  1:1744475‑M1/50‑49 (MQ=255)
                                                 accccgccccgttcgcgccaatcagcgtgacaatctcgccctgattgata  >  1:1444501‑M1/1‑1 (MQ=255)
                                                 aaGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACAATCTCGcc  >  1:1331213/1‑50 (MQ=255)
                                                 aaGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACAATCTCGcc  >  1:3317151/1‑50 (MQ=255)
                                                 |                                                 
CACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACAATCTCGCC  >  NC_000913/3593249‑3593347

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