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–1979269	1979278	10–776	47 [0]	[0] 47	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)	45 (1.000)	39/266	0.1	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

GTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTA  >  NC_000913/1978363‑1978638                                                                                                                                            |                                                                                                                                         gTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                                                                       >  1:62308‑M1/1‑140 (MQ=255)  tCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtata                                                                                                                                      >  1:472499‑M1/1‑139 (MQ=255)   cTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                     <  2:611986‑M1/143‑6 (MQ=255)    tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATCTTCAATAATGCGTGATGCAGATCACACAAACCACTCAATTACTTAACATCAAAGtataag                                                                                                                                    >  2:528064‑M1/1‑137 (MQ=255)     ttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagt                                                                                                                                   >  1:515552‑M1/1‑136 (MQ=255)     ttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagt                                                                                                                                   >  2:217351‑M1/1‑136 (MQ=255)           ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatact                                                                                                                             >  2:777293‑M1/1‑130 (MQ=255)                tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgt                                                                                                                        <  1:197885‑M1/143‑19 (MQ=255)                tAAGATTTTTCCTAAAGCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgt                                                                                                                        <  1:152497‑M1/143‑19 (MQ=255)                 aaGATTTTTCCTAAATCGACGCAACTGTACTAGTCACGACACTCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtt                                                                                                                       <  1:596143‑M1/143‑20 (MQ=255)                  aGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGAACGCTGCAAATCACACAACACAATAAATTACTTAACATAAATAtataagccctactatagtat                                                                                                                      >  2:13329‑M1/1‑123 (MQ=255)                  aGATGTTTCCTAAATCGACCCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttt                                                                                                                      <  1:760685‑M1/143‑21 (MQ=255)                   gttttttccTATATCGACTCAAGTTTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGGGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttt                                                                                                                     <  1:671257‑M1/141‑22 (MQ=255)                   gATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataggtcatacttttgtttt                                                                                                                     >  2:288370‑M1/1‑122 (MQ=255)                             aaaTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACTCTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattt                                                                                                           <  2:21348‑M1/143‑32 (MQ=255)                                 cGACGCAACCGTACCCGTCACTGCACGCACATACAACGGAGGGGGGCGGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                       <  2:405050‑M1/143‑36 (MQ=255)                                       aaCTGTAACCGTCACTACACGCACATACAACGGAGGCGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtta                                                                                                 <  1:768726‑M1/143‑42 (MQ=255)                                          tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaa                                                                                              >  2:661390‑M1/1‑99 (MQ=255)                                                           gctcatacAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCCGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta                                                                             <  1:566452‑M1/140‑62 (MQ=255)                                                           gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtattaatcataatgttgtcttgggtctattttaacatgttacaaagtttttctgtacgcta                                                                             >  2:421767‑M1/1‑82 (MQ=255)                                                           gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctt                                                                             >  1:57201‑M1/1‑82 (MQ=255)                                                               aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaa                                                                         <  2:877128‑M1/143‑66 (MQ=255)                                                                tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACGTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttgttcgctacgctagcagg                                                                        >  2:571290‑M1/1‑77 (MQ=255)                                                                      ggAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                  <  2:596358‑M1/143‑73 (MQ=255)                                                                         ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatatcaatgtgttaaaaagtttttcgctacgctagcaagcaaaaatga                                                               >  2:238547‑M1/1‑68 (MQ=255)                                                                               cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaacaatgaaacagg                                                         >  2:372867‑M1/1‑62 (MQ=255)                                                                                   gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaata                                                     >  1:466310‑M1/1‑58 (MQ=255)                                                                                    cttttcAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaacaagtttttcgctacgctagcacgcaaaaatgaaacaggaacaa                                                    >  1:534580‑M1/2‑57 (MQ=255)                                                                                          aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaacatgattcaggcatcatcgaaa                                              >  1:424237‑M1/1‑51 (MQ=255)                                                                                                gCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat                                        <  2:702902‑M1/143‑99 (MQ=255)                                                                                                 cGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatg                                       >  1:77780‑M1/1‑44 (MQ=255)                                                                                                   tGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgtt                                     >  1:590284‑M1/1‑42 (MQ=255)                                                                                                    gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattgcaagctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttg                                    <  1:21680‑M1/143‑103 (MQ=255)                                                                                                       gCAGATCACACAAAACACTCAATTACTTAACATAAATGtatcagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgtggctc                                 >  1:824547‑M1/1‑38 (MQ=255)                                                                                                       gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgc                                 <  2:130943‑M1/143‑106 (MQ=255)                                                                                                       gCAGAGCACACAAAACACGCAATTACTTAACCTAAGTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgc                                 <  1:750766‑M1/143‑106 (MQ=255)                                                                                                              acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtca                          <  1:858852‑M1/143‑113 (MQ=255)                                                                                                                  aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaat                      >  1:697864‑M1/1‑27 (MQ=255)                                                                                                                     acagttaattacttGTCATAAATGtataagtggtactggagatttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataac                   <  2:180965‑M1/137‑120 (MQ=255)                                                                                                                     acacTCAATTACTTAACATAAATGtataaggcatacttttgttttgtgtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataac                   >  2:313426‑M1/1‑24 (MQ=255)                                                                                                                        cTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactca                <  1:218488‑M1/143‑123 (MQ=255)                                                                                                                           aaTTACTTAACATAAATGtataagtcatacttttgttttgtgtgtatttcaatctgtcaaaaagtttttcgatccgctagcaagcaaaaatgaaacaggaataaccgaaatgggatgttgcgcacagtccaaatacctaaccg             >  2:43807‑M1/1‑18 (MQ=255)                                                                                                                           aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccg             >  1:563424‑M1/1‑18 (MQ=255)                                                                                                                                   aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataat     >  1:190280‑M1/1‑10 (MQ=255)                                                                                                                                    aCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatc    <  1:288370‑M1/143‑135 (MQ=255)                                                                                                                                      aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcat  >  1:772670‑M1/1‑7 (MQ=255)                                                                                                                                      aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcat  <  1:393314‑M1/143‑137 (MQ=255)                                                                                                                                            |                                                                                                                                         GTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTA  >  NC_000913/1978363‑1978638

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑