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	51 [0]	[0] 51	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)	49 (1.150)	42/268	0.0	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

TCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGC  >  NC_000913/1978376‑1978641                                                                                                                               |                                                                                                                                            tCTAAGACTTCGCCTACATCGACGCAGCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttt                                                                                                                             <  2:238738‑M1/143‑17 (MQ=255) gctAAGATTTTTCCTAATTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttt                                                                                                                             <  2:639728‑M1/142‑17 (MQ=255) gctAAGATGTTTCCTAAATCGACGCGACTGTACTCGTCACTACAGGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttt                                                                                                                             <  2:467948‑M1/142‑17 (MQ=255)   tAAGATTTTTCCTAGATCGACGCAACTGTTCTCGTCACTACACGCACATACAGCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACTCTCAATTACTTAACATAAATGtataagtcatacttttgt                                                                                                                           <  2:455694‑M1/143‑19 (MQ=255)   tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgt                                                                                                                           >  2:738514‑M1/1‑125 (MQ=255)    aaGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtt                                                                                                                          >  1:380437‑M1/1‑124 (MQ=255)       ttttttcCCAAGTCGACGCAACTGCTCTCGTCACTACACGCACATACAACGGAGGGTGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                                                                       <  2:141420‑M1/142‑23 (MQ=255)       aattttccTATATCGACGCAACTGTACTCGTCACTACACGCACCTACAACGGAGGGGGGCTGCCATTTTCAATAATGCGTGATGCAGATCACACAAAACATTTAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                                                                       <  2:52574‑M1/141‑23 (MQ=255)       aTTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATAAAACGGAGGGGGGCTGCGATTTTCAAGAATACGTGAAGCAGATCACACAAAACACTCAATTACTCAACATAACTGtcgaagacaaacctttcatttc                                                                                                                       >  2:47109‑M1/1‑121 (MQ=255)              cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                                                                                                <  1:500216‑M1/143‑30 (MQ=255)               tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttggttttcggtgtatt                                                                                                               >  2:287146‑M1/1‑113 (MQ=255)                        gCAACTGTGCTCGTCACTACCCGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                                                                      <  2:272559‑M1/143‑40 (MQ=255)                              gTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                <  1:732543‑M1/143‑46 (MQ=255)                                    ctcaccccacgcacctaCAACGGCGCGGGGCTGCCATTTTCAATACTGCGTGATGCAGATCACACAAAACACTCAATTACTTAACACAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                          <  1:186950‑M1/136‑52 (MQ=255)                                     tCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGAGCACACAAAACACTCAATTACTTAACATACATGtatcagtcatacttttgatcagagtgtctttcaatatgttaaacagttattc                                                                                         >  2:419599‑M1/1‑91 (MQ=255)                                      cACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATGATACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaacagtttttcg                                                                                        <  2:470991‑M1/143‑54 (MQ=255)                                               cacaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactattgttttgggtgtatttcaatctgttcaaaagtttttcactacactag                                                                               >  2:43664‑M1/1‑81 (MQ=255)                                                   tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtgtttgtctacgctagcaag                                                                           >  1:20447‑M1/1‑77 (MQ=255)                                                    aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtgaaaaagtttttcgctacgcgagcacgc                                                                          >  2:499289‑M1/1‑76 (MQ=255)                                                      aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcat                                                                        >  1:490681‑M1/1‑74 (MQ=255)                                                      aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaa                                                                        >  1:490806‑M1/1‑74 (MQ=255)                                                         ggAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                     <  1:678453‑M1/143‑73 (MQ=255)                                                           aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatg                                                                   <  1:645643‑M1/143‑75 (MQ=255)                                                              ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                                                                <  1:500396‑M1/143‑78 (MQ=255)                                                               gggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttaaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaac                                                               >  2:479938‑M1/1‑65 (MQ=255)                                                                   tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTTACATAAATGtatcagtcatacttttgttttgggtgtatttcaatctgttaaaaagtatttcgctacgctagcaagcaaaaatgaaacagga                                                           >  1:216601‑M1/1‑61 (MQ=255)                                                                   tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagga                                                           >  1:594452‑M1/1‑61 (MQ=255)                                                                      gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaata                                                        >  2:256365‑M1/1‑58 (MQ=255)                                                                            cAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaa                                                  <  1:459799‑M1/143‑92 (MQ=255)                                                                                aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg                                              >  1:478009‑M1/1‑48 (MQ=255)                                                                                  tGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga                                            <  1:79488‑M1/143‑98 (MQ=255)                                                                                  tGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga                                            >  1:666797‑M1/1‑46 (MQ=255)                                                                                   gCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat                                           >  2:485921‑M1/1‑45 (MQ=255)                                                                                          gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaacaatgaaacaggaataatcgaaatgggatgttgcgc                                    >  1:756044‑M1/1‑38 (MQ=255)                                                                                           cAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgca                                   <  1:9550‑M1/143‑107 (MQ=255)                                                                                             gATCACACAAATCTCTGAATTACTTAACTTAGATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                 <  2:540123‑M1/143‑109 (MQ=255)                                                                                              agcacacaaaCCACTCAATTACTTATCATAAATGtataagtcatacttttgtattgaatgaaattcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacag                                <  1:780704‑M1/141‑110 (MQ=255)                                                                                                cacacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtc                              >  1:453769‑M1/1‑32 (MQ=255)                                                                                                 caacaaaACATTCAATTACTTAGCATATCTCtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtca                             <  2:780120‑M1/141‑113 (MQ=255)                                                                                                 acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaactgggctgttgcgcacagtca                             >  1:695396‑M1/1‑31 (MQ=255)                                                                                                     aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaat                         >  1:654820‑M1/1‑27 (MQ=255)                                                                                                     aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgtt                                        >  2:28993‑M1/1‑27 (MQ=255)                                                                                                     aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgtt                                        <  1:28993‑M1/128‑102 (MQ=255)                                                                                                            tCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcac                  >  1:657978‑M1/1‑20 (MQ=255)                                                                                                              aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatatgttaaaaagtttttcgctacgctagcaagcaaaaatgaagcaggaagaagcgaaatgggatgttgcggacagtcaaaataactaaccg                >  1:418051‑M1/1‑18 (MQ=255)                                                                                                              aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatgtcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccg                <  2:288178‑M1/143‑126 (MQ=255)                                                                                                                  aCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaa            <  1:738514‑M1/143‑130 (MQ=255)                                                                                                                    ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttttgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaacgtaaata          >  1:100664‑M1/1‑12 (MQ=255)                                                                                                                    ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaata          >  2:524119‑M1/1‑12 (MQ=255)                                                                                                                         aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcat     >  1:107515‑M1/1‑7 (MQ=255)                                                                                                                            aaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctg  <  1:521516‑M1/143‑140 (MQ=255)                                                                                                                               |                                                                                                                                            TCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGC  >  NC_000913/1978376‑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: ‑