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	55 [0]	[0] 55	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)	50 (1.330)	38/268	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

CTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGC  >  NC_000913/1978365‑1978641                                                                                                                                          |                                                                                                                                            cTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                        <  2:470875‑M1/143‑6 (MQ=255)            tCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttt                                                                                                                             >  1:645303‑M1/1‑127 (MQ=255)            tCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTCAACATAAATGtataagtcatactttt                                                                                                                             >  1:649087‑M1/1‑127 (MQ=255)              tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACCCGCACATACAACGGAGGGGGGCTGCGATTTTCCATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtatactgcatacttgtgt                                                                                                                           >  2:241137‑M1/1‑125 (MQ=255)               aaGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGCTTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTAGTTAACATAAATGtataagtcatacatctgtt                                                                                                                          >  2:95208‑M1/1‑124 (MQ=255)                   tttttCCTAAATCTGCGCGACTGTACTCGTCACTACTCGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAATCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                      <  2:392055‑M1/143‑24 (MQ=255)                   tttttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                      <  2:356360‑M1/143‑24 (MQ=255)                        ccTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCAAACAAAACACTCAATTACTTAACATAAATGtatacgtcatactattgttttgggtgta                                                                                                                 >  2:697723‑M1/1‑115 (MQ=255)                            aaTCGACGCAACTGTACTCGTCACTACACGCACATACCACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAAAGtataaacaataatgttgatttaggagtataca                                                                                                             >  2:41181‑M1/1‑111 (MQ=255)                               cGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                          <  2:553403‑M1/143‑36 (MQ=255)                                gACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtcttgggtgtatttcaatc                                                                                                         >  2:126672‑M1/1‑107 (MQ=255)                                  cGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCCGATCACACAAAACCCTCATTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctg                                                                                                       <  1:400340‑M1/143‑39 (MQ=255)                                    cAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                                                                     >  1:626794‑M1/1‑103 (MQ=255)                                     aaCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtta                                                                                                    >  1:644206‑M1/1‑102 (MQ=255)                                       cTGTACTCGTCCCTACACGCACATATAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAACACACTCTATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaa                                                                                                  <  1:627502‑M1/143‑44 (MQ=255)                                        tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaa                                                                                                 <  1:53100‑M1/143‑45 (MQ=255)                                             tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCGGATCACACAAAACACTCAATTACTTAACATAAATGtagaagtcatacttttgttttgggtgtatttcaatctgttaaaaagtgt                                                                                            >  1:148627‑M1/1‑94 (MQ=255)                                               gTCACTACACGCAGATACAACGGAGTGAGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                          <  2:372055‑M1/143‑52 (MQ=255)                                                         gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta                                                                                >  1:598834‑M1/1‑82 (MQ=255)                                                                 aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaa                                                                        >  1:277702‑M1/1‑74 (MQ=255)                                                                    ggAGGTGCGCTGCAATTTTCAATAATGTGTGGTGCAGATTACATAAAGCACTCAATTACTTAACATAAATGtattagtcatacttttgtgttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                     <  2:573754‑M1/143‑73 (MQ=255)                                                                    ggAGGGGGGCTGCGAGTTGCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                     <  1:683982‑M1/143‑73 (MQ=255)                                                                     gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaacaat                                                                    >  1:8595‑M1/1‑70 (MQ=255)                                                                     gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcacacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                                                                    >  1:151201‑M1/1‑70 (MQ=255)                                                                               gCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaa                                                          >  2:525793‑M1/1‑60 (MQ=255)                                                                                cGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaat                                                         <  1:66024‑M1/143‑85 (MQ=255)                                                                                   ttttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat                                                      <  1:298480‑M1/143‑88 (MQ=255)                                                                                   ttttAAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat                                                      <  2:212069‑M1/143‑88 (MQ=255)                                                                                     ttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcg                                                    >  2:32606‑M1/1‑54 (MQ=255)                                                                                      tCAATAATGCTTGATGCAGATCACACAAAACACTCAAGTACTTAACATAAATGtataagtcatacttttgttgtgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcga                                                   <  2:476293‑M1/143‑91 (MQ=255)                                                                                        aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataaccgaac                                                 >  1:635244‑M1/1‑51 (MQ=255)                                                                                               cGTGGTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatg                                          >  2:269986‑M1/1‑44 (MQ=255)                                                                                               cGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaacatgaaacaggaatagtcgaaatgggatg                                          >  1:10423‑M1/1‑44 (MQ=255)                                                                                                 tGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgtt                                        >  1:59756‑M1/1‑42 (MQ=255)                                                                                                   aTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcacaaatgaaacaggaataatcgaaatgggatgttgc                                      >  1:128402‑M1/1‑40 (MQ=255)                                                                                                   aTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgagtgtatttcaacctgtacaaaagataatggcttcgttagccagcaaacatgagacaggaataatcgaaatgggatgttgc                                      >  1:157538‑M1/1‑40 (MQ=255)                                                                                                     gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgc                                    >  2:477861‑M1/1‑38 (MQ=255)                                                                                                        gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                 >  2:371739‑M1/1‑35 (MQ=255)                                                                                                        gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgctttgggtgtatttcaatctgttacaaagtttttcgctacgctagcaagcaaaaacgaaacagaaataatcgacatgccatgttccccaca                                 >  2:160539‑M1/1‑35 (MQ=255)                                                                                                         atgacacaaaaCACTCGATTACTTAACATAAATGtataagtcctacttttgtgttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacag                                <  2:26185‑M1/140‑110 (MQ=255)                                                                                                              acaAAACCGTCAAATGCTTAATATTAATGtataaggcatacttgtgtgttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaa                           <  2:553262‑M1/143‑115 (MQ=255)                                                                                                              acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaacatgaaacaggactactcgacttgggctgttgtgcacattcaaa                           >  1:273085‑M1/1‑29 (MQ=255)                                                                                                              acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaa                           >  1:252090‑M1/1‑29 (MQ=255)                                                                                                              acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgctttgggtgtatttcaatctgtttaaaagtttttcgctacgctagccagcaacaatgtaatatgcataatcgaaatggaatgttacgcaaattcata                           >  1:361424‑M1/1‑29 (MQ=255)                                                                                                                   acccacaattaCTGCACATAAATCtataagtcatactgttgctttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataac                      <  2:304350‑M1/138‑120 (MQ=255)                                                                                                                   acacTCAATTACTTAACATAAATGtatatgtcagacttttgcgttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataac                      <  1:351789‑M1/143‑120 (MQ=255)                                                                                                                       tCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcac                  >  1:336577‑M1/1‑20 (MQ=255)                                                                                                                           ttACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgta              >  1:563542‑M1/1‑16 (MQ=255)                                                                                                                               ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaata          >  1:416803‑M1/1‑12 (MQ=255)                                                                                                                                 aaCGTAAGTGtataagtcatacttttgttttgggtgtagttcaatctgttaaaaagtttttcgcgacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataat        <  2:68867‑M1/143‑134 (MQ=255)                                                                                                                                  tcaTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatc       <  1:578079‑M1/142‑135 (MQ=255)                                                                                                                                    aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcat     >  1:142175‑M1/1‑7 (MQ=255)                                                                                                                                    aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatattgcgcacagtcaaaataactcaccgtaaataatcat     >  1:23210‑M1/1‑7 (MQ=255)                                                                                                                                      aaaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaatactcatct   >  1:177650‑M1/1‑5 (MQ=255)                                                                                                                                       aaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctg  >  2:137635‑M1/1‑4 (MQ=255)                                                                                                                                          |                                                                                                                                            CTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGC  >  NC_000913/1978365‑1978641

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