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	46 [0]	[0] 46	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)	43 (1.020)	40/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

CTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGC  >  NC_000913/1978365‑1978644                                                                                                                                          |                                                                                                                                               cTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                           >  1:524467‑M1/1‑138 (MQ=255)  tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataag                                                                                                                                          >  1:325641‑M1/1‑137 (MQ=255)     ccacgTATCTAAGATTTTTCCAAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCTTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                                                                       <  2:695637‑M1/142‑10 (MQ=255)        cTTATCTCAGATTTTTCCTAAATCGACGGACCTGCACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatac                                                                                                                                    <  2:587540‑M1/143‑13 (MQ=255)        cTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatac                                                                                                                                    <  1:659179‑M1/143‑13 (MQ=255)            tCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATACTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATAtataagtcatactgtt                                                                                                                                >  1:603180‑M1/1‑127 (MQ=255)                   tttgtCCTAAATCGACGCAACGGGACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCAATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                         <  2:793580‑M1/143‑24 (MQ=255)                    ttttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                                                                        >  1:790934‑M1/1‑119 (MQ=255)                         cTAAATCGACGCAACTGTACTCGTCACTGCACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                                                                                                   <  2:214543‑M1/143‑30 (MQ=255)                         cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGGGATTTTCAGTAATGCGGGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgggttggaggtaa                                                                                                                   >  1:735800‑M1/1‑114 (MQ=255)                          tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                                                                                                  >  1:667180‑M1/1‑113 (MQ=255)                          tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtaatacttttattttgggtgtatt                                                                                                                  >  2:677421‑M1/1‑113 (MQ=255)                            aaTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacctttgttttgggtgtactat                                                                                                                >  1:176597‑M1/1‑111 (MQ=255)                                        tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaa                                                                                                    <  1:283270‑M1/143‑45 (MQ=255)                                         gTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                   >  1:777215‑M1/1‑98 (MQ=255)                                               gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttacaaagttttt                                                                                             >  2:325045‑M1/1‑92 (MQ=255)                                                    tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                                                                        <  2:249897‑M1/143‑57 (MQ=255)                                                     tctcccacaaacaaTGGGGGGGGACGGCGAGTATCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctac                                                                                       <  1:178376‑M1/138‑58 (MQ=255)                                                      accgcACATACAGCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTGCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacg                                                                                      <  2:62493‑M1/141‑59 (MQ=255)                                                       acGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgccacgc                                                                                     >  1:237418‑M1/1‑84 (MQ=255)                                                         gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta                                                                                   >  2:291746‑M1/1‑82 (MQ=255)                                                            caTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgctaaaaagtttttcgctacgctagca                                                                                <  2:790934‑M1/143‑65 (MQ=255)                                                               aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcgacgctagcaagc                                                                             >  1:241102‑M1/1‑76 (MQ=255)                                                                   cGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttaaatctgtcaaaaaggttttcgctacgctagcaagcaaaa                                                                         >  1:543631‑M1/1‑72 (MQ=255)                                                                    ggAGGCGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttcgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                        <  1:9451‑M1/143‑73 (MQ=255)                                                                            gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                                >  2:429763‑M1/1‑63 (MQ=255)                                                                              tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagga                                                              >  1:16801‑M1/1‑61 (MQ=255)                                                                                cGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaat                                                            >  1:176234‑M1/1‑59 (MQ=255)                                                                                      tCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaagatgaaacaggactaatcga                                                      >  2:650513‑M1/1‑53 (MQ=255)                                                                                        aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgtgtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa                                                    >  2:132622‑M1/1‑51 (MQ=255)                                                                                        aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa                                                    <  2:288095‑M1/143‑93 (MQ=255)                                                                                        aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaaacagataaaaagtttttcgctacgctagcaagcacaaatgaaacaggaataatcgaaa                                                    >  1:586546‑M1/1‑51 (MQ=255)                                                                                                gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgt                                            >  1:538539‑M1/1‑43 (MQ=255)                                                                                                       agctcacACCAAACACTCGAATACGTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcgacgccagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcac                                     <  2:735800‑M1/140‑108 (MQ=255)                                                                                                       aGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcac                                     >  1:537429‑M1/1‑36 (MQ=255)                                                                                                         aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacag                                   >  1:18235‑M1/1‑34 (MQ=255)                                                                                                         aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacag                                   <  2:26026‑M1/143‑110 (MQ=255)                                                                                                              acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaa                              >  1:388598‑M1/1‑29 (MQ=255)                                                                                                                     acTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactc                       <  1:291746‑M1/143‑122 (MQ=255)                                                                                                                       tCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcac                     >  2:578619‑M1/1‑20 (MQ=255)                                                                                                                             aCTTAACATAAATGtataagtcatacttttgttttggttgtatttcactcttataaatagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaa               <  2:524604‑M1/143‑130 (MQ=255)                                                                                                                               ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaacta             >  2:47832‑M1/1‑12 (MQ=255)                                                                                                                                  aCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatc          <  1:698588‑M1/143‑135 (MQ=255)                                                                                                                                     tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatc       >  2:716851‑M1/1‑6 (MQ=255)                                                                                                                                        aTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgc    <  2:176234‑M1/143‑141 (MQ=255)                                                                                                                                          gtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgcta  >  1:384017‑M1/1‑1 (MQ=255)                                                                                                                                          |                                                                                                                                               CTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGC  >  NC_000913/1978365‑1978644

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