New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4255890 =	35 (0.860)	6 (0.160)	6/282	2.6	14.3%	coding (577/2424 nt)	plsB	glycerol‑3‑phosphate O‑acyltransferase
	?	NC_000913	4256077 =	39 (1.010)					coding (390/2424 nt)	plsB	glycerol‑3‑phosphate O‑acyltransferase
Rejected: Frequency below cutoff threshold.

CCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCG  >  NC_000913/4255935‑4256226                                                                                                                                               |                                                                                                                                                      ccGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACAt                                                                                                                                               <  1:426144/151‑1 (MQ=255)    ccATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCGCCTTTTTTACGCCCCCGCGCGCGACCCAACATCa                                                                                                                                             >  1:356431/1‑150 (MQ=255)     cATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCAcc                                                                                                                                           <  2:222745/151‑1 (MQ=255)           gcgcAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGacac                                                                                                                                      <  2:538656/150‑1 (MQ=255)               aGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGGCCAAACATCACCGACACTGGCa                                                                                                                                 >  1:356779/1‑151 (MQ=255)                gCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTTACCTCGCCTTTTTTACGCCCCCGCCCGCGCACACACATTACCCACACTgggc                                                                                                                                 >  1:155885/1‑148 (MQ=255)                 cGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCAc                                                                                                                                <  2:468403/150‑1 (MQ=255)                 cGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCAc                                                                                                                                <  1:640644/150‑1 (MQ=255)                   aaatttttcgctgtactgtggctcggtcgcgacagttttgtgcgtttctcgccgtcagtttcgctgcgccgtatggcggatgaacacggcacggataaaactatcgctcagaaactggcgcgcgACCAGACATCACCGACACTGGCAc                                                                                                                                <  2:509617‑M2/24‑1 (MQ=255)                    aatttttcgctgtactgtggctcggtcgcgacagttttgtgcgtttctcgccgtcagtttcgctgcgccgtatggcggaggaacacggcacggataaaactatcgctcagaaactggcgcgcgACCAAACATCACCGACACTGGCACCATc                                                                                                                            <  1:443870‑M2/28‑1 (MQ=255)                    aaCTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATc                                                                                                                            <  2:283488/151‑1 (MQ=255)                         aCGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGcac                                                                                                                       <  1:427207/151‑1 (MQ=255)                         aCGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGca                                                                                                                        <  2:155882/150‑1 (MQ=255)                                  aaCGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGAtt                                                                                                              <  2:358112/151‑1 (MQ=255)                                       acaAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGt                                                                                                         <  1:461839/151‑1 (MQ=255)                                             cTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTAc                                                                                                   <  1:635277/151‑1 (MQ=255)                                               gTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACCCTGGCACCATCTGCACATCCAGATTTGGGTTGCTAccg                                                                                                 >  1:258229/1‑149 (MQ=255)                                                      ttttgtgcgtttctcgccgtcagtttcgctgcgccgtatggcggatgaacacggcacggataaaactatcgctcagaaactggcgcgcgACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTg                                                                                               <  1:4920‑M2/57‑1 (MQ=255)                                                       cGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATc                                                                                         <  1:187248/151‑1 (MQ=255)                                                         aGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCg                                                                                        <  2:181064/150‑1 (MQ=255)                                                         aGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCCGCGCGCGGCCAAACATCACCCACACTGGGACCATCTGCACATCCCGATTTTGGGTGCTACGGTGCAAATCga                                                                                       >  2:533663/1‑151 (MQ=255)                                                         aGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCGCCTTTTTTACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGGTGCTACGGTGCAAATCga                                                                                       >  1:138897/1‑151 (MQ=255)                                                           ccACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGATTAATTTTTCTTTTTCTCGCCCCGGGGCGCCACCAAACATCACAGACACACACACCCTCTGCCCACCCACATTTTGGGTGCTCCGGAGAAAATCgaga                                                                                     >  2:483388/1‑151 (MQ=255)                                                              caGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGGGCAAATCGAGATAg                                                                                  >  1:140516/1‑151 (MQ=255)                                                              caGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCGCCTTTTTCACGCCCCCGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACCTCCAGATTTGGGTTTCTACGGTGCAAATAGAGATaa                                                                                  >  2:117251/1‑150 (MQ=255)                                                                gTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTTACGCCCCCGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCCGATTTTGGTTGCTACGGTGCAAATCGAGATAGt                                                                                 >  1:39567/1‑150 (MQ=255)                                                                gTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCCAACCTCACCGCCACTGGCACCATCTGCACATCCAGATTTGGGGTGCTACGGTGCAAATCGAGATAGt                                                                                 >  2:230866/1‑150 (MQ=255)                                                                       cgtcagtttcgctgcgccgtatggcggatgaacacggcacggataaaactatcgctcagaaactggcgcgcgACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAAc                                                                         <  2:534509‑M2/79‑1 (MQ=255)                                                                        aaaaaTTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAAc                                                                         <  1:315691/150‑1 (MQ=255)                                                                        aaaaaTTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAAc                                                                         <  2:579116/150‑1 (MQ=255)                                                                          cattttcgctgcgccgtagggcggatgaacccggcacggataaaactatcgctcagaaactggcgcgcgACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGc                                                                      <  2:308943‑M2/82‑1 (MQ=255)                                                                            aTTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGACACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCtt                                                                    <  1:227257/151‑1 (MQ=255)                                                                                 tGTACGCCGTTAAGCATTCGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAg                                                               <  2:176239/151‑1 (MQ=255)                                                                                 tGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAg                                                               <  1:84217/151‑1 (MQ=255)                                                                                 tGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAg                                                               <  1:356163/151‑1 (MQ=255)                                                                                 tGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAg                                                               <  1:183569/151‑1 (MQ=255)                                                                                 tGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAg                                                               <  1:141848/151‑1 (MQ=255)                                                                                 tGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAg                                                               <  1:133883/151‑1 (MQ=255)                                                                                  ctgcgccgtatggcggatgaacccggcacggataaaactatcgctcagaaactggcgcgcgACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAg                                                               <  1:268228‑M2/89‑1 (MQ=255)                                                                                       ccGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACtctt                                                         >  1:587433/1‑151 (MQ=255)                                                                                                        cggcggGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAg                                        <  1:456462/151‑1 (MQ=255)                                                                                                           cggGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTg                                     <  2:74116/151‑1 (MQ=255)                                                                                                                   cTTCGCCTTTTTCACGCCCCGGCGCGAGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACAcg                             <  1:545899/151‑1 (MQ=255)                                                                                                                        ccTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGcc                        <  2:153788/151‑1 (MQ=255)                                                                                                                              tCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGt                  <  2:592576/151‑1 (MQ=255)                                                                                                                                         cgcgcgACAAAACACCCCCGCCACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAAATAGTCGTGGAACAGCTTAATAGCCTCTTCTTCCGGCGTGTAATAGGTGACCACACGCGGCCCGCCGTGAATGAacac        <  2:356426‑M2/144‑1 (MQ=255)                                                                                                                                         cgcgcgACAAAACACCCCCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGAAAATCGAGATAGTCGGGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAacac        <  1:612671‑M2/144‑1 (MQ=255)                                                                                                                                           cgcgACCAAACATCCCCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACAt      <  1:102837‑M2/146‑1 (MQ=255)                                                                                                                                               aCCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAgcg  >  1:347527/1‑150 (MQ=255)                                                                                                                                               |                                                                                                                                                      CCGCCATACGGCGCAGCGAAACTGACGGCGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTTCACGCCCCGGCGCGCGACCAAACATCACCGACACTGGCACCATCTGCACATCCAGATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTTAATAGACTCTTCTTTCGGCGTGTAATAGGTGAACACACGCGGCCCGCCGTGAATGAACACATAGCG  >  NC_000913/4255935‑4256226

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