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

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1978503–1979269	1979278	10–776	36 [0]	[0] 36	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)	32 (0.660)	27/260	NT	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

ATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCT  >  NC_000913/1979135‑1979419                                                                                                                                                 |                                                                                                                                             atgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgtgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATaa                                                                                                                                          <  2:757281‑M2/5‑1 (MQ=255)               ttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgctcatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacctaaatgTATAAGTCATACttttgtt                                                                                                                            >  2:103996‑M2/131‑149 (MQ=255)                 atctaagatttttcctaaatcgacgcaactgtagtcgtcactacacgcacatacaacggaagggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattcattaacatgaatgTATAAGTCATACtttattttt                                                                                                                          >  2:314143‑M2/129‑149 (MQ=255)                   ctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgcgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                        <  2:798132‑M2/23‑1 (MQ=255)                     aagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                                                                      >  1:215190‑M2/125‑149 (MQ=255)                      agatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtg                                                                                                                     >  2:89435‑M2/124‑149 (MQ=255)                             tcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggtgggcggcgtcgagtttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCa                                                                                                              <  1:623873‑M2/33‑1 (MQ=255)                                taaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTtcaagc                                                                                                           >  2:188544‑M2/114‑147 (MQ=255)                                 aaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatccagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTTGGTGTATTTTAATCt                                                                                                          >  2:476390‑M2/113‑149 (MQ=255)                                  aatcgacgcaaccgtacccgtcaccacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTg                                                                                                         <  1:184636‑M2/38‑1 (MQ=255)                                               gtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt                                                                                            >  1:289921‑M2/99‑149 (MQ=255)                                                   tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTtggca                                                                                        >  2:637062‑M2/95‑148 (MQ=255)                                                         ctacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                                                                  >  2:523905‑M2/89‑149 (MQ=255)                                                           acacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagc                                                                                >  1:534145‑M2/87‑149 (MQ=255)                                                                         cggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                                                                  <  1:530921‑M2/77‑1 (MQ=255)                                                                               ggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaa                                                            >  2:315735‑M2/67‑149 (MQ=255)                                                                               ggggctgcgattttcaataatacgtgatgcagagcacacaaaacactcaattacttaacataaatgTATAAGTGATACTTTAGTTTTGGGTGTATTTCAATCTGTTGAAAAGTTTTATGCTACGCTAGCAAGAaaaaaggaaaaaggat                                                            >  2:512013‑M2/67‑148 (MQ=255)                                                                                 ggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaata                                                          <  2:115512‑M2/85‑1 (MQ=255)                                                                                 ggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaata                                                          >  2:118187‑M2/65‑149 (MQ=255)                                                                                      cgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGa                                                     >  1:620482‑M2/60‑149 (MQ=255)                                                                                        attttcaataatgcctgatgcagatcacacaaaacactcaattacttaagataaatgTATAAGTCATACTTTTGTATTGGGTGTATTTCAATCCGTTAAAACGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa                                                   >  1:673297‑M2/58‑149 (MQ=255)                                                                                        attttcaataatgcatgatgcagatcacacaaaacactcaattacttaaaataaatgTATAAGTCATACTTTTGTGTTGGGTGTATTTCAATCTGTTAAAACGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa                                                   >  1:673264‑M2/58‑149 (MQ=255)                                                                                           ttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTCTAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCACTACGCGAGCAATCAAAAATGGAACAGGAATGATcgatatag                                                >  2:129697‑M2/55‑147 (MQ=255)                                                                                              aataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATGGGACAGACTGCTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt                                             <  1:664260‑M2/98‑1 (MQ=255)                                                                                                      gtgatgcagatcacacaaaacacttaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGca                                     <  2:253920‑M2/106‑1 (MQ=255)                                                                                                      gtgatgcagatcacaaaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGca                                     <  1:597272‑M2/106‑1 (MQ=255)                                                                                                       tgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTCTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATTTTGCGcac                                    >  2:592260‑M2/43‑149 (MQ=255)                                                                                                          tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGt                                 >  1:767552‑M2/40‑149 (MQ=255)                                                                                                          tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGt                                 >  1:720264‑M2/40‑149 (MQ=255)                                                                                                            cagatcacacaaaacactccattccttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCa                               <  1:197187‑M2/112‑1 (MQ=255)                                                                                                               atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaaa                            <  2:363352‑M2/115‑1 (MQ=255)                                                                                                               atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaaa                            >  2:488366‑M2/35‑149 (MQ=255)                                                                                                               atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaaa                            >  1:753176‑M2/35‑149 (MQ=255)                                                                                                                                    cttaacataaatgTATAAGTCATACTTTTCTTTTGGGTGTATTTCAATCTGTTTAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAATCAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAATTCACCGTAAATAatcat       <  2:673264‑M2/136‑1 (MQ=255)                                                                                                                                    cttaacataaatgTATAAGTCATACTTTTCTTTTGGGTGTATTTCAATCTGTTCAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAGCAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatcat       <  2:673297‑M2/136‑1 (MQ=255)                                                                                                                                         aataaatgTATAGGTCCTACTTTTGTTTTGGGTGTATTTCAATGTGTTAAAAAGTTTGTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCt  <  2:656817‑M2/141‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                             ATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCT  >  NC_000913/1979135‑1979419

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