New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4035244	19 (0.180)	7 (0.070)	7/168	NT	28.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

CTGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCG  >  NC_000913/4035145‑4035343                                                                                                    |                                                                                                    cTGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGAC                                                                                                     >  2:1469762‑M1/1‑100 (MQ=255) cTGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGAC                                                                                                     >  1:373573‑M1/1‑100 (MQ=255) cTGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGAC                                                                                                     <  2:76014‑M1/100‑1 (MQ=255) cTGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGAC                                                                                                     >  1:185628‑M1/1‑100 (MQ=255)   gAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGacac                                                                                                   >  1:416448/1‑100 (MQ=255)   gAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGacac                                                                                                   <  1:866685/100‑1 (MQ=255)   gAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGacac                                                                                                   <  1:2295364/100‑1 (MQ=255)    aaaTAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACg                                                                                                  <  2:2346468/100‑1 (MQ=255)      aTAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGa                                                                                                <  1:1314668/100‑1 (MQ=255)        aaGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaac                                                                                              >  2:2381516/1‑100 (MQ=255)          gCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGAAACGGaacaa                                                                                            >  2:813993/1‑100 (MQ=255)           cATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaacaac                                                                                           >  1:482665/1‑100 (MQ=40)           cATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaacaac                                                                                           <  2:2110073/100‑1 (MQ=40)           cATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaacaac                                                                                           >  1:2787305/1‑100 (MQ=40)           cATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaacaac                                                                                           <  1:637269/100‑1 (MQ=40)            aTAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACg                                                                                          >  2:1586848/1‑100 (MQ=40)             tAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACgg                                                                                         <  2:885341/100‑1 (MQ=38)             tAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACg                                                                                          >  1:1375793/1‑99 (MQ=38)               aaGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCa                                                                                       <  2:1249992/100‑1 (MQ=35)               aaGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCa                                                                                       >  2:1641153/1‑100 (MQ=35)                aGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCaa                                                                                      <  1:867976/100‑1 (MQ=35)                  aaTAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAac                                                                                    <  2:2700916/100‑1 (MQ=32)                  aaTAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAac                                                                                    <  2:1952160/100‑1 (MQ=32)                  aaTAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAac                                                                                    <  1:1402123/100‑1 (MQ=32)                  aaTAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAac                                                                                    <  2:1948311/100‑1 (MQ=32)                   aTAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAaca                                                                                   <  1:458623/100‑1 (MQ=31)                    tAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAacac                                                                                  <  1:2308722/100‑1 (MQ=11)                    tAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAacac                                                                                  >  2:1766309/1‑100 (MQ=11)                               cgGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACccggcggcgtgtttgccgttgttccgtgtc                                                                       <  2:2596188‑M1/100‑31 (MQ=255)                                 gTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACccggcggcgtgtttgccgttgttccgtgtcag                                                                     <  1:660963‑M1/100‑33 (MQ=255)                                  tCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACccggcggcgtgtttgccgttgttccgtgtcagt                                                                    <  1:495963‑M1/100‑34 (MQ=255)                                     gAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACccggcggcgtgtttgccgttgttccgtgtcagtggt                                                                 <  1:767360‑M1/100‑37 (MQ=255)                                      aaaaTTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACccggcggcgtgtttgccgttgttccgtgtcagtggtg                                                                <  2:261071‑M1/100‑38 (MQ=255)                                         attattTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACccggcggcgtgtttgccgttgttccgtgtcagtggtggcg                                                             <  1:2190022‑M1/100‑41 (MQ=255)                                                  aaTTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagg                                                    <  2:913463‑M1/100‑50 (MQ=255)                                                                                                cTGACccggcggcgtatttgccgttgtaccgtgtcagtagtggcgcattatagagagttattccggcctgacaagaggaaatttaaaataattttctg        <  1:2536384‑M1/98‑94 (MQ=255)                                                                                                 ttgcccgggggcgtttttgccgtggttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaactttaaaataattttctgac      <  1:1696317‑M1/96‑96 (MQ=255)                                                                                                  gACccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggacatttaaaataattttctgaccg    >  1:2274174‑M1/1‑3 (MQ=255)                                                                                                  gACccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccg    >  1:962631‑M1/1‑3 (MQ=255)                                                                                                  gACccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccg    <  1:2449096‑M1/100‑98 (MQ=255)                                                                                                    cccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcg  >  1:2871655‑M1/1‑1 (MQ=255)                                                                                                    cccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcg  <  2:2195249‑M1/100‑100 (MQ=255)                                                                                                    cccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcg  <  2:2276111‑M1/100‑100 (MQ=255)                                                                                                    cccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcg  >  1:1717297‑M1/1‑1 (MQ=255)                                                                                                    cccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcg  >  1:151393‑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 < 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.