New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1530096 =	36 (0.470)	3 (0.040)	3/248	NT	8.2%	coding (175/483 nt)	ydcD	putative immunity protein for RhsE
	?	NC_000913	1530130 =	33 (0.450)					coding (209/483 nt)	ydcD	putative immunity protein for RhsE

TCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAGCGTGCGATTAAAAAACATTAGTAAACTGGA  >  NC_000913/1529950‑1530244                                                                                                                                                   |                                                                                                                                                     ccgttgttttgtagcgcagaactccatattctccagtttactaatgttttttaatcgcacgctataaagaatatcattgaaaaaattaagatttatatcaaccactccacagtgcgtcctgtaattattaaccattatacttataaAAt                                                                                                                                                    <  2:523715‑M1/3‑1 (MQ=255)   aTTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGa                                                                                                                                                  <  2:572628/149‑1 (MQ=255)   aTTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGa                                                                                                                                                  <  2:140224/149‑1 (MQ=255)           aaTTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCaaaaaa                                                                                                                                          <  3:452364/149‑1 (MQ=255)                  cTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagtag                                                                                                                                   >  4:28015/1‑149 (MQ=255)                   tttGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTaga                                                                                                                                  >  1:616486/1‑149 (MQ=255)                    ttGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTagag                                                                                                                                 <  1:131227/149‑1 (MQ=255)                        ttATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGaa                                                                                                                             <  1:122795/149‑1 (MQ=255)                                tccagtttactaatgttttttaatcgcacgctataaagaatatcattgaaaaaattaagatttatatcaaccactccacagtgcgtcctgtaattattaaccattatacttataaAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt                                                                                                                     >  1:338892‑M1/116‑149 (MQ=255)                                 gttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAttt                                                                                                                    >  2:649094/1‑149 (MQ=255)                                   ttACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTta                                                                                                                  <  1:673144/149‑1 (MQ=255)                                      ccTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAGAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATaa                                                                                                               <  2:503072/149‑1 (MQ=255)                                       cTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                                                                                                              >  2:235735/1‑149 (MQ=255)                                        ttGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGt                                                                                                             <  4:136796/149‑1 (MQ=255)                                             tttGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAAt                                                                                                        <  3:374390/149‑1 (MQ=255)                                               tGATAACAATTCAAATTCTGACAACGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATgg                                                                                                      <  1:219685/149‑1 (MQ=255)                                                gATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGt                                                                                                     <  2:225477/149‑1 (MQ=255)                                                       aTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGt                                                                                                     >  3:526631/1‑142 (MQ=255)                                                       aTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGt                                                                                                     <  4:526631/142‑1 (MQ=255)                                                                 tGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGc                                                                                    >  2:337722/1‑149 (MQ=255)                                                                      tCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACtgt                                                                               >  2:494775/1‑149 (MQ=255)                                                                      tCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACtgt                                                                               >  2:252787/1‑149 (MQ=255)                                                                       cGTATTTTATGGACATAAAACACCAAAGAGCTTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACtgtg                                                                              <  3:582504/149‑1 (MQ=255)                                                                       cGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACtgtg                                                                              <  4:247876/149‑1 (MQ=255)                                                                          aTTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAg                                                                           <  4:307323/149‑1 (MQ=255)                                                                             ttgaaaaaattaagatttatatcaaccactccacagtgcgtcctgtaattattaaccattatacttataaAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATgg                                                                                                      <  4:329136‑M1/49‑1 (MQ=255)                                                                             ttgaaaaaattaagatttatatcaaccactccacagtgcgtcctgtaattattaaccattatacttataaAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATgg                                                                                                      >  3:329136‑M1/71‑119 (MQ=255)                                                                              tATGGACATAAAACACCAAAGATCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGt                                                                       <  4:561275/149‑1 (MQ=255)                                                                              tATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGt                                                                       <  4:612344/149‑1 (MQ=255)                                                                                tGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTg                                                                     >  4:622235/1‑149 (MQ=255)                                                                                      tAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGat                                                                   <  3:281202/145‑1 (MQ=255)                                                                                      tAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGat                                                                   >  4:281202/1‑145 (MQ=255)                                                                                       aaaaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATaaa                                                              <  3:622235/149‑1 (MQ=255)                                                                                           cacCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCtt                                                          >  3:48401/1‑149 (MQ=255)                                                                                                               tatCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCt                                      <  1:179013/149‑1 (MQ=255)                                                                                                                           aaaaaTATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAGCGTGCg                          >  4:133738/1‑149 (MQ=255)                                                                                                                                 attattTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAGCGTGCGATTaaa                    <  1:193536/149‑1 (MQ=255)                                                                                                                                    atttatAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAGCGTGCGATTaaaaaa                 <  2:606919/149‑1 (MQ=255)                                                                                                                                        ataAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAGCGTGCGATTAAAAAACAtt             >  4:588780/1‑149 (MQ=255)                                                                                                                                        ataAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATAT‑AATCTTAATTTTTTCAATGATATTCTTTATAGCGTGCGATTAAAAAACATTa            <  2:616486/149‑1 (MQ=255)                                                                                                                                         taAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAGCGTGCGATTAAAAAACATTa            <  1:649094/149‑1 (MQ=255)                                                                                                                                         taAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAGCGTGCGATTAAAAAACATTa            <  1:244922/149‑1 (MQ=255)                                                                                                                                              taattaatGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATAtt                                        <  3:306299/116‑1 (MQ=255)                                                                                                                                              taattaatGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATAtt                                        >  4:306299/1‑116 (MQ=255)                                                                                                                                                   AATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAGCGTGCGATTAAAAAACATTAGTAAACTGga  >  3:500453‑M1/1‑149 (MQ=255)                                                                                                                                                   AATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTCATATAAATCTTAATTTTTTCAATGATATTCTTTATAGCGTGCGATTAAAAAACATTAGTAAACTGga  >  4:48401‑M1/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     TCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAGCGTGCGATTAAAAAACATTAGTAAACTGGA  >  NC_000913/1529950‑1530244

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