New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 345392719 (0.170)76 (0.740) 46/92 0.5 86.5% coding (1/420 nt) gspB putative general secretion pathway protein B
?NC_000913 3469854 = 6 (0.060)intergenic (‑1/+291) chiA/tufA endochitinase/translation elongation factor Tu 1

ACACCAGCCCCATACCAATCATAGATTTAGTAAATATATTTAATTTCATTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTACGCAGGAA  >  NC_000913/3469805‑3469903
                                                 |                                                 
ggtgacctgtttctttctgctcacgcgcggcgatatagaattcaaacatt                                                   <  1:105828‑M2/1‑1 (MQ=255)
ggtgacctgtttctttctgctcacgcgcggcgatatagaattcaaacatt                                                   <  2:2630528‑M2/1‑1 (MQ=255)
 gtgacctgtttctttctgctcacgcgcggcgatatagaattcaaacatta                                                  <  1:666194‑M2/2‑1 (MQ=255)
 gtgacctgtttctttctgctcacgcgcggcgatatagaattcaaacatta                                                  >  1:3391744‑M2/49‑50 (MQ=255)
  tgacctgtttctttctgctcacgcgcggcgatatagaattcaaacattat                                                 <  2:1912412‑M2/3‑1 (MQ=255)
  tgacctgtttctttctgctcacgcgcggcgatatagaattcaaacatta                                                  >  2:2820142‑M2/48‑49 (MQ=255)
   gacctgtttctttctgctcacgcgcggcgatatagaattcaaacattat                                                 >  2:2547185‑M2/47‑49 (MQ=255)
    acctgtttctttctgctcacgcgcggcgatatagaattcaaacattata                                                <  2:2302764‑M2/4‑1 (MQ=255)
    acctgtttctttctgctcacgcgcggcgatatagaattcaaacattata                                                <  1:944946‑M2/4‑1 (MQ=255)
    acctgtttctttctgctcacgcgcggcgatatagaattcaaacattata                                                <  1:4875518‑M2/4‑1 (MQ=255)
    acctgtttctttctgctcacgcgcggcgatatagaattcaaacatTATaa                                               <  1:4685729‑M2/5‑1 (MQ=255)
    acctgtttctttctgctcacgcgcggcgatatagaattcaaacatTATaa                                               <  1:4361653‑M2/5‑1 (MQ=255)
    acctgtttctttctgctcacgcgcggcgatatagaattaaaacattata                                                <  1:4204668‑M2/4‑1 (MQ=255)
     cctgtttctttctgctcacgcgcggcgatatagaattcaaacatTATAAt                                              <  2:1265941‑M2/6‑1 (MQ=255)
      ctgtttctttctgctcacgcgcggcgatatagaattcaaacatTATATTc                                             >  1:4719847‑M2/44‑50 (MQ=255)
         tttctttctgctcacgcgcggcgatatagaattcaaacatTATAATccc                                           <  2:1918404‑M2/9‑1 (MQ=255)
         tttctttctgctcacgcgcggcgatatagaattcaaacatTATAATCCCt                                          >  1:4011803‑M2/41‑50 (MQ=255)
         tttctttctgctcacgcgcggcgatatagaattcaaacatTATAATCCCt                                          >  1:2849297‑M2/41‑50 (MQ=255)
         tttctttctgctcacgcgcggcgatatagaattcaaacatTATAATCCCt                                          <  2:528530‑M2/10‑1 (MQ=255)
          ttctttctgctcacgcgcggcgatatagaattcaaacatTATAATCCCtt                                         >  1:892637‑M2/40‑50 (MQ=255)
           tctttctgctcacgcgcggcgatatagaattcaaacatTATAATCCCtt                                         <  1:2422440‑M2/11‑1 (MQ=255)
           tctttctgctcacgcgcggcgatatagaattcaaacatTATAATCCCtt                                         <  2:4901613‑M2/11‑1 (MQ=255)
           tctttctgctcacgcgcggcgatatagaattcaaacatTATAATCCCtt                                         <  2:1212645‑M2/11‑1 (MQ=255)
           tctttctgctcacgcgcggcgatatagaattcaaacatTATAATCCCTtg                                        >  2:3181907‑M2/39‑50 (MQ=255)
           tctttctgctcacgcgcggcgatatagaattcaaacatTATAATCCCTtg                                        >  2:1232810‑M2/39‑50 (MQ=255)
             tttctgctcacgcgcggcgatatagaattcaaacatTATAATCCCTtgt                                       >  2:1160826‑M2/37‑49 (MQ=255)
              ttctgctcacgcgcggcgatatagaattcaaacatTATAATCCCTTGTGa                                     >  2:101242‑M2/36‑50 (MQ=255)
              ttctgctcacgcgcggcgatatagaattcaaacatTATAATCCCTTGTGa                                     >  1:2448521‑M2/36‑50 (MQ=255)
               tctgctcacgcgcggcgatatagaattcaaacatTATAATCCCTtgtg                                      <  1:738408‑M2/14‑1 (MQ=255)
               tctgctcacgcgcggcgatatagaattcaaacatTATAATCCCTtgtg                                      <  2:3295227‑M2/14‑1 (MQ=255)
               tctgctcacgcgcggcgatatagaattcaaacatTATAATCCCTtgt                                       <  1:3564608‑M2/13‑1 (MQ=255)
               tctgctcacgcgcggcgatatagaattcaaacatTATAATCCCTtgt                                       <  1:3300497‑M2/13‑1 (MQ=255)
                ctgctcacgcgcggcgatatagaattcaaacatTATAATCCCTTGTGAc                                    <  1:3719104‑M2/16‑1 (MQ=255)
                aaTCATAGATTTAGTAAATATATTTAATTTCATTATAATCCCTTGTGACg                                   >  1:1111794/1‑50 (MQ=255)
                 tgctcacgcgcggcgatatagaattcaaacatTATAATCCCTTGTGACg                                   >  1:3791346‑M2/33‑49 (MQ=255)
                  gctcacgcgcggcgatatagaattcaaacatTATAATCCCTTGTGACGt                                  >  1:5259300‑M2/32‑49 (MQ=255)
                   ctcacgcgcggcgatatagaattcaaacatTATAATCCCTTGTGACGt                                  <  2:1292284‑M2/18‑1 (MQ=255)
                      acgcgcggcgatatagaattcaaacatTATAATCCCTTGTGACGTaaaa                              <  1:4625305‑M2/22‑1 (MQ=255)
                      acgcgcggcgatatagaattcaaacatTATAATCCCTTGTGACGTaaaa                              <  2:4513118‑M2/22‑1 (MQ=255)
                      acgcgcggcgatatagaattcaaacatTATAATCCCTTGTGACGTaaaa                              <  2:1993009‑M2/22‑1 (MQ=255)
                      acgcgcggcgatatagaattcaaacatTATAAACCCTTGTGACGTaaaa                              <  2:3240135‑M2/22‑1 (MQ=255)
                       cgcgcggcgatatagaattcaaacatTATAATCTCTTGTGACGTAAAAAc                            <  2:270749‑M2/24‑1 (MQ=255)
                       cgcgcggcgatatagaattcaaacatTATAATCCCTTGTGACGTAAAAAc                            <  1:4605947‑M2/24‑1 (MQ=255)
                       cgcgcggcgatatagaattcaaacatTATAATCCCTTGTGACGTAAAAAc                            <  2:764842‑M2/24‑1 (MQ=255)
                        gcgcggcgatatagaattcaaacatTATAATCCCTTGTGACGTAAAAACt                           <  1:4234898‑M2/25‑1 (MQ=255)
                          gcggcgatatagaattcaaacatTATAATCCCTTGTGACGTAAAAACTg                          >  2:825135‑M2/24‑49 (MQ=255)
                           cggcgatatagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGc                         >  2:5071413‑M2/23‑49 (MQ=255)
                           cggcgatatagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGc                         >  1:2448928‑M2/23‑49 (MQ=255)
                            ggcgatatagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGca                        >  1:1526594‑M2/22‑49 (MQ=255)
                            aGTAAATATATTTAATTTCATTATAATCCCTTGTGACGTAAAAACTGca                        >  1:4431685/1‑49 (MQ=255)
                             gcgatatagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGcaaa                      <  2:2284690‑M2/30‑1 (MQ=255)
                              tAAATATATTTAATTTCATTATAATCCCTTGTGACGTAAAAACTGcaaaa                     <  2:2538978/50‑1 (MQ=255)
                               gatatagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGcaaaac                    <  1:1473466‑M2/32‑1 (MQ=255)
                               gatatagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGcaaaac                    <  2:4157676‑M2/32‑1 (MQ=255)
                               gatatagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGcaaaac                    <  1:3161148‑M2/32‑1 (MQ=255)
                                atatagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGcaaaaca                   <  1:1160993‑M2/33‑1 (MQ=255)
                                 tatagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGcaaaacaa                  <  1:3983249‑M2/34‑1 (MQ=255)
                                  atagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGcaaaacaa                  >  2:2964237‑M2/16‑49 (MQ=255)
                                  atagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGcaaaacaa                  >  2:3053963‑M2/16‑49 (MQ=255)
                                  atagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGcaaaacaa                  >  2:2435986‑M2/16‑49 (MQ=255)
                                   tagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGcaaaacaaa                 <  2:3949857‑M2/35‑1 (MQ=255)
                                   tagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGcaaaacaaa                 <  2:2526057‑M2/35‑1 (MQ=255)
                                      tttAATTTCATTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGt              <  1:90906/49‑1 (MQ=255)
                                      aattcaaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTc             >  2:665349‑M2/12‑50 (MQ=255)
                                       ttAATTTCATTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTcc            >  1:1941815/1‑50 (MQ=255)
                                       attcaaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTcc            >  1:3295637‑M2/11‑50 (MQ=255)
                                       attcaaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTcc            <  1:2060791‑M2/40‑1 (MQ=255)
                                        ttcaaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTcc            >  1:1405571‑M2/10‑49 (MQ=255)
                                        ttcaaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCt           <  1:2222676‑M2/41‑1 (MQ=255)
                                        ttcaaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCt           <  1:1816043‑M2/41‑1 (MQ=255)
                                         tcaaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTa          >  2:4889553‑M2/9‑50 (MQ=255)
                                          caaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTa          <  1:5263563‑M2/42‑1 (MQ=255)
                                          caaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTa          <  1:1458197‑M2/42‑1 (MQ=255)
                                          caaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTAc         >  1:4375126‑M2/8‑50 (MQ=255)
                                          aTTTCATTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTAc         <  1:4172751/50‑1 (MQ=255)
                                           aaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTAc         <  2:2117509‑M2/43‑1 (MQ=255)
                                           aaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTAc         <  1:348486‑M2/43‑1 (MQ=255)
                                           aaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTAc         <  2:1456757‑M2/43‑1 (MQ=255)
                                           aaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTAc         <  1:2174874‑M2/43‑1 (MQ=255)
                                           aaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTAc         <  1:3319889‑M2/43‑1 (MQ=255)
                                           aaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTAc         <  1:2813109‑M2/43‑1 (MQ=255)
                                             acatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTACg        >  2:5111636‑M2/5‑48 (MQ=255)
                                              catTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTACGCAg     <  2:1708199‑M2/47‑1 (MQ=255)
                                                tTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTACGCAgg    >  1:1044258‑M2/2‑49 (MQ=255)
                                                tTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTACGCAGGa   <  1:480116‑M2/49‑1 (MQ=255)
                                                 tataATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTACGCAGGaa  >  1:4062597/1‑50 (MQ=255)
                                                 |                                                 
ACACCAGCCCCATACCAATCATAGATTTAGTAAATATATTTAATTTCATTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTACGCAGGAA  >  NC_000913/3469805‑3469903

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