New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	526081 =	NA (NA)	11 (0.100)	9/166	NT	100%	coding (2821/4281 nt)	rhsD	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
	?	NC_000913	526087 =	0 (0.000)					coding (2827/4281 nt)	rhsD	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

TCCCGTATGCCACGGACCCGGCGGGCAACCGGCTGCCGGACCCGGAGCTGCACCCGGACAGTACACTCACAGTGTGGCCGGATAACCGCATCGCGGAGGATGCGCACTATGTCTACCGCCACGATGAATACGGCAGGCTGACGGAGAAGACGGACCGCATCCCGGCGGGTGTGATACGGACGGACGACGAGCGGACCC  >  NC_000913/525989‑526186                                                                                                   |                                                                                                    tcgtcgtccgtccgtatcacacccgccgggatgcggtccgtcttctccgtcagcctgccgtattcatcgtggcggtagacatagagcgcatcctccgcGa                                                                                                    >  2:3302492‑M2/99‑100 (MQ=255)  cgtcgtccgtccgtatcacacccgccgggatgcggtccgtcttctccgtcagcctgccgtattcatcgtggcggtagacatagtgcgcatcctccgcGAt                                                                                                   >  1:2596828‑M2/98‑100 (MQ=255)  cgtcgtccgtccgtatcacacccgccgggatgcggtccgtcttctccgtcagcctgccgtattcatcgtggcggtagacatagtgcgcatcctccgcGAt                                                                                                   >  2:117310‑M2/98‑100 (MQ=255)   gtcgtccgtccgtatcacacccgccgggatgcggtccgtcttctccgtcagcctgccgtattcatcgtggcggtagacatagtgcgcatcctccgcGATg                                                                                                  <  2:1783537‑M2/4‑1 (MQ=255)   gtcgtccgtccgtatcacacccgccgggatgcggtccgtcttctccgtcagcctgccgtattcatcgtggcggtagacatagtgcgcatcctccgcGATg                                                                                                  >  2:848107‑M2/97‑100 (MQ=255)   gtcgtccgtccgtatcacacccgccgggatgcggtccgtcttctccgtcagcctgccgtattcatcgtggcggtagacatagtgcgcatcctccgcGATg                                                                                                  >  2:847966‑M2/97‑100 (MQ=255)    tcgtccgtccgtatcacacccgccgggatgcggtccgtcttctccgtcagcctgccgtattcatcgtggcggtagacatagtgcgcatcctccgcGATgc                                                                                                 <  1:684384‑M2/5‑1 (MQ=255)    tcgtccgtccgtatcacacccgccgggatgcggtccgtcttctccgtcagcctgccgtattcatcgtggcggtagacatagtgcgcatcctccgcGATgc                                                                                                 <  2:2719771‑M2/5‑1 (MQ=255)     cgtccgtccgtatcacacccgccgggatgcggtccgtcttctccgtcagcctgccgtattcatcgtggcggtagacatagtgcgcatcctccgcGATgcg                                                                                                <  2:1289110‑M2/6‑1 (MQ=255)     cgtccgtccgtatcacacccgccgggatgcggtccgtcttctccgtcagcctgccgtattcatcgtggcggtagacatagtgcgcatcctccgcGATgcg                                                                                                >  1:1196733‑M2/95‑100 (MQ=255)     cgtccgtccgtatcacacccgccgggatgcggtccgtcttctccgtcagcctgccgtattcatcgtggcggtagacatagtgcgcatcctccgcGATgcg                                                                                                >  1:344584‑M2/95‑100 (MQ=255)     cgtccgtccgtatcacacccgccgggatgcggtccgtcttctccgtcagcctgccgtattcatcgtggcggtagacatagtgcgcatcctccgcGATgc                                                                                                 <  2:18013‑M2/5‑1 (MQ=255)                                                                 catcgtggcggtagacatagtgcgcatcctccgcGATGCGCACTATGTCTACCGCCACGATGAATACGGCAGGCTGACGGAGAAGACGGACCGCATCCcg                                    >  1:236921‑M2/35‑100 (MQ=255)                                                                          ggtagacatagtgcgcatcctccgcGATGCGCACTATGTCTACCGCCACGATGAATACGGCAGGCTGACGGAGAAGACGGACCGCATCCCGGCGGGtgtg                           <  2:236922‑M2/75‑1 (MQ=255)                                                                              gacatagtgcgcatcctccgcGATGCGCACTATGTCTACCGCCACGATGAATACGGCAGGCTGACGGAGAAGACGGACCGCATCCCGGCGGGTGTGATac                       >  1:2829313‑M2/22‑100 (MQ=255)                                                                                catagtgcgcatcctccgcGATGCGCACTATGTCTACCGCCACGATGAATACGGCAGGCTGACGGAGAAGACGGACCGCATCCCGGCGGGTGTGATacgg                     >  2:379641‑M2/20‑100 (MQ=255)                                                                                catagtgcgcatcctccgcGATGCGCACTATGTCTACCGCCACGATGAATACGGCAGGCTGACGGAGAAGACGGACCGCATCCCGGCGGGTGTGATacgg                     >  2:324205‑M2/20‑100 (MQ=255)                                                                                 atagtgcgcatcctccgcGATGCGCACTATGTCTACCGCCACGATGAATACGGCAGGCTGACGGAGAAGACGGACCGCATCCCGGCGGGTGTGATacgga                    >  2:2401944‑M2/19‑100 (MQ=255)                                                                                 atagtgcccatcctccgcGATGCGCACTATGTCTACCGCCACGATGAATACGGCAGGCTGACGGAGAAGACGGACCGCATCCCGGCGGGTGTGATacgga                    >  2:1760981‑M2/19‑100 (MQ=255)                                                                                    gtgcgcatcctccgcGATGCGCACTATGTCTACCGCCACGATGAATACGGCAGGCTGACGGAGAAGACGGACCGCATCCCGGCGGGTGTGATacggacgg                 >  2:1311318‑M2/16‑100 (MQ=255)                                                                                     tgcgcatcctccgcGATGCGCACTATGTCTACCGCCACGATGAATACGGCAGGCTGACGGAGAAGACGGACCGCATCCCGGCGGGTGTGATacggacgg                 >  1:1407850‑M2/15‑99 (MQ=255)                                                                                        gcatcctccgcGATGCGCACTATGTCTACCGCCACGATGAATACGGCAGGCTGACGGAGAAGACGGACCGCATCCCGGCGGGTGTGATACGGACGgacga             >  2:255686‑M2/12‑100 (MQ=255)                                                                                           tcctccgcGATGCGCACTATGTCTACCGCCACGATGAATACGGCAGGCTGACGGAGAAGACGGACCGCATCCCGGCGGGTGTGATACGGACGGacgacg           >  2:3145166‑M2/9‑99 (MQ=255)                                                                                                   GATGCGCACTATGTCTACCGCCACGATGAATACGGCAGGCTGACGGAGAAGACGGACCGCATCCCGGCGGGTGTGATACGGACGGACGACGAGCGGAccc  <  1:2720803‑M2/100‑1 (MQ=255)                                                                                                   |                                                                                                    TCCCGTATGCCACGGACCCGGCGGGCAACCGGCTGCCGGACCCGGAGCTGCACCCGGACAGTACACTCACAGTGTGGCCGGATAACCGCATCGCGGAGGATGCGCACTATGTCTACCGCCACGATGAATACGGCAGGCTGACGGAGAAGACGGACCGCATCCCGGCGGGTGTGATACGGACGGACGACGAGCGGACCC  >  NC_000913/525989‑526186

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