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	1979278	776	32 [0]	[0] 32	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)	28 (0.680)	20/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

ATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAA  >  NC_000913/1978359‑1978648                                                                                                                                                |                                                                                                                                                   aTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                               <  2:454808‑M1/149‑6 (MQ=255) aTGGGCCTTTACACTTATCTAACATTTTTCCTAAATTGACGCAACTTTACACGTCACTACACTCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCATTTACTTAACATAAAGGtataa                                                                                                                                               <  2:424160‑M1/149‑6 (MQ=255)      tCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCAAACACAACACTCAATTACTTAACATAAATGtataagtcat                                                                                                                                          >  1:299680‑M1/1‑139 (MQ=255)        tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatac                                                                                                                                        <  2:502721‑M1/149‑13 (MQ=255)                     aaGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                                                                           >  1:708243‑M1/1‑124 (MQ=255)                               cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                                 >  2:763531‑M1/1‑114 (MQ=255)                               cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                                 >  2:178606‑M1/1‑114 (MQ=255)                                tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                                                                >  1:232955‑M1/1‑113 (MQ=255)                                       aCGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggggtatttcaatctgttaaa                                                                                                         >  1:118767‑M1/1‑106 (MQ=255)                                              tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttt                                                                                                  >  1:722156‑M1/1‑99 (MQ=255)                                                     tgcacTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCTATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctac                                                                                           <  1:108331‑M1/147‑58 (MQ=255)                                                       cACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgc                                                                                         <  1:446072‑M1/149‑60 (MQ=255)                                                         cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtgaaaaaatttttctctacgcta                                                                                       >  2:638507‑M1/1‑88 (MQ=255)                                                          tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctag                                                                                      >  2:97343‑M1/1‑87 (MQ=255)                                                          tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctaagctag                                                                                      >  2:275364‑M1/1‑87 (MQ=255)                                                             acGCACATACAACGGAGTGGTGCTGCGATTTTCAATATTGCTTGATGCAGATCACACAAAACACTCAATTACTTAACATAATTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaactttttcgctacgctagcta                                                                                   <  2:404504‑M1/149‑66 (MQ=255)                                                                            gggggggcTGCTATTTTCAATAATGCGTGATGCAGATCACACCACACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtgtttcaatctgttaaaaaggttttcgctccgctagcaagcaaaaatgaaacag                                                                    <  2:718014‑M1/148‑81 (MQ=255)                                                                                           ttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcgagcaaaaatgaaacaggaataatcgaaatgg                                                     >  2:624356‑M1/1‑54 (MQ=255)                                                                                           ttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg                                                     <  1:5840‑M1/149‑96 (MQ=255)                                                                                                         aTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacag                                       <  2:708243‑M1/149‑110 (MQ=255)                                                                                                               aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcggacagtcaaaa                                 >  1:629482‑M1/1‑34 (MQ=255)                                                                                                               aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                                 >  1:466440‑M1/1‑34 (MQ=255)                                                                                                                      aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctccgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataaccca                          >  1:367608‑M1/1‑27 (MQ=255)                                                                                                                      aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtggatttcaatctgttaaaaagtttttcttttcttttgcaagcaaaaatgaaacaggaataatcgaaatgggatgttgggcccagccaaaataactcc                          >  2:419178‑M1/1‑27 (MQ=255)                                                                                                                       aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttttcttttttttttagcaaaaaagaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcac                         >  2:458849‑M1/1‑26 (MQ=255)                                                                                                                       aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcac                         >  1:325498‑M1/1‑26 (MQ=255)                                                                                                                               aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaaccaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaata                 >  2:76687‑M1/1‑18 (MQ=255)                                                                                                                               aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaaccaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaata                 >  2:73523‑M1/1‑18 (MQ=255)                                                                                                                                  tACTTAACATAAATGtataattcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaacatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatc              <  2:333636‑M1/149‑135 (MQ=255)                                                                                                                                  tACTTAACATAAACGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatc              <  2:463584‑M1/149‑135 (MQ=255)                                                                                                                                            aaaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctata    >  2:336240‑M1/1‑5 (MQ=255)                                                                                                                                              aTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaa  >  1:57332‑M1/1‑3 (MQ=255)                                                                                                                                                |                                                                                                                                                   ATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAA  >  NC_000913/1978359‑1978648

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