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–1979270	1979278	9–776	59 [0]	[0] 58	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)	51 (1.000)	44/266	0.1	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

TCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGC  >  NC_000913/1978364‑1978644                                                                                                                                           |                                                                                                                                               tcgttacaCTTACCTACGATTTGTCCTAAAGCGACGCAACTGTACTCGTCACTACACGCACCTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtata                                                                                                                                            <  1:185198‑M1/140‑5 (MQ=255) tCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCACTAATGCGTGATGCGGATCACACAAAACACTCCATTACTTAACCTCAATGcata                                                                                                                                            >  2:347068‑M1/1‑139 (MQ=255)        acTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcata                                                                                                                                     >  2:806927‑M1/1‑132 (MQ=255)         cgtATTTAAGATTTTTTCTAACTCGGCGCTACTGTTCTCTTCACTACAGGCACTTCCAGCGATTGGGGTTTGTGATTTTCAAGTCTGCGTGATGCAGATCACACCAAACTCTCACTTGCTTAACACAAATgtataagtcatac                                                                                                                                    <  1:828614/141‑13 (MQ=255)               tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgt                                                                                                                              >  2:753092‑M1/1‑125 (MQ=255)               tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtaaaactcttgt                                                                                                                              >  1:618185‑M1/1‑125 (MQ=255)                   aTTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                                                                          >  2:345454‑M1/1‑121 (MQ=255)                    tttttCCTAAATCGAGGCAACTGTACTGGTCTCTACGCGCTCAGACCACCGATGGAGGCTGCGATTTTCCATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                         <  1:771732‑M1/143‑24 (MQ=255)                           tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAGCATAAATGtataagtcatacttttgttttgggtgtatt                                                                                                                  >  1:281818‑M1/1‑113 (MQ=255)                                     cAAGTGTAGTCGTTACTACACGCACATACGACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                                                                        <  1:646019‑M1/143‑41 (MQ=255)                                      aaCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtaataattttgttttgggtgtatatcaatctgttc                                                                                                       >  1:359360‑M1/1‑102 (MQ=255)                                        cTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaa                                                                                                     >  2:4134‑M1/1‑100 (MQ=255)                                        cTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGAGTGATGCAGATCACACAAGACACTGAATTAATTAAAATAAATGtataagacatactttcgatttgggtgtatttcaatgtataaaa                                                                                                     >  1:592361‑M1/1‑100 (MQ=255)                                              tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttacgggtgtatttcaatctgttaaaaagttt                                                                                               >  2:899488‑M1/1‑94 (MQ=255)                                                gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaacgttctt                                                                                             >  2:734407‑M1/1‑92 (MQ=255)                                                     tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                                                                        <  2:637580‑M1/143‑57 (MQ=255)                                                     tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttaaatccgttaaaaagtttttcgcta                                                                                        >  1:290006‑M1/1‑87 (MQ=255)                                                        acGCACAGACAACGGAAGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACTCTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatgtcaatctgttaaaaagtttttcgctacgc                                                                                     <  1:142729‑M1/143‑60 (MQ=255)                                                          gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttatgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta                                                                                   >  2:545738‑M1/1‑82 (MQ=255)                                                          gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATCtataagtcacaatcttgcgtgtgctttatttcaatttgttaaaaagtttttcgctacgcta                                                                                   >  1:674992‑M1/1‑82 (MQ=255)                                                              aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaa                                                                               >  2:299643‑M1/1‑78 (MQ=255)                                                              aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaa                                                                               <  2:704317‑M1/143‑66 (MQ=255)                                                                 cAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctaggctagcaagca                                                                            >  2:669376‑M1/1‑75 (MQ=255)                                                                 cAACGGAGGGGGGCAGCGAGTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttaaatctgttaaaaagtttgtcactaagctagcaagca                                                                            >  1:477855‑M1/1‑75 (MQ=255)                                                                  aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACCTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctccgctcgcaaccaa                                                                           >  1:632908‑M1/1‑74 (MQ=255)                                                                     ggAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                        >  1:137981‑M1/1‑71 (MQ=255)                                                                       aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatg                                                                      <  1:299643‑M1/143‑75 (MQ=255)                                                                          ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                                                                   <  1:112793‑M1/143‑78 (MQ=255)                                                                          ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                                                                   >  2:752247‑M1/1‑66 (MQ=255)                                                                               tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagga                                                              >  1:875700‑M1/1‑61 (MQ=255)                                                                                     ttgcaatAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatc                                                        <  1:753976‑M1/140‑89 (MQ=255)                                                                                          ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat                                                   >  2:876108‑M1/1‑50 (MQ=255)                                                                                            aatTCGTGAAGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg                                                 <  1:25477‑M1/143‑96 (MQ=255)                                                                                             atGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg                                                <  1:753092‑M1/143‑97 (MQ=255)                                                                                               gCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttcttcgctacccccgcaagcaaaaatgaaacaggaataatcgaaatgggat                                              >  1:625118‑M1/1‑45 (MQ=255)                                                                                                cgggatgCCGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatg                                             <  2:485918‑M1/140‑100 (MQ=255)                                                                                                      gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgc                                       >  2:818399‑M1/1‑38 (MQ=255)                                                                                                       cAGCTCACACAAAACTCTCAATCACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctggtaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgca                                      <  2:616553‑M1/143‑107 (MQ=255)                                                                                                        aGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcac                                     >  2:262994‑M1/1‑36 (MQ=255)                                                                                                          aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacag                                   >  2:675628‑M1/1‑34 (MQ=255)                                                                                                          aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacag                                   <  1:544292‑M1/143‑110 (MQ=255)                                                                                                              cacaAAACACTCAATTACTTAACATAAATGgataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaa                               >  1:118534‑M1/1‑30 (MQ=255)                                                                                                                caAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                             <  2:249280‑M1/143‑116 (MQ=255)                                                                                                                        tCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaattttttcgctacgctcgcaaccaaaaatgaaacaagaatcatcgaaattggattttgcacacagtcaaaacaactcac                     >  1:340022‑M1/1‑20 (MQ=255)                                                                                                                          aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaaccggaataatcgaaatgggatgttgcgcacagtcaaaataactacacg                   >  2:153339‑M1/1‑18 (MQ=255)                                                                                                                          aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccg                   >  2:352456‑M1/1‑18 (MQ=255)                                                                                                                          aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataacccactg                   >  2:210389‑M1/1‑18 (MQ=255)                                                                                                                           aTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaattatgtaacaagaacaatcgaaatgggatgttgcgcacagtcaaaataactcaccgt                  >  2:652573‑M1/1‑17 (MQ=255)                                                                                                                            ttACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgtaatggtatgttgcgcacaatcaaaataactcaccgta                 >  2:156407‑M1/1‑16 (MQ=255)                                                                                                                               cTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaatgactcaccgtaaat              >  2:797584‑M1/1‑13 (MQ=255)                                                                                                                                ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtgaata             >  2:872089‑M1/1‑12 (MQ=255)                                                                                                                                  aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttagctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtaaaaataactcaccgtgaataat           >  2:393990‑M1/1‑10 (MQ=255)                                                                                                                                    cATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatca         <  2:618185‑M1/143‑136 (MQ=255)                                                                                                                                     aTAACTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaaccaaaaatgaaacacgaataatcgaaatgggatgttgctcacagtcaaaataactcaccgtaaataatcat        <  1:502187‑M1/143‑137 (MQ=255)                                                                                                                                     aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcat        >  2:4971‑M1/1‑7 (MQ=255)                                                                                                                                     aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcacagtaaataatcat        >  1:455956‑M1/1‑7 (MQ=255)                                                                                                                                      tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttagctacgctagcacgcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaaagtcaaaagaactcaccgacaaaaatcaac       >  1:845915‑M1/1‑6 (MQ=255)                                                                                                                                        aaTGtataagtcatacttttgttttgggtgtagttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgatgcgcttagtcaaaataactcaccgtaaataatcatctg     <  2:750841‑M1/143‑140 (MQ=255)                                                                                                                                           gtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgcta  <  1:876108‑M1/143‑143 (MQ=255)                                                                                                                                           |                                                                                                                                               TCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGC  >  NC_000913/1978364‑1978644

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: ‑