New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4035244	20 (0.230)	7 (0.090)	5/168	NT	27.1%	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

TGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCG  >  NC_000913/4035146‑4035343                                                                                                   |                                                                                                    tGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGaca                                                                                                    <  1:824154/100‑1 (MQ=255) tGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGaca                                                                                                    >  2:1082312/1‑100 (MQ=255)  gAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCGCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGacac                                                                                                   <  2:20156/100‑1 (MQ=255)   aaaTAAGCATAAAGAATAAACAATGCGCGTTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGAC                                                                                                     <  2:1180398‑M1/97‑1 (MQ=255)    aaTAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACgg                                                                                                 >  2:1587904/1‑100 (MQ=255)     aTAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGa                                                                                                <  2:763515/100‑1 (MQ=255)      tAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaa                                                                                               <  2:2153932/100‑1 (MQ=255)      tAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTACCTATAATGCGCCACCACTGACACGGaa                                                                                               >  2:578554/1‑100 (MQ=255)       aaGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaac                                                                                              <  2:309520/100‑1 (MQ=255)        aGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaaca                                                                                             <  2:447704/100‑1 (MQ=255)        aGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaaca                                                                                             >  1:1879485/1‑100 (MQ=255)        aGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCAACCACTGACACGGaac                                                                                              >  2:1416676/1‑99 (MQ=255)          cATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACccggcggcgt                                                                                           >  2:1612222‑M1/1‑90 (MQ=255)          cATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaacaac                                                                                           >  1:2320833/1‑100 (MQ=40)          cATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaacaac                                                                                           >  2:113521/1‑100 (MQ=40)           aTAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGTAACAACg                                                                                          >  2:1186238/1‑100 (MQ=255)           aTAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACg                                                                                          <  2:2360042/100‑1 (MQ=40)            tAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACgg                                                                                         <  2:1290000/100‑1 (MQ=38)             aaaGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGc                                                                                        <  1:1924297/100‑1 (MQ=38)             aaaGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGc                                                                                        >  2:97080/1‑100 (MQ=38)              aaGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCa                                                                                       >  1:1349184/1‑100 (MQ=35)              aaGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCa                                                                                       >  2:260003/1‑100 (MQ=35)                 aaTAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAac                                                                                    <  1:1919769/100‑1 (MQ=32)                 aaTAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAac                                                                                    <  2:183494/100‑1 (MQ=32)                 aaTAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAac                                                                                    <  1:2295975/100‑1 (MQ=32)                  aTAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAaca                                                                                   >  1:113307/1‑100 (MQ=31)                   tAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAacac                                                                                  <  1:2219419/100‑1 (MQ=11)                                 tCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACccggcggcgtgtttgccgttgttccgtgtcagt                                                                    <  1:160951‑M1/100‑34 (MQ=255)                                  cAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACccggcggcgtgtttgccgttgttccgtgtcagcg                                                                   <  1:1612234‑M1/100‑35 (MQ=255)                                   aGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACccggcggcgtgtttgccgttgttccgtgtcagtgg                                                                  <  1:1510142‑M1/100‑36 (MQ=255)                                      aaaTTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACccggcggcgtgtttgcctttgttccgtgtcagtggtgg                                                               <  1:2245051‑M1/100‑39 (MQ=255)                                        attattTTAAATTTCCTCTTGTCAGGCCGTAATAACTCCCTATAATGCGCCACCACTGACccggcggcgtgtttgccgttgttccgtgtcagtggtggcg                                                             <  1:1990766‑M1/100‑41 (MQ=255)                                        attattTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACccggcggcgtgtttgccgttgttccgtgtcagtggtggcg                                                             <  1:309413‑M1/100‑41 (MQ=255)                                                                                              gatggcccagcggcgtgtttgccgttgtttcgtgttagtggtggcgaattatagggagttattccggcctgacaagaggaaatttaaaataattttctg        <  2:1608770‑M1/94‑94 (MQ=255)                                                                                               cTGCCccggcggcgtgtttgccgttgtttcgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgac      <  1:1735702‑M1/100‑96 (MQ=255)                                                                                               cTGACccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgac      <  1:786438‑M1/100‑96 (MQ=255)                                                                                                tGACccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgacc     <  1:479403‑M1/100‑97 (MQ=255)                                                                                                 gACccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccg    <  1:348795‑M1/100‑98 (MQ=255)                                                                                                 gACccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccg    >  1:427770‑M1/1‑3 (MQ=255)                                                                                                 gACccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccg    >  1:703838‑M1/1‑3 (MQ=255)                                                                                                 gACccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccg    >  2:937219‑M1/1‑3 (MQ=255)                                                                                                 gACccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccg    >  1:862440‑M1/1‑3 (MQ=255)                                                                                                   cccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcg  >  1:1760287‑M1/1‑1 (MQ=255)                                                                                                   cccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcg  >  1:70394‑M1/1‑1 (MQ=255)                                                                                                   cccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcg  <  2:1041122‑M1/100‑100 (MQ=255)                                                                                                   cccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcg  >  2:1005905‑M1/1‑1 (MQ=255)                                                                                                   |                                                                                                    TGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCG  >  NC_000913/4035146‑4035343

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