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–1979269	1979278	10–776	38 [0]	[0] 38	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)	37 (0.940)	33/274	0.2	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

ATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAA  >  NC_000913_3_bme_pgi/1979135‑1979424                                                                                                                                                 |                                                                                                                                                  atgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacgcaaaacactcaattacttaacataaatgTATaa                                                                                                                                               >  2:383266‑M2/145‑149 (MQ=255)  tgggtctttacacttatctaagatttttcctaaatcgacccaaataaactcgtcactacacgcacatacaacgaagcggcgctacgattttcaataatacgtgattcagatcacacaaaacactcaattcctttccataaatgTAtaaa                                                                                                                                              >  2:42322‑M2/144‑148 (MQ=255)       ctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATa                                                                                                                                         >  1:588186‑M2/139‑149 (MQ=255)             acttatctaggatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaaccgaggggggctgcgattttcaataatgcgtcatgcagagaacacaaagaagtccattacttaacataaatgTATAACACGTCCttttg                                                                                                                                   >  2:2799‑M2/133‑149 (MQ=255)                 atcaaagagtgttccgaaatcgacgcaactgtactcatcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaactctcaatttcttaacagaaatgTATAAGTCATACttttgtttt                                                                                                                               <  2:123775‑M2/21‑1 (MQ=255)                   ctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                             <  1:58037‑M2/23‑1 (MQ=255)                                    tcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTTGGTgtatttaaacttttt                                                                                                            >  2:384305‑M2/110‑143 (MQ=255)                                             ctgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttt                                                                                                   >  2:290491‑M2/101‑149 (MQ=255)                                             ctgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttt                                                                                                   >  2:659145‑M2/101‑149 (MQ=255)                                              tgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtttt                                                                                                  >  1:591693‑M2/100‑149 (MQ=255)                                                    cgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgcta                                                                                            <  1:226742‑M2/56‑1 (MQ=255)                                                     gtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCATTCTGTTAAAAAGTTTTTcgctac                                                                                           >  1:400868‑M2/93‑149 (MQ=255)                                                           acacgcagatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagc                                                                                     <  2:668676‑M2/63‑1 (MQ=255)                                                              cgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACATTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaag                                                                                  <  1:481427‑M2/66‑1 (MQ=255)                                                                        acggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCCAAAATGaa                                                                        >  1:649588‑M2/74‑149 (MQ=255)                                                                        acggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaa                                                                        >  2:392474‑M2/74‑149 (MQ=255)                                                                                   ctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat                                                             >  1:519986‑M2/63‑149 (MQ=255)                                                                                   ctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat                                                             <  2:69505‑M2/87‑1 (MQ=255)                                                                                   ctgcgattttcaacaatgcgtgatgcagaacacacaaaacactcaattacgtaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATGTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGgaataag                                                             <  1:679398‑M2/87‑2 (MQ=255)                                                                                   cgtcgatttttaataatgcgggaggcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTTTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCGAGCAAGCAAAAATGAAACAGGaataat                                                             <  2:47989‑M2/87‑1 (MQ=255)                                                                                    tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTTTTTTTGTTTTATTTCAATCTGTTAAAAAGTTTTTCGTTACGCTAGTAAGCAAACGTGAAACAGgagtaata                                                            >  2:93740‑M2/62‑148 (MQ=255)                                                                                       gattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAttt                                                                                                                     >  2:279785‑M2/59‑89 (MQ=255)                                                                                       gattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAttt                                                                                                                     <  1:279785‑M2/31‑1 (MQ=255)                                                                                        attttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCATACAAAAATGAAACAGGAAtaatctcaa                                                        >  2:85510‑M2/58‑145 (MQ=255)                                                                                         ttttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTGTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTGGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt                                                       >  2:384979‑M2/57‑149 (MQ=255)                                                                                            tcaataagtcgtgatgaagatcatacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATggg                                                    <  1:383266‑M2/96‑1 (MQ=255)                                                                                               ataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg                                                 >  2:241925‑M2/51‑149 (MQ=255)                                                                                               ataatgactgatccagaacaaacaaaaccctcaattacttaacacaaaggTATAAGTCCTACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTCCGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg                                                 <  1:42322‑M2/99‑1 (MQ=255)                                                                                                    gcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTgcg                                            <  1:275472‑M2/104‑1 (MQ=255)                                                                                                              gatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaa                                  <  1:91890‑M2/114‑1 (MQ=255)                                                                                                               atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaaa                                 >  1:712875‑M2/35‑149 (MQ=255)                                                                                                               atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaaa                                 <  2:712875‑M2/115‑1 (MQ=255)                                                                                                                  acaaaaaccacgcaattccttaccataaatgTATAAGTCATACTTTTGTTTGGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATaa                              <  2:505603‑M2/118‑1 (MQ=255)                                                                                                                    acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACt                            >  1:703599‑M2/30‑149 (MQ=255)                                                                                                                     caaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTc                           >  1:471023‑M2/29‑149 (MQ=255)                                                                                                                             tcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTaaa                   >  1:307585‑M2/21‑149 (MQ=255)                                                                                                                                  tacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatc              <  2:588186‑M2/134‑1 (MQ=255)                                                                                                                                              atgTATAAGGTAGACTTTTGTTTTGCGTGTAGTTCAATCGGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATACTCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTataaa  <  2:620897‑M2/146‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                  ATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAA  >  NC_000913_3_bme_pgi/1979135‑1979424

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