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

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1978503	1979278	776	40 [0]	[1] 41	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)	35 (0.960)	29/268	0.2	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

TGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGC  >  NC_000913/1979146‑1979418                                                                                                                                      |                                                                                                                                            catcccatttcgatttttcctaaatcgacgcaactgtactcgtccctacactcacatacaacggaggggggctgcgattttcaataatcggtgatgcagaccacacaaaacactcaattacttaacctaaataTATAAGTCAt                                                                                                                                    >  1:108691‑M2/134‑143 (MQ=255)      gtgagttatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaagactcaattacttaacataaatgTTTCAGTCATACttt                                                                                                                               >  1:356500‑M2/129‑143 (MQ=255)              tttttcctaaatcgacgcatgtgtacccgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacataaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                       <  1:255830‑M2/23‑1 (MQ=255)               ttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                                                                                                      <  2:448809‑M2/24‑1 (MQ=255)                 ttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtg                                                                                                                    >  2:474204‑M2/118‑143 (MQ=255)                    ctaaatcgacgcaactgtactcgtcactacacgcacataccactgaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAt                                                                                                                 <  2:248188‑M2/29‑1 (MQ=255)                        atcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCa                                                                                                             <  1:65015‑M2/33‑1 (MQ=255)                         gcgacgcagctgtactcgtctctacactcacatacaacgggggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacacttaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCaa                                                                                                            <  1:105886‑M2/34‑1 (MQ=255)                            acgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCt                                                                                                         >  1:573723‑M2/107‑143 (MQ=255)                               gagctctactcctcactccactcacctacaacggtggggggctgcgattttcactaatgcgtgatgcagatcccacaaaaccctcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGtt                                                                                                      <  1:7924‑M2/40‑1 (MQ=255)                                        tcgtctctaaatgcacatacaacgtaggagtgctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttt                                                                                             <  1:422674‑M2/49‑1 (MQ=255)                                         cgtcactacacgcacatacaacggaggggggatgacatttcgaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtttt                                                                                            <  1:473701‑M2/50‑1 (MQ=255)                                              ctacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgct                                                                                       >  1:527336‑M2/89‑143 (MQ=255)                                                    gcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTCGGGTGTCTTATAGTCTGTTAAAAAGTTTTTcgctacgcta                                                                                 >  2:634443‑M2/83‑143 (MQ=255)                                                    gcacagacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                                                                 <  1:254364‑M2/61‑1 (MQ=255)                                                              cggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaa                                                                       >  2:247341‑M2/73‑143 (MQ=255)                                                                 aggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATg                                                                    >  2:90339‑M2/70‑143 (MQ=255)                                                                 aggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAGTTTAGTCTGCTATAGAGCGTGTCGCTACGCGAGCAGGAAACAATg                                                                    >  2:279155‑M2/70‑143 (MQ=255)                                                                       gctgctatttttaataatgcattatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATATGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAg                                                              <  1:208809‑M2/80‑1 (MQ=255)                                                                         tgcgattttctataatgcgtgatgcagatcacacaaaacactcaattacttaacataactcTATAAGTCATACTTTTGTTCTGGGTGTCTTTCAATCTGTTAAAAAGTTTTTCTCTACTTTAGCAAGCAAAAACTAAAcagta                                                            >  1:255240‑M2/62‑141 (MQ=255)                                                                           cgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaat                                                          >  1:262998‑M2/60‑143 (MQ=255)                                                                           cgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaat                                                          >  1:246710‑M2/60‑143 (MQ=255)                                                                                ttaaatactccgtgatgcagatcacacaaaacacccaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg                                                     <  2:71994‑M2/89‑1 (MQ=255)                                                                                      aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATgg                                               <  2:573723‑M2/95‑1 (MQ=255)                                                                                          cgttatgcacatcacacaaaccacgcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGGATTTCAATCTGTTAAAAAGTGTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg                                           <  2:108691‑M2/99‑1 (MQ=255)                                                                                              tttgagattacataaaacactcaattacttaatattaaagTATAAGACATACGTTGGGTTTGGGGGGATTTCAATCTGTTAAAAAGTTTTTCGCTACGCGAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTgc                                       <  2:391734‑M2/103‑1 (MQ=255)                                                                                                  agatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTTGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCACAAATGAAACAGTAATAACCGAACTGGGATAttgcgccc                                   >  1:224536‑M2/37‑141 (MQ=255)                                                                                                       acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCa                              >  1:176914‑M2/32‑143 (MQ=255)                                                                                                             accactcagttacttaacctcaacgGATAAGTCATACTTTTGTTTTGGGTGGATATCAATCCGTGAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAAGCGAAATGGGATGTTGCGCACAGTCAAAATaa                        <  2:382777‑M2/118‑1 (MQ=255)                                                                                                             aacactcaattactgaacataaatgTATAAGTCATACATTTGTGTTGGGTGTATTTCACTCTGTTAAAAAGTTCTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATaa                        <  2:458248‑M2/118‑1 (MQ=255)                                                                                                                  tccgttacgtaacatcaaggTATAAGTCATACTATTGTTTTGGGTGTATTTCAATCTGTCAAAAAGGTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCAc                   <  2:644638‑M2/123‑1 (MQ=255)                                                                                                                  tcaaagacttaacataaaagTATAAGTCATACTTTAGTTTTGGGTGAAATTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCAc                   <  2:356500‑M2/123‑1 (MQ=255)                                                                                                                       ttgtgcacATTAATGTATGAGTTATGTTTTTTTTTTGGGGGTATTTCAATCTGTGAAAAAGGTTTTCGCTATGCCAGCAATCAAAAATGAAACAGGAATACTCGAAATGGGATGTTGCGCACAGTCAAAATACCTCACCGTaa              <  1:279155/136‑1 (MQ=255)                                                                                                                       tacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTaa              <  1:90339‑M2/128‑1 (MQ=255)                                                                                                                        acttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTaaa             >  2:94897‑M2/15‑143 (MQ=255)                                                                                                                          ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGCTGTTGCGAACAGTCAAAATAACTCACCGTAaata           >  2:169406‑M2/13‑143 (MQ=255)                                                                                                                               ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAGATAatcat      >  2:452648‑M2/8‑143 (MQ=255)                                                                                                                               ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACAgtaaataattaa      >  1:72529‑M2/8‑140 (MQ=255)                                                                                                                                   atgTATAAGTCATACTTTTGTTTTGGGTGTATTTCGATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGt                                          <  2:513938‑M2/100‑1 (MQ=255)                                                                                                                                   atgTATAAGTCATACTTTTGTTTTGGGTGTATTTCGATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGt                                          >  1:513938‑M2/4‑103 (MQ=255)                                                                                                                                   atgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTACATAATCATCTGc  >  1:199953‑M2/4‑143 (MQ=255)                                                                                                                                      |                                                                                                                                            TGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGC  >  NC_000913/1979146‑1979418

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