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	1979278	776	43 [0]	[0] 43	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)	40 (1.010)	35/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

GCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATC  >  NC_000913/1979142‑1979415                                                                                                                                          |                                                                                                                                         tttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagataacacaagacactcaattacttaacataaatgTATAAg                                                                                                                                     >  1:707882‑M2/138‑143 (MQ=255)      acttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCatc                                                                                                                                >  2:543693‑M2/133‑142 (MQ=255)             taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgt                                                                                                                         >  2:569606‑M2/126‑143 (MQ=255)             taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgt                                                                                                                         >  1:149350‑M2/126‑143 (MQ=255)                 attgttcctaaatcgacgcacctgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttttg                                                                                                                     <  1:11797‑M2/22‑1 (MQ=255)                  tttttcctaagtcgacgcgactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                    <  1:452709‑M2/23‑1 (MQ=255)                   ttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                                                                                                   <  2:624689‑M2/24‑1 (MQ=255)                    tttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgactttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACATTTGTTTTGGgt                                                                                                                  >  1:455011‑M2/119‑143 (MQ=255)                     ttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtg                                                                                                                 >  1:111155‑M2/118‑143 (MQ=255)                            atcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacgcaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCa                                                                                                          >  1:573301‑M2/111‑143 (MQ=255)                              cgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                                                        <  2:328515‑M2/35‑1 (MQ=255)                               gacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                                                                                                       >  2:586929‑M2/108‑143 (MQ=255)                                  gcaactgtactcgtcactacacgcacatacaacggaggggggctgcgactttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGt                                                                                                    <  2:447200‑M2/39‑1 (MQ=255)                                        gtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaaa                                                                                              <  2:37707‑M2/45‑1 (MQ=255)                                           ctcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgagtgctgcagatcacacaaaacactcaattacttaacataaatgTATAAATCATACTTTTGTTTTGGGTGTATTTCAATCTTTTAAAAAGtt                                                                                           >  1:46457‑M2/96‑143 (MQ=255)                                            tcgtcactacacgcacatacaacggaggggggctgcgcttttcaataatgcgtgatgcagatcacacaaaacactcacttacttaacataaatgTATAAGTCATACTTTTGTTTTGCGTGTATTTCACTctcttcaaacgctt                                                                                          >  1:240386‑M2/95‑138 (MQ=255)                                            tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttt                                                                                          >  1:113696‑M2/95‑143 (MQ=255)                                              gtcactacacgcacatacaacggaggggggctgcgatttccaatgatgcgtgatgcagatcacacaaaactctcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt                                                                                        <  1:548831‑M2/51‑1 (MQ=255)                                                   tacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgcta                                                                                   >  2:479139‑M2/88‑143 (MQ=255)                                                   tacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGCTTTGGTTGTGTTTCGATGTGTTAACAAGGTTTtcgctg                                                                                   >  2:24700‑M2/88‑142 (MQ=255)                                                      acgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTCTTTcgctacgc                                                                                >  1:302890‑M2/85‑143 (MQ=255)                                                        gcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                                                              <  2:149350‑M2/61‑1 (MQ=255)                                                                aacggaggggtgctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTGGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAgcaagcaa                                                                      <  1:531392‑M2/69‑1 (MQ=255)                                                                aacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAgcaagcaa                                                                      >  2:324740‑M2/75‑143 (MQ=255)                                                                 acggaggggggctgcgattttcaataatgcgtgatgaagatcacacaaaacactgaattacttaacataaatgTATAAGTGATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGATTTTCTCTACGCTAGCAAGCaaa                                                                     >  1:19724‑M2/74‑143 (MQ=255)                                                                           gctgcgattttcaagaatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAg                                                           <  2:622576‑M2/80‑1 (MQ=255)                                                                                   ttttgataatgcgtgatgcagatcatacgaaacactcaatttcttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGGTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGgaataatg                                                   <  1:5353‑M2/88‑2 (MQ=255)                                                                                    ttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg                                                  >  1:411598‑M2/55‑143 (MQ=255)                                                                                          aatgcgtgatgcagagtacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATgg                                            <  1:104276‑M2/95‑1 (MQ=255)                                                                                           atgcctgatgcagatcacccaaaacacgcaatgacttaacataaatgTATAAGTCATACTTGTGTGTTGCATGACTTTACATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATggg                                           <  1:185017‑M2/96‑1 (MQ=255)                                                                                            tgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGa                                          <  2:72478‑M2/97‑1 (MQ=255)                                                                                                    gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTccgc                                  >  1:406095‑M2/39‑143 (MQ=255)                                                                                                          cacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTc                            >  2:332637‑M2/33‑143 (MQ=255)                                                                                                             acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTGGCTACGCGATCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaa                         >  2:64860‑M2/30‑143 (MQ=255)                                                                                                                aaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTAAAGATa                      >  2:451935‑M2/27‑143 (MQ=255)                                                                                                                 aacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATaa                     >  2:215168‑M2/26‑143 (MQ=255)                                                                                                                        aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCg              >  2:366207‑M2/19‑143 (MQ=255)                                                                                                                         attacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGTTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGCTGTTGCGAACAGTCAAAATAACTCACAGt             >  1:46513‑M2/18‑143 (MQ=255)                                                                                                                              ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaata        >  1:559887‑M2/13‑143 (MQ=255)                                                                                                                               taatataaacgTATCAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTGCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaataa       <  2:227518‑M2/132‑1 (MQ=255)                                                                                                                                  cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatca    <  2:302890‑M2/135‑1 (MQ=255)                                                                                                                                   ataaatgTATAAGTCAGACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatcat   <  1:469702‑M2/136‑1 (MQ=255)                                                                                                                                    taaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatcatc  >  2:427140‑M2/7‑143 (MQ=255)                                                                                                                                          |                                                                                                                                         GCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATC  >  NC_000913/1979142‑1979415

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