Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC JC	NC_000913	1,979,486	IS5 (+) +4 bp	100%	intergenic (‑271/‑264)	insA‑5 ← / → uspC	IS1 protein InsA/universal stress protein C

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1979486 =	0 (0.000)	23 (1.410)	20/264	NT	100%	intergenic (‑271/‑267)	insA‑5/uspC	IS1 protein InsA/universal stress protein C
	?	NC_000913	2101749 =	NA (NA)					noncoding (1195/1195 nt)	IS5	repeat region
*	?	NC_000913	= 1979489	0 (0.000)	9 (0.550)	9/264	NT	100%	intergenic (‑274/‑264)	insA‑5/uspC	IS1 protein InsA/universal stress protein C
	?	NC_000913	= 2290113	NA (NA)					noncoding (1/1195 nt)	IS5	repeat region

TCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGC  >  NC_000913/2289981‑2290243                                                                                                                                     |                                                                                                                                   tctTTTCTGGTCTGACGGCGGTTACTGCTTAAGTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGt                                                                                                                          >  2:39198/1‑143 (MQ=255)   ttttCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATCACGCTGTAg                                                                                                                        >  1:137246/1‑143 (MQ=37)       cTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCATTTTATGTAATGTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAGACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGccac                                                                                                                    <  2:311048/143‑1 (MQ=14)        tGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGCAGccacc                                                                                                                   >  2:199962/1‑143 (MQ=255)            cTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGACCCCTGTTCTATGGCGCCAGATGACAACCATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGaa                                                                                                               <  1:272471/143‑1 (MQ=255)                            ctgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACg                                                                                               >  1:312716/1‑143 (MQ=255)                                        tCGGTGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAAc                                                                                   <  1:301315/143‑1 (MQ=255)                                           gtgAAGGTAATTTTATGACTCATTATGTACCCTGCTCTATCGCTCCAGATGACAAACATGCTCTCATAGCAGAGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCCGAACGACAAAGAGAAACAGa                                                                                <  2:248498/141‑1 (MQ=255)                                            cGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAACCACCAGAACAGATATTGCGGGAAAATAAAGAGAAACAGaa                                                                               >  2:49044/1‑143 (MQ=18)                                                                      aaCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCCGGGACTTGTTCGCACCACCCATACCGCCGTAGCCACCAGAACAGAAATTGCGAAACCACAAACAGAACCAGAACCAGATTGATGCAATGAGCTTTCATc                                                     >  1:81263/1‑143 (MQ=21)                                                                                                     tGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTgg                      <  2:198428/143‑1 (MQ=255)                                                                                                                     aCTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCGGATGAAATTAAGTGCTGTTAAAAGCATTGGGTACAGAAAATACCCa      >  1:104216/1‑143 (MQ=255)                                                                                                                        tGTTCGCACCTTCCATAACGCTGTAGCCAGCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAg   <  1:303255/143‑1 (MQ=255)                                                                                                                         gTTCGCACCTTCCATAACGCTGTAGCCACCATAACAGATACTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGc  <  1:199962/143‑1 (MQ=255)                                                                                                                                     |                                                                                                                                   TCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGC  >  NC_000913/2289981‑2290243

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