New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	= 2000783	9 (0.620)	3 (0.230)	3/154	1.3	24.8%	coding (1147/1488 nt)	amyA	alpha‑amylase
	?	NZ_CP009273	= 2000836	10 (0.760)					coding (1200/1488 nt)	amyA	alpha‑amylase

GCGCATTACGAAGATGTCGGTGGTGACGGGCAAACCTATCCGATAGATATGCCAATAATCGAACAGCTTGATGAGTTAATTCTCGCCCGTCAGCGTTTCGCCCACGGTGTACAGACGTTATTTTTCGACCATCCGAACTGCATTG  >  NZ_CP009273/2000753‑2000897                                                                                    |                                                              gcgcATTACGAAGATGTCGGTGGTGACGGGCAAACCTATCCGATAGATATGCCAATAATCGAACAGCTTGATGAGTTAATTCTCgcccgtc                                                        <  2:55344‑M2/91‑8 (MQ=255)   gcATTACGAAGATGTCGGTGGTGACGGGCAAACCTATCCGATAGATATGCCAATAATCGAACAGCTTGATGAGTTAATTCTCgcccgtcac                                                      >  1:159967‑M2/1‑82 (MQ=255)               tGTCGGTGGTGACGGGCAAACCTATCCGATAGATATGCCAATAATCGAACAGCTTGATGAGTTAATTCTCGCCCGTCAGCGTTTCGCCCAc                                          <  1:55344/91‑1 (MQ=255)                      ggtgACGGGCAAACCTATCCGATAGATATGCCAATAATCGAACAGCTTGATGAGTTAATTCTCgcccgtcaccaccgacatcttcgtaatg                                   >  2:173522‑M2/1‑63 (MQ=255)                       gtgACGGGCAAACCTATCCGATAGATATGCCAATAATCGAACAGCTTGATGAGTTAATTCTCGCCCGTCAGCGTTTCGCCCACGGTGTACa                                  <  1:129356/91‑1 (MQ=255)                             ggCAAACCTATCCGATAGATATGCCAATAATCGAACAGCTTGATGAGTTAATTCTCGCCCGTCAGCGTTTCGCCCACGGTGTACAGACGtt                            <  1:182601/91‑1 (MQ=255)                             ggCAAACCTATCCGATAGATATGCCAATAATCGAACAGCTTGATGAGTTAATTCTCGCCCGTCAGCGTTTCGCCCACGGTGTACAGACGtt                            <  2:334806/91‑1 (MQ=255)                                     tATCCGATAGATATGCCAATAATCGAACAGCTTGATGAGTTAATTCTCgcccgtcaccaccgacatcttcgtaatgcgcaccgtagaggtc                    >  2:306580‑M2/1‑48 (MQ=255)                                      aTCCGATAGATATGCCAATAATCGAACAGCTTGATGAGTTAATTCTCGCCCGTCAGCGTTTCGCCCACGGTGTACAGCCGTTATTTTTCGa                   <  2:273416/91‑1 (MQ=255)                                      aTCCGATAGATATGCCAATAATCGAACAGCTTGATGAGTTAATTCTCGCCCGTCAGCGTTTCGCCCACGGTGTACAGACGTTATTTTTCGa                   >  1:214143/1‑91 (MQ=255)                                                tatGCCAATAATCGAACAGCTTGATGAGTTAATTCTCGCCCGTCAGCGTTTCGCCCACGGTGTACAGACGTTATTTTTCGACCATCCGAAc         <  2:214143/91‑1 (MQ=255)                                                       ataatCGAACAGCTTGATGAGTTAATTCTCGCCCGTCAGCGTTTCGCCCACGGTGTACAGACGTTATTTTTCGACCATCCGAACTGCATTg  >  1:233682/1‑91 (MQ=255)                                                                            ttAATTCTCGCCCGTCAGCGTTTCGCCCACGGTGTAc                                   >  1:157504/1‑37 (MQ=255)                                                                            ttAATTCTCGCCCGTCAGCGTTTCGCCCACGGTGTAc                                   <  2:157504/37‑1 (MQ=255)                                                                                    |                                                              GCGCATTACGAAGATGTCGGTGGTGACGGGCAAACCTATCCGATAGATATGCCAATAATCGAACAGCTTGATGAGTTAATTCTCGCCCGTCAGCGTTTCGCCCACGGTGTACAGACGTTATTTTTCGACCATCCGAACTGCATTG  >  NZ_CP009273/2000753‑2000897

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