New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_pgi	730954 =	15 (0.870)	4 (0.240)	4/260	1.8	21.8%	coding (1372/4194 nt)	rhsC	Rhs protein with putative toxin domain; putative neighboring cell growth inhibitor
	?	NC_000913_3_pae_pgi	3620588 =	NA (NA)					coding (1397/4236 nt)	rhsB	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

GTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCG  >  NC_000913_3_pae_pgi/730813‑731102                                                                                                                                              |                                                                                                                                                     gtgAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGAt                                                                                                                                              >  2:186481/1‑150 (MQ=255)                     cGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGtggtg                                                                                                                          <  3:50167/150‑5 (MQ=255)                      ggCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTAc                                                                                                                         <  4:77979/150‑1 (MQ=255)                                       tCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCa                                                                                                        <  1:186481/150‑1 (MQ=255)                                         aGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGt                                                                                                      >  3:13362/1‑150 (MQ=255)                                              gACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACg                                                                                                 <  4:116514/150‑1 (MQ=255)                                                                     tcatattcccggcgcaattccagcccgtcaggcccggtggctgacgttaactggttgtggtggttatagtaaaACGCCGGATGGCAGGGCATCGGCGttt                                                                                                                            <  4:63710‑M1/27‑1 (MQ=255)                                                                     tcatattcccggcgcaattccagcccgtcaggcccggtggctgacgttaactggttgtggtggttatagtaaaACGCCGGATGGCAGGGCATCGGCGttt                                                                                                                            >  3:63710‑M1/74‑100 (MQ=255)                                                                      caGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTg                                                                         <  1:65399/150‑1 (MQ=255)                                                                                    aattccagcccgtcaggcccggtggctgacgttaactggttgtggtggttatagtaaaACGCCGGATGGCAGGGCATCGGCGtt                                                                                                                             <  1:47480‑M1/26‑1 (MQ=255)                                                                                    aattccagcccgtcaggcccggtggctgacgttaactggttgtggtggttatagtaaaACGCCGGATGGCAGGGCATCGGCGtt                                                                                                                             >  2:47480‑M1/59‑84 (MQ=255)                                                                                              acaGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGaa                                                 >  4:90091/1‑150 (MQ=255)                                                                                                  aaTACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAATggg                                             >  2:177874/1‑150 (MQ=255)                                                                                                       aGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGt                                        >  3:16229/1‑150 (MQ=255)                                                                                                           cGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGa                                    >  4:11019/1‑150 (MQ=255)                                                                                                                                          ccccacGCCGGATGGCAGGGCCTCGGCGTTTTACTATACCCACCACAGCCAGTTAACGTCGGCCACCGGGCCTGAGGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCAc     <  4:13362/148‑1 (MQ=255)                                                                                                                                           accacGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCAcc    >  1:183180/1‑150 (MQ=255)                                                                                                                                            ccacGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCAccc   <  3:11019/150‑1 (MQ=255)                                                                                                                                             cacGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCg  >  3:86125/1‑150 (MQ=255)                                                                                                                                              |                                                                                                                                                     GTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCG  >  NC_000913_3_pae_pgi/730813‑731102

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