New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3428910 =	20 (0.200)	9 (0.100)	5/168	NT	25.4%	intergenic (‑148/‑326)	rrsD/yrdA	16S ribosomal RNA of rrnD operon/bacterial transferase hexapeptide domain protein
	?	NC_000913	3428937 =	34 (0.360)					intergenic (‑175/‑299)	rrsD/yrdA	16S ribosomal RNA of rrnD operon/bacterial transferase hexapeptide domain protein

ACTCTTGAGACTTGGTATTCATTTTTCGTCTTGCGACGTTAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTG  >  NC_000913/3428811‑3429006                                                                                                    |                                                                                                 acgtgaaacacttcacaggatggtcggaacaacgaagagaaaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcg                                                                                                  <  2:357668‑M1/1‑1 (MQ=255)   tgttcagaCTTGGTATTCATTGTCCGTCTTGCGACGTTAAGAACCCGTAGCTTCGATTGCCCACACAGATTGTCTTATAAATTGTTAAAGAGCAGTTGCg                                                                                                <  2:919412/95‑1 (MQ=1)   gtgaaacacttcacaggatggtcggaacaacgaagagaaaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcgcg                                                                                                >  2:2542682‑M1/98‑100 (MQ=255)   gtgaaacacttcacaggatggtcggaacaacgaagagaaaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcgcg                                                                                                <  2:1776170‑M1/3‑1 (MQ=255)     gaaacacttcacaggatggtcggaacaacgaagagaaaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcGCGAc                                                                                              >  1:467807‑M1/96‑100 (MQ=255)                           gaacaacgaagagaaaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcGCGACGCGGCTTTCAGCTCACTGTcgc                                                                        >  1:457389‑M1/74‑100 (MQ=255)                            aacaacgaagagaaaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcGCGACGCGGCTTTCAGCTCACTGTcgcg                                                                       >  1:268928‑M1/73‑100 (MQ=255)                              tcggttgaagagaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcGCGACGCGGCTTTCAGCTCACTGTCGCGAg                                                                     >  1:1091258‑M1/71‑100 (MQ=255)                              tcggttgaagagaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcGCGACGCGGCTTTCAGCTCACTGTCGCGAg                                                                     >  1:2045238‑M1/71‑100 (MQ=255)                              caacgaagagaaaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcGCGACGCGGCTTTCAGCTCACTGTCGCGAg                                                                     >  1:1591089‑M1/71‑100 (MQ=255)                               cggttgaagagaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcGCGACGCGGCTTTCAGCTCACTGTCGCGAgg                                                                    >  1:297694‑M1/70‑100 (MQ=255)                               aacgaagagaaaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcGCGACGCGGCTTTCAGCTCACTGTCGCGAgg                                                                    >  2:2057634‑M1/70‑100 (MQ=255)                               aacgaagagaaaaaaatcctgaaattcagggttgactctgaaagaggaaagcgtaatatacgccacctcGCGACGCGGCTTTCAGCTCACTGTCGCGAgg                                                                    >  1:2666461‑M1/70‑100 (MQ=255)                                                    tcgtCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGTCGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTCTcc                                                <  2:1537371/96‑1 (MQ=0)                                                                          gtttGATAAATTGTTAAAGAGCAGTTGCTAAGCAGCTGTCAGCTCACTGTCGCGAGGTCGCGTAAATTACGCTTGCCTCTTTCAGAGTCAACCCTGAAtt                         <  2:774099/97‑1 (MQ=1)                                                                             gaaagagtaatatacgccacctcGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCa                      >  1:1283744‑M1/24‑100 (MQ=255)                                                                                      gTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGAtttttt             <  2:1993849/100‑1 (MQ=11)                                                                                      gTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGAtttttt             >  2:287933/1‑100 (MQ=11)                                                                                      gTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGAtttttt             <  2:350652/100‑1 (MQ=11)                                                                                      gTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGAtttttt             >  2:980171/1‑100 (MQ=11)                                                                                      gTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGAtttttt             <  1:520225/100‑1 (MQ=11)                                                                                       ttAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGAttttttt            >  1:1460523/1‑100 (MQ=31)                                                                                       ttAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGAttttttt            >  1:2247179/1‑100 (MQ=31)                                                                                        tAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTtc           <  1:766663/100‑1 (MQ=32)                                                                                        tAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTtc           >  1:1496907/1‑100 (MQ=32)                                                                                         aaaGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTtct          >  1:2674166/1‑100 (MQ=33)                                                                                         aaaGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTtct          <  2:1054102/100‑1 (MQ=32)                                                                                          aaGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTtctc         <  2:1102382/100‑1 (MQ=33)                                                                                          aaGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTtctc         <  1:2227705/100‑1 (MQ=33)                                                                                          aaGAGCAGTTGCGACGCGGCTTTCAGCTCAATGTCGCGAGGTGGCGTATCTTACGCTTTCCTCTTTCAGAGTCAACCATGAATTTCAGGATTTTTTtctc         >  1:1923292/1‑100 (MQ=11)                                                                                            gagCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCtt       <  2:1062024/100‑1 (MQ=35)                                                                                             agCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGCGTCAACCCTGAATTTCAGGATTTTTTTCTCTTc      >  1:2365893/1‑100 (MQ=255)                                                                                              gCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACTCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCg     <  2:90873/100‑1 (MQ=255)                                                                                              gCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCg     <  1:2014893/100‑1 (MQ=38)                                                                                               cAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCgt    >  1:533454/1‑100 (MQ=40)                                                                                               cAGTTCCGACGCGGCTTTCAGCTTTCTGTCGCGAGGTTGCGTACATTACTCTTTCCTCTTTCAGAGTCAACCCTGAATTTCATGATTTTTTTCTCTTCg     <  2:409997/99‑1 (MQ=14)                                                                                                aGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCgtt   <  1:2423319/100‑1 (MQ=255)                                                                                                 gTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCgttg  <  2:1324517/100‑1 (MQ=255)                                                                                                 gTTGCGACGCGGCATTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCgttg  >  1:1848207/1‑100 (MQ=255)                                                                                                    |                                                                                                 ACTCTTGAGACTTGGTATTCATTTTTCGTCTTGCGACGTTAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTG  >  NC_000913/3428811‑3429006

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.