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	63 [0]	[0] 64	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)	58 (1.160)	43/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

GTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTAT  >  NC_000913/1979137‑1979421                                                                                                                                               |                                                                                                                                               gggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcggatcacacttcaccctcatttacttagcataaatgt                                                                                                                                                >  1:250503‑M2/143‑143 (MQ=255)   gtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtat                                                                                                                                              <  2:845771‑M2/3‑1 (MQ=255)   gtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtat                                                                                                                                              >  1:837091‑M2/141‑143 (MQ=255)           acttatctaagatttttcctacatcgacgcaaccgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacctagatgTCTCAGTCatc                                                                                                                                      >  2:180299‑M2/133‑142 (MQ=255)           acttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacgacggaggggggctgcggttttcaataatgcgggatgcagatcacacagaacactcaattactaaacatgaatgTATAAGTCATa                                                                                                                                      >  2:528137‑M2/133‑143 (MQ=255)           acttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATa                                                                                                                                      >  1:77869‑M2/133‑143 (MQ=255)             ttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataactgTATAAGTCATACt                                                                                                                                    >  1:157832‑M2/131‑143 (MQ=255)                tctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttcacataaatgTATAAGTCATACtttt                                                                                                                                 >  1:735367‑M2/128‑143 (MQ=255)                  taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgt                                                                                                                               >  1:683866‑M2/126‑143 (MQ=255)                  taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgt                                                                                                                               >  1:729716‑M2/126‑143 (MQ=255)                  taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgt                                                                                                                               >  1:901080‑M2/126‑143 (MQ=255)                   aagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacgagacactcaattaattaacataaatgTATAAGTCATACttttgtt                                                                                                                              >  2:9050‑M2/125‑143 (MQ=255)                   aagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATactcttgat                                                                                                                              >  2:342210‑M2/125‑141 (MQ=255)                       tttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggcgggggtctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATActttctttttgc                                                                                                                          >  2:518569‑M2/121‑142 (MQ=255)                       tttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaacacacccaatttcgttacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                          <  1:654218‑M2/23‑1 (MQ=255)                       tttttcccaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                          <  1:408631‑M2/23‑1 (MQ=255)                       ttgttcttagatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                          <  2:485094‑M2/23‑1 (MQ=255)                        ttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATCAGTCATACTTCTGTTTTggg                                                                                                                         >  1:465084‑M2/120‑143 (MQ=255)                             ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAt                                                                                                                    >  1:236230‑M2/115‑143 (MQ=255)                             ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACATTTGTTTTGGGTGTAt                                                                                                                    >  2:940108‑M2/115‑143 (MQ=255)                               aaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAttt                                                                                                                  <  1:720711‑M2/31‑1 (MQ=255)                                   cgacgctactgtactcgtcactacatgcacctacaacggaggggggctgggattttcaataatgcgtgatgcagatcacagaaaacaatcaactatgtatcatgagggTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                                                              <  1:404188‑M2/35‑1 (MQ=255)                                           ctgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaccaaccaatgacttaacatagatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaa                                                                                                      >  1:130571‑M2/101‑143 (MQ=255)                                               actcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGt                                                                                                  >  2:741213‑M2/97‑143 (MQ=255)                                                 tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGCATTGGTTGTATTTCAATCTGTTAAAAAGttt                                                                                                >  2:8452‑M2/95‑143 (MQ=255)                                                 tcggcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAtttt                                                                                                <  1:315532‑M2/49‑1 (MQ=255)                                                  cgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtttt                                                                                               >  1:526386‑M2/94‑143 (MQ=255)                                                             gcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgttacgcta                                                                                    >  1:272166‑M2/83‑143 (MQ=255)                                                               acatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagc                                                                                  >  1:215289‑M2/81‑143 (MQ=255)                                                                 atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaa                                                                                <  1:698149‑M2/65‑1 (MQ=255)                                                                  tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGTTagcgag                                                                               >  1:936627‑M2/78‑143 (MQ=255)                                                                       cggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaa                                                                          <  1:103718‑M2/71‑1 (MQ=255)                                                                        ggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataactgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaaa                                                                         <  2:658691‑M2/72‑1 (MQ=255)                                                                          aggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATg                                                                       <  2:662260‑M2/74‑1 (MQ=255)                                                                                gctgcgattttcactaatgcgtgatgcagatcaaacaaaactatcaattacttatcattattgTATAAGTTATTTTTTTGTGTTGGGTTTATTTCAATCTGTTAAAAAGTTGTTCGCTACGCTAGCAAGCAAAAATGAAACAg                                                                 <  1:384898‑M2/80‑1 (MQ=255)                                                                                gctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAg                                                                 >  1:535842‑M2/64‑143 (MQ=255)                                                                                  tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGa                                                               >  1:110690‑M2/62‑143 (MQ=255)                                                                                  tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGa                                                               >  1:190686‑M2/62‑143 (MQ=255)                                                                                     gattttcaataatgcgtgatgcaaatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaata                                                            <  2:405561‑M2/85‑1 (MQ=255)                                                                                         ttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGATACGCTAGCAAGCAAAAATGAAACAGGAATAATCg                                                        >  2:824427‑M2/55‑143 (MQ=255)                                                                                            aataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa                                                     <  2:735367‑M2/92‑1 (MQ=255)                                                                                             ataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt                                                    <  1:518569‑M2/93‑1 (MQ=255)                                                                                               aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATCAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATgg                                                  >  1:86350‑M2/49‑143 (MQ=255)                                                                                               aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGCAATgg                                                  >  1:506184‑M2/49‑143 (MQ=255)                                                                                                       atgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGCATAATCGAAATGGTATGTTgc                                          >  1:589065‑M2/41‑143 (MQ=255)                                                                                                       atgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGCAATGGGATGTTgc                                          >  2:227041‑M2/41‑143 (MQ=255)                                                                                                            gatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGcaca                                     >  2:212219‑M2/36‑143 (MQ=255)                                                                                                            gatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGcaca                                     <  2:837091‑M2/108‑1 (MQ=255)                                                                                                             atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAATATGAAACAGGAATAATCGAAATGGGATGTTGCGAACAg                                    >  2:716978‑M2/35‑143 (MQ=255)                                                                                                             atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGACACAGGAATAATCGAAATGGGATGTTGCGCACAg                                    >  2:51921‑M2/35‑143 (MQ=255)                                                                                                             atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAg                                    >  2:464957‑M2/35‑143 (MQ=255)                                                                                                             atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAg                                    >  1:67590‑M2/35‑143 (MQ=255)                                                                                                             atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATAAAGCAGGAATAATCGCAATGGGAtgttgcgccacg                                    >  2:57183‑M2/35‑139 (MQ=255)                                                                                                                acacaaaacactccattactgaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAGAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCa                                 <  2:175368‑M2/112‑1 (MQ=255)                                                                                                                    aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAAt                             <  2:77869‑M2/116‑1 (MQ=255)                                                                                                                                tacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGTTGAAGGGAAACAGGAATAATCGAAATGGGATGTGGCGCACAGTCAAAATAAATCAGCGTaa                 >  1:743047‑M2/16‑143 (MQ=255)                                                                                                                                   ttaacataaatgTATAAGTCATACTTTAGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTAACCGTAaata              >  1:930337‑M2/13‑143 (MQ=255)                                                                                                                                    taacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaataa             >  2:434629‑M2/12‑143 (MQ=255)                                                                                                                                           aatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTg      >  2:616176‑M2/5‑143 (MQ=255)                                                                                                                                           aatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTg      >  2:375721‑M2/5‑143 (MQ=255)                                                                                                                                           aatgTATAAATTATATTTTTGTTTTGGGTTTATTTCAATCTGTTAAAAAGTTTTTCGCTACTCAAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTg      <  2:766189‑M2/139‑1 (MQ=255)                                                                                                                                            atgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTaaaaaataattagc     >  1:640184‑M2/4‑139 (MQ=255)                                                                                                                                             tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCt    <  2:361798‑M2/141‑1 (MQ=255)                                                                                                                                               tataAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCtat  >  1:641869/1‑143 (MQ=255)                                                                                                                                               |                                                                                                                                               GTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTAT  >  NC_000913/1979137‑1979421

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