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	57 [0]	[1] 59	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)	54 (1.110)	44/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

GGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCG  >  NC_000913/1978361‑1978643                                                                                                                                              |                                                                                                                                              gggTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACCAACCACTCAATTCCTTAACATAAATGt                                                                                                                                              >  1:792233‑M1/1‑142 (MQ=255)      tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATCAATGtataag                                                                                                                                         >  1:770143‑M1/1‑137 (MQ=255)      tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataag                                                                                                                                         >  1:180635‑M1/1‑137 (MQ=255)      tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataag                                                                                                                                         >  2:631420‑M1/1‑137 (MQ=255)           acTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACCCTCACATACAACGGAGGGGGGCTGCGATTTTCAAAACTGCGCCATACAGATCGCACAAGACACTCAATTACTTAACATAAATGcataactccac                                                                                                                                    >  2:592349‑M1/1‑132 (MQ=255)            cTTATCTAAGAATTTTCCTAAAGCGACGCAACTGTACGCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatac                                                                                                                                   <  1:670223‑M1/143‑13 (MQ=255)                tCTAAGTTTCTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGTGCTGCGATTTTCAATAATTCGTGATGCAGATCACACAAAACTCTCAATTACTTAACATAAATGtataagtcatactttt                                                                                                                               <  2:479111‑M1/143‑17 (MQ=255)                  tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttga                                                                                                                             >  2:37508‑M1/1‑125 (MQ=255)                   aagtcttttcgcaaATCGACGCAACGGGAGTCGTCACCACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACGCTCAATTACTTAACATAAATGtataagtcatacttttgtt                                                                                                                            <  2:476286‑M1/138‑20 (MQ=255)                     gATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttt                                                                                                                          <  1:52947‑M1/143‑22 (MQ=255)                      aTTTTGTCTAAGTCGACGCAACTGTACTCGTCACTACATGCACATACAACTGAGGGTGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                                                                         <  2:247399‑M1/143‑23 (MQ=255)                       tttttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                        <  2:437300‑M1/143‑24 (MQ=255)                         tttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAACTGtataattcatacttttgttttggga                                                                                                                      >  1:752136‑M1/1‑118 (MQ=255)                             cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                                                                                                  >  1:830833‑M1/1‑114 (MQ=255)                              tAAATCGACGCAACTGTACTCGTCACTACCCGCACATACGACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACGCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                                                                                                 >  1:38738‑M1/1‑113 (MQ=255)                              tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACAAAAGGATGAGAAGTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                                                                                                 >  1:659607‑M1/1‑113 (MQ=255)                               aaaTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattt                                                                                                                <  2:827375‑M1/143‑32 (MQ=255)                               aaaTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattt                                                                                                                >  1:340552‑M1/1‑112 (MQ=255)                                 aTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttca                                                                                                              >  1:890834‑M1/1‑110 (MQ=255)                                   cGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                            <  1:156527‑M1/143‑36 (MQ=255)                                    gACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                                                           >  1:624808‑M1/1‑107 (MQ=255)                                     aCGCAACTGTACTCGTCACTACACGCACATACCACGGCGCGACCCTGCGATTTTCAATAATACGTGATGCAGATCACACAAAACACTCAATTACTTAACATAACTGtataagtcctacttctgctttcggtctatttcaatct                                                                                                          >  2:404781‑M1/1‑106 (MQ=255)                                        cAACTGTACCCGGCACTACACGCACACACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                                                                       <  1:93863‑M1/143‑41 (MQ=255)                                         aaCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACCCTCAATTCCTTAACATAAATGtataagtcatacttttgtttttgttctatttccatctcttc                                                                                                      >  1:454541‑M1/1‑102 (MQ=255)                                            tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaa                                                                                                   >  1:304315‑M1/1‑99 (MQ=255)                                             gTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                  <  2:1510‑M1/143‑46 (MQ=255)                                                 tCGTCACGACACGAACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttt                                                                                              <  1:86543‑M1/143‑50 (MQ=255)                                                        tccacgcACATAGCACTGAGGGTGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                                                                       <  2:771060‑M1/141‑57 (MQ=255)                                                        tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTCCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctc                                                                                       >  2:269955‑M1/1‑87 (MQ=255)                                                        tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                                                                       >  1:65340‑M1/1‑87 (MQ=255)                                                           acGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgc                                                                                    <  2:431983‑M1/143‑60 (MQ=255)                                                               acaTACAGCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACGCTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagc                                                                                <  2:289857‑M1/143‑64 (MQ=255)                                                                   aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtaatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagc                                                                            <  2:850917‑M1/143‑68 (MQ=255)                                                                     aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaa                                                                          >  1:526829‑M1/1‑74 (MQ=255)                                                                         gAGGGGGGCTGTGATTTTCAATAATGCGTGATGCAGATCACACAAACCACTCCATTACTTATCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                                                                      <  2:121088‑M1/143‑74 (MQ=255)                                                                         gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACCCTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                                                                      <  1:512862‑M1/143‑74 (MQ=255)                                                                            gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaatagtttttcgcttcgctcgcccgcatgaattca                                                                   >  1:288615‑M1/1‑67 (MQ=255)                                                                                    cGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaat                                                           >  1:158780‑M1/1‑59 (MQ=255)                                                                                    cGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaat                                                           >  2:130510‑M1/1‑59 (MQ=255)                                                                                      aTTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataa                                                         >  1:274349‑M1/1‑57 (MQ=255)                                                                                       ttttCAATAATGCGTGATGCAGATCACACTAATCACTCCAGTACCTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat                                                        <  1:165425‑M1/143‑88 (MQ=255)                                                                                          tCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcga                                                     <  1:695457‑M1/143‑91 (MQ=255)                                                                                           cAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaa                                                    >  1:346399‑M1/1‑52 (MQ=255)                                                                                                 ttggtgaTGGAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga                                              <  1:727726‑M1/140‑98 (MQ=255)                                                                                                 tGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgtgtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaacaatgaaacaggaataatcgaaatggga                                              >  2:439666‑M1/1‑46 (MQ=255)                                                                                                             aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacag                                  >  1:152866‑M1/1‑34 (MQ=255)                                                                                                              ccacacaaAACACTCAATTACTTAACATAAACGtataagtcatacttttgtttggggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagt                                 <  1:244337‑M1/142‑111 (MQ=255)                                                                                                                      aaCACTCAATTACTTAGCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                         >  2:744372‑M1/1‑25 (MQ=255)                                                                                                                             aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgagcacagtcaaaataactcaccg                  >  2:3935‑M1/1‑18 (MQ=255)                                                                                                                                   ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaata            >  1:437141‑M1/1‑12 (MQ=255)                                                                                                                                   ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatattgcgcacagtcaaaataagtcaccgtcaata            >  1:461129‑M1/1‑12 (MQ=255)                                                                                                                                   ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataaccaaaatgggatcgcgcacacagtcaaaataactcaccgtaaata            >  2:434275‑M1/1‑12 (MQ=255)                                                                                                                                     aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataat          >  1:15189‑M1/1‑10 (MQ=255)                                                                                                                                      aCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatc         >  1:185482‑M1/1‑9 (MQ=255)                                                                                                                                        aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgtgatgttgcgcacagtcaaaataactcaccgtccataatcat       >  2:152189‑M1/1‑7 (MQ=255)                                                                                                                                           aaTGtataagtcatccttgtgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggattttgcgcacagtcaaaataactcaccgtaaataatcatctg    <  1:305980‑M1/143‑140 (MQ=255)                                                                                                                                             tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgct  <  1:137680‑M1/143‑142 (MQ=255)                                                                                                                                              |                                                                                                                                              GGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCG  >  NC_000913/1978361‑1978643

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