Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,978,503	Δ776 bp		insB1–insA	insB1, insA

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1978503–1979269	1979278	10–776	51 [0]	[0] 52	insB1–insA	insB1,insA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1978502	0 (0.000)	49 (1.250)	42/266	0.0	100%	intergenic (‑305/+16)	flhD/insB1	flagellar class II regulon transcriptional activator, with FlhC/IS1 transposase B
	?	NC_000913	1979279 =	0 (0.000)					intergenic (‑64/‑474)	insA/uspC	IS1 repressor TnpA/universal stress protein

TTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCC  >  NC_000913/1978366‑1978642                                                                                                                                         |                                                                                                                                             tttACACTTATCTAAGTTTTGTCCTGAATCGACGCATCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGCTTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataag                                                                                                                                        <  2:575591‑M1/143‑7 (MQ=255) tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataag                                                                                                                                        >  1:685575‑M1/1‑137 (MQ=255)    acacTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACCTAAATGtataagtca                                                                                                                                     >  2:185166‑M1/1‑134 (MQ=255)           tCTAAGAGTTTTCCTAAATCGCCGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttt                                                                                                                              <  1:743709‑M1/143‑17 (MQ=255)           tCTAACATTTCTCCGCAAGCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCTATTTTCAATAATGCTTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttt                                                                                                                              <  2:378945‑M1/143‑17 (MQ=255)             tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAACAAGGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgt                                                                                                                            >  1:557563‑M1/1‑125 (MQ=255)                gATTTTTCCTAAATCGACGCGACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttt                                                                                                                         <  2:227515‑M1/143‑22 (MQ=255)                  ttgttcctAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                       <  1:651965‑M1/140‑24 (MQ=255)                         tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                                                                                                >  2:383344‑M1/1‑113 (MQ=255)                         tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttaatttgggtgtatt                                                                                                                >  2:522155‑M1/1‑113 (MQ=255)                         tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAACAAGGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                                                                                                >  2:704518‑M1/1‑113 (MQ=255)                            aTCGACGCAACTGTACTCGTGACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttca                                                                                                             <  1:286724‑M1/143‑34 (MQ=255)                               gACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtcttgggtgtctttcactc                                                                                                          >  1:607439‑M1/1‑107 (MQ=255)                                aCGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatct                                                                                                         >  1:253524‑M1/1‑106 (MQ=255)                                 cGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctg                                                                                                        <  1:542573‑M1/143‑39 (MQ=255)                                    aaCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggctgtatttcaatctgtta                                                                                                     >  2:221919‑M1/1‑102 (MQ=255)                                       tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaa                                                                                                  >  2:603719‑M1/1‑99 (MQ=255)                                        gTACTCGTCACTACACGCTCGTACAACGGGGGGGGGCTGCGATTTTCAATAATTCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                 <  2:281706‑M1/143‑46 (MQ=255)                                        gTACTCGTCACTACACGCACCTACACCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                 <  1:93585‑M1/143‑46 (MQ=255)                                        gTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGGTTTTCAATAATGCGTGATGCAGATCACACAAAACACGCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                 <  2:740328‑M1/143‑46 (MQ=255)                                        gTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAAAAGCCCGGTGCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                 >  1:174034‑M1/1‑98 (MQ=255)                                            tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtctaagtcatacttttgttttgggtgtatttcaaactgttcaaaagtcg                                                                                             >  2:83344‑M1/1‑94 (MQ=255)                                            tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactcttgttttggttctatttcaatctgtgtaaaagttt                                                                                             >  2:150171‑M1/1‑94 (MQ=255)                                              gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTGGATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                           >  1:321698‑M1/1‑92 (MQ=255)                                              gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggcatatttcaatctctaaaaaagttttt                                                                                           >  2:602866‑M1/1‑92 (MQ=255)                                               acaCTACACACACAAACAACGGAGGGGGGCCGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttc                                                                                          <  1:499601‑M1/142‑53 (MQ=255)                                                cACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcg                                                                                         <  1:54357‑M1/143‑54 (MQ=255)                                                   tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaagaagtgtttcgcta                                                                                      >  1:470134‑M1/1‑87 (MQ=255)                                                    acacGCACATACGACGGAGGGGGGCGGCGATTCTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTATCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctac                                                                                     <  1:308054‑M1/143‑58 (MQ=255)                                                     cacacacaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttgtgggtgtatttcaatctgttaaaaagtttttcgctacg                                                                                    >  2:250630‑M1/1‑85 (MQ=255)                                                       cGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgct                                                                                  <  1:651108‑M1/143‑61 (MQ=255)                                                          acaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtcttgggtgtatttcaatctgttaaaaagtttttcgctacgctagc                                                                               >  1:633703‑M1/1‑80 (MQ=255)                                                               cAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACCCTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagca                                                                          >  2:626669‑M1/1‑75 (MQ=255)                                                                 aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATACATGtataagtcatacttttgtttttggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaa                                                                        >  2:618199‑M1/1‑73 (MQ=255)                                                                   ggAGGGGTGCTGCGAATTTCAAGAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                      <  2:3396‑M1/143‑73 (MQ=255)                                                                    gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                                                                     >  1:574470‑M1/1‑70 (MQ=255)                                                                    gAGGGGGGCTGCGAGTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATATATGtataagtcatacttttgtgttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                                                                     <  2:367650‑M1/143‑74 (MQ=255)                                                                          ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaatcaaaaataaaaca                                                               >  2:192817‑M1/1‑64 (MQ=255)                                                                           gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                              >  2:227419‑M1/1‑63 (MQ=255)                                                                           gCTGCGATTTTCAATAACGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                              <  1:183906‑M1/143‑81 (MQ=255)                                                                                          aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg                                               <  1:536019‑M1/143‑96 (MQ=255)                                                                                                        aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgagaccggaataatcgaaatgggatgttgcgcaccg                                 >  1:538404‑M1/1‑34 (MQ=255)                                                                                                             acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaa                            >  2:552422‑M1/1‑29 (MQ=255)                                                                                                               aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaat                          >  2:721383‑M1/1‑27 (MQ=255)                                                                                                                aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaata                         >  1:404741‑M1/1‑26 (MQ=255)                                                                                                                 aaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                        <  1:543831‑M1/143‑119 (MQ=255)                                                                                                                      tCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcac                   >  1:511502‑M1/1‑20 (MQ=255)                                                                                                                           tACTTAACATAAATGtataagtcatacttttattttgggtgtatttcaatatgataaaaattttttcgctaagctagcaagcaaaaatgaaacaggaataatcgaactgagatcttacgcacagtaaaaaaaactcagcctaa              >  2:343606‑M1/1‑15 (MQ=255)                                                                                                                            aCTTAACATAATTGtataagtcatccttttgtgttgggtgtatttcaatctgttaaaaaggttctcgcgacgctcgcaagcaaaaatgaaacaggaataatcgcaatgggatgttgcgcacagtcaaaataactcaccgtaaa             <  2:194093‑M1/143‑130 (MQ=255)                                                                                                                                   aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcat      >  2:69062‑M1/1‑7 (MQ=255)                                                                                                                                       aTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaccgtcaaaataactcaccgtaaataatcatctgc  <  2:398584‑M1/143‑141 (MQ=255)                                                                                                                                         |                                                                                                                                             TTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCC  >  NC_000913/1978366‑1978642

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑