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	1979278	776	34 [0]	[0] 34	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)	32 (0.840)	28/268	0.3	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

CTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGC  >  NC_000913/1978365‑1978641                                                                                                                                          |                                                                                                                                            ttttaGACTTATCTAAGTGGTTTCCTAAATCGACGCAACTGTACTCGTGACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGAGTGATGCAGATCACACAAAATACTCAATTACTTAACATAAATGtataa                                                                                                                                        <  1:262863‑M1/142‑6 (MQ=255) cTTTACCCTTATCGAAGATTTTTTCTAAATCGACGCAACTGTACTCTTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                        <  1:93678‑M1/143‑6 (MQ=255)  tttACACTTATCTAAGATTTTTCCGAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataaa                                                                                                                                       >  2:102385‑M1/1‑137 (MQ=255)   ttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagt                                                                                                                                      <  2:579106‑M1/143‑8 (MQ=255)         ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATACAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacc                                                                                                                                >  1:284627‑M1/1‑130 (MQ=255)              tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgt                                                                                                                           >  2:199133‑M1/1‑125 (MQ=255)              tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgt                                                                                                                           >  1:334333‑M1/1‑125 (MQ=255)                  atctttcctAAAGCGACGCAACTGTACGCGTAACTACACGCACATACGACGGAGGGGGGCTGCGATTTTCATTAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                                                                       <  1:391811‑M1/140‑23 (MQ=255)                           aaaTCGAGGCAGTTGTACTCGTCCCTACACGCGCATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattt                                                                                                              <  1:627990‑M1/143‑32 (MQ=255)                              tCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaa                                                                                                           >  2:446551‑M1/1‑109 (MQ=255)                                         gaactCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAACACTCTCAATTACTTACCATAAATGtataagtcatacttttgttttgggtgtattttaatctgttaaaaa                                                                                                <  1:138674‑M1/141‑46 (MQ=255)                                            cTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttaatgtagggtatatatcaaaaagttaaaaagtt                                                                                             >  2:310184‑M1/1‑95 (MQ=255)                                             tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttt                                                                                            >  1:29080‑M1/1‑94 (MQ=255)                                               gTCACTACACGCACATACAACTGAGGGGGGCTGCGATTTTCAATTATGCGTGATGCAGATCACACAAAACTCTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                          <  1:298734‑M1/143‑52 (MQ=255)                                                 cACTACACGCACATAAAACGGAGGTGTGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAAGGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcg                                                                                        <  2:48726‑M1/143‑54 (MQ=255)                                                       acGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgc                                                                                  <  1:199133‑M1/143‑60 (MQ=255)                                                             agacaACGGAGGGGGGCGGCGAGTTTCGATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaa                                                                            <  2:55564‑M1/141‑66 (MQ=255)                                                               aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattt                                                                                                              <  1:523242‑M1/107‑32 (MQ=255)                                                               aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattt                                                                                                              >  2:523242‑M1/1‑76 (MQ=255)                                                                    ggAGGGGTGCTGCGATTTTCAATAATGCGTTATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                     <  1:305463‑M1/143‑73 (MQ=255)                                                                    ggAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcactacgctagcaggcaaaaa                                                                     >  1:55025‑M1/1‑71 (MQ=255)                                                                      aGGGGGGCTGCGATTTTCAATAATGCGTGTTGCTGATTATACATAATACTTAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatg                                                                   <  2:562836‑M1/143‑75 (MQ=255)                                                                         ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtctaagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaaacaaaaatgaaa                                                                >  1:145086‑M1/1‑66 (MQ=255)                                                                         ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                                                                >  2:632451‑M1/1‑66 (MQ=255)                                                                             cTGCGGTTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgcaacagg                                                            >  2:552081‑M1/1‑62 (MQ=255)                                                                              tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacaatagcaagcaaaaatgaaacagga                                                           >  2:656351‑M1/1‑61 (MQ=255)                                                                               gCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaa                                                          >  2:715503‑M1/1‑60 (MQ=255)                                                                                                                aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttctctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaat                         >  2:559273‑M1/1‑27 (MQ=255)                                                                                                                  aaCACTTAATTACTTAACATAAATGtataagtcatacttttgtttttgttgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                       <  2:356651‑M1/143‑119 (MQ=255)                                                                                                                      cTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggaggttgcgcacaggccaaataactca                   >  2:424391‑M1/1‑21 (MQ=255)                                                                                                                       tCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctaagctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttacgcacaataaaaataacacac                  >  1:411346‑M1/1‑20 (MQ=255)                                                                                                                              cTTAACATAATTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaat           <  1:102385‑M1/143‑131 (MQ=255)                                                                                                                                   cATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatca      >  2:664172‑M1/1‑8 (MQ=255)                                                                                                                                       aacgtataagtcatacttttgttttgggtgtatctcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctg  <  1:145920‑M1/140‑140 (MQ=255)                                                                                                                                          |                                                                                                                                            CTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGC  >  NC_000913/1978365‑1978641

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