Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	458,790	IS186 (–) +6 bp :: Δ1 bp	intergenic (+90/‑93)	clpX → / → lon	ATPase and specificity subunit of ClpX‑ClpP ATP‑dependent serine protease/DNA‑binding ATP‑dependent protease La

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	458790 =	0 (0.000)	30 (1.140)	28/278	0.1	100%	intergenic (+90/‑98)	clpX/lon	ATPase and specificity subunit of ClpX‑ClpP ATP‑dependent serine protease/DNA‑binding ATP‑dependent protease La
	?	NC_000913	608008 =	NA (NA)					noncoding (2/1345 nt)	IS186	repeat region
*	?	NC_000913	= 458795	0 (0.000)	26 (1.000)	25/276	0.1	100%	intergenic (+95/‑93)	clpX/lon	ATPase and specificity subunit of ClpX‑ClpP ATP‑dependent serine protease/DNA‑binding ATP‑dependent protease La
	?	NC_000913	= 609351	NA (NA)					noncoding (1345/1345 nt)	IS186	repeat region

TGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGC  >  NC_000913/607870‑608147                                                                                                                                           |                                                                                                                                            tGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCAt                                                                                                                                         <  1:87686/143‑1 (MQ=255)   cTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATaa                                                                                                                                       <  2:459757/143‑1 (MQ=255)   cTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATaa                                                                                                                                       >  1:440788/1‑143 (MQ=255)     gtgtTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAgc                                                                                                                                     >  1:228632/1‑143 (MQ=255)               gTGCTGGGGTGTGCGCTTCTGGTCGGAGAGTCGCTGGGGGAGTGTACGGTGAAGGATCGTAATATTGAGTGTAGGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTaa                                                                                                                           <  2:55200/143‑1 (MQ=255)                        tttACGCTTCTGGGCGGGGAGTCGCTGTGTTAGTTGACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGTCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTgg                                                                                                                  <  1:74888/143‑1 (MQ=39)                             gCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAGGTAGCCGTTGTGCGGGGAGTAGCCCCCTTACCCCTCACTTCAGTGTTGTGGTAttc                                                                                                             >  1:29298/1‑142 (MQ=255)                                      cGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGACGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAAGGTTGTGGTATTACGCCTGata                                                                                                    >  1:430988/1‑143 (MQ=255)                                          ggcTCGCTGGGTGATGTTTCGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAt                                                                                                <  2:322490/141‑1 (MQ=255)                                           tcTCGTTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGGCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAtt                                                                                               <  2:8011/142‑1 (MQ=255)                                                           tACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAt                                                                               <  1:420889/143‑1 (MQ=255)                                                              ggTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAc                                                                            >  1:267126/1‑143 (MQ=255)                                                                tGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActc                                                                          >  2:6120/1‑143 (MQ=255)                                                                       aCGTAATATTGAGCTTAAGGCTGTTCTCGCTGACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGAt                                                                   <  2:12500/143‑1 (MQ=255)                                                                                tGAGTGTAAGGCTGTTGTCGCTTACGCACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTCACTTAAGGCTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCa                                                          <  2:6871/143‑1 (MQ=21)                                                                                 gAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAATCTCACGATAACTGGTCAg                                                         >  2:139622/1‑143 (MQ=255)                                                                                           cTGTTCTCCCGTACGATCCGAAGATGGAGCCGGTGTGCGGGGAGTAATCCCATAAGCGCTAATTTGAGGGTCGCGGTCTCACCTGTGTTTTGATTTTACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc                                               <  2:319605/143‑1 (MQ=17)                                                                                                  cGCTTACGAACCGAAGAAGTAGCCGTTGTTCGGGGAGTAATCCCATAAGCACTAACTTAAGGGATGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATAtt                                        >  1:369640/1‑143 (MQ=255)                                                                                                    cTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCCTATTgg                                      >  1:276545/1‑143 (MQ=255)                                                                                                         gAACCGAAGAAGGAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATCACTGGTCAGCAATTCTGGCCCATATTGGTAAGc                                 <  2:440110/143‑1 (MQ=255)                                                                                                              gaagaaGTAGCCTTGGTGCGGGCATTAATCCCATGAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGGTAACCGGTCAGCAATTCTGGCCCATATTGGTAAGCCCgaa                            <  2:200769/143‑1 (MQ=14)                                                                                                                   agTAGCCGTTGTGCGGGGAGTAAGCCCATAAGCGCTAACTTAAGGGTTGTGGGATTACGCCTGATATGATGTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACt                       <  1:138532/143‑1 (MQ=255)                                                                                                                     tAGCCGGTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACGCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTgg                     <  1:127062/143‑1 (MQ=255)                                                                                                                      aGCCGTTGAGCGGAGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGa                    <  1:3976/143‑1 (MQ=255)                                                                                                                        ccGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCCATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATa                  >  2:306973/1‑143 (MQ=255)                                                                                                                         cGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATAc                 >  1:160447/1‑143 (MQ=255)                                                                                                                                  gggAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAGGCCCGACGAACTGGATACTTAGGCaag        >  2:16601/1‑141 (MQ=255)                                                                                                                                       tAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATg   >  1:246890/1‑143 (MQ=32)                                                                                                                                        aaTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCGCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGc  <  1:269235/143‑1 (MQ=14)                                                                                                                                           |                                                                                                                                            TGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGC  >  NC_000913/607870‑608147

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.