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

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913_3_bme_pgi	1978503–1979268	1979278	11–776	62 [0]	[0] 62	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)	55 (1.250)	39/256	NT	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

TGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAA  >  NC_000913_3_bme_pgi/1979133‑1979427                                                                                                                                                   |                                                                                                                                                     aaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtat                                                                                                                                                    <  1:631660‑M2/3‑1 (MQ=255)   atgggtctgtggacttatcgaagttgtttcctatatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacataacactcaattacttaacataaatgTATaa                                                                                                                                                  <  2:93063‑M2/5‑1 (MQ=255)          tttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcggacaaaacactcaattacttaacataaatgTATAAGTCATAc                                                                                                                                           >  1:744489‑M2/138‑149 (MQ=255)          tttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATAc                                                                                                                                           >  2:866935‑M2/138‑149 (MQ=255)          tttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATAc                                                                                                                                           >  1:450347‑M2/138‑149 (MQ=255)          tttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggagggctgcgattttcaatactcccccatgcacatcacacaaaacacacaattacttaacataaatgTATAAGTCATAc                                                                                                                                           >  2:317639‑M2/138‑149 (MQ=255)                     ctaatatttttcttaaatcggcgcaactgtactcgtcactacacgcacatacaaagaaggatggatgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATActtttgttttgt                                                                                                                                <  2:680457‑M2/23‑2 (MQ=255)                     ctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacagaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                                <  2:199684‑M2/23‑1 (MQ=255)                      taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggcggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactaaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                                                                                                               >  1:612784‑M2/126‑149 (MQ=255)                      taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTAATACTTTTGTTTTggg                                                                                                                               >  1:463454‑M2/126‑149 (MQ=255)                      taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgagtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTtgg                                                                                                                               >  2:232872‑M2/126‑149 (MQ=255)                          atttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgagtgatgcaaatcacacaaaaaactcaattacttaacataaatgAATAAGCCATAATTTTGTTTCGGGtata                                                                                                                           >  2:655016‑M2/122‑149 (MQ=255)                          atttttcctaaatcgacgcaactgcactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTa                                                                                                                           >  1:779255‑M2/122‑149 (MQ=255)                          atttttcctaaatcaacgcaactgcactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTa                                                                                                                           >  1:779004‑M2/122‑149 (MQ=255)                                  taaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                                                                                                                   <  1:441556‑M2/36‑1 (MQ=255)                                  taaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGATTTGGGTGTATTTCAATc                                                                                                                   >  1:632816‑M2/114‑149 (MQ=255)                                      tcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGtt                                                                                                               >  2:754777‑M2/110‑149 (MQ=255)                                         acgcaactgtactcgtcactacactcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaa                                                                                                            <  1:674125‑M2/43‑1 (MQ=255)                                           tcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTTACACTTTCTTTTTGGGTGTATTTCAATCTGTTaaaaa                                                                                                          <  1:892238‑M2/45‑1 (MQ=255)                                           ggaactgtactcgacactaaacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaaa                                                                                                          <  2:626734‑M2/45‑1 (MQ=255)                                             aactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcagacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTATTAGATGTATTTCAATATGATAAAAAGt                                                                                                        >  2:284672‑M2/103‑149 (MQ=255)                                             aactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGt                                                                                                        >  1:756477‑M2/103‑149 (MQ=255)                                             aactgtactcgtcactacacgcacatacaacggaggagggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTATATTTCAAtctgttaccaacg                                                                                                        >  2:444072‑M2/103‑143 (MQ=255)                                                 gtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt                                                                                                    <  1:819274‑M2/51‑1 (MQ=255)                                                 gtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt                                                                                                    <  1:218656‑M2/51‑1 (MQ=255)                                                    ctcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgc                                                                                                 <  2:868084‑M2/54‑1 (MQ=255)                                                     tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTCTGGCTGTATTACAATATATTAAAAAGCTtttctca                                                                                                >  2:806669‑M2/95‑146 (MQ=255)                                                     tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagataacacaaaagcctcaattacttaacataaatgTATAAATCATACGTTTATTTTGGGTGTATGTCAATCTGTTAAAAAATTTTTcgct                                                                                                >  1:525097‑M2/95‑149 (MQ=255)                                                        tccctacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaactctcaattacttaacattaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacg                                                                                             <  1:270452‑M2/58‑1 (MQ=255)                                                                      tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATATGTTAAAAAGTTTTTCGCTCCGCTAGCAGGCacaaa                                                                               >  2:510132‑M2/78‑149 (MQ=255)                                                                      tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGCTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaaa                                                                               >  2:679316‑M2/78‑149 (MQ=255)                                                                       acaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcg                                                                                                  <  1:621065‑M2/53‑1 (MQ=255)                                                                       acaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcg                                                                                                  >  2:621065‑M2/77‑129 (MQ=255)                                                                          acggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaa                                                                           >  2:124175‑M2/74‑149 (MQ=255)                                                                           cggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                                                                          >  1:446426‑M2/73‑149 (MQ=255)                                                                              aggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAg                                                                       <  1:113132‑M2/80‑1 (MQ=255)                                                                                        cgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGATTTTCGCTACTCTAGCCAGCCAAAATGGAAAAGCAATAatcca                                                             >  2:226552‑M2/60‑147 (MQ=255)                                                                                          attttcaataatgcgtgatgcagatcaaaaaaatcattcaattacttaacataaatgTATAAGTTATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa                                                           <  2:721474‑M2/92‑1 (MQ=255)                                                                                            tttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTAAATATGTTAAAAAGTTTTTAGCTACGATAGCAAGCAAAAATGAAACAGGAATAATAGAAATg                                                         >  1:742661‑M2/56‑149 (MQ=255)                                                                                               caataatgcatgatgaagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGTGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAACGAAACAGGAATAATCGAAATGGGa                                                      <  2:553761‑M2/97‑1 (MQ=255)                                                                                                        gtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGATGTATTTTAATCTGTTAAAAACTTTTTCGCTACGCTAGCAAGCAAAAATAACACAGGAATAATCGAAATGGTATGttgcgac                                             >  2:337074‑M2/44‑147 (MQ=255)                                                                                                          gatgcagattatacaaaatactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAACAGTTTCTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGcaca                                           <  1:849432‑M2/108‑1 (MQ=255)                                                                                                          gatgcagatcacacaaaatacttaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAATTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGcaca                                           <  1:393741‑M2/108‑1 (MQ=255)                                                                                                          gatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGcaca                                           <  2:549031‑M2/108‑1 (MQ=255)                                                                                                              cagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCa                                       <  1:820522‑M2/112‑1 (MQ=255)                                                                                                                gatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTaaaa                                     >  1:86839‑M2/36‑149 (MQ=255)                                                                                                                 atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAAGTTACGCACAGTCaaaa                                    >  2:890852‑M2/35‑149 (MQ=255)                                                                                                                  ttacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAAt                                   <  1:32511‑M2/116‑1 (MQ=255)                                                                                                                      acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACt                               >  2:785851‑M2/30‑149 (MQ=255)                                                                                                                      aaaaaacactcaattacttaacataaatgTATAAGTTATACTTTTGTTTTAGGTGTATTTCAATCTGTTAAAAAGAATTTCGCGACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACt                               <  1:123227‑M2/120‑1 (MQ=255)                                                                                                                         aaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTCCGATAGGAAGCAAAAATTAAAAAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATaacttca                            >  2:802244‑M2/27‑146 (MQ=255)                                                                                                                               gcaattacataaaataaatgAAATAGTCATAGTTTTGTTTTCGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTAGGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTaaa                      <  2:93153‑M2/129‑1 (MQ=255)                                                                                                                                 aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACTCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAAATCACCGTAaata                    >  1:514366‑M2/19‑149 (MQ=255)                                                                                                                                  attacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaataa                   >  1:483245‑M2/18‑149 (MQ=255)                                                                                                                                   ttacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTgc                                                <  2:168327‑M2/103‑1 (MQ=255)                                                                                                                                   ttacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTgc                                                >  1:168327‑M2/17‑119 (MQ=255)                                                                                                                                     acttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAAAAatca                >  1:762385‑M2/15‑149 (MQ=255)                                                                                                                                      ttgacaataaatgAATAAGTTATACTTTTGTTTTGGTTGTATTTCTATCTTTTAAAAATTTTTTTGCTACGCTAGAAAGCCAAAATGAAACAGGAAAAATCGAAAGGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatcat               <  1:450200‑M2/136‑1 (MQ=255)                                                                                                                                      cttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTACCCGTAAACAaacat               >  2:668448‑M2/14‑149 (MQ=255)                                                                                                                                      cttaacataaatgTAAAAGTTACACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCATCGTAAATAatcat               <  2:446426‑M2/136‑1 (MQ=255)                                                                                                                                            ataaatgTATAAGTCATACTTTTTTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGAAAAAATGAAAAAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAAATTACAGTAAATAATCATCTGctc         >  2:534165‑M2/8‑148 (MQ=255)                                                                                                                                                   tataAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATATGCTataaataa  >  1:863391/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     TGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAA  >  NC_000913_3_bme_pgi/1979133‑1979427

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑