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	54 [0]	[0] 54	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.350)	37/268	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

GGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTG  >  NC_000913/1979141‑1979417                                                                                                                                           |                                                                                                                                           ctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATaa                                                                                                                                        <  2:470875‑M2/5‑1 (MQ=255)            tctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACtttt                                                                                                                             >  1:645303‑M2/128‑143 (MQ=255)            tctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattactcaacataaatgTATAAGTCATACtttt                                                                                                                             >  1:649087‑M2/128‑143 (MQ=255)              taagatttttcctaaatcgacgcaactgtactcgtcactacccgcacatacaacggaggggggctgcgattttccataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATACTGCATACttgtgt                                                                                                                           >  2:241137‑M2/126‑143 (MQ=255)               aagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgcttttcaataatgcgtgatgcagatcacacaaaacactcaattagttaacataaatgTATAAGTCATACatctgtt                                                                                                                          >  2:95208‑M2/125‑143 (MQ=255)                   tttttcctaaatctgcgcgactgtactcgtcactactcgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaatcactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                      <  2:392055‑M2/23‑1 (MQ=255)                   tttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                      <  2:356360‑M2/23‑1 (MQ=255)                        cctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcaaacaaaacactcaattacttaacataaatgTATACGTCATACTATTGTTTTGGGTGTa                                                                                                                 >  2:697723‑M2/116‑143 (MQ=255)                               cgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                                                          <  2:553403‑M2/35‑1 (MQ=255)                                gacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTCTTGGGTGTATTTCAATc                                                                                                         >  2:126672‑M2/108‑143 (MQ=255)                                  cgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgccgatcacacaaaaccctcatttacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTg                                                                                                       <  1:400340‑M2/38‑1 (MQ=255)                                    caactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGtt                                                                                                     >  1:626794‑M2/104‑143 (MQ=255)                                     aactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTa                                                                                                    >  1:644206‑M2/103‑143 (MQ=255)                                       ctgtactcgtccctacacgcacatataacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaacacactctattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaa                                                                                                  <  1:627502‑M2/43‑1 (MQ=255)                                        tgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaa                                                                                                 <  1:53100‑M2/44‑1 (MQ=255)                                             tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcggatcacacaaaacactcaattacttaacataaatgTAGAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAagtgt                                                                                            >  1:148627‑M2/95‑141 (MQ=255)                                               gtcactacacgcagatacaacggagtgaggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt                                                                                          <  2:372055‑M2/51‑1 (MQ=255)                                                         gcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                                                                >  1:598834‑M2/83‑143 (MQ=255)                                                                 aacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAgcaagcaa                                                                        >  1:277702‑M2/75‑143 (MQ=255)                                                                    ggaggtgcgctgcaattttcaataatgtgtggtgcagattacataaagcactcaattacttaacataaatgTATTAGTCATACTTTTGTGTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaaa                                                                     <  2:573754‑M2/72‑1 (MQ=255)                                                                    ggaggggggctgcgagttgcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaaa                                                                     <  1:683982‑M2/72‑1 (MQ=255)                                                                     gaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAACAAt                                                                    >  1:8595‑M2/71‑143 (MQ=255)                                                                     gaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCACACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAAt                                                                    >  1:151201‑M2/71‑143 (MQ=255)                                                                               gcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaa                                                          >  2:525793‑M2/61‑143 (MQ=255)                                                                                cgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaat                                                         <  1:66024‑M2/84‑1 (MQ=255)                                                                                   ttttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat                                                      <  1:298480‑M2/87‑1 (MQ=255)                                                                                   ttttaaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat                                                      <  2:212069‑M2/87‑1 (MQ=255)                                                                                     ttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg                                                    >  2:32606‑M2/55‑143 (MQ=255)                                                                                      tcaataatgcttgatgcagatcacacaaaacactcaagtacttaacataaatgTATAAGTCATACTTTTGTTGTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGa                                                   <  2:476293‑M2/90‑1 (MQ=255)                                                                                        aataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAACCgaac                                                 >  1:635244‑M2/52‑142 (MQ=255)                                                                                               cgtggtgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg                                          >  2:269986‑M2/45‑143 (MQ=255)                                                                                               cgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAACATGAAACAGGAATAGTCGAAATGGGATg                                          >  1:10423‑M2/45‑143 (MQ=255)                                                                                                 tgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGtt                                        >  1:59756‑M2/43‑143 (MQ=255)                                                                                                   atgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCACAAATGAAACAGGAATAATCGAAATGGGATGTTgc                                      >  1:128402‑M2/41‑143 (MQ=255)                                                                                                   atgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGAGTGTATTTCAACCTGTACAAAAGATAATGGCTTCGTTAGCCAGCAAACATGAGACAGGAATAATCGAAATGGGATGTTgc                                      >  1:157538‑M2/41‑143 (MQ=255)                                                                                                     gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTgcgc                                    >  2:477861‑M2/39‑143 (MQ=255)                                                                                                        gatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGcaca                                 >  2:371739‑M2/36‑143 (MQ=255)                                                                                                        gatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGCTTTGGGTGTATTTCAATCTGTTACAAAGTTTTTCGCTACGCTAGCAAGCAAAAACGAAACAGAAATAATCGACATGCCATGTTCCCcaca                                 >  2:160539‑M2/36‑143 (MQ=255)                                                                                                         atgacacaaaacactcgattacttaacataaatgTATAAGTCCTACTTTTGTGTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAg                                <  2:26185‑M2/109‑1 (MQ=255)                                                                                                              acaaaaccgtcaaatgcttaatattaatgTATAAGGCATACTTGTGTGTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaa                           <  2:553262‑M2/114‑1 (MQ=255)                                                                                                              acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAACATGAAACAGGACTACTCGACTTGGGCTGTTGTGCACATTCaaa                           >  1:273085‑M2/30‑143 (MQ=255)                                                                                                              acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaa                           >  1:252090‑M2/30‑143 (MQ=255)                                                                                                              acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGCTTTGGGTGTATTTCAATCTGTTTAAAAGTTTTTCGCTACGCTAGCCAGCAACAATGTAATATGCATAATCGAAATGGAATGTTACGCAAATtcata                           >  1:361424‑M2/30‑141 (MQ=255)                                                                                                                   acccacaattactgcacataaatcTATAAGTCATACTGTTGCTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAAc                      <  2:304350‑M2/119‑1 (MQ=255)                                                                                                                   acactcaattacttaacataaatgTATATGTCAGACTTTTGCGTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAAc                      <  1:351789‑M2/119‑1 (MQ=255)                                                                                                                       tcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCAc                  >  1:336577‑M2/21‑143 (MQ=255)                                                                                                                           ttacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTa              >  1:563542‑M2/17‑143 (MQ=255)                                                                                                                               ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaata          >  1:416803‑M2/13‑143 (MQ=255)                                                                                                                                 aacgtaagtgTATAAGTCATACTTTTGTTTTGGGTGTAGTTCAATCTGTTAAAAAGTTTTTCGCGACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaataat        <  2:68867‑M2/133‑1 (MQ=255)                                                                                                                                  tcataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatc       <  1:578079‑M2/134‑1 (MQ=255)                                                                                                                                    ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatcat     >  1:142175‑M2/8‑143 (MQ=255)                                                                                                                                    ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATATTGCGCACAGTCAAAATAACTCACCGTAAATAatcat     >  1:23210‑M2/8‑143 (MQ=255)                                                                                                                                      aaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATACTCATCt   >  1:177650‑M2/6‑143 (MQ=255)                                                                                                                                       aatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTg  >  2:137635‑M2/5‑143 (MQ=255)                                                                                                                                           |                                                                                                                                           GGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTG  >  NC_000913/1979141‑1979417

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