Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	2,016,437	Δ2 bp :: IS186 (–) +8 bp :: Δ1 bp	coding (570‑577/687 nt)	fliH →	flagellar biosynthesis protein FliH

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2016437 =	1 (0.020)	75 (1.460)	69/482	0.0	98.6%	coding (570/687 nt)	fliH	flagellar biosynthesis protein FliH
	?	NC_000913	2514274 =	NA (NA)					noncoding (2/1345 nt)	IS186	repeat region
*	?	NC_000913	= 2016444	2 (0.040)	61 (1.190)	59/484	0.1	96.7%	coding (577/687 nt)	fliH	flagellar biosynthesis protein FliH
	?	NC_000913	= 2515615	NA (NA)					noncoding (1343/1345 nt)	IS186	repeat region

GCAGAAAATTGCTTCCACGCTCTCTCCGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTG‑CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTC‑TGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATGTCATTACGTGAAGTCACTGCATGGGC  >  NC_000913/2514025‑2514521                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                         gCAGAAAATTGCTTCCACGCTCTCTCCGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTG‑CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTc                                                                                                                                                                                                                                                          <  2:218961/250‑1 (MQ=255)                  cGCTCTCTCCGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGCAGCTGGTTTACGGCTCTACCCTGGTGAGTGTG‑CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTa                                                                                                                                                                                                                                         >  2:243084/1‑250 (MQ=255)                   gCTCTCTCCGCGTGCTGCAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTG‑CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTa                                                                                                                                                                                                                                         >  1:172674/1‑249 (MQ=255)                        ctcCGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTG‑CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAggg                                                                                                                                                                                                                                     <  1:422401/248‑1 (MQ=255)                                   gAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTG‑CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc                                                                                                                                                                                                                         <  1:274/249‑1 (MQ=255)                                      ggCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTACTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTG‑CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                                                                                                                                                                                                                    >  1:188779/1‑251 (MQ=255)                                                      tCCTGGTTTAATTCGCTACCTTCTATTCAATCGGGATTATCGCAGGCGCAGTTAAAGGACTGAATGAAGAGCAAGGAACCGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTG‑CTGTCTGCGTCAATCGCAGCCCTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGGTGTGGTATTACGCCTGATATGATTTACCGTGc                                                                                                                                                                                                    >  2:502809/1‑251 (MQ=255)                                                       ccTGGTTTCCGTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTG‑CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGAGATGATTTAACGTGcc                                                                                                                                                                                                   <  2:277086/251‑1 (MQ=255)                                                               ccTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAACGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTG‑CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAtt                                                                                                                                                                                           >  2:100615/1‑251 (MQ=255)                                                                   cGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTG‑CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGGATTACGCCTGATATGATTTAACGTGCCGATGAATTAct                                                                                                                                                                                        <  2:442587/250‑1 (MQ=255)                                                                          ttCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTG‑CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGAt                                                                                                                                                                                >  2:9118/1‑251 (MQ=255)                                                                                 caaacGGGATTATCGCAGGGGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGTGGAGTGTG‑CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCTCTAACTTAAGGGTTGTGGTTTTACGCCTTATATGATTTAACGTGCCTATGAGTTACTCTCACGATAACTGGt                                                                                                                                                                         <  2:535322/247‑1 (MQ=255)                                                                                   aTCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTG‑CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCa                                                                                                                                                                       <  1:9118/251‑1 (MQ=255)                                                                                      gggATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTG‑CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGc                                                                                                                                                                     >  2:173370/1‑250 (MQ=255)                                                                                           tATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTACCGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTG‑CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAAc                                                                                                                                                                                                                                            >  1:514641/1‑174 (MQ=255)                                                                                           tATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGTTCTGAAGCGGGTTTACGGCTCTACCCTGGTGAGTGTG‑CTGTCTTCGTCAATCGCAGCACTTTTGCTGTAAGACCATACATAAAAAAGCCGGGTATAATTCCCATTAGCGCTAAc                                                                                                                                                                                                                                            <  2:514641/174‑1 (MQ=255)                                                                                                       aGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTG‑CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTCAGGGTTGTGGTATTACGCCTGCTATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTC‑TGGCCCATATTgg                                                                                                                                                  >  2:357169/1‑251 (MQ=255)                                                                                                             aGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTG‑CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTC‑TGGCCCATATTGGTAAGcc                                                                                                                                            >  1:196474/1‑251 (MQ=255)                                                                                                                 cTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTG‑CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTC‑TGGCCCATATTGGTAAGCCCga                                                                                                                                         >  2:488711/1‑250 (MQ=255)                                                                                                                           agCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTG‑CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTC‑TGGCCCATATTGGTAAGCCCGAAGAACTGGAt                                                                                                                               >  1:533350/1‑250 (MQ=255)                                                                                                                            gCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTG‑CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTC‑TGGCCCATATTGGTAAGCCCGAAGAACTGGATAc                                                                                                                             <  1:357169/251‑1 (MQ=255)                                                                                                                                 gTAACGTGGTTTCTCGCGTCGGTCTGAAGCTGGTTTACTGCTCTACCCTGGTGATGGTG‑CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGGATTACTCCTGATATGATTTAACGTGCCTATGAATTACTCTCACGATAACTGTGCAGGAATTC‑TGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCgg                                                                                                                        <  1:51025/251‑1 (MQ=25)                                                                                                                                        ggTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTG‑CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTC‑TGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTa                                                                                                                  <  1:488711/250‑1 (MQ=255)                                                                                                                                         gTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTG‑CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTC‑TGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAAt                                                                                                                <  1:167789/251‑1 (MQ=255)                                                                                                                                                   tCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTG‑GTTTGGCTGTCTGCGTCAATCGCAGCACTGGTTCTGGAAGTCCATACATAAAAAAGCCGGGTATAATTCCCATAAGCGCTACCTTAAGTGTTTTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCTGCAATTC‑TGGCGCATATCGGTAAGCCCGAAGACCTGGATACTTCGGCACGTAGTGCCGGGGct                                                                                                       <  2:488722/250‑1 (MQ=18)                                                                                                                                                          gAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTG‑CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATAGGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTC‑TGGCCCATATTGGTAAGCTCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAAcc                                                                                                 <  1:373088/249‑1 (MQ=255)                                                                                                                                                            aGCTGGTTTACGGCTCTACCCTGGTGAGTGTG‑CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTC‑TGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACccgcc                                                                                             <  2:172674/251‑1 (MQ=255)                                                                                                                                                                     aCGGCTCTACCCTGGTGAGTGTG‑CTGGCTGCGTCAATCGCATCACTGGTGCTGTAAGACCATACATAAAAAAGGCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTTTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGTTCAGCAATCCGTGGCC‑ATATTGGTAAGCCCGCAGAAGTGGATACTTCGGCACGTAATGCCGGGATTCTGACCCGCCGCCGCGaa                                                                                     <  2:65588/250‑1 (MQ=18)                                                                                                                                                                        gCTCTACCCTGGTGAGTGTG‑CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTC‑TGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTc                                                                                 >  1:22918/1‑251 (MQ=255)                                                                                                                                                                                  ggTGAGTGTG‑CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTC‑TGGCCCATATTGGTAAGCCCGAAGAACTGGTTACTTCGGCACGTAAGGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGAtgct                                                                         <  2:42183/249‑1 (MQ=255)                                                                                                                                                                                       gtgtG‑CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTC‑TGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAActc                                                                  >  2:341792/1‑251 (MQ=255)                                                                                                                                                                                          tg‑CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTC‑TGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATCCGTGCTGCTGCACCTCTGc                                                               >  2:101116/1‑251 (MQ=255)                                                                                                                                                                                                   gCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTC‑TGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTc                                                         <  1:178442/249‑1 (MQ=255)                                                                                                                                                                                                                       gCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTC‑TGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCgg                                   >  1:99059/1‑251 (MQ=255)                                                                                                                                                                                                                                         aaaaaGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTC‑TGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATGTCATTACGTg                  >  1:116188/1‑250 (MQ=37)                                                                                                                                                                                                                                           aaaGCCGGGGATAATGCCCATAAGCGCTCTCTTAAGGGTTGTGGTATTACGCCTGTTATGTTTTAACGTGCTGTGGAATTGCTCGCACGATAACTGGTAAGCAATGC‑TGTCTCAGATTGGTAAGCCCGAAGACCTGGATTCTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGGCGCGAAATTCGTGAGGCTGCAACTCTGCTACGTTGGGGGCTGGCTTACGGCCCCGGGGGGATGTCATTACGTGAAg               <  2:394035/251‑1 (MQ=11)                                                                                                                                                                                                                                                        aTTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTC‑TGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATGTCATTACGTGAAGTCACTGCATGGGc  >  1:4157/1‑251 (MQ=11)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                         GCAGAAAATTGCTTCCACGCTCTCTCCGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTG‑CTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTC‑TGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATGTCATTACGTGAAGTCACTGCATGGGC  >  NC_000913/2514025‑2514521

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