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	1979278	776	41 [2]	[1] 42	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	2 (0.040)	36 (0.860)	32/270	0.3	97.4%	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

CCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAA  >  NC_000913_3_bme_pgi/1979139‑1979427                                                                                                                                             |                                                                                                                                                     gtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacatcaatgTATAAGTCa                                                                                                                                              >  2:422120‑M2/141‑149 (MQ=255)        cacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACtttt                                                                                                                                       <  1:494976‑M2/16‑1 (MQ=255)               ctaagattgttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaatcacttaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                                <  2:303636‑M2/23‑1 (MQ=255)                taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTAGATTTggg                                                                                                                               >  1:559341‑M2/126‑149 (MQ=255)                         tcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCa                                                                                                                      <  2:699485‑M2/33‑1 (MQ=255)                           ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaatcacgcgtgaggcagagcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                                                                    >  2:510034‑M2/115‑149 (MQ=255)                           ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagctcaaacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                                                                    >  1:372055‑M2/115‑149 (MQ=255)                           ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgatttacaataatgcgtgatgcagatcacacaaaacactcaattacttaaaataaatgTATAAGTCATAGTTTTGTTTTGGTGGTATTTCAAt                                                                                                                    >  2:320260‑M2/115‑149 (MQ=255)                            taaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                                                                                                                   >  2:635684‑M2/114‑149 (MQ=255)                             aaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactaaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCt                                                                                                                  <  1:411896‑M2/37‑1 (MQ=255)                               ctcgacgcaactctgctctccacaacacgcacatacaactgaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGt                                                                                                                <  2:411545‑M2/39‑1 (MQ=255)                                tcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAATCATACTTTTGTTTTGGGTGTATCTCAATCTGtt                                                                                                               >  2:727059‑M2/110‑149 (MQ=255)                                            tactcgtcactacacgcacatacaacggaggggggctgcgagtttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTc                                                                                                   >  2:709219‑M2/98‑149 (MQ=255)                                               tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgct                                                                                                >  1:554377‑M2/95‑149 (MQ=255)                                                cgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgcta                                                                                               <  1:24803‑M2/56‑1 (MQ=255)                                                 gtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctac                                                                                              >  1:732335‑M2/93‑149 (MQ=255)                                                  tcactacacgcccatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTAAATCTGTTCAAACGTTTTTcgctccg                                                                                             >  2:271989‑M2/92‑149 (MQ=255)                                                      ttcacgcacatacaacggaggcgggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattactgaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTag                                                                                         <  2:207723‑M2/62‑1 (MQ=255)                                                      tacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCACTCTGTTGTAATGTTTTTCGCTACGCTag                                                                                         >  1:561036‑M2/88‑149 (MQ=255)                                                             acatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAgcaagcaa                                                                                  <  2:732335‑M2/69‑1 (MQ=255)                                                                tacagcggaggggggctgctatttacaataatgcgtgctgcagatcacataaaacactttattacttaacataaatgTATAAGTCATATTTTTGTTTTGAGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaaa                                                                               <  2:118075‑M2/72‑1 (MQ=255)                                                                    acggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaa                                                                           >  2:721164‑M2/74‑149 (MQ=255)                                                                                        tcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATggg                                                       >  1:268570‑M2/54‑149 (MQ=255)                                                                                            taatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGt                                                   >  1:633877‑M2/50‑149 (MQ=255)                                                                                              atgcgtgaagcagatcacaaaaaacgctcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTg                                                 <  1:422120‑M2/102‑1 (MQ=255)                                                                                                   tgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTCTTTCAATCTGCTGAAAAGTTGTTCGCTACGCTAGCAAGCAAAAATAAAACAGTAATAATCGAAATGGGAtgttgccgac                                            >  2:277675‑M2/43‑145 (MQ=255)                                                                                                    gatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGcaca                                           <  1:267839‑M2/108‑1 (MQ=255)                                                                                                      tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGt                                         >  2:108923‑M2/40‑149 (MQ=255)                                                                                                           atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaaa                                    >  1:46961‑M2/35‑149 (MQ=255)                                                                                                           atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaaa                                    >  2:686564‑M2/35‑149 (MQ=255)                                                                                                               cacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTTAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAAc                                >  2:741470‑M2/31‑149 (MQ=255)                                                                                                               cacaaaacactcaatgacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAAc                                <  2:414041‑M2/119‑1 (MQ=255)                                                                                                                       actcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATTTTGCGCACAGTCAAAATAactccccttc                        >  2:308938‑M2/23‑146 (MQ=255)                                                                                                                           aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAGTCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATGGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaata                    >  1:156913‑M2/19‑149 (MQ=255)                                                                                                                           aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAACATAACTCACCGTAaata                    >  1:425628‑M2/19‑149 (MQ=255)                                                                                                                                 ttaacATAAATGAAAAAGTCATACTGTTGTTTAAACAATCATACAAGCTGTTAAAAGGTTTTACGCTACGCTAGCAAGCAAAAAAGAAACAGGAATAATCGAAATGGGATGTTGCGCACAATCAAAATAACCCACCGTAAATAatcttc              <  2:629849/145‑1 (MQ=255)                                                                                                                                   aacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTAACCGTAAATAATCATCTg            >  2:237358‑M2/11‑149 (MQ=255)                                                                                                                                      ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATg                                                         >  2:388923‑M2/8‑101 (MQ=255)                                                                                                                                      ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATg                                                         <  1:388923‑M2/94‑1 (MQ=255)                                                                                                                                           tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAACAAGGAATAATCAAAATGGGATGTTGCGCACAGTCAAAATGACTCACCGTAAATAATCAtctgctactaat    >  2:519595‑M2/3‑144 (MQ=255)                                                                                                                                             tataAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTataaataa  <  1:108923/149‑1 (MQ=255)                                                                                                                                             tataAGTCATACTTTTGTTGTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTataaataa  <  1:746418/149‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                     CCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAA  >  NC_000913_3_bme_pgi/1979139‑1979427

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