New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 31106790 (0.000)49 (0.450)
+GA
36/94 0.9 100% coding (448/537 nt) yghD putative type II secretion system M‑type protein
?NC_000913 3136382 = 0 (0.000)intergenic (‑11/+281) pitB/gss metal phosphate:H(+) symporter PitB/fused glutathionylspermidine amidase/glutathionylspermidine synthetase

CAGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATATTTTTCATCCAGTGCGTTTAACCATTTCAGCAGGTCATTAAACACCACA  >  NC_000913/3110631‑3110728
                                                |                                                 
cAGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAg                                                  <  1:5039711‑M1/50‑2 (MQ=255)
cAGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAg                                                  <  1:2371577‑M1/50‑2 (MQ=255)
  ggCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaa                                                <  1:4242171‑M1/50‑4 (MQ=255)
  ggCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaa                                                <  1:7123461‑M1/50‑4 (MQ=255)
  ggCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaa                                                >  1:2426914‑M1/1‑47 (MQ=255)
   gCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaac                                               <  1:2022399‑M1/50‑5 (MQ=255)
    cTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaacg                                              <  1:7359841‑M1/50‑6 (MQ=255)
     ttCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgc                                             >  1:7101253‑M1/1‑44 (MQ=255)
     ttCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgc                                             >  1:7282794‑M1/1‑44 (MQ=255)
       ctcAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgc                                             <  1:8304345‑M1/48‑7 (MQ=255)
         cagcagCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccatt                                         >  1:1697906‑M1/1‑40 (MQ=255)
           gcagcaCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccatttt                                       <  1:6492635‑M1/50‑13 (MQ=255)
           gcagcaCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccatttt                                       <  1:3054664‑M1/50‑13 (MQ=255)
            cagcaCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttg                                      <  1:2050623‑M1/50‑14 (MQ=255)
             agcaCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttga                                     <  1:5490072‑M1/50‑15 (MQ=255)
             agcaCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttga                                     <  1:3542044‑M1/50‑15 (MQ=255)
              gcaCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaa                                    <  1:5733946‑M1/50‑16 (MQ=255)
                 cTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatt                                  >  1:6963717‑M1/1‑32 (MQ=255)
                  tGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaattta                                >  1:5234788‑M1/1‑31 (MQ=255)
                    aCATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttacg                              <  1:257434‑M1/50‑22 (MQ=255)
                     cATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttacga                             <  1:4051527‑M1/50‑23 (MQ=255)
                       tCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaac                           <  1:7725607‑M1/50‑25 (MQ=255)
                       tCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaac                           <  1:1320749‑M1/50‑25 (MQ=255)
                       tCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaac                           <  1:2627537‑M1/50‑25 (MQ=255)
                         aaTTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaacgg                         <  1:597636‑M1/50‑27 (MQ=255)
                         aaTTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaacgg                         <  1:2074027‑M1/50‑27 (MQ=255)
                         aaTTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaacg                          >  1:1275965‑M1/1‑24 (MQ=255)
                         aaTTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaacg                          >  1:2090122‑M1/1‑24 (MQ=255)
                          aTTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaacgga                        >  1:6652295‑M1/1‑23 (MQ=255)
                           tttGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaacggac                       >  1:6977032‑M1/1‑22 (MQ=255)
                           tttGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaacggac                       <  1:1042451‑M1/50‑29 (MQ=255)
                            ttGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaacggacg                      >  1:2820864‑M1/1‑21 (MQ=255)
                                 tCACCCGCAGCGCATAgaacgccattttgaatttacgaacggacgcatta                 <  1:76370‑M1/50‑35 (MQ=255)
                                  cACCCGCAGCGCATAgaacgtcattttgaatttacgaacggacgcattat                >  1:2670641‑M1/1‑15 (MQ=255)
                                  cACCCGCAGCGCATAgaacgccattttgaatttacgaacggacgcattat                >  1:4473979‑M1/1‑15 (MQ=255)
                                   aCCCGCAGCGCATAgaacgccattttgaatttacgaacggacgcattat                <  1:8005410‑M1/49‑36 (MQ=255)
                                   aCCCGCAGCGCATAgaacgccattttgaatttacgaacggacgcattat                <  1:5213407‑M1/49‑36 (MQ=255)
                                   aCCCGCAGCGCATAgaacgccattttgaatttacgaacggacgcattat                >  1:6769902‑M1/1‑14 (MQ=255)
                                     ccGCAGCGCATAgaacgccattttgaatttacgaacggacgcattatcag             >  1:2895095‑M1/1‑12 (MQ=255)
                                      cGCAGCGCATAgaacgccattttgaatttacgaacggacgcattatcagt            >  1:7314335‑M1/1‑11 (MQ=255)
                                       gCAGCGCATAgaacgccattttgaatttacgaacggacgcattatcagtg           <  1:5433848‑M1/50‑41 (MQ=255)
                                       gCAGCGCATAgaacgccattttgaatttacgaacggacgcattatcagtg           <  1:6840972‑M1/50‑41 (MQ=255)
                                       gCAGCGCATAgaacgccattttgaatttacgaacggacgcattatcagt            >  1:7011644‑M1/1‑10 (MQ=255)
                                       gCAGCGCATAgaacgccatgttgaatttacgaacggacgcattatcagtg           <  1:3317244‑M1/50‑41 (MQ=255)
                                        cAGCGCATAgaacgccattttgaatttacgaacggacgcattatcagtga          <  1:8374081‑M1/50‑42 (MQ=255)
                                          gcgcATAgaacgccattttgaatttacgaacggacgcattatcagtgact        <  1:7391388‑M1/50‑44 (MQ=255)
                                            gcATAgaacgccattttgaatttacgaacggacgcattatcagtgacttt      <  1:3824924‑M1/50‑46 (MQ=255)
                                            gcATAgaacgccattttgaatttacgaacggacgcattatcagtgacttt      <  1:1149658‑M1/50‑46 (MQ=255)
                                            gcATAgaacgccattttgaatttacgaacggacgcattatcagtgactt       >  1:57922‑M1/1‑5 (MQ=255)
                                               tAgaacgccattttgaatttacgaacggacgcattatcagtgactttaac   <  1:6831809‑M1/50‑49 (MQ=255)
                                                agaacgccattttgaatttacgaacggacgcattatcagtgactttaacg  >  1:5448706‑M1/1‑1 (MQ=255)
                                                |                                                 
CAGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATATTTTTCATCCAGTGCGTTTAACCATTTCAGCAGGTCATTAAACACCACA  >  NC_000913/3110631‑3110728

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