Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	944,912	Δ1 bp :: IS186 (+) +6 bp :: Δ1 bp	intergenic (+16/+14)	dmsC → / ← ycaC	dimethyl sulfoxide reductase, anaerobic, subunit C/putative isochorismatase family hydrolase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	608008 =	NA (NA)	38 (0.960)	36/292	0.3	100%	noncoding (2/1345 nt)	IS186	repeat region
	?	NC_000913	= 944917	0 (0.000)					intergenic (+21/+14)	dmsC/ycaC	dimethyl sulfoxide reductase, anaerobic, subunit C/putative isochorismatase family hydrolase
*	?	NC_000913	= 609350	NA (NA)	46 (1.170)	45/292	0.1	100%	noncoding (1344/1345 nt)	IS186	repeat region
	?	NC_000913	944912 =	0 (0.000)					intergenic (+16/+19)	dmsC/ycaC	dimethyl sulfoxide reductase, anaerobic, subunit C/putative isochorismatase family hydrolase

CATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTATCTTTTCTTTTAT  >  NC_000913/609205‑609500                                                                                                                                                  |                                                                                                                                                       cATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTa                                                                                                                                                    >  1:111890/1‑150 (MQ=31) cATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTa                                                                                                                                                    <  2:111890/150‑1 (MQ=31)  aTCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAAt                                                                                                                                                  >  1:381617/1‑151 (MQ=33)       gAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCg                                                                                                                                             <  2:100665/151‑1 (MQ=35)            tACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCatat                                                                                                                                        >  1:612287/1‑151 (MQ=39)             aCAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGACCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCatat                                                                                                                                        <  1:437401/150‑1 (MQ=255)              cAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATcc                                                                                                                                      <  2:230296/151‑1 (MQ=40)              cAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATcc                                                                                                                                      <  1:424116/151‑1 (MQ=40)               aGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCg                                                                                                                                     >  1:288867/1‑151 (MQ=40)               aGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCg                                                                                                                                     >  1:561819/1‑151 (MQ=40)               aGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCg                                                                                                                                     >  2:619422/1‑151 (MQ=40)                   ggCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTg                                                                                                                                 <  2:366488/151‑1 (MQ=255)                     ccATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTgtc                                                                                                                               <  2:326012/151‑1 (MQ=255)                      cagtCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTgtca                                                                                                                              <  2:252681/148‑1 (MQ=255)                       aTTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTgtc                                                                                                                               >  1:38206/1‑149 (MQ=255)                          cGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTgtcaggtc                                                                                                                          <  2:605072/151‑1 (MQ=255)                                  aGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGt                                                                                                                  <  1:600461/151‑1 (MQ=255)                                   gTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGt                                                                                                                  >  2:371485/1‑150 (MQ=255)                                    tATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTa                                                                                                                 >  1:424040/1‑150 (MQ=255)                                       cccTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAgg                                                                                                              <  1:350401/150‑1 (MQ=255)                                       cccTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGc                                                                                                             >  2:97805/1‑151 (MQ=255)                                                 cTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCttt                                                                                                    >  2:616369/1‑150 (MQ=255)                                                   aCAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTaca                                                                                                 <  1:488115/151‑1 (MQ=255)                                                         cGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCg                                                                                           >  2:215575/1‑151 (MQ=255)                                                           cAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGcc                                                                                         <  2:609993/151‑1 (MQ=255)                                                           cAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGcc                                                                                         >  1:509618/1‑151 (MQ=255)                                                            aaTTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGccc                                                                                        >  2:157958/1‑151 (MQ=255)                                                                    aaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTc                                                                                  >  2:524778/1‑149 (MQ=255)                                                                    aaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAg                                                                                <  1:331228/151‑1 (MQ=255)                                                                                     tAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCg                                                                 <  2:385801/149‑1 (MQ=255)                                                                                                 gAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGAtg                                                   <  1:483233/151‑1 (MQ=255)                                                                                                        cGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCa                                             >  1:508686/1‑150 (MQ=255)                                                                                                           aTTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTa                                         <  1:283936/151‑1 (MQ=255)                                                                                                                 cAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCgtg                                   <  1:611649/151‑1 (MQ=255)                                                                                                                        ttCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGc                              >  2:386596/1‑149 (MQ=255)                                                                                                                          cAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGcaca                           <  2:509607/150‑1 (MQ=255)                                                                                                                            aCCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCg                        <  1:599944/151‑1 (MQ=255)                                                                                                                              ccTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGtt                      >  1:267351/1‑151 (MQ=255)                                                                                                                               cTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGtt                      >  2:563664/1‑150 (MQ=255)                                                                                                                                 tAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTAt                    <  2:446813/150‑1 (MQ=255)                                                                                                                                 tAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTAt                    <  1:157784/150‑1 (MQ=255)                                                                                                                                 tAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATg                   >  2:3218/1‑151 (MQ=255)                                                                                                                                 tAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATg                   >  1:567759/1‑151 (MQ=255)                                                                                                                                       aGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCGCAGCGTTATGGCttt              <  2:40701/150‑3 (MQ=255)                                                                                                                                       aGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTAtc             <  1:97805/151‑1 (MQ=255)                                                                                                                                          gcTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTAtcttt          <  1:103311/151‑1 (MQ=255)                                                                                                                                             tATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTAtctttt         <  2:293253/149‑1 (MQ=255)                                                                                                                                                gggAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTAtcttttcttt     <  2:390794/150‑1 (MQ=255)                                                                                                                                                  gAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTATCTTTTCTTTTAt  <  2:246123/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                       CATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTATCTTTTCTTTTAT  >  NC_000913/609205‑609500

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