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

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1978503–1979270	1979278	9–776	56 [2]	[0] 54	insB‑5–insA‑5	insB‑5,insA‑5

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1978502	2 (0.040)	52 (1.210)	38/262	0.1	98.2%	intergenic (‑305/+16)	flhD/insB‑5	DNA‑binding transcriptional dual regulator FlhD/IS1 protein InsB
	?	NC_000913	1979279 =	0 (0.000)					intergenic (‑64/‑474)	insA‑5/uspC	IS1 protein InsA/universal stress protein C

TCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATG  >  NC_000913/1978364‑1978640                                                                                                                                           |                                                                                                                                           tCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGt                                                                                                                                           <  1:703697‑M1/140‑2 (MQ=255)           tATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcata                                                                                                                                 >  1:464576‑M1/1‑129 (MQ=255)             tCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGta                                                                                                                                          <  1:255921‑M1/129‑3 (MQ=255)             tCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGta                                                                                                                                          >  2:255921‑M1/1‑127 (MQ=255)               tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttt                                                                                                                             >  1:323442‑M1/1‑125 (MQ=255)                  gATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTAAATAAGGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgt                                                                                                                          >  2:673834‑M1/1‑122 (MQ=255)                      tttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttt                                                                                                                       <  2:366713‑M1/139‑22 (MQ=255)                      tttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttt                                                                                                                       >  1:366713‑M1/1‑118 (MQ=255)                         ccTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                                                                   <  1:566574‑M1/140‑26 (MQ=255)                         ccTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                                                                   <  1:524382‑M1/140‑26 (MQ=255)                         ccTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                                                                   <  1:481991‑M1/140‑26 (MQ=255)                         ccTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                                                                   <  2:243847‑M1/140‑26 (MQ=255)                                   cGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                         <  1:725095‑M1/140‑36 (MQ=255)                                           tACTCGTCACTACACGCACATACAACGGATGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaa                                                                                                 <  1:80828‑M1/140‑44 (MQ=255)                                           tACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaa                                                                                                 <  1:764889‑M1/140‑44 (MQ=255)                                                gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtt                                                                                            <  2:248091‑M1/140‑49 (MQ=255)                                                  cACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttt                                                                                          >  1:277427‑M1/1‑90 (MQ=255)                                                     tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATCACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcg                                                                                       <  2:190853‑M1/140‑54 (MQ=255)                                                       cacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatttgttaaaaagtttttcgct                                                                                     >  1:375772‑M1/1‑85 (MQ=255)                                                        acGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                                                                    >  2:47085‑M1/1‑84 (MQ=255)                                                               tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagc                                                                             >  1:121141‑M1/1‑77 (MQ=255)                                                               tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagc                                                                             >  2:463141‑M1/1‑77 (MQ=255)                                                               tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagc                                                                             >  2:31728‑M1/1‑77 (MQ=255)                                                                    cGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagca                                                                        >  2:93029‑M1/1‑72 (MQ=255)                                                                        ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                    >  2:14736‑M1/1‑68 (MQ=255)                                                                         gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                                                                   >  1:676046‑M1/1‑67 (MQ=255)                                                                         gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                                                                   >  2:272413‑M1/1‑67 (MQ=255)                                                                          ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatg                                                                  >  2:659277‑M1/1‑66 (MQ=255)                                                                            ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaa                                                                >  2:505077‑M1/1‑64 (MQ=255)                                                                            ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaa                                                                <  2:755908‑M1/140‑77 (MQ=255)                                                                                     tttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaata                                                       >  1:2021‑M1/1‑55 (MQ=255)                                                                                              tGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg                                              >  2:465182‑M1/1‑46 (MQ=255)                                                                                                cGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg                                            <  1:465182‑M1/140‑97 (MQ=255)                                                                                                cGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg                                            <  2:676046‑M1/140‑97 (MQ=255)                                                                                                     tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgtt                                       >  1:415993‑M1/1‑39 (MQ=255)                                                                                                         gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgc                                   <  1:680169‑M1/140‑106 (MQ=255)                                                                                                         gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgc                                   >  1:646303‑M1/1‑35 (MQ=255)                                                                                                         gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgc                                   <  1:416144‑M1/140‑106 (MQ=255)                                                                                                         gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgc                                   <  2:310164‑M1/140‑106 (MQ=255)                                                                                                         gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgc                                   <  2:415993‑M1/140‑106 (MQ=255)                                                                                                          aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgca                                  <  1:80774‑M1/140‑107 (MQ=255)                                                                                                           tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcac                                 <  1:659277‑M1/140‑108 (MQ=255)                                                                                                            cacacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                >  2:451607‑M1/1‑32 (MQ=255)                                                                                                            cacacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                <  1:272413‑M1/140‑109 (MQ=255)                                                                                                                        tCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataact                    >  2:447674‑M1/1‑20 (MQ=255)                                                                                                                        tCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataact                    >  1:684781‑M1/1‑20 (MQ=255)                                                                                                                         cAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactc                   <  1:93029‑M1/140‑122 (MQ=255)                                                                                                                            ttACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatttgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcacc                <  2:375772‑M1/140‑125 (MQ=255)                                                                                                                               cTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgta             >  2:683013‑M1/1‑13 (MQ=255)                                                                                                                                  aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaat          >  1:570883‑M1/1‑10 (MQ=255)                                                                                                                                  aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactc                   <  2:406273‑M1/131‑122 (MQ=255)                                                                                                                                  aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactc                   >  1:406273‑M1/1‑10 (MQ=255)                                                                                                                                      tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatc      >  1:12892‑M1/1‑6 (MQ=255)                                                                                                                                       aaaTGtataagtcatacttttgttttgggtgtattcaatctgttaaaaagtttgtcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcat     <  1:215375‑M1/140‑136 (MQ=255)                                                                                                                                        agcgggTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTa    >  1:754463/4‑140 (MQ=11)                                                                                                                                          cgggTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATg  >  2:691705/2‑140 (MQ=11)                                                                                                                                           |                                                                                                                                           TCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATG  >  NC_000913/1978364‑1978640

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑