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–1979270	1979278	9–776	59 [0]	[0] 58	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)	56 (1.160)	45/268	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

TGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGC  >  NC_000913/1978360‑1978641                                                                                                                                               |                                                                                                                                            tGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                                                                                             >  1:678891‑M1/1‑143 (MQ=255)   ggTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGta                                                                                                                                           >  2:18134‑M1/1‑141 (MQ=255)       tttACACTTATCTAAGATTTTTCCTAAATCGACGCGACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataag                                                                                                                                       >  1:719572‑M1/1‑137 (MQ=255)              ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCGCGCAAAACACTCAATTACGTAACATAGATGtatgagtgatact                                                                                                                                >  1:429699‑M1/1‑130 (MQ=255)                    aaGGTTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtt                                                                                                                          <  2:503447‑M1/143‑20 (MQ=255)                     aGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTGAACATAAATGtataagtcatacttttgttt                                                                                                                         >  1:206255‑M1/1‑123 (MQ=255)                        tttttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACACAGAACACAAGTACCTGACACAAAGGtttaagtcattcttttgttttgg                                                                                                                      >  1:125237‑M1/1‑120 (MQ=255)                        tttttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACCTAAATGtataagtcatacttttgttttgg                                                                                                                      >  1:124384‑M1/1‑120 (MQ=255)                        tttttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                      <  2:137867‑M1/143‑24 (MQ=255)                        tttttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                      <  2:652045‑M1/143‑24 (MQ=255)                              cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgatttgggtgtat                                                                                                                >  2:89313‑M1/1‑114 (MQ=255)                                         cAACTGTACTCGTCACTACACGCACATACAACGGAGGTGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                                                                     <  2:491245‑M1/143‑41 (MQ=255)                                         cAACTGGACTCGTCACTACACGCACAGACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                                                                     <  2:485146‑M1/143‑41 (MQ=255)                                             tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaa                                                                                                 <  2:130338‑M1/143‑45 (MQ=255)                                              gttctcgTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                <  2:571027‑M1/140‑46 (MQ=255)                                              gTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                <  1:847762‑M1/143‑46 (MQ=255)                                                  tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtcttcgctgtatttccatctgttacaaagttt                                                                                            >  2:145991‑M1/1‑94 (MQ=255)                                                    gTCACGACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTATCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                          <  2:232443‑M1/143‑52 (MQ=255)                                                        cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgct                                                                                      <  2:577495‑M1/143‑56 (MQ=255)                                                         tACACCCACATACACCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                                                                     <  1:558345‑M1/143‑57 (MQ=255)                                                          acacGCACATACAACGTAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctac                                                                                    <  1:332033‑M1/143‑58 (MQ=255)                                                           cacGCACATACAACACAGGGTGGCGGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacg                                                                                   <  1:457895‑M1/143‑59 (MQ=255)                                                            acGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaggtttttcgctacgc                                                                                  >  2:590636‑M1/1‑84 (MQ=255)                                                            acGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgc                                                                                  >  1:725422‑M1/1‑84 (MQ=255)                                                            acGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgc                                                                                  <  2:540031‑M1/143‑60 (MQ=255)                                                                acaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagc                                                                              <  1:577130‑M1/143‑64 (MQ=255)                                                                  aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACATATATCACTCAATTACTTAACATAAATGtataagtcatacttttgtttagggtgtatttcaatctgttaaaaagtttttcgctacgctagcaa                                                                            <  1:1083‑M1/143‑66 (MQ=255)                                                                            ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatga                                                                  <  2:682914‑M1/143‑76 (MQ=255)                                                                             gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaa                                                                 >  1:114794‑M1/1‑67 (MQ=255)                                                                             gggggCTGCGATTGTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATGAATGtataagtcatacttttggtttgggtgtattgcaagctgttaaaaagtttttcgctacgctagcaagcaaaaatgaa                                                                 <  2:421562‑M1/143‑77 (MQ=255)                                                                              ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                                                                <  1:518864‑M1/143‑78 (MQ=255)                                                                              gcggctGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                                                                <  1:185066‑M1/141‑78 (MQ=255)                                                                               gggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaaggaaac                                                               >  1:564240‑M1/1‑65 (MQ=255)                                                                                 gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                             >  1:640963‑M1/1‑63 (MQ=255)                                                                                 gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                             <  2:875311‑M1/143‑81 (MQ=255)                                                                                 gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                             <  1:780118‑M1/143‑81 (MQ=255)                                                                                   tGCGATTTTCAAGAATGCGTAATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagga                                                           <  2:278848‑M1/143‑83 (MQ=255)                                                                                      gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCCATTACTTAACATAAATGtataagtcatacttttgttttgggtgtctttcaatctgttaaaaagttgttcgctccgctagcaagcaagaatgaaacaggcagc                                                        >  1:431315‑M1/1‑58 (MQ=255)                                                                                       aTTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataa                                                       >  1:297687‑M1/1‑57 (MQ=255)                                                                                       aTTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataa                                                       <  2:321312‑M1/143‑87 (MQ=255)                                                                                         tttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatc                                                     >  1:215733‑M1/1‑55 (MQ=255)                                                                                             aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa                                                 >  1:311933‑M1/1‑51 (MQ=255)                                                                                                  tGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatattgttgttttgtgtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga                                            <  1:907378‑M1/143‑98 (MQ=255)                                                                                                  cgcGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga                                            <  1:323430‑M1/142‑98 (MQ=255)                                                                                                   gCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat                                           <  2:108905‑M1/143‑99 (MQ=255)                                                                                                         tGCAGATCACACAAAACACTCAATTACTTAACCTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcg                                     <  1:348712‑M1/143‑105 (MQ=255)                                                                                                             gATCACACAAACGCCTCAAGTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                 <  2:432349‑M1/143‑109 (MQ=255)                                                                                                              attacacaaaaCACACAATAACTTCACATACCTGtaaatgttacacttttgtcttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacag                                <  1:630389‑M1/140‑110 (MQ=255)                                                                                                              aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacag                                <  1:33470‑M1/143‑110 (MQ=255)                                                                                                                 acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatggtgcaaacagtac                             >  2:185345‑M1/1‑31 (MQ=255)                                                                                                                    caAAACTCTCAATTACTTAACATCAGTGtttgagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                          <  2:8488‑M1/143‑116 (MQ=255)                                                                                                                          acTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtatatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactc                    >  1:871477‑M1/1‑22 (MQ=255)                                                                                                                            tCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcac                  >  1:255805‑M1/1‑20 (MQ=255)                                                                                                                            tCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcac                  >  2:808795‑M1/1‑20 (MQ=255)                                                                                                                            tCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcactacgctagcaagcaaaaatgaaacaggaataatcgaactggggcgttgcgcacagtgaaaataaaaaaa                  >  1:568529‑M1/1‑20 (MQ=255)                                                                                                                                    ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgagatgttgcgcacagtcaaaataactcaccctaaata          >  2:620512‑M1/1‑12 (MQ=255)                                                                                                                                     tAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataa         >  1:32670‑M1/1‑11 (MQ=255)                                                                                                                                           aaaTGgataagtcgtaattttgtgtggagtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatct   <  1:613585‑M1/143‑139 (MQ=255)                                                                                                                                            aaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctg  >  2:219735‑M1/1‑4 (MQ=255)                                                                                                                                               |                                                                                                                                            TGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGC  >  NC_000913/1978360‑1978641

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