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

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913_3_bme_pgi	1978503–1979270	1979278	9–776	63 [0]	[0] 61	insB1–insA	insB1,insA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	= 1978502	0 (0.000)	59 (1.040)	47/282	0.1	100%	intergenic (‑305/+16)	flhD/insB1	flagellar class II regulon transcriptional activator, with FlhC/IS1 transposase B
	?	NC_000913_3_bme_pgi	1979279 =	0 (0.000)					intergenic (‑64/‑474)	insA/uspC	IS1 repressor TnpA/universal stress protein

CTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAAT  >  NC_000913_3_bme_pgi/1979132‑1979425                                                                                                                                                    |                                                                                                                                                   taaatggttctttacacttatctaatatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacgtaggggggctgcgattttcaataatgcgtgatacagatcacacaaaacactcaattacttaacataaatgtat                                                                                                                                                  <  1:531946‑M2/3‑1 (MQ=255)         gtctttacacttatctaagattttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATa                                                                                                                                          <  1:243429‑M2/11‑1 (MQ=255)               cacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttg                                                                                                                                    <  1:912927‑M2/17‑1 (MQ=255)                    atcttagatttttcctaaatcgacgcaactgtactagtcactacacgcacaaacaacggcgggtggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttttg                                                                                                                               <  1:811625‑M2/22‑1 (MQ=255)                    atctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttttg                                                                                                                               <  1:826694‑M2/22‑1 (MQ=255)                      ctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                                                                                                             >  2:306076‑M2/127‑150 (MQ=255)                      ctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                                                                                                             >  1:225667‑M2/127‑150 (MQ=255)                       taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                                                                            >  1:753621‑M2/126‑150 (MQ=255)                          gatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaataacttaacataaatgTATAATTCATACTTTTGTTTTGGGTGTa                                                                                                                         >  2:573418‑M2/123‑150 (MQ=255)                            tttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAtt                                                                                                                       <  1:27970‑M2/30‑1 (MQ=255)                              tttcctaactcgacgcaactggactcgtcactacacgcacatagaacggaggggggctgggattttcaataatgcgggatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTc                                                                                                                     <  2:31013‑M2/32‑1 (MQ=255)                               ttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCa                                                                                                                    >  1:731162‑M2/118‑150 (MQ=255)                               ttcctaaatcgacgcaactgtactcggcactacacgcacatacaacggaggggggctgagattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCa                                                                                                                    <  1:200107‑M2/33‑1 (MQ=255)                                   taaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCt                                                                                                                <  1:452911‑M2/37‑1 (MQ=255)                                    aaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTg                                                                                                               <  2:796373‑M2/38‑1 (MQ=255)                                     aatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGt                                                                                                              >  2:160135‑M2/112‑150 (MQ=255)                                              aactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtt                                                                                                     <  2:489683‑M2/48‑1 (MQ=255)                                                ctgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtttt                                                                                                   >  1:318895‑M2/101‑150 (MQ=255)                                                  ttactcgtcactacacgcacatacaacggatgggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcacttacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTc                                                                                                 <  2:154030‑M2/52‑1 (MQ=255)                                                      tcgtcactacacgcacatacaacggaggggggctgagattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgcta                                                                                             <  1:13876‑M2/56‑1 (MQ=255)                                                         tcactacacgcacatacaacggaggggggctgcgattatcaataatccgtgatgcagatcacacaagacactcaattacttaacataaatgTATAAGTCATACTGTTGTTTAGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgc                                                                                          >  2:134816‑M2/92‑150 (MQ=255)                                                           aatacacgaacatacaacggaggggggtggcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaagataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                                                                        <  2:537307‑M2/61‑1 (MQ=255)                                                            ctacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTag                                                                                       <  1:3308‑M2/62‑1 (MQ=255)                                                            ctacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTag                                                                                       <  1:206637‑M2/62‑1 (MQ=255)                                                                 cgcacatacaacggagggtggctgcgattttcaataatgcgtgatgcagatcacacaaaagactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGATACGCTagcaagc                                                                                  <  2:865431‑M2/67‑1 (MQ=255)                                                                 cgcacatacaacggaggggggctgcgattttcaataatgcttgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaagc                                                                                  <  2:68899‑M2/67‑1 (MQ=255)                                                                   cacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAgcaagcaa                                                                                >  1:362863‑M2/82‑150 (MQ=255)                                                                    acatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaa                                                                               <  2:225667‑M2/70‑1 (MQ=255)                                                                        acaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacacccaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATg                                                                           <  2:471815‑M2/74‑1 (MQ=255)                                                                          aacggaggggggctgcgattttcaataatgcgtgatccagatcaaacaaaacactcaattacttaacataaagcTATAACTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGAATTTCGCTACGAAAGcaagcaaaattcaa                                                                         >  1:914563‑M2/75‑145 (MQ=255)                                                                           acggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACTCTAGCAAGCAAAAATGaaa                                                                        >  1:217442‑M2/74‑150 (MQ=255)                                                                                  ggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaat                                                                 >  1:404317‑M2/67‑150 (MQ=255)                                                                                      ctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc                                                             <  1:315510‑M2/88‑1 (MQ=255)                                                                                      ctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAACCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc                                                             <  1:97262‑M2/88‑1 (MQ=255)                                                                                       tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg                                                            <  1:160135‑M2/89‑1 (MQ=255)                                                                                         cgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaa                                                          <  1:129774‑M2/91‑1 (MQ=255)                                                                                          gattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa                                                         >  1:278352‑M2/59‑150 (MQ=255)                                                                                             tttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATgg                                                      >  1:595055‑M2/56‑150 (MQ=255)                                                                                              ttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATggg                                                     >  1:1098‑M2/55‑150 (MQ=255)                                                                                               tcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGa                                                    <  2:731162‑M2/97‑1 (MQ=255)                                                                                                    aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTg                                               >  1:783662‑M2/49‑150 (MQ=255)                                                                                                    aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTg                                               >  2:625869‑M2/49‑150 (MQ=255)                                                                                                      tgcgtgatgcagatcacacaaaacactcacttacataacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTgcg                                             <  2:788963‑M2/104‑1 (MQ=255)                                                                                                              gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTTCGCACAGTCa                                     >  2:872576‑M2/39‑150 (MQ=255)                                                                                                              gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCa                                     >  2:376361‑M2/39‑150 (MQ=255)                                                                                                              gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCa                                     >  2:628028‑M2/39‑150 (MQ=255)                                                                                                                 gatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCCATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaaa                                  >  2:683342‑M2/36‑150 (MQ=255)                                                                                                                 gatcacacaaaacactcaattacttaacataaatgTATAAGTCATAATTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAAACGAAATGGGATGTTGCGCACAGTCaaaa                                  >  2:713422‑M2/36‑150 (MQ=255)                                                                                                                   tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATa                                <  1:347595‑M2/117‑1 (MQ=255)                                                                                                                           aacactctattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGAGCACAGTCAAAATAACTCACCg                        <  2:616021‑M2/125‑1 (MQ=255)                                                                                                                            acactcatttacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGt                       >  2:312600‑M2/25‑150 (MQ=255)                                                                                                                               ctcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTaaa                    <  1:207577‑M2/129‑1 (MQ=255)                                                                                                                                 caattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaata                  >  1:677415‑M2/20‑150 (MQ=255)                                                                                                                                 caattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaata                  >  1:714112‑M2/20‑150 (MQ=255)                                                                                                                                   attacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaataat                >  1:165358‑M2/18‑150 (MQ=255)                                                                                                                                        ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCt           <  1:336849‑M2/138‑1 (MQ=255)                                                                                                                                        ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCt           <  1:547706‑M2/138‑1 (MQ=255)                                                                                                                                        ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCt           <  1:306076‑M2/138‑1 (MQ=255)                                                                                                                                          aacataaatgTATAAGTCATAATTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGc         <  2:248661‑M2/140‑1 (MQ=255)                                                                                                                                            cattaatgTATAAGTCATACGTTTGTTTTGGGTGTAGTTCACTCTGTTAAAAAGTTTTTCGCTAAGCTAGCAAGCAAAATTGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAGCTCACCGTATATAATCATCTGCta       >  1:800080‑M2/9‑150 (MQ=255)                                                                                                                                                 atgTATAAGTCATACTTTTGTTTTGGGTGTATTTCTATGTGTTAAAAAGTTTTTATCTACTCTAGCAAGCAAAAAAGACAAAGGAATAATCGAAATTGTATGTTGCGCACAGTCAAATTAAATCACAGTAAATAATCATCTGCTataaat  <  2:914563‑M2/147‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAAT  >  NC_000913_3_bme_pgi/1979132‑1979425

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑