New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4035244	18 (0.210)	3 (0.040)	3/168	NT	15.0%	intergenic (+91/‑287)	hemG/rrsA	protoporphyrin oxidase, flavoprotein/16S ribosomal RNA of rrnA operon
	?	NC_000913	= 4035271	NA (NA)					intergenic (+118/‑260)	hemG/rrsA	protoporphyrin oxidase, flavoprotein/16S ribosomal RNA of rrnA operon

CTGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCG  >  NC_000913/4035145‑4035343                                                                                                    |                                                                                                    cTGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGAC                                                                                                     >  2:1491498‑M1/1‑100 (MQ=255) cTGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGAC                                                                                                     >  1:1364244‑M1/1‑100 (MQ=255) cTGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGAC                                                                                                     >  2:1536183‑M1/1‑100 (MQ=255) cTGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATGATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGAC                                                                                                     <  2:1653148‑M1/100‑1 (MQ=255)  tGAAATAAGCATAAAGAATAAAAAATTCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGaca                                                                                                    <  2:1114897/100‑1 (MQ=255)  tGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGaca                                                                                                    >  2:108646/1‑100 (MQ=255)   gAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGacac                                                                                                   >  2:568091/1‑100 (MQ=255)     aaTAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACgg                                                                                                 <  2:1979867/100‑1 (MQ=255)      aTAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTCACACGGa                                                                                                <  1:1237701/100‑1 (MQ=255)       tAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaa                                                                                               <  1:111061/100‑1 (MQ=255)         aGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaaca                                                                                             >  1:547212/1‑100 (MQ=255)         aGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaaca                                                                                             >  1:123475/1‑100 (MQ=255)           cATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaacaac                                                                                           >  2:655882/1‑100 (MQ=40)           cATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaacaac                                                                                           >  2:1972278/1‑100 (MQ=40)            aTAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACg                                                                                          <  1:800685/100‑1 (MQ=40)             tAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACgg                                                                                         >  2:632973/1‑100 (MQ=38)              aaaGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGc                                                                                        <  1:644904/100‑1 (MQ=38)              aaaGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGc                                                                                        >  2:963899/1‑100 (MQ=38)               aaGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCa                                                                                       <  2:1966574/100‑1 (MQ=35)               aaGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCa                                                                                       <  2:1031757/100‑1 (MQ=35)               aaGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCa                                                                                       >  1:452747/1‑100 (MQ=35)               aaGAATAAAAAATGCGCGGTCACAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACgg                                                                                         <  2:1943593/98‑1 (MQ=25)                aGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCTCCACCACTGACACGGAACAACGGCaa                                                                                      <  1:1437857/100‑1 (MQ=25)                 gAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCaaa                                                                                     <  1:1622156/100‑1 (MQ=33)                  aaTAAAAAATGCGTGGTCAAAAAATTTTTTTAATTTTCCTCACGTAAGTCCGGAATATGACACTATTGTGCGCCACCACTGACACGGAACAACGGCAAac                                                                                    <  2:2351347/100‑1 (MQ=2)                  aaTAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAac                                                                                    <  2:212742/100‑1 (MQ=32)                   aTAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAaca                                                                                   >  2:2047795/1‑100 (MQ=31)                                 gTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACccggcggcgtgtttgccgttgttccgtgtcag                                                                     <  1:1082880‑M1/100‑33 (MQ=255)                                  tCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACccggcggcgtgtttgccgttgttccgtgtcagt                                                                    <  2:2086823‑M1/100‑34 (MQ=255)                                       aaaTTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACccggcggcgtgtttgccgttgttccgtgtcagtggtgg                                                               <  2:979968‑M1/100‑39 (MQ=255)                                                                                                cTGACccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgac      <  2:1353826‑M1/100‑96 (MQ=255)                                                                                                 tGACccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgacc     <  2:137617‑M1/100‑97 (MQ=255)                                                                                                 tGACccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgacc     <  2:137520‑M1/100‑97 (MQ=255)                                                                                                  gACccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccg    <  2:1115873‑M1/100‑98 (MQ=255)                                                                                                  gACccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccg    >  1:1205541‑M1/1‑3 (MQ=255)                                                                                                    cccggcggcgtgtttgccgttgttccgtgtccgtggtggcgccttatagggagttattccggcctgacacgaggaaatttaaaataattttctgacagcg  >  2:883718‑M1/1‑1 (MQ=255)                                                                                                    cccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcg  <  1:770520‑M1/100‑100 (MQ=255)                                                                                                    cccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcg  >  1:518451‑M1/1‑1 (MQ=255)                                                                                                    cccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcg  <  2:913489‑M1/100‑100 (MQ=255)                                                                                                    cccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcg  >  1:1187485‑M1/1‑1 (MQ=255)                                                                                                    |                                                                                                    CTGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCG  >  NC_000913/4035145‑4035343

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