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

CGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTA  >  NC_000913/1979140‑1979420                                                                                                                                            |                                                                                                                                              tctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcactaatgcgtgatgcggatcacacaaaacactccattacttaacctcaatgcata                                                                                                                                            >  2:347068‑M2/140‑143 (MQ=255) tcgttacacttacctacgatttgtcctaaagcgacgcaactgtactcgtcactacacgcacctacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtata                                                                                                                                            <  1:185198‑M2/4‑1 (MQ=255)        acttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATa                                                                                                                                     >  2:806927‑M2/133‑143 (MQ=255)               taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgt                                                                                                                              >  2:753092‑M2/126‑143 (MQ=255)               taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTAAAACtcttgt                                                                                                                              >  1:618185‑M2/126‑143 (MQ=255)                   atttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttttg                                                                                                                          >  2:345454‑M2/122‑143 (MQ=255)                    tttttcctaaatcgaggcaactgtactggtctctacgcgctcagaccaccgatggaggctgcgattttccataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                         <  1:771732‑M2/23‑1 (MQ=255)                           taaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttagcataaatgTATAAGTCATACTTTTGTTTTGGGTGTAtt                                                                                                                  >  1:281818‑M2/114‑143 (MQ=255)                                     caagtgtagtcgttactacacgcacatacgacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGtt                                                                                                        <  1:646019‑M2/40‑1 (MQ=255)                                      aactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTAATAATTTTGTTTTGGGTGTATATCAATCTgttc                                                                                                       >  1:359360‑M2/103‑142 (MQ=255)                                        ctgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaa                                                                                                     >  2:4134‑M2/101‑143 (MQ=255)                                        ctgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgagtgatgcagatcacacaagacactgaattaattaaaataaatgTATAAGACATACTTTCGATTTGGGTGtatttcaatgtataaaa                                                                                                     >  1:592361‑M2/101‑135 (MQ=255)                                              tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTACGGGTGTATTTCAATCTGTTAAAAAGttt                                                                                               >  2:899488‑M2/95‑143 (MQ=255)                                                gtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAACGttctt                                                                                             >  2:734407‑M2/93‑143 (MQ=255)                                                     tacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgcta                                                                                        <  2:637580‑M2/56‑1 (MQ=255)                                                     tacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTAAATCCGTTAAAAAGTTTTTcgcta                                                                                        >  1:290006‑M2/88‑143 (MQ=255)                                                        acgcacagacaacggaagggggctgcgattttcaataatgcgtgatgcagatcacacaaaactctcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATGTCAATCTGTTAAAAAGTTTTTcgctacgc                                                                                     <  1:142729‑M2/59‑1 (MQ=255)                                                          gcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTATGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                                                                   >  2:545738‑M2/83‑143 (MQ=255)                                                          gcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatcTATAAGTCACAATCTTGCGTGTGCTTTATTTCAATTTGTTAAAAAGTTTTTcgctacgcta                                                                                   >  1:674992‑M2/83‑143 (MQ=255)                                                              atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaa                                                                               >  2:299643‑M2/79‑143 (MQ=255)                                                              atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaa                                                                               <  2:704317‑M2/65‑1 (MQ=255)                                                                 caacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTAGGCTagcaagca                                                                            >  2:669376‑M2/76‑143 (MQ=255)                                                                 caacggaggggggcagcgagtttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTAAATCTGTTAAAAAGTTTGTCACTAAGCTagcaagca                                                                            >  1:477855‑M2/76‑143 (MQ=255)                                                                  aacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacctaaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTCCGCTCgcaaccaa                                                                           >  1:632908‑M2/75‑143 (MQ=255)                                                                     ggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaaa                                                                        >  1:137981‑M2/72‑143 (MQ=255)                                                                       aggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATg                                                                      <  1:299643‑M2/74‑1 (MQ=255)                                                                          ggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                                                                   <  1:112793‑M2/77‑1 (MQ=255)                                                                          ggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                                                                   >  2:752247‑M2/67‑143 (MQ=255)                                                                               tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGa                                                              >  1:875700‑M2/62‑143 (MQ=255)                                                                                     ttgcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc                                                        <  1:753976‑M2/88‑1 (MQ=255)                                                                                          ataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt                                                   >  2:876108‑M2/51‑143 (MQ=255)                                                                                            aattcgtgaagcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATgg                                                 <  1:25477‑M2/95‑1 (MQ=255)                                                                                             atgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATggg                                                <  1:753092‑M2/96‑1 (MQ=255)                                                                                               gcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTCTTCGCTACCCCCGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt                                              >  1:625118‑M2/46‑143 (MQ=255)                                                                                                cgggatgccgatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg                                             <  2:485918‑M2/99‑1 (MQ=255)                                                                                                      gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTgcgc                                       >  2:818399‑M2/39‑143 (MQ=255)                                                                                                       cagctcacacaaaactctcaatcacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGGTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGca                                      <  2:616553‑M2/106‑1 (MQ=255)                                                                                                        agatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGcac                                     >  2:262994‑M2/37‑143 (MQ=255)                                                                                                          atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAg                                   >  2:675628‑M2/35‑143 (MQ=255)                                                                                                          atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAg                                   <  1:544292‑M2/109‑1 (MQ=255)                                                                                                              cacaaaacactcaattacttaacataaatgGATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaa                               >  1:118534‑M2/31‑143 (MQ=255)                                                                                                                caaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaaa                             <  2:249280‑M2/115‑1 (MQ=255)                                                                                                                        tcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAATTTTTTCGCTACGCTCGCAACCAAAAATGAAACAAGAATCATCGAAATTGGATTTTGCACACAGTCAAAACAACTCAc                     >  1:340022‑M2/21‑143 (MQ=255)                                                                                                                          aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACCGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAataactacacg                   >  2:153339‑M2/19‑138 (MQ=255)                                                                                                                          aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCg                   >  2:352456‑M2/19‑143 (MQ=255)                                                                                                                          aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACCcactg                   >  2:210389‑M2/19‑141 (MQ=255)                                                                                                                           attacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAATTATGTAACAAGAACAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGt                  >  2:652573‑M2/18‑143 (MQ=255)                                                                                                                            ttacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGTAATGGTATGTTGCGCACAATCAAAATAACTCACCGTa                 >  2:156407‑M2/17‑143 (MQ=255)                                                                                                                               cttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATGACTCACCGTAaat              >  2:797584‑M2/14‑143 (MQ=255)                                                                                                                                ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTGaata             >  2:872089‑M2/13‑143 (MQ=255)                                                                                                                                  aacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTAGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTAAAAATAACTCACCGTGaataat           >  2:393990‑M2/11‑143 (MQ=255)                                                                                                                                    cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatca         <  2:618185‑M2/135‑1 (MQ=255)                                                                                                                                     ataactgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAACCAAAAATGAAACACGAATAATCGAAATGGGATGTTGCTCACAGTCAAAATAACTCACCGTAAATAatcat        <  1:502187‑M2/136‑1 (MQ=255)                                                                                                                                     ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatcat        >  2:4971‑M2/8‑143 (MQ=255)                                                                                                                                     ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACAGTAAATAatcat        >  1:455956‑M2/8‑143 (MQ=255)                                                                                                                                      taaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTAGCTACGCTAGCACGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCAAAGTCAAAAGAACTCACCGACAAaaatcaac       >  1:845915‑M2/7‑141 (MQ=255)                                                                                                                                        aatgTATAAGTCATACTTTTGTTTTGGGTGTAGTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGATGCGCTTAGTCAAAATAACTCACCGTAAATAATCATCTg     <  2:750841‑M2/139‑1 (MQ=255)                                                                                                                                           gTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTa  <  1:876108‑M2/142‑1 (MQ=255)                                                                                                                                            |                                                                                                                                              CGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTA  >  NC_000913/1979140‑1979420

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