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 =	0 (0.000)	60 (1.320)	55/478	0.0	100%	coding (570/687 nt)	fliH	flagellar biosynthesis protein FliH
	?	NC_000913	2514274 =	NA (NA)					noncoding (2/1345 nt)	IS186	repeat region
*	?	NC_000913	= 2016444	0 (0.000)	57 (1.250)	52/480	0.1	100%	coding (577/687 nt)	fliH	flagellar biosynthesis protein FliH
	?	NC_000913	= 2515615	NA (NA)					noncoding (1343/1345 nt)	IS186	repeat region

GACATAATCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGCC  >  NC_000913/2515381‑2515865                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                           gacATAATCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGttt                                                                                                                                                                                                                                            <  2:103360/251‑1 (MQ=34)                     ggATTTCCCCCCCAGAAGTGACGGAGCCGTAAAGAAGAACTACCTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTGACAGGGGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACTCAATTGAAGAACGCGAGGCATCGTCTGAACGAGGGACCGAGGCGTCTCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGc                                                                                                                                                                                                                        <  1:469439/251‑1 (MQ=14)                       aTTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCt                                                                                                                                                                                                                       >  2:457889/1‑250 (MQ=255)                                               ccGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAAt                                                                                                                                                                                              >  1:104864/1‑251 (MQ=255)                                                    aagaagAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATa                                                                                                                                                                                                                                                        >  2:61298/1‑188 (MQ=11)                                                    aagaagAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATa                                                                                                                                                                                                                                                        <  1:61298/188‑1 (MQ=11)                                                           actaacaCGTTGTGGAGAATCAGAAACATGCCCATCTGGAGCTTACCCGTGGCCATTCGTGGGACAGCATCCCTGACAGCCTACCAACCGCAATAGAAGAACGCGAGGCATCGTATTAACGAGGCACCGAGGCGTCGCATTCTTCAGACGGTTCAACCCTTCAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGAGGGCGACCAATCAAATACCCTTGc                                                                                                                                                                                      >  2:404535/1‑247 (MQ=18)                                                                                              cTGGAGCTTACAGGTGGCAATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCctgctgc                                                                                                                                               >  2:380148/1‑251 (MQ=255)                                                                                               tGGAGCTTACAGGTGGCCATTCGTGGGACCGTATCCCTGACAGCCTACAAAACCCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCAACGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCctgctgct                                                                                                                                              >  1:194157/1‑251 (MQ=255)                                                                                                ggAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGCTGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCctgctgctg                                                                                                                                             >  1:197244/1‑251 (MQ=255)                                                                                                                     gTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCcgc                                                                                                                        >  1:425230/1‑251 (MQ=255)                                                                                                                                 ctcTGACAGCCTAAAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAAGCGTCGCATTCTTCAGATGGCTAAACCCTTAAGTTAGCGCTTATGGGATAATTCCCAGGTTTTTTTAAGCCTGTTAATCAGCTAATGGCTGCGGGCGCCCAATCAAATACACTTGCAGATATTGCACCCCGTGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCCCCTCCACACTCa                                                                                                             >  2:369310/3‑250 (MQ=25)                                                                                                                                 cccTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCa                                                                                                             <  2:210721/250‑1 (MQ=255)                                                                                                                                     gACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGAc                                                                                                        >  1:350252/1‑251 (MQ=255)                                                                                                                                           cTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCt                                                                                                  >  1:127952/1‑251 (MQ=255)                                                                                                                                              cAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTc                                                                                               <  1:152769/251‑1 (MQ=255)                                                                                                                                               aaaaCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCg                                                                                              <  2:345242/251‑1 (MQ=255)                                                                                                                                                 aaCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCTACTACACTCAATGACTTCTCTTTCGcc                                                                                            <  1:392876/251‑1 (MQ=255)                                                                                                                                                      aattgaccaaCGCGAGACACAGTCTTAACGAGGCGCCGAGGCGTCGCCTTCTTCAGATGGTTCACCACTTCAGTTAGCGCTTATGGGATAACTCACCGGTTTTTTCACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGAACCCCGCGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGcctaacg                                                                                       >  2:288099/9‑246 (MQ=255)                                                                                                                                                      aaTTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGGTCAACCCTTAAGTTAGCGCGTATGGGATAATTCCCCGGTTTTTTTACGCCTGGTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATGGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGGGGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGc                                                                                             <  1:288099/245‑1 (MQ=255)                                                                                                                                                            agaaCGCGAGGCATCGTCTTACCGAGGAAACGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAAt                                                                                 >  2:468248/1‑251 (MQ=255)                                                                                                                                                                    aGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGAtctc                                                                          >  1:102741/1‑250 (MQ=255)                                                                                                                                                                              ttAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGc                                                               >  1:297823/1‑251 (MQ=255)                                                                                                                                                                                aaCGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCAt                                                             >  1:106660/1‑251 (MQ=255)                                                                                                                                                                                  cGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTACGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATcc                                                           >  1:461283/1‑251 (MQ=255)                                                                                                                                                                                                gTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATAcc                                                                                                                                                                                           >  1:416724/1‑109 (MQ=255)                                                                                                                                                                                                gTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATAcc                                                                                                                                                                                           <  2:416724/109‑1 (MQ=255)                                                                                                                                                                                                          tcAGATGGTTCAACCCTTCCGTTAGCGCTTATGGGATACTTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTc                                   >  1:130955/1‑251 (MQ=255)                                                                                                                                                                                                           cAGATGGGTCAACCCTTAAGTCAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCGAATGGCTGCGGGCGACCAATCAAATACCCTTTCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCTACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCTTCCAGCGTGTTCGTGACAATATCTTTc                                   <  2:252028/250‑1 (MQ=255)                                                                                                                                                                                                                      aCCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCa                       <  2:106660/251‑1 (MQ=255)                                                                                                                                                                                                                        ccTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATc                     <  2:106207/251‑1 (MQ=255)                                                                                                                                                                                                                                        aTGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCa     <  1:468248/251‑1 (MQ=255)                                                                                                                                                                                                                                           ggATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGAGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGCAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAACGACGCCATCAATTTTGATGATATCAGcc  >  1:164923/1‑251 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                           GACATAATCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCAATTTTGATGATATCAGCC  >  NC_000913/2515381‑2515865

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.