Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
MC JC	NC_000913	1,978,503	Δ776 bp	100%		insB1–insA	insB1, insA

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1978503–1979269	1979278	10–776	56 [0]	[0] 56	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)	50 (0.720)	33/260	NT	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:1092825‑M1/1‑146 (MQ=255)    tGGGTCTTTACCCACATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataag                                                                                                                                           >  2:479280‑M1/1‑143 (MQ=255)          tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCCGATCACACAAAACCCTCAATTACTTAACATAAATGtataagtgctac                                                                                                                                     >  1:197988‑M1/1‑137 (MQ=255)          tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatac                                                                                                                                     >  1:29773‑M1/1‑137 (MQ=255)          tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataaatcatac                                                                                                                                     >  1:681313‑M1/1‑137 (MQ=255)              cacTTTTCTAAGATTTTTCCGATATCGACGCAACTGTACTCGTCACTAAACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactctt                                                                                                                                 <  2:857658‑M1/149‑17 (MQ=255)               aattatTTAATATTCTTCCTAAATCGACGCAACTGTCCTCGTCACTACACGCCCATGAAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcctacttttg                                                                                                                                >  2:364786‑M1/3‑132 (MQ=255)                     ctcagctttttcccaacTCGACGCACCTGTACTCGTCACTACACGCACATACAAGGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                          <  2:290265‑M1/143‑24 (MQ=255)                     cTAAGATTTTTCCTAAATCTACGCAACTGGCCTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATACTGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                          <  2:320599‑M1/149‑24 (MQ=255)                     cTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                          <  2:586902‑M1/149‑24 (MQ=255)                      tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGTTTAAGATGACACAAAACACTCAATTCCTTAACATAAATGtataagtcatacttttgttttggg                                                                                                                         >  1:946188‑M1/1‑125 (MQ=255)                      tAAGATTTTTCCTAAATCGACGAAACTGTACTCGCCCCTCCAGGCGCGTCCAACGCAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                                                                         <  1:297710‑M1/149‑25 (MQ=255)                          aTTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTCCTTAGCATAAAGGtataagtcatacttttgttttgggtgta                                                                                                                     >  1:1001076‑M1/1‑121 (MQ=255)                          aTTTTTCCTAAATCGACGCAACTGTACTCGTCACGACACGCACATACAACTTAGGAGGGCTTCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataaggcatactttcttttaaggtgta                                                                                                                     >  1:44364‑M1/1‑121 (MQ=255)                                         aCGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaa                                                                                                      <  1:904890‑M1/149‑44 (MQ=255)                                                   aCTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcg                                                                                            >  1:218805‑M1/1‑96 (MQ=255)                                                       gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCAAACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttattcgctac                                                                                        >  1:984234‑M1/1‑92 (MQ=255)                                                          aCTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgct                                                                                     >  1:94992‑M1/1‑89 (MQ=255)                                                             acacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagc                                                                                  <  1:1192168‑M1/149‑64 (MQ=255)                                                             acacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcgacgctagc                                                                                  <  2:300937‑M1/149‑64 (MQ=255)                                                                ctctcatacaACTGATGGTGGCAGCGATTTTCAATAACGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtatatggtatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag                                                                               <  1:268970‑M1/145‑67 (MQ=255)                                                                 gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaggc                                                                              >  1:1171991‑M1/1‑82 (MQ=255)                                                                 gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagc                                                                              >  2:1054750‑M1/1‑82 (MQ=255)                                                                 gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGAGGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagc                                                                              >  2:1283699‑M1/1‑82 (MQ=255)                                                                 gCACACACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagc                                                                              <  1:154000‑M1/149‑68 (MQ=255)                                                                       aCAACGGAGGGGGGCTGCGATTTTTAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                                                                        <  1:191190‑M1/149‑74 (MQ=255)                                                                       aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGAACACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                                                                        <  2:671193‑M1/149‑74 (MQ=255)                                                                       aCAACGGAGGGGGGCTGCGATTTTACATAATGCGAGATGCAGATCACAAAGATCACTCATTTACTTAACGTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                                                                        <  2:74139‑M1/149‑74 (MQ=255)                                                                              aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggttgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                                 <  2:842104‑M1/149‑81 (MQ=255)                                                                                 ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACATAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaa                                                              <  1:31893‑M1/149‑84 (MQ=255)                                                                                  gggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaat                                                             >  2:740486‑M1/1‑65 (MQ=255)                                                                                  gggCTGCGATTTTCAATAAGGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaat                                                             <  1:811480‑M1/149‑85 (MQ=255)                                                                                             ttCAATAATGCGTGATGCAGATCAAACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg                                                  >  2:345776‑M1/1‑54 (MQ=255)                                                                                                  taatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgt                                             <  1:522462‑M1/149‑101 (MQ=255)                                                                                                  taatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgt                                             <  1:521884‑M1/149‑101 (MQ=255)                                                                                                     tGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgc                                          <  1:522965‑M1/149‑104 (MQ=255)                                                                                                      gCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgttacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcg                                         <  1:892254‑M1/149‑105 (MQ=255)                                                                                                        ggcatgcAGAACACACAAAGCGCTCAATTACTTAACATAAATGtataagtcatacttttgttttgcgtgtatttcaatctgttaaaaagtttttcgttacaaaaaaaaaaaaaaatgaaacaggaataatcgaaatgggatgttgcgca                                       <  2:1263438‑M1/146‑107 (MQ=255)                                                                                                                gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaa                               >  1:68973‑M1/1‑35 (MQ=255)                                                                                                                 aTCACACAAAACACTCAATTACTTAACATTAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagggtgtcgctacgctagcaggcaaagatgaaacgggagtagtcgaaatgggatgtggcgcagagtcaaat                              >  1:912215‑M1/1‑34 (MQ=255)                                                                                                                 aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                              >  2:988519‑M1/1‑34 (MQ=255)                                                                                                                 aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                              >  2:1040799‑M1/1‑34 (MQ=255)                                                                                                                    acacaAAACACACAATTACTAATAATAAATGtataagtcatacttttgttttgggtgttttgttatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                           <  1:623460‑M1/149‑119 (MQ=255)                                                                                                                      acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataact                         >  1:78778‑M1/1‑29 (MQ=255)                                                                                                                      acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataact                         >  1:375‑M1/1‑29 (MQ=255)                                                                                                                           acacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccg                    <  1:79163‑M1/149‑126 (MQ=255)                                                                                                                           acacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataacgcaccg                    <  2:910380‑M1/149‑126 (MQ=255)                                                                                                                                cAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaat               <  2:1001076‑M1/149‑131 (MQ=255)                                                                                                                                  aTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgctcagtcaaaataactcaccgtaagtga             >  1:36369‑M1/1‑17 (MQ=255)                                                                                                                                  aTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataa             >  1:1270626‑M1/1‑17 (MQ=255)                                                                                                                                     aCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatca          >  1:398057‑M1/1‑14 (MQ=255)                                                                                                                                           cATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgct    <  2:510587‑M1/149‑142 (MQ=255)                                                                                                                                            atcaatggagaagtcatacattggttttgggtgtatttcaaactgttaaaaagtttttcgctacgaaaaaaaaaaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgcta   <  2:1219658‑M1/146‑143 (MQ=255)                                                                                                                                            aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgcta   >  2:495439‑M1/1‑7 (MQ=255)                                                                                                                                            aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgagt   >  2:186851‑M1/1‑7 (MQ=255)                                                                                                                                             tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctat  >  1:1016052‑M1/1‑6 (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: ‑