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

ATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCG  >  NC_000913/1978359‑1978643                                                                                                                                                |                                                                                                                                              aTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGTGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                          <  2:757281‑M1/149‑6 (MQ=255)               ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCTCATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACCTAAATGtataagtcatacttttgtt                                                                                                                            >  2:103996‑M1/1‑130 (MQ=255)                 aTCTAAGATTTTTCCTAAATCGACGCAACTGTAGTCGTCACTACACGCACATACAACGGAAGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTCATTAACATGAATGtataagtcatactttattttt                                                                                                                          >  2:314143‑M1/1‑128 (MQ=255)                   cTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGCGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                        <  2:798132‑M1/149‑24 (MQ=255)                     aaGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                                                                      >  1:215190‑M1/1‑124 (MQ=255)                      aGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtg                                                                                                                     >  2:89435‑M1/1‑123 (MQ=255)                             tCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGTGGGCGGCGTCGAGTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttca                                                                                                              <  1:623873‑M1/149‑34 (MQ=255)                                tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaagc                                                                                                           >  2:188544‑M1/1‑113 (MQ=255)                                 aaaTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATCCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttttggtgtattttaatct                                                                                                          >  2:476390‑M1/1‑112 (MQ=255)                                  aaTCGACGCAACCGTACCCGTCACCACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctg                                                                                                         <  1:184636‑M1/149‑39 (MQ=255)                                               gTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                            >  1:289921‑M1/1‑98 (MQ=255)                                                   tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttggca                                                                                        >  2:637062‑M1/1‑94 (MQ=255)                                                         cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta                                                                                  >  2:523905‑M1/1‑88 (MQ=255)                                                           acacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagc                                                                                >  1:534145‑M1/1‑86 (MQ=255)                                                                         cGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                                                                  <  1:530921‑M1/149‑78 (MQ=255)                                                                               ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaa                                                            >  2:315735‑M1/1‑66 (MQ=255)                                                                               ggggCTGCGATTTTCAATAATACGTGATGCAGAGCACACAAAACACTCAATTACTTAACATAAATGtataagtgatactttagttttgggtgtatttcaatctgttgaaaagttttatgctacgctagcaagaaaaaaggaaaaaggat                                                            >  2:512013‑M1/1‑66 (MQ=255)                                                                                 ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaata                                                          <  2:115512‑M1/149‑86 (MQ=255)                                                                                 ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaata                                                          >  2:118187‑M1/1‑64 (MQ=255)                                                                                      cGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcga                                                     >  1:620482‑M1/1‑59 (MQ=255)                                                                                        aTTTTCAATAATGCCTGATGCAGATCACACAAAACACTCAATTACTTAAGATAAATGtataagtcatacttttgtattgggtgtatttcaatccgttaaaacgtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa                                                   >  1:673297‑M1/1‑57 (MQ=255)                                                                                        aTTTTCAATAATGCATGATGCAGATCACACAAAACACTCAATTACTTAAAATAAATGtataagtcatacttttgtgttgggtgtatttcaatctgttaaaacgtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa                                                   >  1:673264‑M1/1‑57 (MQ=255)                                                                                           ttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtctaagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcactacgcgagcaatcaaaaatggaacaggaatgatcgatatag                                                >  2:129697‑M1/1‑54 (MQ=255)                                                                                              aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtatgggacagactgctgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat                                             <  1:664260‑M1/149‑99 (MQ=255)                                                                                                      gTGATGCAGATCACACAAAACACTTAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgca                                     <  2:253920‑M1/149‑107 (MQ=255)                                                                                                      gTGATGCAGATCACAAAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgca                                     <  1:597272‑M1/149‑107 (MQ=255)                                                                                                       tGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttctgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggattttgcgcac                                    >  2:592260‑M1/1‑42 (MQ=255)                                                                                                          tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagt                                 >  1:767552‑M1/1‑39 (MQ=255)                                                                                                          tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagt                                 >  1:720264‑M1/1‑39 (MQ=255)                                                                                                            cAGATCACACAAAACACTCCATTCCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtca                               <  1:197187‑M1/149‑113 (MQ=255)                                                                                                               aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                            <  2:363352‑M1/149‑116 (MQ=255)                                                                                                               aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                            >  2:488366‑M1/1‑34 (MQ=255)                                                                                                               aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                            >  1:753176‑M1/1‑34 (MQ=255)                                                                                                                                    cTTAACATAAATGtataagtcatacttttcttttgggtgtatttcaatctgtttaaaagtttttcgctacgctagcaagcaaaaatgaatcaggaataatcgaaatgggatgttgcgcacagtcaaaataattcaccgtaaataatcat       <  2:673264‑M1/149‑137 (MQ=255)                                                                                                                                    cTTAACATAAATGtataagtcatacttttcttttgggtgtatttcaatctgttcaaaagtttttcgctacgctagcaagcaaaaatgaagcaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcat       <  2:673297‑M1/149‑137 (MQ=255)                                                                                                                                         aatAAATGtataggtcctacttttgttttgggtgtatttcaatgtgttaaaaagtttgtcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgct  <  2:656817‑M1/148‑142 (MQ=255)                                                                                                                                                |                                                                                                                                              ATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCG  >  NC_000913/1978359‑1978643

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