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)	46 (1.010)	46/482	0.1	97.8%	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)	62 (1.350)	58/484	0.0	98.4%	coding (577/687 nt)	fliH	flagellar biosynthesis protein FliH
	?	NC_000913	= 2515615	NA (NA)					noncoding (1343/1345 nt)	IS186	repeat region

TTCCACGCTCTCTCCGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAA‑GCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATGTCATTACGTGAAGTCACTGCAT  >  NC_000913/2514037‑2514517                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     ttCCACGCTCTCTCCGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTa                                                                                                                                                                                                                                          >  1:119456/1‑250 (MQ=255)              ccGCGTGCTGCAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGCAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTtgtg                                                                                                                                                                                                                            >  2:182303/1‑251 (MQ=255)              ccGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGtt                                                                                                                                                                                                                               <  2:152151/248‑1 (MQ=255)                       gAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACGTAAGTGGTGTGGTATTACGc                                                                                                                                                                                                                   <  2:463663/251‑1 (MQ=255)                          ggCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                                                                                                                                                                                                                <  1:6821/251‑1 (MQ=255)                                   tctGTGATCCTGGTTTCCTTCGCTAACTTCTCTTAAATCGGGATTATAGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTCAGCCTGATATGAttt                                                                                                                                                                                                       >  2:401925/1‑251 (MQ=255)                                    ctGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCACGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTTCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAttt                                                                                                                                                                                                       >  1:27688/1‑250 (MQ=255)                                     tgtgTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAACGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTa                                                                                                                                                                                                      >  1:33189/1‑250 (MQ=255)                                              ggTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt                                                                                                                                                                                            <  1:80167/251‑1 (MQ=255)                                                      tCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa                                                                                                                                                                                      <  1:212050/249‑1 (MQ=255)                                                       cGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActc                                                                                                                                                                                   >  1:277792/1‑251 (MQ=255)                                                                               tATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCGGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTc                                                                                                                                                            <  2:453117/250‑1 (MQ=255)                                                                                aTCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTg                                                                                                                                                          <  1:401865/251‑1 (MQ=255)                                                                                     aGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTATCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGccc                                                                                                                                                      >  1:451383/1‑250 (MQ=255)                                                                                       gTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCatt                                                                                                                                                   >  1:156227/1‑250 (MQ=255)                                                                                                               agCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGCGAGTGTGCTGTCTGCGTCACTCGCAGCACTGGTGCTGTACGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAA‑GCCCGACGAACTGGAt                                                                                                                           >  2:454411/1‑250 (MQ=255)                                                                                                                       aaCGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTAAGGCTCTACCCTGGCGAGTGTGCTGTCTGCGTAAATCGCAGCACTGGTGCTGTAAGACCATACCTAAAAAAGCCGGGGATAATTCCCATACGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAC‑GCTCGCAGAACTGGATACTTCGGcc                                                                                                                  >  2:407151/1‑250 (MQ=25)                                                                                                                               ttCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAA‑GCCCGAAGAACTGGATACTTCGGCACGTAATg                                                                                                           <  1:49977/250‑1 (MQ=255)                                                                                                                                                aGCTGGTTGACGGCTCTACCCTGGGGAGTGTGCTGTCAGCGACAATCGCAGCACTGGTGCTGGAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGAGATGATTTAACGTGCCGATGAATGAGTCGCACGATAACTGGACAGCAAGTCTGGCCCATAG‑GGTAAAGCCCGAAGAACTGGATACTTCGGCACGTGACGGCGGGGGTCTAACccgcc                                                                                         <  2:451383/251‑1 (MQ=18)                                                                                                                                                                   ccTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAA‑GCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGAtgc                                                                      >  1:102717/1‑251 (MQ=255)                                                                                                                                                                                  gTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAA‑GCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACg                                                       >  1:390451/1‑251 (MQ=255)                                                                                                                                                                                      gCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCGGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAA‑GCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTg                                                   <  2:277792/251‑1 (MQ=255)                                                                                                                                                                                        gTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAA‑GCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTggg                                                 >  1:157934/1‑251 (MQ=255)                                                                                                                                                                                          cAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAA‑GCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTgggg                                                <  1:112890/250‑1 (MQ=255)                                                                                                                                                                                                  gcaCTGGTGCTGTAAGCACATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTTTGGTATTACGCCTGATATGATTTACCGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAA‑GCCCTAAGAACTGGATACTTCGGCGCGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGGGATGCTGCAACTCTGCTACGTCGGGGGCTGGCTTAc                                       <  2:239781/251‑1 (MQ=14)                                                                                                                                                                                                      tGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAA‑GCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGcc                                   >  2:32087/1‑251 (MQ=255)                                                                                                                                                                                                               aaGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAA‑GCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCgggggg                           >  1:344022/1‑250 (MQ=255)                                                                                                                                                                                                                                   cGGGGATAATTCCCAGAAGCGTTAACTTAAGGGTTGTGTTATTACGCCTGATATGATTAAACGGGCCGAGGAATTACGCTCACGAGAACTGGTCAGCAATTCAGGCCCATATTGGTAA‑GCCCGAAGAACTGGATACTTCGGCACGGAATGCCGGGGCTCTAACCCGCCGGCGCGAGATGCGTGATGCTGCAACTCTGCTACTTCTGGGGCTGGCTTACGGCCCCGGGGGGATGTCATTACGTGAAGACACt      <  2:448302/251‑1 (MQ=11)                                                                                                                                                                                                                                       gATAATTCCCATAATGGCTAACTTAAGGGTGGTGGGTTTAGCCCTGATATGATTTAACGTGCCGAAGAATTTCTCTCACGATAACGGGACAGCAATTCTGGCACATATTGGTAA‑GCCCGAAGAACTGGATAATTCGGCACGAAATGCCGGGGCGCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGGCTGGGGCAGGCATCCGGCCCCGGGGGGATGTCATTTCGTGAAGTCACTGCAt  <  2:27688/251‑1 (MQ=11)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     TTCCACGCTCTCTCCGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAA‑GCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATGTCATTACGTGAAGTCACTGCAT  >  NC_000913/2514037‑2514517

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