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	50 [0]	[0] 50	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)	48 (0.970)	43/266	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

AAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCT  >  NC_000913/1978357‑1978645                                                                                                                                                  |                                                                                                                                                aaaTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                                                                              <  1:125211‑M1/149‑4 (MQ=255)   aTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                            <  1:875938‑M1/149‑6 (MQ=255)          tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACCACGGAGGGGGGCTGCGATTTTCAATAATCCGTGAGGCAGATCACACAAAACACTCAATTACTTACCATAAATGtacaagtcatac                                                                                                                                     >  2:800689‑M1/1‑137 (MQ=255)           ttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAAATCACACAAAACACTCAATTACTTACCATAAATGtttaagtcatact                                                                                                                                    >  2:31674‑M1/1‑136 (MQ=255)                     cTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCATTAATGCGTGATTCAGAGCACACAAAACGCTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                          <  1:183526‑M1/149‑24 (MQ=255)                      tAAGAGTTTTCCTAACTCGACGGAACTGTATTCGTCAGTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                                                                         <  1:68615‑M1/149‑25 (MQ=255)                       aaGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                                                                        >  1:395942‑M1/1‑124 (MQ=255)                       aaGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAAGGCGTGATGCAGATCACAAAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                                                                        >  1:846791‑M1/1‑124 (MQ=255)                            ttttCCTAAATCGACGCAACCGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGCGATGCAGATCACACAAAACACTCAATTACGTAACATAAATGtataagtcatacttttgttttgggtgtatt                                                                                                                   <  2:353627‑M1/149‑31 (MQ=255)                             tttCCTAAAGCGACGCAACTGGACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattt                                                                                                                  <  2:178900‑M1/149‑32 (MQ=255)                              ttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                                                                                                 >  1:256487‑M1/1‑117 (MQ=255)                                 cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                              <  1:483989‑M1/149‑36 (MQ=255)                                     aTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatatttttgttttgggtgtatttcaatctct                                                                                                          >  1:756170‑M1/1‑110 (MQ=255)                                     aTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                                                                          >  2:525056‑M1/1‑110 (MQ=255)                                        gACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaa                                                                                                       >  2:183402‑M1/1‑107 (MQ=255)                                           ggaactGTACTCGTTACCACACGTAAATACAACGGAGGGGGGCTCACATTTTCAATAATGCGTGATGCAGATCACACAAATCTCTCAAGTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                    <  1:356077‑M1/147‑46 (MQ=255)                                             aaCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggttgtatttcattctgttaaaaagt                                                                                                  >  2:224378‑M1/1‑102 (MQ=255)                                              aCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtt                                                                                                 <  1:525056‑M1/149‑49 (MQ=255)                                                  tACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatagtttttttttgggtgtatttcaatctgttaaaaagtttttt                                                                                             >  1:248877‑M1/1‑97 (MQ=255)                                                       gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctac                                                                                        >  1:846842‑M1/1‑92 (MQ=255)                                                        tCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcctactttagttttgggtgtatttcactctgttaaacagtttttcgctacg                                                                                       >  1:857601‑M1/1‑91 (MQ=255)                                                        tCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacg                                                                                       >  1:816942‑M1/1‑91 (MQ=255)                                                         cACTACCCGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgc                                                                                      <  2:757530‑M1/149‑60 (MQ=255)                                                           cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta                                                                                    <  2:627711‑M1/149‑62 (MQ=255)                                                              cacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagca                                                                                 <  1:409925‑M1/149‑65 (MQ=255)                                                                 gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAACCCCTGCCGTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttacaaagtttttcgctacgctagcaagc                                                                              >  1:901960‑M1/1‑82 (MQ=255)                                                                 gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagc                                                                              >  1:415718‑M1/1‑82 (MQ=255)                                                                   acaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaa                                                                            >  1:469481‑M1/1‑80 (MQ=255)                                                                       aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                                                                        >  2:136231‑M1/1‑76 (MQ=255)                                                                         cgagggggggggcggcgaGTTTCAATGATGAGTGATGCAGATCACACAAAGCGCTCAATTACTTAACATAAATGtataagtcagacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatga                                                                      <  2:340029‑M1/143‑76 (MQ=255)                                                                           cGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                                                                    <  1:224378‑M1/149‑78 (MQ=255)                                                                              aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                                 <  2:388115‑M1/149‑81 (MQ=255)                                                                                gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaggtttttcgctacgctagcaagcaaaaatgaaacagga                                                               >  2:848460‑M1/1‑67 (MQ=255)                                                                                  gggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaat                                                             >  1:331123‑M1/1‑65 (MQ=255)                                                                                      tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatc                                                         >  1:738790‑M1/1‑61 (MQ=255)                                                                                       gCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaa                                                                     <  2:92380‑M1/136‑77 (MQ=255)                                                                                       gCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaa                                                                     >  1:92380‑M1/1‑60 (MQ=255)                                                                                            tttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaa                                                                                >  2:740642‑M1/1‑55 (MQ=255)                                                                                            tttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaa                                                                                <  1:740642‑M1/120‑66 (MQ=255)                                                                                                aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataaatcatactcttgttttgggtgtatttcaatattttaaaaactttttctctacgctagcaagcaaaaatgaacaaggaataatcgaaataggat                                               >  2:508392‑M1/1‑51 (MQ=255)                                                                                                 ataatGCGTGATGCAGATCACACAAAACACTTAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatg                                              <  2:816942‑M1/149‑100 (MQ=255)                                                                                                 ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgttacgcttgcaagaaaaaataaaaaaggaataaacgacatatgctg                                              >  1:65842‑M1/1‑50 (MQ=255)                                                                                                            tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctaaaaaaaaaaaaaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagt                                   <  2:395942‑M1/149‑111 (MQ=255)                                                                                                                         aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcac                      <  2:368281‑M1/149‑124 (MQ=255)                                                                                                                             acTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaattgggatgttgcgcacagtcacaataactcaccgta                  >  1:304813‑M1/1‑22 (MQ=255)                                                                                                                               tCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtctatgcaggtgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaa                <  1:634801‑M1/149‑130 (MQ=255)                                                                                                                                  aTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataa             <  2:256487‑M1/149‑133 (MQ=255)                                                                                                                                      cTTAACATAAATGtattagtcatacttttgttttgggtgtatttcaatctgttaaaaagatttacgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaaaaatcac         <  2:846791‑M1/149‑137 (MQ=255)                                                                                                                                            aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgcta   >  2:60‑M1/1‑7 (MQ=255)                                                                                                                                             tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctat  <  1:89425‑M1/149‑144 (MQ=255)                                                                                                                                                  |                                                                                                                                                AAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCT  >  NC_000913/1978357‑1978645

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