New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 225685	86 (1.320)	12 (0.200)	10/186	6.9	12.9%	intergenic (+110/‑74)	alaV/rrlH	tRNA‑Ala/23S ribosomal RNA of rrnH operon
	?	NC_000913	3426842 =	NA (NA)					intergenic (‑59/+116)	rrlD/alaU	23S ribosomal RNA of rrnD operon/tRNA‑Ala
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below cutoff threshold.

CAGGGCATCCACCGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAAGAGCAAAACTTCGCAGTGAAC  >  NC_000913/3426749‑3426910                                                                                              |                                                                     cAGGGCATCCACCGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCgt                                                                >  1:1000888/1‑100 (MQ=31)   gggCATCCACCGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCgttg                                                              >  2:246646/1‑100 (MQ=33)   gggCATCCACCGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCgttg                                                              <  2:1349834/100‑1 (MQ=32)   gagCATCCACCGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCgttg                                                              >  1:1008242/3‑100 (MQ=17)      cATCCACCGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTc                                                           <  2:1110913/100‑1 (MQ=33)         ccACCGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAgt                                                        >  2:1585833/1‑100 (MQ=32)          cACCGTGTACGCTTAGTCGCTTAACCTCCCAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAgtg                                                       <  2:1331831/100‑1 (MQ=18)          cACCGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAgtg                                                       <  1:1104153/100‑1 (MQ=33)          cACCGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAgtg                                                       <  2:539219/100‑1 (MQ=33)            ccGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGtt                                                     >  2:56657/1‑100 (MQ=32)             cGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGttt                                                    >  1:618907/1‑100 (MQ=33)              gtgtACGCTTAGTCGCTTAACCTCACAACCCGACGATGTGTCACTTCGGTGTTGCGAACATTTGAGAGACTCACGACCATCTTTCGTTGTTCAGTGttt                                                    >  1:760222/1‑99 (MQ=11)              gtgtACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTc                                                   >  1:402227/1‑100 (MQ=32)                  aCGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAACATTTGAGAGACTCACTACCAACTTTCGTTGTTCAGCGCTTCAAtt                                               >  1:674446/1‑100 (MQ=11)                   cGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAAttt                                              >  2:581895/1‑100 (MQ=33)                     cTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGTACAACTTTCGTTGTTCACTGTTTCAATTTTc                                            <  2:1075031/100‑1 (MQ=16)                     cTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTc                                            <  2:341659/100‑1 (MQ=33)                       tAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAg                                          >  1:753211/1‑100 (MQ=32)                        aGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGc                                         >  2:1501407/1‑100 (MQ=33)                         gTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCt                                        >  2:274006/1‑100 (MQ=33)                           cGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTg                                      >  2:462300/1‑100 (MQ=33)                           cGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTg                                      >  2:1073142/1‑100 (MQ=33)                            gCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGa                                     >  2:1555550/1‑100 (MQ=33)                             cTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGAt                                    >  1:1522498/1‑100 (MQ=33)                              ttAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATc                                   >  1:75525/1‑100 (MQ=32)                                 aCCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAg                                <  2:1008252/100‑1 (MQ=33)                                 aCCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAg                                >  2:454991/1‑100 (MQ=33)                                   cTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGAt                              >  2:387241/1‑100 (MQ=32)                                    tCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGAtt                             >  1:70506/1‑100 (MQ=33)                                      acaACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGAtttt                           <  1:1169510/100‑1 (MQ=33)                                        gatctgAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTa                         <  2:772854/95‑1 (MQ=18)                                          cccGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTaaa                       <  1:697895/100‑1 (MQ=32)                                           ccGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAag                      <  1:1525578/100‑1 (MQ=33)                                             gAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTaagatc                    >  1:772847/1‑96 (MQ=18)                                             gAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAagag                    <  2:860355/100‑1 (MQ=33)                                             gAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAagag                    <  1:28563/100‑1 (MQ=33)                                              aaGATTTTTCACTTCTGAGTTGGGAAAATTTGTTTTTCTCACGAACCACTTCCGTTGTTCAGTGTTTCACTTTTCACCTTGATCCCGATTTTTAAaga                     <  1:273766/98‑1 (MQ=2)                                                   gTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAAGAGCAAAAc              <  2:1261800/100‑1 (MQ=35)                                                     ttCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAAGAGCAAAACtt            >  1:786425/1‑100 (MQ=38)                                                       cACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAAGAGCAAAACTTCg          >  1:811798/1‑100 (MQ=40)                                                               agTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAAGAGCAAAACTTCGCAGTGAAc  <  1:1006240/100‑1 (MQ=11)                                                               agTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAAGAGCAAAACTTCGCAGTGAAc  <  1:1006190/100‑1 (MQ=11)                                                                                              |                                                                     CAGGGCATCCACCGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAAGAGCAAAACTTCGCAGTGAAC  >  NC_000913/3426749‑3426910

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.