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)	63 (1.210)	59/480	0.1	98.4%	coding (570/687 nt)	fliH	flagellar biosynthesis protein FliH
	?	NC_000913	2514274 =	NA (NA)					noncoding (2/1345 nt)	IS186	repeat region
*	?	NC_000913	= 2016444	1 (0.020)	74 (1.410)	63/482	0.0	98.6%	coding (577/687 nt)	fliH	flagellar biosynthesis protein FliH
	?	NC_000913	= 2515615	NA (NA)					noncoding (1343/1345 nt)	IS186	repeat region

CAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCA‑TTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGAC  >  NC_000913/2515390‑2515840                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                  cAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCAGCTGGAGCTTACAGGTGGCCA‑TTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGcc                                                                                                                                                                                                          <  2:367854/251‑1 (MQ=39)        cGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCA‑TTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATGCCCCGGTTTTTTTACGCCTGTTaa                                                                                                                                                                                                    <  1:350192/250‑1 (MQ=255)                                  ggATCCGCAAAGAAGAACTAACTCGTTGTGGAGACTAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCA‑TTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCggg                                                                                                                                                                                  >  2:324405/1‑242 (MQ=255)                                  ggATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCA‑TTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATc                                                                                                                                                                         <  1:346755/251‑1 (MQ=255)                                  ggATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCA‑TTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATGCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCggg                                                                                                                                                                                  <  1:324405/242‑1 (MQ=255)                                   gATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCA‑TTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCa                                                                                                                                                                        <  1:220021/251‑1 (MQ=255)                                                  actaactCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCA‑TTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAAGGGCTGCGGGCGACCAATCAAATACCCTTGCAGat                                                                                                                                                         <  2:196880/251‑1 (MQ=255)                                                     aactCGTTGTGGAGAATAACAAAAATGGGCATCTGGAGCTTACAGGTGGCCA‑TTCGTGGGACAGTATCCCTGACAGCCTACAAAGCTCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATACTTCGGATGGTTCGACCCGTAGGTGAGCGCTTATTGGATAATACCCCGGTGTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATAtt                                                                                                                                                      <  1:308495/251‑1 (MQ=18)                                                     aactCGTTGTGGAGAATAACAAAAATGGACATCTGGAGCTTACAGGTGGCCA‑TTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATAtt                                                                                                                                                      >  1:117930/1‑251 (MQ=255)                                                                     tAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCA‑TTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGCGCTTATg                                                                                                                                      >  1:74707/1‑251 (MQ=255)                                                                     tAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCA‑TTCGTGGGAAAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCATTACGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGGGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATg                                                                                                                                      >  2:236976/1‑251 (MQ=255)                                                                          aaaaTGGCAATCTGGCGTTTACAGGTGCACA‑TTCGTGGGAAAGAATCCCTGACCGCCTACCAACAGCAATTGAAGACCGCGAGGCATCGTCTTAACGAGGAACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGAGCTTATGGGATAATTCCCTGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAACTACCCTTGCAGATCTTGCACCCCGAGCTTCTGCAAt                                                                                                                                  >  2:284101/1‑250 (MQ=21)                                                                            aaTGGTCATCTGGAGCTTACAGGTGGCCA‑TTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAgc                                                                                                                               <  1:269548/251‑1 (MQ=255)                                                                                   aTCTGGTGCTTACAGGTG‑CCACTTCATGGGTCAGTATACCTGACCGCCCACAAACCGCCATTGAAGTACGCGCAGCATCGTCTTAACGAGGCACCGCGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTACTCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCctgct                                                                                                                        >  2:458676/1‑251 (MQ=25)                                                                                       ggAGCTTACAGGCGGCCA‑TTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCctgctgctg                                                                                                                    >  1:363050/1‑251 (MQ=255)                                                                                              aCAGGTGGCCA‑TTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGtc                                                                                                             >  1:329932/1‑251 (MQ=255)                                                                                                 ggTGGCCA‑TTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGATTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAACAGCGCCTGTTGCTGTTGAGTCTCg                                                                                                          >  2:8804/1‑251 (MQ=255)                                                                                                  gTGGCCA‑TTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCg                                                                                                          >  2:155581/1‑250 (MQ=255)                                                                                                  gTGGCCA‑TTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGa                                                                                                         <  2:329932/251‑1 (MQ=255)                                                                                                                       aTCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGGTAATCAGCTAATGGCTGCTGGCGACCAATCAAATGCCCTTGCAGATATTGCGCCCCGAGCTTATGCGATAGCGCCTTCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACt                                                                                      <  1:458676/250‑1 (MQ=255)                                                                                                                             ggcAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGAATTCCCCGGTTTTTGTACGCCTGGTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTTCAAATAGTGCTCCCCGAGCTTATGCAAGATCGCCTGCTACTTTTGCGTCTCGACAAACTCCGCGACCACACTCAATGAc                                                                               <  1:284101/249‑1 (MQ=25)                                                                                                                             gACAGCCTACAAAACGCAATTGAAGACCGAGAGGCTTCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGAGCCTGCTGCTGTTGCGTCTTGACAAACTCCGCGACCACACTCAATGAc                                                                               >  2:197603/1‑251 (MQ=255)                                                                                                                              aCAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACt                                                                              <  1:155581/251‑1 (MQ=255)                                                                                                                                                    agaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAAt                                                        >  2:47315/1‑251 (MQ=255)                                                                                                                                                      aaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTACGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTg                                                      >  2:293626/1‑251 (MQ=255)                                                                                                                                                      aaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTg                                                      >  2:338189/1‑251 (MQ=255)                                                                                                                                                           gAGGCATCGTCTTAACGCGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTCTGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGAtctc                                                 >  2:224336/1‑251 (MQ=255)                                                                                                                                                                  cGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAatc                                           <  1:325328/250‑1 (MQ=255)                                                                                                                                                                       tAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCCTAATCATCGc                                      >  1:359288/1‑250 (MQ=255)                                                                                                                                                                          cGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTTCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATc                                                              >  2:114354/1‑223 (MQ=255)                                                                                                                                                                          cGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGt                                                                                                              <  1:23690/175‑1 (MQ=255)                                                                                                                                                                          cGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGt                                                                                                              >  2:23690/1‑175 (MQ=255)                                                                                                                                                                          cGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATc                                                              <  1:114354/223‑1 (MQ=255)                                                                                                                                                                                  cGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGtt                          <  1:123523/251‑1 (MQ=255)                                                                                                                                                                                   gAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAg                                                                                                                                                           >  2:130378/1‑121 (MQ=255)                                                                                                                                                                                   gAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAg                                                                                                                                                           <  1:130378/121‑1 (MQ=255)                                                                                                                                                                                        gTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGAc                    <  1:131421/251‑1 (MQ=255)                                                                                                                                                                                          cGCATTCTTCAGATGGTTCAAACCTTAAGTTAGCGCTTATGGGATCATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGCCCACACTCAATGACTTCGCTTTCGCCAGATCGGTACTTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACa                   >  2:292284/1‑250 (MQ=255)                                                                                                                                                                                                          ttCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTATGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGAc  >  2:132237/1‑251 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                  CAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCA‑TTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGAC  >  NC_000913/2515390‑2515840

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