New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 22362312 (0.120)7 (0.080) 5/168 NT 21.1% intergenic (+215/‑148) gmhB/rrsH D,D‑heptose 1,7‑bisphosphate phosphatase/16S ribosomal RNA of rrnH operon
?NC_000913 3428937 = 41 (0.450)intergenic (‑175/‑299) rrsD/yrdA 16S ribosomal RNA of rrnD operon/bacterial transferase hexapeptide domain protein

GGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGT  >  NC_000913/223524‑223722
                                                                                                   |                                                                                                   
ggttcggttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGC                                                                                                     >  2:2483429‑M1/1‑100 (MQ=255)
 gttcggttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCa                                                                                                    >  2:1978047/1‑100 (MQ=255)
    cggttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACt                                                                                                 <  1:2495952/100‑1 (MQ=40)
    cggttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACt                                                                                                 >  1:2366186/1‑100 (MQ=255)
     ggttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTg                                                                                                >  1:1270221/1‑100 (MQ=40)
     ggttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTg                                                                                                >  1:819558/1‑100 (MQ=40)
     ggttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTg                                                                                                >  1:984175/1‑100 (MQ=40)
        tGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGctc                                                                                             >  2:2068686/1‑100 (MQ=38)
        tGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGctc                                                                                             <  1:600659/100‑1 (MQ=35)
        aagagagAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGctc                                                                                             >  2:1976550/4‑100 (MQ=1)
         gAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGctct                                                                                            >  2:466323/1‑100 (MQ=35)
         gAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGctct                                                                                            >  2:948599/1‑100 (MQ=35)
           agagAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGTGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCttt                                                                                          >  2:1415024/1‑100 (MQ=17)
           agagAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCttt                                                                                          >  2:671861/1‑100 (MQ=32)
           agagAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCttt                                                                                          <  2:905572/100‑1 (MQ=33)
            gagaAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCTTTa                                                                                         <  1:432586/100‑1 (MQ=32)
            gagaAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCTTTa                                                                                         >  2:2321368/1‑100 (MQ=32)
              gaAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCTTTAAc                                                                                       <  1:1386088/100‑1 (MQ=11)
                                         ctGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCgaggtggcgtatattacgctttcctctttcagagtcaaccctg                                                          >  1:1215376‑M1/1‑57 (MQ=255)
                                                   aaaGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCgaggtggcgtatattacgctttcctctttcagagtcaaccctgaaattcag                                                  >  1:351084‑M1/1‑49 (MQ=255)
                                                                  ccACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCgaggtggcgtatattacgctttcctctttcagagtcaaccctgaatttcaggatttttttctcttc                                   >  1:2131265‑M1/1‑34 (MQ=255)
                                                                   cACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCgaggtggcgtatattacgctttcctctttcagagtcaaccctgaatttcaggatttttttctcttcg                                  >  1:1183991‑M1/1‑33 (MQ=255)
                                                                     ccTCGCGACAGTGAGCTGAAAGCCGCGTCGCgaggtggcgtatattacgctttcctctttcagagtcaaccctgaatttcaggatttttctcttcaaccg                                >  1:1396276‑M1/1‑31 (MQ=255)
                                                                         gcgACAGTGAGCTGAAAGCCGCGTCGCgaggtggcgtatattacgctttcctctttcagagtcaaccctgaatttcaggatttttttctcttcgttgttc                            >  1:2633760‑M1/1‑27 (MQ=255)
                                                                         gcgACAGTGAGCTGAAAGCCGCGTCGCgaggtggcgtatattacgctttcctctttcagagtcaaccctgaatttcaggatttttttctcttcgttgttc                            >  1:2636886‑M1/1‑27 (MQ=255)
                                                                                    cttaaaGCCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAAt                  <  2:2277630/96‑1 (MQ=1)
                                                                                               gTCGCgaggtggcgtatattacgctttcctctttcagagtcaaccctgaatttcaggatttttttctcttcgttgttccgaccatcctgtgaagtgtttc      >  1:2266393‑M1/1‑5 (MQ=255)
                                                                                               gTCGCgaggtggcgtatattacgctttcctctttcagagtcaaccctgaatttcaggatttttttctcttcgttgttccgaccatcctgtgaagtgtttc      <  2:1641691‑M1/100‑96 (MQ=255)
                                                                                               gTCGCgaggtggcgtatattacgctttcctctttcagagtcaaccctgaatttcaggatttttttctcttcgttgttccgaccatcctgtgaagtgtttc      <  1:336405‑M1/100‑96 (MQ=255)
                                                                                               gTCGCgaggtggcgtatattacgctttcctctttcagactcaaccctgaatttcaggatttttttctcttcgtttttccgaccatcctgagaagtgattc      >  1:1746444‑M1/1‑5 (MQ=255)
                                                                                                 cgcgaggtggcgtatattacgctttcctctttcagagtcaaccctgaatttcaggatttttttctcttcgttgttccgaccatcctgtgaagtgtttcac    >  1:175467‑M1/1‑3 (MQ=255)
                                                                                                 cgcgaggtggcgtatattacgctttcctctttcagagtcaaccctgaattacaggatttttatatcttcgttgttccgaccatcctgtgaagtgtttaa     >  1:1418873‑M1/1‑3 (MQ=255)
                                                                                                  gcgaggtggcgtatattacgctttcctctttcagagtcaaccctgaatttcaggatttttttctcttcgttgttccgaccatcctgtgaagtgtttcacg   <  1:1544948‑M1/100‑99 (MQ=255)
                                                                                                  gcgaggtggcgtatattacgctttcctctttcagagtcaaccctgaatttcaggatttttttctcttcgttgttccgaccatcctgtgaagtgtttcacg   >  2:814540‑M1/1‑2 (MQ=255)
                                                                                                  gcgaggtggcgtatattacgctttcctctttcagagtcaaccctgaatttcaggatttttttctcttcgttgttccaaccatcctgtgaagtgtttcacg   <  2:1183980‑M1/100‑99 (MQ=255)
                                                                                                   cgaggtggcgtatattacgctttcctctttcagagtcaaccctgaatttcaggatttttttctcttcgttgttccgaccatcctgtgaagtgtttcacgt  <  1:1696077‑M1/100‑100 (MQ=255)
                                                                                                   |                                                                                                   
GGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGT  >  NC_000913/223524‑223722

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