Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	575,786	Δ991 bp	IS5‑mediated	nmpC ←	putative outer membrane porin NmpC

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	574779–575790	576776	987–1998	10 [8]	[0] 39	[insH2]–[nmpC]	[insH2],[nmpC]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 275149	NA (NA)	39 (0.590)	37/488	0.7	100%	noncoding (1/1195 nt)	IS5	repeat region
	?	NC_000913	576777 =	0 (0.000)					pseudogene (49/1040 nt)	nmpC	putative outer membrane porin NmpC

TTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCAAGATGTAAGTTTTGGTTTCTCGTTTGATTTGTGAAGTAGATCTCTATTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAAT  >  NC_000913/576542‑577026                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                          ggataaggtcgccggccattaccagccttggggtaaaacggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGc                                                                                                                                                                                                                                            <  1:162913‑M2/16‑1 (MQ=255)     aaggtcgccggccattaccagccttggggtaaaacggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTAc                                                                                                                                                                                                                                         >  1:431195‑M2/232‑250 (MQ=255)       ggtcgccggccattaccagccttggggtaaaacggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattaactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTAcagc                                                                                                                                                                                                                                      >  1:466092‑M2/230‑251 (MQ=255)         tcgccggccattaccagccttggggtaaaacggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTAcagc                                                                                                                                                                                                                                      >  2:104699‑M2/228‑249 (MQ=255)              ggccattaccagccttggggtaaaacggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAAtt                                                                                                                                                                                                                               <  2:18386‑M2/29‑1 (MQ=255)                  attaccagccttggggtaaaacggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCa                                                                                                                                                                                                                           <  1:97527‑M2/33‑1 (MQ=255)                     accagccttggggtaaaacggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTg                                                                                                                                                                                                                        <  1:143152‑M2/36‑1 (MQ=255)                                ggtaaaacggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgagctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCATAAATTGCCACTGGTAATTTTTTc                                                                                                                                                                                                             <  2:444548‑M2/47‑1 (MQ=255)                                        ggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAg                                                                                                                                                                                                     <  1:338604‑M2/55‑1 (MQ=255)                                            cggtgacctcccccatgttttgccatggcagaatctgctccatgggggacaagcaaatctcttttttggtctgacggcgcttactgctgaattcactgtcggcgtaggtacgttgatgactcatgatgaaccctgttctatggctccagatgacaaacaatatctcatatccgggactttttcgccccttccTTAGTACTGATGCAGCTACAGCAGAAATTTCCCCTGTTAATTTTTTTCTCGGGAgcccc                                                                                                                                                                                                 >  2:564808‑M2/193‑250 (MQ=255)                                             gatgacttccaccatgttttgccatggcagaatctgcgccatgcgggacaagaaaatctcttttctgttctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatcgcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCGACAGGTGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCtt                                                                                                                                                                                                <  1:623721‑M2/60‑1 (MQ=255)                                                gacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCttttt                                                                                                                                                                                             <  1:337798‑M2/63‑1 (MQ=255)                                                     ccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgacccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGACATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGa                                                                                                                                                                                         >  1:222951‑M2/184‑250 (MQ=255)                                                        ccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACta                                                                                                                                                                                     <  1:599712‑M2/71‑1 (MQ=255)                                                           tgttttgccatggcagaatctgctccatgagggacaagaaaatctcttttctggtctgacggcgcttactgctgcattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTACTTTTTTCATCGTGAGCCCTTTTTTTGAActattt                                                                                                                                                                                  >  2:328528‑M2/178‑250 (MQ=255)                                                            gttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACtattat                                                                                                                                                                                 >  2:277891‑M2/177‑251 (MQ=255)                                                                       gcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTGTTTCAGCGTGAGCCCTTTTTTTGAACTAGTAGTAAAAAatgat                                                                                                                                                                      <  2:287022‑M2/86‑1 (MQ=255)                                                                        cagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAatgatg                                                                                                                                                                     <  2:235278‑M2/87‑1 (MQ=255)                                                                         agaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGt                                                                                                                                                                    >  1:448256‑M2/164‑251 (MQ=255)                                                                                gctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTgcg                                                                                                                                                             >  1:92700‑M2/157‑251 (MQ=255)                                                                                           gacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATAtt                                                                                                                                                  >  2:523053‑M2/146‑251 (MQ=255)                                                                                            acaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTc                                                                                                                                                 >  1:501767‑M2/145‑251 (MQ=255)                                                                                                 aaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCTCGATAAATATTCATCTa                                                                                                                                            <  2:623257‑M2/112‑1 (MQ=255)                                                                                                 aaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCt                                                                                                                                             >  2:25727‑M2/140‑250 (MQ=255)                                                                                                                gtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTTAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGattatt                                                                                                                             <  1:52435‑M2/127‑1 (MQ=255)                                                                                                                   tgacggcgcttactgctgaattcactgtcggcgaaggtacgttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTc                                                                                                                           >  1:669473‑M2/122‑250 (MQ=255)                                                                                                                    gacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCa                                                                                                                          >  1:129575‑M2/121‑250 (MQ=255)                                                                                                                    gacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCa                                                                                                                          <  1:523053‑M2/130‑1 (MQ=255)                                                                                                                       ggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATCGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCa                                                                                                                          >  1:498015‑M2/118‑247 (MQ=255)                                                                                                                         cgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCAAGATGt                                                                                                                    <  2:575869‑M2/136‑1 (MQ=255)                                                                                                                                                                           gaaccctgttcgatggctccagatgacaaacatgatctcatatcagggacgtgttcgcacctttcTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAAGGTGATTATTTCAAGATGTAAGTTTTGGTTTCTCGTTTGATGTGGGAAGTAGATCTCTATTTTTATctt                                                                  <  2:334810‑M2/186‑2 (MQ=255)                                                                                                                                                                                             ccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCAAGATGTAAGTTTTGGTTTCTCGTTTGATTTGTGAAGTAGATCTCTATTTTTATCTGAACTTTTTTCTATCGAAt                                                <  2:466092‑M2/204‑1 (MQ=255)                                                                                                                                                                                              cagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCAAGATGTAAGTTTTGGTTTCTCGTTTGATTTGTGAAGTAGATCTCTATTTTTATCTGAACTTTTTTCTATCGAATc                                               >  2:387391‑M2/47‑251 (MQ=255)                                                                                                                                                                                                   gacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCAAGATGTAAGTTTTGGTTTCTCGTTTGATTTGTGAAGTAGATCTCTATTTTTATCTGAACTTTTTTCTATCGAATCCTAtt                                          <  1:131943‑M2/210‑1 (MQ=255)                                                                                                                                                                                                   gacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCAAGATGTAAGTTTTGGTTTCTCGTTTGATTTGTGAAGTAGATCTCTATTTTTATCTGAACTTTTTTCTATCGAATCCTAt                                           <  2:669473‑M2/209‑1 (MQ=255)                                                                                                                                                                                                                         ggactagttcgcacattccTTAGAACTGATGCAGCTACAGCAGAAATTGCGACTGTTAATTCTTTAATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGCTAAATATTCATCTAATCAATGTGATTATTTCAAGATGTAAGTTTTGGTTTCTCGTTTGATTTGTGAAGTAGATCTCTATTTTTATCTGAACTTTTTTCAATCGAATCCTATTCATGGCTCTTGGCTGaataaa                     >  2:173130‑M2/20‑250 (MQ=255)                                                                                                                                                                                                                                    caccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCAAGATGTAAGTTTTGGTTTCTCGTTTGATTTGTGAAGTAGATCTCTATTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTAtt         <  1:24216‑M2/243‑1 (MQ=255)                                                                                                                                                                                                                                         tccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCAAGATGTAAGTTTTGGTTTCTCGTTTGATTTGTGAAGTAGATCTCTATTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGcc     >  1:471035‑M2/4‑250 (MQ=255)                                                                                                                                                                                                                                           cTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCAAGATGTAAGTTTTGGTTTCTCGTTTGATTTGTGAAGTAGATCTCTATTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAAt  >  2:229736‑M2/2‑251 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                          TTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCAAGATGTAAGTTTTGGTTTCTCGTTTGATTTGTGAAGTAGATCTCTATTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAAT  >  NC_000913/576542‑577026

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