Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,978,503	Δ776 bp		insB5–insA5	insB5, insA5

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1978503–1979270	1979278	9–776	38 [0]	[1] 39	insB5–insA5	insB5,insA5

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1978502	0 (0.000)	38 (0.350)	35/466	1.1	100%	intergenic (‑305/+16)	flhD/insB5	DNA‑binding transcriptional dual regulator FlhD/IS1 family protein InsB
	?	NC_000913	1979279 =	0 (0.000)					intergenic (‑64/‑474)	insA5/uspC	IS1 family protein InsA/universal stress protein C

CCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCAT  >  NC_000913/1978269‑1978731                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                      ccGGTAAAAAATAATTAGCATTAGAATAGTTGCGATACGCTGCAATAAGAAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatccttttg                                                                                                                                                                                                                      >  2:500242‑M1/1‑234 (MQ=255)      aaaaaaTAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCGCACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                                                                                                                                                                 >  2:187181‑M1/1‑229 (MQ=255)      aaaaaaTAATTAGCATAAGAATAGTTGCGATAAGCTGCTATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                                                                                                                                                                 >  2:597680‑M1/1‑229 (MQ=255)              atTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                                                                                                                                                                                         >  2:16139‑M1/1‑221 (MQ=255)                      tAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatct                                                                                                                                                                                                  <  2:58746‑M1/250‑38 (MQ=255)                       aGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCATATCTCACAAAACACTCAATTACTTAACATAAATGtataaggcatacttttgttttgggtggatttcaatctgt                                                                                                                                                                                                <  2:430600‑M1/251‑40 (MQ=255)                          aTAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtta                                                                                                                                                                                              >  1:595142‑M1/1‑209 (MQ=255)                                     aaGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcg                                                                                                                                                                                  >  2:224237‑M1/1‑198 (MQ=255)                                     aaGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcg                                                                                                                                                                                  >  1:772304‑M1/1‑198 (MQ=255)                                          gCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACGACACGCACATACAACGGAGGGGGGCTTCGATTTTCAATAATGCGTGATGCATATCACACAAATCTCTCAATTACTTAACATAAATGtataattcattcttttgttttgggtgtattgcaatctgttaaaaattttttcgctacg                                                                                                                                                                             <  2:924408‑M1/251‑59 (MQ=255)                                                  cAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATCGAACTGGGTCTTTAACCTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGAACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggagtatttcaatctgataaaaagtttttcgctacgctagca                                                                                                                                                                       >  1:75876‑M1/1‑185 (MQ=255)                                                                                               gTCTTTACACTTAACTAAGATTTTTCCTAAATCGACGCAAATGTACTCGTCACTACACGAACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagaaaaaatgaaacaggaataatcgaactgggatgttgcgcaccgt                                                                                                                         >  2:537187‑M1/1‑140 (MQ=255)                                                                                                   ttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaa                                                                                                                     >  2:308931‑M1/1‑136 (MQ=255)                                                                                                     acacTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGATGGGGGCTGCGATTTTCAATAATGCGTGAGGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaata                                                                                                                  <  1:500242‑M1/251‑118 (MQ=255)                                                                                                                             aTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgct                                                                                          <  2:595142‑M1/251‑142 (MQ=255)                                                                                                                                   gCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaat                                                                                    >  2:29915‑M1/1‑104 (MQ=255)                                                                                                                                      aCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataat                                                                                 <  2:838496‑M1/251‑151 (MQ=255)                                                                                                                                           aCTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcctacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgcaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactt                                                                            >  1:905884‑M1/1‑96 (MQ=255)                                                                                                                                              cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttca                                                                         <  1:16139‑M1/251‑159 (MQ=255)                                                                                                                                               gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATACTGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcat                                                                        >  1:174282‑M1/1‑92 (MQ=255)                                                                                                                                                      cacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCGCACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaatac                                                                 <  1:187181‑M1/251‑167 (MQ=255)                                                                                                                                                      cacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgcta                                                                                         <  1:831709‑M1/227‑143 (MQ=255)                                                                                                                                                      cacGCACATAAAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTCCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgcta                                                                                         >  2:831709‑M1/1‑85 (MQ=255)                                                                                                                                                                cAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagct                                                       <  1:308931‑M1/251‑177 (MQ=255)                                                                                                                                                                   cGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATGACTTAACATAAATGtataagtcatactttggtttttggtgtatttcaatctgttaaaaagtttttcgctactctagcaagcaagaatgaaacaggaataggcgaaatgggatggtgcgctcagtcaaaatagctcaccgtaaaggatcagcggctgaaaataatcactttcatgcaataccagataagctatt                                                    <  1:52112‑M1/251‑180 (MQ=255)                                                                                                                                                                      aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctattttt                                                 >  2:542772‑M1/1‑69 (MQ=255)                                                                                                                                                                      aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcat                                                                        >  2:786033‑M1/1‑69 (MQ=255)                                                                                                                                                                      aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcat                                                                        <  1:786033‑M1/228‑160 (MQ=255)                                                                                                                                                                        gggggCTGCGATGTTCAAGAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcattcttttgttttgggtgtatttcaatctgttaaaaagtttttcgcgacgctagcaagcaaaaatgaaacaggaagaatcgaaatggtatgttgcgcacagtcaaaataactcaccgtaaataatcatctgttagaaataagcgcagtcattcaataccagataagctatgtttaa                                               <  2:905884‑M1/251‑185 (MQ=255)                                                                                                                                                                               gCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacaga                                          <  1:542772‑M1/249‑190 (MQ=255)                                                                                                                                                                                cGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacac                                       <  1:29915‑M1/251‑193 (MQ=255)                                                                                                                                                                                         ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacacttaccgcac                              <  1:224237‑M1/251‑202 (MQ=255)                                                                                                                                                                                                     gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacacttaccgcacaacaaactaatt                  >  2:234958‑M1/1‑38 (MQ=255)                                                                                                                                                                                                         aTCACACAAACCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacacttaccgcacaacaaactaattaaat              >  1:299923‑M1/1‑34 (MQ=255)                                                                                                                                                                                                           cacacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacacttaccgcacaacaaactaattaaatac            >  2:414034‑M1/1‑32 (MQ=255)                                                                                                                                                                                                             cacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacacttaccgcacaacaaactaattaaatacaa          <  2:772304‑M1/251‑222 (MQ=255)                                                                                                                                                                                                                   acacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacacttaccgcacaacaaactaattaaatacaaagaa      >  1:205427‑M1/1‑24 (MQ=255)                                                                                                                                                                                                                     acTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacacttaccgcacaacaaactaattaaatacaaagaaaggc  <  2:299923‑M1/251‑230 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                      CCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCAT  >  NC_000913/1978269‑1978731

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑