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 ← / → uspC	IS1 repressor TnpA/universal stress protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1979486 =	0 (0.000)	75 (0.800)	55/178	NT	100%	intergenic (‑271/‑267)	insA/uspC	IS1 repressor TnpA/universal stress protein
	?	NC_000913	2101749 =	NA (NA)					noncoding (1195/1195 nt)	IS5	repeat region
*	?	NC_000913	= 1979489	0 (0.000)	71 (0.760)	49/178	NT	100%	intergenic (‑274/‑264)	insA/uspC	IS1 repressor TnpA/universal stress protein
	?	NC_000913	= 2290113	NA (NA)					noncoding (1/1195 nt)	IS5	repeat region

CTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTT  >  NC_000913/2290017‑2290212                                                                                                 |                                                                                                    cTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATa                                                                                                  <  1:630513/100‑1 (MQ=32)   gTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAAc                                                                                                <  1:358035/100‑1 (MQ=33)   gTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAAc                                                                                                <  1:365141/100‑1 (MQ=33)     cGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGc                                                                                              >  2:1623364/1‑100 (MQ=38)     cGGCGAAGGTAAGCTGATGACTCATGAGGAACCCGGTTCTATGGCTCCAGTTGACAAACATTATCTCATATCAGGGACTTGTACGCACCTTCCATAACGc                                                                                              <  2:946488/100‑1 (MQ=11)      ggCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCt                                                                                             >  1:2504898/1‑100 (MQ=38)      ggCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCt                                                                                             >  1:1332481/1‑100 (MQ=38)      ggCGAAGGTAAGTTGAGGACTCATGATGAACCCTGTTCTATGGCACCAGAGGACAAACAAAAACTCATATCAGGGAATTGTTCGAACCTTCTATAAcac                                                                                              >  2:1207674/1‑97 (MQ=11)        cGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTg                                                                                            <  2:2243806/99‑1 (MQ=40)         gAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTa                                                                                          <  1:1824345/100‑1 (MQ=255)         gAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTa                                                                                          <  1:1413403/100‑1 (MQ=255)          aaGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCGTATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAg                                                                                         <  2:1388937/100‑1 (MQ=255)          aaGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAg                                                                                         <  2:1632742/100‑1 (MQ=255)          aaGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAg                                                                                         <  1:2156779/100‑1 (MQ=255)          aaGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAg                                                                                         <  2:436973/100‑1 (MQ=255)            ggTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGcc                                                                                       <  2:269290/100‑1 (MQ=255)             gTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGcca                                                                                      <  2:1753131/100‑1 (MQ=255)                aGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGccacca                                                                                   <  1:984838/100‑1 (MQ=255)                 gTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAg                                                                                  >  1:914604/1‑100 (MQ=255)                 gTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAg                                                                                  <  1:38247/100‑1 (MQ=255)                 gTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAg                                                                                  >  1:1793196/1‑100 (MQ=255)                    gatgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAAc                                                                               <  2:1577226/100‑1 (MQ=255)                     atgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACa                                                                              >  1:1988856/1‑100 (MQ=255)                      tgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAg                                                                             <  1:349834/100‑1 (MQ=255)                      tgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAg                                                                             <  1:420310/100‑1 (MQ=255)                      tgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAg                                                                             <  1:223090/100‑1 (MQ=255)                       gaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGa                                                                            <  2:801636/100‑1 (MQ=255)                       gaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGa                                                                            <  1:752067/100‑1 (MQ=255)                         cTCATGATGAACCCTGTTCTATGGCTCCAGATGCCAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAACGCTGTAGCCACCAGAACAGat                                                                           >  2:1717364/1‑99 (MQ=255)                           cATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATAtt                                                                        >  1:846526/1‑100 (MQ=255)                           cATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATAtt                                                                        >  2:1182832/1‑100 (MQ=255)                           cATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATAtt                                                                        >  2:2055971/1‑100 (MQ=255)                             tgatgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGc                                                                      >  2:849235/1‑100 (MQ=255)                              gatgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTACATAACGCTGTAGCCACCAGAACAGATATTGCg                                                                     >  2:186828/1‑100 (MQ=255)                               atgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCgg                                                                    >  1:888485/1‑100 (MQ=255)                                  aaCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAAc                                                                 <  1:2330195/100‑1 (MQ=255)                                  aaCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAAc                                                                 <  1:675489/100‑1 (MQ=255)                                  aaCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAAc                                                                 >  2:664418/1‑100 (MQ=255)                                    cccTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGa                                                               >  1:1090277/1‑100 (MQ=255)                                     ccTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGAc                                                              <  2:833118/100‑1 (MQ=255)                                     ccTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGAc                                                              <  1:143245/100‑1 (MQ=255)                                             agggCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGaa                                                      <  2:1580057/98‑1 (MQ=255)                                                gCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACa                                                   >  1:258859/1‑100 (MQ=255)                                                      gATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCa                                             >  1:225079/1‑100 (MQ=255)                                                        tGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGa                                           >  2:67066/1‑100 (MQ=255)                                                        tGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGa                                           >  1:1130204/1‑100 (MQ=255)                                                          aCAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGAtt                                         <  1:861180/100‑1 (MQ=255)                                                          aCAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGAtt                                         >  1:1490347/1‑100 (MQ=255)                                                           cAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTg                                        >  1:2147259/1‑100 (MQ=255)                                                           cAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATAGCGCAACGACAAAGAGAAACAGAACCAGact                                         >  2:590103/1‑97 (MQ=255)                                                           cAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGACATTGCGGAACGACAAAGAGAAACAGAACCAGATTg                                        >  1:592102/1‑100 (MQ=255)                                                            aaaCATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGa                                       >  1:219160/1‑100 (MQ=255)                                                             aaCATGATCTCATATCAGGGGCTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGAt                                      >  2:398438/1‑100 (MQ=255)                                                                   aTCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTg                                >  2:671016/1‑100 (MQ=255)                                                                    actcATAGCAGTGACGTGTTCGCACCTGCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGCCCAAGAGAAACCGAACCAGATTGCTGCATTGa                               >  2:1126306/2‑100 (MQ=255)                                                                       cATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCt                            >  2:2505008/1‑100 (MQ=255)                                                                       cATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCt                            >  2:2365580/1‑100 (MQ=255)                                                                         tatCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCttt                          <  2:75329/100‑1 (MQ=255)                                                                         tatCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCttt                          >  1:2558973/1‑100 (MQ=255)                                                                         tatCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCttt                          <  1:1567706/100‑1 (MQ=255)                                                                          atCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTc                         >  1:2219569/1‑100 (MQ=255)                                                                          agcAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTc                         <  2:2266653/98‑1 (MQ=255)                                                                           tCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCa                        >  2:962938/1‑100 (MQ=255)                                                                           tCAGGGACTTGTTCGCACCTTCCATAACGCTGTAACCACCAGAACAGATATTTCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCa                        >  2:1260604/1‑100 (MQ=255)                                                                              gggACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATcc                     <  2:2329448/100‑1 (MQ=255)                                                                                gACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTa                   <  2:1113508/100‑1 (MQ=255)                                                                                 aCTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTAt                  <  2:1402962/100‑1 (MQ=255)                                                                                  cTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATg                 >  1:1015236/1‑100 (MQ=255)                                                                                  cTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATg                 >  1:1995710/1‑100 (MQ=255)                                                                                  cTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATg                 >  2:91374/1‑100 (MQ=255)                                                                                    tggtCGCACATGCCATAGCGCTGTAGGCACCATAACAGATATTGCGGAACGTCAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTAt                  <  2:1424113/94‑1 (MQ=255)                                                                                     gTTCGCACCTTCCCTAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGCAACAGAACCAGACTGATGCATTGAGCTGTCATCCTATgca               >  2:1993337/1‑97 (MQ=255)                                                                                     gTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGaaa              <  2:1130199/100‑1 (MQ=255)                                                                                     gTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGaaa              >  2:1836856/1‑100 (MQ=255)                                                                                       tCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAaatt            >  2:927654/1‑100 (MQ=255)                                                                                           aCCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAaattaatt        <  1:2312298/100‑1 (MQ=255)                                                                                           aCCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAaattaatt        <  2:1415457/100‑1 (MQ=255)                                                                                              ttCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCt     >  1:997106/1‑100 (MQ=255)                                                                                               tCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTg    >  1:2157544/1‑100 (MQ=255)                                                                                                 cATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGtt  >  1:387657/1‑100 (MQ=255)                                                                                                 cATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGtt  <  1:2211014/100‑1 (MQ=255)                                                                                                 |                                                                                                    CTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTT  >  NC_000913/2290017‑2290212

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