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	39 [0]	[1] 40	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)	34 (0.970)	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

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