New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3620675 =	NA (NA)	3 (0.040)	3/248	NT	100%	coding (1484/4236 nt)	rhsB	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
	?	NC_000913	3620695 =	0 (0.000)					coding (1504/4236 nt)	rhsB	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

TAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAG  >  NC_000913/3620558‑3620735                                                                                                                      |                                                             gtaaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATc                               >  1:324549‑M2/25‑149 (MQ=255)  taaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCa                              >  1:542207‑M2/24‑149 (MQ=255)   aaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCAc                             >  1:217656‑M2/23‑149 (MQ=255)     aacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCAc                             <  4:234861‑M2/127‑1 (MQ=255)     aacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCAc                             >  3:234861‑M2/21‑147 (MQ=255)     aacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGTCCTGACGGGCTGGAATTGCGCCGGGAATATGCTGAATTGGGACGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCAccc                           >  3:580970‑M2/21‑149 (MQ=255)     aacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCAccc                           >  1:384815‑M2/21‑149 (MQ=255)     aacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCAccc                           >  1:178502‑M2/21‑149 (MQ=255)     aacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTAAGGAAACGGCCCCTGACGGCGATATAAccc                           >  4:276711‑M2/21‑149 (MQ=255)      acgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCg                          >  3:645687‑M2/20‑149 (MQ=255)      acgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCg                          >  2:645235‑M2/20‑149 (MQ=255)                cctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTAt                <  1:645235‑M2/140‑1 (MQ=255)                         CGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCca       <  3:655449‑M2/149‑1 (MQ=255)                         CGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCca       >  3:611036‑M2/1‑149 (MQ=255)                          GGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCcac      >  1:488410‑M2/1‑149 (MQ=255)                          GGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCcac      <  2:246641‑M2/149‑1 (MQ=255)                           GCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATGTGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCGACCGTTATGATAATCcaca     >  4:72670‑M2/1‑149 (MQ=255)                           GCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCg                                    >  2:95506‑M2/1‑118 (MQ=255)                           GCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCg                                    >  2:251498‑M2/1‑118 (MQ=255)                           GCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCg                                    <  1:95506‑M2/118‑1 (MQ=255)                           GCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCg                                    <  1:251498‑M2/118‑1 (MQ=255)                           GCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCcaca     >  3:571280‑M2/1‑149 (MQ=255)                           GCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCcaca     >  2:35663‑M2/1‑149 (MQ=255)                           GCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCcaca     >  3:686121‑M2/1‑149 (MQ=255)                           GCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCcaca     <  4:666032‑M2/149‑1 (MQ=255)                           GCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTCACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCctca     <  4:228170‑M2/149‑1 (MQ=255)                             GTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCcacaca   >  1:682992‑M2/1‑149 (MQ=255)                              TTTTACTATAACCACCACAACCAGTTAACGTCAGCCTCCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAg  >  4:652882‑M2/1‑149 (MQ=255)                              TTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGGGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGa                                         <  2:555174‑M2/110‑1 (MQ=255)                              TTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGGGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGa                                         >  1:555174‑M2/1‑110 (MQ=255)                              TTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAg  >  1:614204‑M2/1‑149 (MQ=255)                                                                                                                      |                                                             TAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAG  >  NC_000913/3620558‑3620735

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