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)	61 (0.980)	54/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)	56 (0.900)	53/292	0.3	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

GTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGG  >  NC_000913/607863‑608151                                                                                                                                                  |                                                                                                                                                gTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATa                                                                                                                                            <  2:295856/151‑1 (MQ=255)   cATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATaa                                                                                                                                           <  2:160973/150‑1 (MQ=255)      aGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAgcgc                                                                                                                                       >  1:731369/1‑151 (MQ=255)      aGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAgcg                                                                                                                                        <  1:518016/150‑1 (MQ=255)                 tAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAgg                                                                                                                             <  1:500191/150‑1 (MQ=255)                  aaCGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGt                                                                                                                           <  2:1022153/151‑1 (MQ=255)                   aCGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGtt                                                                                                                          <  1:947759/151‑1 (MQ=255)                       tGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTgg                                                                                                                      >  1:908295/1‑151 (MQ=255)                        gatGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTgg                                                                                                                      >  2:780607/3‑150 (MQ=255)                        gCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTgg                                                                                                                      <  1:221749/150‑1 (MQ=255)                        gCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTgg                                                                                                                      >  2:213141/1‑150 (MQ=255)                                  aCGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                                                                                                           <  2:640795/151‑1 (MQ=255)                                  aCGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                                                                                                           >  1:665650/1‑151 (MQ=255)                                  aCGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                                                                                                           <  1:92844/151‑1 (MQ=255)                                   cGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                                                                                                           >  2:534804/1‑150 (MQ=255)                                        cTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGatat                                                                                                       <  1:417460/149‑1 (MQ=255)                                         tGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAt                                                                                                    >  1:554960/1‑151 (MQ=255)                                             cGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTa                                                                                                 <  2:575430/150‑1 (MQ=255)                                             cGGAGACTCGCTGTGTGAGTGTACGGTGAAGGAACGTAATATTGAGATTAAGGCTGATCTCGCTTACGAACCGAAGAAGGAGCCGTTGTGCGTGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTa                                                                                                 <  2:583224/150‑1 (MQ=255)                                                      gCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt                                                                                       <  2:289860/151‑1 (MQ=255)                                                      gCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt                                                                                       >  1:467400/1‑151 (MQ=255)                                                         gtgtgAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTTTGCGGGGAGTAATCCCTTAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGaa                                                                                    <  2:843671/151‑1 (MQ=255)                                                         gtgtgAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATg                                                                                      <  1:374936/149‑1 (MQ=255)                                                         gtgtgAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGaa                                                                                    <  1:72429/151‑1 (MQ=255)                                                          tgtgAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAt                                                                                   >  1:1023218/1‑151 (MQ=255)                                                                 ttACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGCTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACtctc                                                                            >  2:316470/1‑151 (MQ=255)                                                                  tACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACtctc                                                                            >  2:178448/1‑150 (MQ=255)                                                                   aCGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCAc                                                                          >  2:468786/1‑151 (MQ=255)                                                                      gTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGAt                                                                       >  1:153290/1‑151 (MQ=255)                                                                       tGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGAt                                                                       >  1:809430/1‑150 (MQ=255)                                                                                gTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTc                                                               >  1:134551/1‑149 (MQ=255)                                                                                  aaTATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGc                                                            <  1:175586/150‑1 (MQ=255)                                                                                  aaTATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGc                                                            >  1:132849/1‑150 (MQ=255)                                                                                   atatTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCaa                                                          <  1:345130/151‑1 (MQ=255)                                                                                   atatTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCaa                                                          <  2:626327/151‑1 (MQ=255)                                                                                     atTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAAtt                                                        >  1:297884/1‑151 (MQ=255)                                                                                     atTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAAtt                                                        >  1:775848/1‑151 (MQ=255)                                                                                     atTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAAt                                                         <  2:211220/150‑1 (MQ=255)                                                                                     atTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAAt                                                         <  2:390020/150‑1 (MQ=255)                                                                                        gAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCt                                                      <  1:139536/150‑1 (MQ=255)                                                                                        gAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCt                                                      <  1:703208/150‑1 (MQ=255)                                                                                             tAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGccc                                                 <  1:782689/150‑1 (MQ=255)                                                                                              aaGGCTGTTCTCGCTTACGATCCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATTAATTACTCTCACGATAACTGGTCAGCAATTCTGGCACat                                               <  2:146463/151‑1 (MQ=255)                                                                                               aGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCata                                              <  2:653697/151‑1 (MQ=255)                                                                                                 gCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATAtt                                            <  2:82530/151‑1 (MQ=255)                                                                                                      tcgcGCTTACGAACCGAAGAAGTAGCCGTTGCGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATAAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGt                                         >  2:988546/4‑149 (MQ=255)                                                                                                         cGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTaa                                       <  1:614430/148‑1 (MQ=255)                                                                                                         cGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAg                                      <  1:987820/149‑1 (MQ=255)                                                                                                         cGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGca                                    <  2:667856/151‑2 (MQ=255)                                                                                                           cTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCg                                  >  1:21833/1‑151 (MQ=255)                                                                                                             tACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGccca                                  <  2:554929/149‑2 (MQ=255)                                                                                                                gAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCAGATTGGTAAGCCCgaagaa                             <  2:106244/151‑1 (MQ=255)                                                                                                                  aCCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACt                           <  1:625639/151‑1 (MQ=255)                                                                                                                         aagTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACt                    <  1:770477/151‑1 (MQ=255)                                                                                                                          agTAGCCGTTGTGCGGTGAGTAATCCCATAAGCGCTAACTTAAGGGGTGGGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGACCCATATTGGTAAGCCCGAAGAACTGGATGCtt                   <  2:773118/151‑1 (MQ=255)                                                                                                                          agTAGCCGTTGCGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACtt                   >  1:307817/1‑151 (MQ=255)                                                                                                                                  ttGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGt           >  2:222776/1‑151 (MQ=39)                                                                                                                                      gCGGGGAGTAATCCCATAAGCGATAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAATGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCTCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAAt        >  2:197576/1‑150 (MQ=255)                                                                                                                                           gAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCggg  <  1:1009376/151‑1 (MQ=33)                                                                                                                                                  |                                                                                                                                                GTCATAGTGCTGTGTTTAACGGTGCTGGGATTTACGCTTCTGGTCGGAGACTCGCTGTGTGAGTTTACGGTGAAGGAACGTAATATTGAGTTTAAGGCTGTTCTCGCTTACGAACCGAAGAAGTAGCCGTTGTGCGGGGAGTAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGG  >  NC_000913/607863‑608151

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