Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC JC	NC_000913	3,805,650	IS1 (+) +8 bp	23.1%	coding (294‑301/798 nt)	waaP ←	lipopolysaccharide core heptose (I) kinase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3583428 =	NA (NA)	17 (0.190)	16/176	NT	21.0%	noncoding (1/768 nt)	IS1	repeat region
	?	NC_000913	= 3805657	58 (0.700)					coding (294/798 nt)	waaP	lipopolysaccharide core heptose (I) kinase
*	?	NC_000913	= 3584195	NA (NA)	17 (0.190)	16/170	NT	21.9%	noncoding (768/768 nt)	IS1	repeat region
	?	NC_000913	3805650 =	55 (0.690)					coding (301/798 nt)	waaP	lipopolysaccharide core heptose (I) kinase

GTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3583523‑3583428 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGATCTCACTCCCACAATTAGCCTTGAAGATTATTGTGCCGAT  <  NC_000913/3805657‑3805567                                                                                                                                                                                             GTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCGA                                                                                          >  2:480267/1‑99            GGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCGAAAAAGGGTTAA                                                                               >  1:1047983/1‑99             GAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCGAAAAAGGGTTAAA                                                                              <  2:1975989/99‑1                    AGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCGAAAAAGGGTTAAATCCATTA                                                                       >  1:1597338/1‑99                       TGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCGAAAAAGGGTTAAATCCATTAACT                                                                    <  2:1608472/99‑1                          TAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCGAAAAAGGGTTAAATCCATTAACTCGC                                                                 >  1:2190568/1‑99                              AACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCGAAAAAGGGTTAAATCCATTAACTCGCGCAT                                                             <  2:68854/99‑1                                CAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCGAAAAAGGGTTAAATCCATTAACTCGCGCATCA                                                           >  1:1481611/1‑99                                  GAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATT                                                         <  2:144101/99‑1                                           GGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTAC                                                >  1:1324608/1‑99                                                 CCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGA                                          >  1:2131124/1‑99                                                     AAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGATCTC                                      <  2:713012/99‑1                                                                 ATACACTAAATCAGTAAGTTGGCAGCATCACCCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGATCTCACTCCCACAATT                          >  1:1175485/1‑99                                                                      CTAAATCAGTAAGTTGGCAGCATCACCCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGATCTCACTCCCACAATTAGCCT                     >  2:288884/1‑99                                                                              GTAAGTTGGCAGCATCACCCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGATCTCACTCCCACAATTAGCCTTGAAGATT             >  1:1721951/1‑99                                                                                  GTTGGCAGCATCACCCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGATCTCACTCCCACAATTAGCCTTGAAGATTATTG         >  2:1391052/1‑99                                                                                         GCATCACCCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGATCTCACTCCCACAATTAGCCTTGAAGATTATTGTGCCGAT  <  1:2212489/99‑1                                                                                         GCATCACCCGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGATCTCACTCCCACAATTAGCCTTGAAGATTATTGTGCCGAT  <  2:857640/99‑1                                                                                                                                                                                             GTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3583523‑3583428 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CGAAAAAGGGTTAAATCCATTAACTCGCGCATCATTTATTATTACCGAAGATCTCACTCCCACAATTAGCCTTGAAGATTATTGTGCCGAT  <  NC_000913/3805657‑3805567

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