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

GTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAGCGTGCGATTAAAAAACATTAGTAAACTGGAGAATATGGAGTTCTGCGCTACAAAACAACGG  >  NC_000913/1529982‑1530275                                                                                                                                                     |                                                                                                                                                  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                                                                                                  >  4:281202/1‑145 (MQ=255)                                                      tAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGat                                                                                                  <  3:281202/145‑1 (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)                                                                                                    ccattatacttataaaatgaccattccctctactaattttttggtcattTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCaa                                                                              <  4:329136‑M2/70‑1 (MQ=255)                                                                                                    ccattatacttataaaatgaccattccctctactaattttttggtcattTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCaa                                                                              >  3:329136‑M2/50‑119 (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:244922/149‑1 (MQ=255)                                                                                                         taAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAGCGTGCGATTAAAAAACATTa                                           <  1:649094/149‑1 (MQ=255)                                                                                                              taattaatGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATAtt                                                                       <  3:306299/116‑1 (MQ=255)                                                                                                              taattaatGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATAtt                                                                       >  4:306299/1‑116 (MQ=255)                                                                                                                   aatgaccattccctctactaattttttggtcattTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAGCGTGCGATTAAAAAACATTAGTAAACTGga                                 >  1:338892‑M2/35‑149 (MQ=255)                                                                                                                   AATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAGCGTGCGATTAAAAAACATTAGTAAACTGga                                 >  3:500453‑M1/1‑149 (MQ=255)                                                                                                                   AATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTCATATAAATCTTAATTTTTTCAATGATATTCTTTATAGCGTGCGATTAAAAAACATTAGTAAACTGga                                 >  4:48401‑M1/1‑149 (MQ=255)                                                                                                                    ATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAGCGTGCGATTAAAAAACATTAGTAAACTGgag                                <  2:492094‑M1/149‑1 (MQ=255)                                                                                                                       ACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAGCGTGCGATTAAAAAACATTAGTAAACTGGAGAat                             <  4:280277‑M1/149‑1 (MQ=255)                                                                                                                          AAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAGCGTGCGATTAAAAAACATTAGTAAACTGGAGAATATg                          >  3:568167‑M1/1‑149 (MQ=255)                                                                                                                             aaaTTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAGCGTGCGATTAAAAAACATTAGTAAACTGGAGAATATGGAg                       <  1:536173/149‑1 (MQ=255)                                                                                                                             aaaTTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAGCGTGCGATTAAAAAACATTAGTAAACTGGAGAATATGGAg                       <  1:119809/149‑1 (MQ=255)                                                                                                                                         ggAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAGCGTGCGATTAAAAAACATTAGTAAACTGGAGAATATGGAGTTCTGCGCTACa           >  2:214172/1‑149 (MQ=255)                                                                                                                                         ggAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAGCGTGCGATTAAAAAACATTAGTAAACTGGAGAATATGGAGTTCTGCGCTACa           >  4:95761/1‑149 (MQ=255)                                                                                                                                                  attTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAGCGTGCGATTAAAAAACATTAGTAAACTGGAGAATATGGAGTTCTGCGCTACAAAACAACgg  <  2:523715‑M2/146‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAGCGTGCGATTAAAAAACATTAGTAAACTGGAGAATATGGAGTTCTGCGCTACAAAACAACGG  >  NC_000913/1529982‑1530275

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.