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

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1978503	1979278	776	63 [0]	[0] 62	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)	56 (0.770)	42/260	NT	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

TGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATA  >  NC_000913/1979136‑1979426                                                                                                                                                |                                                                                                                                                    ttggtccttacccttatctaatatttttcctaaatcgacgccactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAg                                                                                                                                                <  2:446809‑M2/6‑1 (MQ=255)      ctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATa                                                                                                                                           >  1:61511‑M2/139‑149 (MQ=255)       tttacacttatctaagatttttcctaaatcgacgcaactgtactcgccactacacgcccatacaacggaggggggctgcgattttcaataatgcgtgctgcagatcacacacaacactccattccttaacctaaatgCACAAgtcagcc                                                                                                                                          >  2:151791‑M2/138‑146 (MQ=255)        ttacacttatctaaggtttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACt                                                                                                                                         >  1:488388‑M2/137‑149 (MQ=255)        ttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACt                                                                                                                                         >  2:1347577‑M2/137‑149 (MQ=255)           cacttatctaagattttgcctaaaccgacgcaactgtactcgtcactacaagcacatacaacggagggcggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACtttt                                                                                                                                      <  1:26378‑M2/16‑1 (MQ=255)               tatctaagatttttcccaaaccgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttt                                                                                                                                  <  1:779670‑M2/20‑1 (MQ=255)                atctaagaattttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataattgTATAAGTCATACttttatttt                                                                                                                                 <  2:919574‑M2/21‑1 (MQ=255)                acctaagatttttcctaaatcgacgcaactgtactcggcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtttt                                                                                                                                 <  2:330889‑M2/21‑1 (MQ=255)                    aagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                                                                             >  1:141909‑M2/125‑149 (MQ=255)                         ttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAtt                                                                                                                        >  1:614766‑M2/120‑149 (MQ=255)                              ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                                                                   >  2:1066307‑M2/115‑149 (MQ=255)                               taaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggccgcggttttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCTATc                                                                                                                  >  2:590414‑M2/114‑149 (MQ=255)                                  atcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGt                                                                                                               >  2:1085202‑M2/111‑149 (MQ=255)                                   tcgaggcaactgtactcgtcactacacgcacatacaacggaggggggcttcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGtt                                                                                                              <  2:925766‑M2/40‑1 (MQ=255)                                        gcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacgtaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaaa                                                                                                         <  1:859245‑M2/45‑1 (MQ=255)                                        gcaactgtactcgtcactacacgcacatacaacggaggggggcagcaattttcaataatgcgtgatgcagatcacacaaaacactcatttacttaacataaatgTATACGTCATACTTTTGTTTTGGGTGTATTTCAATCTTTTaaaaa                                                                                                         <  2:649558‑M2/45‑1 (MQ=255)                                         caactgtactcgtcaccacacgcacagacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAg                                                                                                        <  1:731400‑M2/46‑1 (MQ=255)                                             tgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactccattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtttt                                                                                                    <  1:1061438‑M2/50‑1 (MQ=255)                                                   cgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgcta                                                                                              >  1:329753‑M2/94‑149 (MQ=255)                                                     tcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAGAAAGTTTTTcgctacg                                                                                            >  1:298672‑M2/92‑149 (MQ=255)                                                      cactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgc                                                                                           <  2:870381‑M2/59‑1 (MQ=255)                                                       actacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgct                                                                                          >  2:1129444‑M2/90‑149 (MQ=255)                                                        ctccccgcacagacaacggaggggggcggcgattttcaataatgcgtgatgcagatcacacaaaaccctcaattacttaacataagtgTATAAGTCATACTTTTGTTGTGGGTGTCTTGCAATCTGTTAAAAAGTTTTtcgctacgctg                                                                                         <  2:837265‑M2/61‑2 (MQ=255)                                                        ctacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                                                                         >  1:791381‑M2/89‑149 (MQ=255)                                                        ctacacgcacatacaacggaggggggctgcgattttcaataatgcgtgacgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                                                                         >  2:348738‑M2/89‑149 (MQ=255)                                                               cacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTAGCTACGCTagcaagca                                                                                  >  1:1235532‑M2/82‑149 (MQ=255)                                                                  atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaatactcagttacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTGGCCAGCaaaa                                                                               >  1:266513‑M2/79‑149 (MQ=255)                                                                  atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaa                                                                               >  2:1323903‑M2/79‑149 (MQ=255)                                                                   tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaaa                                                                              >  1:569677‑M2/78‑149 (MQ=255)                                                                   tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaaa                                                                              <  2:61511‑M2/72‑1 (MQ=255)                                                                   tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaaa                                                                              <  1:681678‑M2/72‑1 (MQ=255)                                                                      aacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGa                                                                           >  1:1212359‑M2/75‑149 (MQ=255)                                                                      aacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGa                                                                           <  1:479851‑M2/75‑1 (MQ=255)                                                                        cggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                                                                         >  1:547114‑M2/73‑149 (MQ=255)                                                                                   tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc                                                              >  1:858198‑M2/62‑149 (MQ=255)                                                                                    gcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg                                                             >  1:887279‑M2/61‑149 (MQ=255)                                                                                       attttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa                                                          <  2:657985‑M2/92‑1 (MQ=255)                                                                                       attttcaataatgcgggatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa                                                          <  1:1347577‑M2/92‑1 (MQ=255)                                                                                        ttttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTCTTTTATGCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAAGGAAACAGGAATAATCGAAAt                                                         >  1:1327785‑M2/57‑149 (MQ=255)                                                                                        ttttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt                                                         >  2:341213‑M2/57‑149 (MQ=255)                                                                                        ttttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGGTAAAAAGTTTTTGGCTTCGCTATCAAGCAAAAGGGAAACAGGAATAATCGAAAt                                                         >  1:617525‑M2/57‑149 (MQ=255)                                                                                            caataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGa                                                     <  2:170933‑M2/97‑1 (MQ=255)                                                                                              ataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg                                                   >  1:700029‑M2/51‑149 (MQ=255)                                                                                                aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGtt                                                 >  2:1346373‑M2/49‑149 (MQ=255)                                                                                                        atgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGGTGTACCGCACAg                                         >  1:227304‑M2/41‑149 (MQ=255)                                                                                                        atgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAg                                         >  1:226681‑M2/41‑149 (MQ=255)                                                                                                         tgcagatcacacaaaacacttaaatacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATACTCGAAATGGGATGTTGCGCACAGt                                        <  1:1005694‑M2/110‑1 (MQ=255)                                                                                                           cagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCa                                      >  2:810330‑M2/38‑149 (MQ=255)                                                                                                           cagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGAAAAAATCAAACAGGGATAATCGACATGGGATGTTGGGCACATtcg                                      >  2:1348942‑M2/38‑148 (MQ=255)                                                                                                                 acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCCAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATaa                                <  1:719454‑M2/118‑1 (MQ=255)                                                                                                                  cacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAAc                               >  1:350087‑M2/31‑149 (MQ=255)                                                                                                                           ctcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTaa                      <  1:1170892‑M2/128‑1 (MQ=255)                                                                                                                           ctcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTaa                      <  2:454891‑M2/128‑1 (MQ=255)                                                                                                                              aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAACGGGATGTTGCGCACAGGCGAAATAAGACCCCGTGaata                   >  1:1232117‑M2/19‑149 (MQ=255)                                                                                                                              aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGAGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAGTGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaata                   >  1:515166‑M2/19‑149 (MQ=255)                                                                                                                                    ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatcatc             >  1:505551‑M2/13‑149 (MQ=255)                                                                                                                                       acataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGc          <  2:107315‑M2/140‑1 (MQ=255)                                                                                                                                         ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCta        <  2:700029‑M2/142‑1 (MQ=255)                                                                                                                                               gTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAACTGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATACCTCACCGTAAATAATCATCTGCTataaata  >  1:264271‑M2/2‑149 (MQ=255)                                                                                                                                               gTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAACGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTataaata  <  1:1066307‑M2/148‑1 (MQ=255)                                                                                                                                               gTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAACTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAAGTCACCGTAAATAATCATCTGCTataaata  >  2:845340‑M2/2‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                    TGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATA  >  NC_000913/1979136‑1979426

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