Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC	NZ_CP009273	1,193,954	(AATGATGA)6→5	intergenic (‑261/‑190)	BW25113_RS05920 ← / → lit	hypothetical protein/cell death peptidase Lit

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	1193918 =	1 (0.020)	26 (0.830)	19/158	0.4	96.9%	intergenic (‑225/‑233)	BW25113_RS05920/lit	hypothetical protein/cell death peptidase Lit
	?	NZ_CP009273	= 1193909	1 (0.030)					intergenic (‑216/‑242)	BW25113_RS05920/lit	hypothetical protein/cell death peptidase Lit

TGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTGAGAATCTGAAATTGATAGAAATGAAATAATGAAATAATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAACATGT  >  NZ_CP009273/1193783‑1194049                                                                                                                                        |                                                                                                                                    tttgtatctgttcaggctgacctagcttcgctgacagacaatattgtgatcagtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataatg                                                                                                                                   >  2:597933‑M1/136‑138 (MQ=255)           ttcaggctgacctagcttcgctgacagacaatattgtgatcagtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataatgaaatgatgaaa                                                                                                                        <  1:902846‑M1/14‑1 (MQ=255)                gctgacctagcttcgctgacagacaatattgtgatcagtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataatgaaatgatgaaatgat                                                                                                                    <  2:680186‑M1/18‑1 (MQ=255)                   gacctagcttcgctgacagacaatattgtgatcagtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataatgaaatgatgaaatgatgaaa                                                                                                                <  2:301303‑M1/22‑1 (MQ=255)                               ctgacagacaatattgtgatcagtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataatgaaatgatgaaatgatgaaatgatgaaatgat                                                                                                    <  2:841642‑M1/34‑1 (MQ=255)                                         aTCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTGAGAATCTGAAATTGATAGAAATGAAATAATGAAATAATGAaatgatgaaatgatgaaatgatgaaatgatgaaatgatgaaatgatga                                                                                          >  1:869038/1‑139 (MQ=255)                                          tattgtgatcagtagcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAc                                                                                         >  2:618517‑M1/95‑139 (MQ=255)                                                        gcacgtatcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTg                                                                            >  1:105454‑M1/81‑138 (MQ=255)                                                              atcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACg                                                                      <  2:549638‑M1/64‑1 (MQ=255)                                                               tcgaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGa                                                                    >  1:600811‑M1/74‑139 (MQ=255)                                                                 gaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATg                                                                  >  1:616709‑M1/72‑139 (MQ=255)                                                                 gaggaggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATg                                                                  >  1:680186‑M1/72‑139 (MQ=255)                                                                      ggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTc                                                             <  2:474198‑M1/73‑1 (MQ=255)                                                                      ggagtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTc                                                             >  2:440567‑M1/67‑139 (MQ=255)                                                                        agtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTcc                                                                                <  2:575079‑M1/54‑1 (MQ=255)                                                                        agtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTcc                                                                                >  1:575079‑M1/65‑118 (MQ=255)                                                                        agtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATg                                                                  >  2:499300‑M1/65‑132 (MQ=255)                                                                        agtagcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATg                                                                  <  1:499300‑M1/68‑1 (MQ=255)                                                                            gcgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGa                                                       >  1:598645‑M1/61‑139 (MQ=255)                                                                             cgctacaaatttgacgctgggtgagaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGAc                                                      >  1:449541‑M1/60‑139 (MQ=255)                                                                                                     gaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTAATCGTTCCCAgt                              >  1:549638‑M1/36‑139 (MQ=255)                                                                                                     gaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTAATCGTTCCCAgt                              >  1:166880‑M1/36‑139 (MQ=255)                                                                                                     gaatctgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTAATCGTTCCCAgt                              >  1:739046‑M1/36‑139 (MQ=255)                                                                                                          tgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCa                                                                        >  1:744095‑M1/31‑92 (MQ=255)                                                                                                          tgaaattgatagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCa                                                                        <  2:744095‑M1/62‑1 (MQ=255)                                                                                                                  atagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAg                 >  2:902846‑M1/23‑139 (MQ=255)                                                                                                                  atagaaatgaaataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAg                 >  1:474198‑M1/23‑139 (MQ=255)                                                                                                                             ataatgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTaaa       <  1:712694‑M1/127‑1 (MQ=255)                                                                                                                                atgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATg                                                                  <  2:136245‑M1/68‑1 (MQ=255)                                                                                                                                atgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATg                                                                  >  1:136245‑M1/9‑76 (MQ=255)                                                                                                                                 tgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAACATGt  <  2:709348‑M1/132‑1 (MQ=255)                                                                                                                                 tgaaataATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAACATGt  >  1:708548‑M1/8‑139 (MQ=255)                                                                                                                                        |                                                                                                                                    TGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTGAGAATCTGAAATTGATAGAAATGAAATAATGAAATAATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTAATCGTTCCCAGTGTGCCTATAGCAGATATTCTAAACATGT  >  NZ_CP009273/1193783‑1194049

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