New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 38884311 (0.010)32 (0.390) 22/96 2.1 98.5% intergenic (+239/‑4) mnmE/tnaC 5‑carboxymethylaminomethyluridine‑tRNA synthase GTPase subunit/tnaAB operon leader peptide
?NC_000913 3906599 = 0 (0.000)intergenic (‑32/+254) bglG/phoU transcriptional antiterminator BglG/negative regulator of the pho regulon

AAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCATTATGAATATCTTACATATATGTGTGACCTCAAAATGGTTCA  >  NC_000913/3888384‑3888474
                                               |                                           
aaaCTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCag                                           <  1:3707926/50‑2 (MQ=255)
    tctGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgactt                                       >  1:4493296‑M1/1‑44 (MQ=255)
    tctGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgactt                                       >  1:538666‑M1/1‑44 (MQ=255)
    tctGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgactt                                       >  1:4855047‑M1/1‑44 (MQ=255)
      tGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttca                                     >  1:3215640‑M1/1‑42 (MQ=255)
      tGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgactt                                       <  1:5276437‑M1/48‑7 (MQ=255)
        aGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcacc                                   <  1:3980679‑M1/50‑11 (MQ=255)
        aGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcacc                                   <  1:2802856‑M1/50‑11 (MQ=255)
        aGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcac                                    <  1:5144535‑M1/49‑10 (MQ=255)
        aGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcac                                    <  1:4438928‑M1/49‑10 (MQ=255)
        aGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcac                                    <  1:6021024‑M1/49‑10 (MQ=255)
        aGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttca                                     <  1:2642385‑M1/48‑9 (MQ=255)
           gtAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccag                                 >  1:5070692‑M1/1‑37 (MQ=255)
             aataatGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtat                              >  1:1541704‑M1/1‑35 (MQ=255)
             aataatGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtat                              >  1:6284988‑M1/1‑35 (MQ=255)
              ataatGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtatt                             >  1:5443886‑M1/1‑34 (MQ=255)
               taatGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattc                            >  1:1390341‑M1/1‑33 (MQ=255)
               taatGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattc                            >  1:4612308‑M1/1‑33 (MQ=255)
                 atGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctc                          >  1:3399547‑M1/1‑31 (MQ=255)
                 atGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattct                           <  1:4873928‑M1/49‑19 (MQ=255)
                    tAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctgg                       >  1:4526682‑M1/1‑28 (MQ=255)
                       ccTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctggtta                    >  1:6357883‑M1/1‑25 (MQ=255)
                        cTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttat                   >  1:4754394‑M1/1‑24 (MQ=255)
                           gtgtCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttatgt                 <  1:3030215‑M1/49‑29 (MQ=255)
                           gtgtCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttatgt                 <  1:5239942‑M1/49‑29 (MQ=255)
                           gtgtCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttatg                  <  1:2556130‑M1/48‑28 (MQ=255)
                              tCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttatgtcagg             <  1:3279115‑M1/50‑33 (MQ=255)
                                    gAGGATAAGTGCcgacttcaccagtattctctggttatgtcaggttttgc       >  1:2366512‑M1/1‑12 (MQ=255)
                                      ggATAAGTGCcgacttcaccagtattctctggttatgtcaggttttgcct     <  1:239083‑M1/50‑41 (MQ=255)
                                        aTAAGTGCcgacttcaccagtattctctggttatgtcaggttttgcctgc   <  1:2336481‑M1/50‑43 (MQ=255)
                                        aTAAGTGCcgacttcaccagtattctctggttatgtcaggttttgcctgc   <  1:853724‑M1/50‑43 (MQ=255)
                                         tAAGTGCcgacttcaccagtattctctggttatgtcaggttttgcctgcg  <  1:3181115‑M1/50‑44 (MQ=255)
                                         tAAGTGCcgacttcaccagtattctctggttatgtcaggttttgcctgcg  <  1:1036936‑M1/50‑44 (MQ=255)
                                               |                                           
AAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCATTATGAATATCTTACATATATGTGTGACCTCAAAATGGTTCA  >  NC_000913/3888384‑3888474

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