New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3428910 =16 (0.190)8 (0.100) 5/168 NT 34.6% intergenic (‑148/‑326) rrsD/yrdA 16S ribosomal RNA of rrnD operon/bacterial transferase hexapeptide domain protein
?NC_000913 = 4207974 NA (NA)intergenic (‑442/‑173) purH/rrsE IMP cyclohydrolase and phosphoribosylaminoimidazolecarboxamide formyltransferase/16S ribosomal RNA of rrnE operon

ACTCTTGAGACTTGGTATTCATTTTTCGTCTTGCGACGTTAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTT  >  NC_000913/3428811‑3429008
                                                                                                   |                                                                                                  
atgtgaatcacttcacacaaacagccggttcggttgaagagaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcg                                                                                                    >  1:1181063‑M1/100‑100 (MQ=255)
 tgtgaatcacttcacacaaacagccggttcggttgaagagaaaaatcctgaaattcagggttgcctctgaaagaggaaagcgtaatatacgccacctcgc                                                                                                   >  2:1034838‑M1/99‑100 (MQ=255)
  gtgaatcacttcacacaaacagccggttcggttgaagagaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcgcg                                                                                                  <  1:1163201‑M1/3‑1 (MQ=255)
  gtgaatcacttcacacaaacagccggttcggttgaagagaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcgcg                                                                                                  <  1:1910258‑M1/3‑1 (MQ=255)
  gtgaatcacttcacacaaacagccggttcggttgaagagaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcgcg                                                                                                  >  1:280961‑M1/98‑100 (MQ=255)
  gtgaatcacttcacacaaacagccggttcggttgaagagaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcg                                                                                                    >  2:331897‑M1/98‑98 (MQ=255)
   tgaatcacttcacacaaacagccggttcggttgaagagaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcGCGa                                                                                                 <  2:615741‑M1/4‑1 (MQ=255)
   tgaatcacttcacacaaacagccggttcggttgaagagaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcGCGa                                                                                                 <  1:545627‑M1/4‑1 (MQ=255)
    gaatcacttcacacaaacagccggttcggttgaagagaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcGCGa                                                                                                 <  2:2083077‑M1/4‑1 (MQ=255)
    gaatcacttcacacaaacagccggttcggttgaagagaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcGCGa                                                                                                 <  2:1419084‑M1/4‑1 (MQ=255)
    gaatcacttcacacaaacagccggttcggttgaagagaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcGCGAc                                                                                                <  2:1312312‑M1/5‑1 (MQ=255)
    gaatcacttcacacaaacagccggttcggttgaagagaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcGCGAc                                                                                                <  1:1747312‑M1/5‑1 (MQ=255)
     aatcacttcacacaaacagccggttcggttgaagagaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcGCGAc                                                                                                >  2:2381781‑M1/95‑99 (MQ=255)
                 ttaaTTTTTCGTCTTGCGACGTTAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCt                                                                                   <  1:514396/97‑1 (MQ=1)
                        ccggttcggttgaagagaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcGCGACGCGCTTTAGCGCACTGTcgc                                                                          >  1:866695‑M1/76‑100 (MQ=255)
                           gttcggttgaagagaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcGCGACGCGGCTTTCAGCTCACTGTcgcg                                                                         >  2:1641080‑M1/73‑100 (MQ=255)
                           gttcggttgaagagaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcGCGACGCGGCTTTCAGCTCACTGTcgcg                                                                         >  1:1638205‑M1/73‑100 (MQ=255)
                           aacaacgaagagaaaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcGCGACGCGGCTTTCAGCTCACTGTcgcg                                                                         >  2:2098798‑M1/73‑100 (MQ=255)
                             tcggttgaagagaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcGCGTCGCGGCTTTCAGCTCACTGTCGCGAg                                                                       >  1:1448245‑M1/71‑100 (MQ=255)
                             tcggttgaagagaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcGCGACGCGGCTTTCAGCTCACTGTCGCGAg                                                                       >  1:1157986‑M1/71‑100 (MQ=255)
                              cggttgaagagaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcGCGACGCGGCTTTCAGCTCACTGTCGCGAgg                                                                      >  1:139640‑M1/70‑100 (MQ=255)
                                                     ttagAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCctct                                               >  2:525376/4‑100 (MQ=1)
                                                            gttgactctgaaagaggaaagcgtaatatacgccacctcGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCagag                                        >  1:1750222‑M1/40‑100 (MQ=255)
                                                                                     gTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGCATTTCAGGAtttttt               >  1:2352138/1‑100 (MQ=11)
                                                                                     gTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGAtttttt               >  2:1288262/1‑100 (MQ=11)
                                                                                     gTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGAtttttt               <  2:1730071/100‑1 (MQ=11)
                                                                                     gTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGAtttttt               >  1:708137/1‑100 (MQ=11)
                                                                                      ttAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGAttttttt              >  1:977465/1‑100 (MQ=31)
                                                                                       tAAAGAGGAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGGGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTtc             <  2:1048367/100‑1 (MQ=12)
                                                                                       tAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTtc             <  2:1295806/100‑1 (MQ=32)
                                                                                       tAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTtc             <  1:2130954/100‑1 (MQ=32)
                                                                                       tAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTtc             >  2:210377/1‑100 (MQ=32)
                                                                                        aaaGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTtct            >  1:1831110/1‑100 (MQ=33)
                                                                                         aaGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTtctc           <  2:1346078/100‑1 (MQ=33)
                                                                                          agagCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTctct          >  1:285716/1‑100 (MQ=35)
                                                                                            agCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTc        <  1:1547515/100‑1 (MQ=38)
                                                                                             gCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCg       >  1:227875/1‑100 (MQ=40)
                                                                                              cAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCgt      <  2:166395/100‑1 (MQ=40)
                                                                                               aGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCgtt     >  2:985961/1‑100 (MQ=255)
                                                                                                gtggCGACGCGGCTTTGAGCCCAGGGGCTTGAGCTTGCGTATATTACGCTTTCCTCTTTCAGATTCAACCCGGTACTTCAGGATTTTTTTCTCGTCgtt     <  1:1666345/96‑1 (MQ=2)
                                                                                                gTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCgttg    <  2:1050800/100‑1 (MQ=255)
                                                                                                gTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCgttg    >  2:549843/1‑100 (MQ=255)
                                                                                                 ttGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCgttgt   <  1:1196364/100‑1 (MQ=255)
                                                                                                   GCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTTgttgtt  <  1:368538‑M1/99‑1 (MQ=255)
                                                                                                   |                                                                                                  
ACTCTTGAGACTTGGTATTCATTTTTCGTCTTGCGACGTTAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTT  >  NC_000913/3428811‑3429008

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.