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	574812–575788	576776	989–1965	5 [3]	[0] 25	[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)	25 (0.740)	25/486	0.4	100%	noncoding (1/1195 nt)	IS5	repeat region
	?	NC_000913	576777 =	0 (0.000)					pseudogene (49/1040 nt)	nmpC	putative outer membrane porin NmpC

CTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCAAGATGTAAGTTTTGGTTTCTCGTTTGATTTGTGAAGTAGATCTCTATTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATT  >  NC_000913/576573‑577019                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   ggtaaaacggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATGGCCACTGTTAATTTTTTc                                                                                                                                                                                                      <  2:48610‑M2/47‑1 (MQ=255)       acggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGc                                                                                                                                                                                                                      >  1:312927‑M2/199‑229 (MQ=255)       acggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcagatcagggacttgttcgcaccttccGTAGTACTGATGCAGCTACAGCAGAAATTGc                                                                                                                                                                                                                      <  2:312927‑M2/31‑1 (MQ=255)         ggctcgatgacttccaccatggtttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAg                                                                                                                                                                                              <  1:287547‑M2/55‑1 (MQ=255)              gatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTg                                                                                                                                                                                                                 >  1:267075‑M2/192‑227 (MQ=255)              gatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTg                                                                                                                                                                                                                 <  2:267075‑M2/36‑1 (MQ=255)                                       ggcagcatatgctccatgaggaacaagaaaatctcttttctggtctgacggcgcttcctgctgaattcactgtcggcgaaggtaagttgatgcctcatgatgcaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAACatg                                                                                                                                                                 >  2:161407‑M2/167‑250 (MQ=255)                                          agaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGt                                                                                                                                                             >  1:77583‑M2/164‑251 (MQ=255)                                               ctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccaggtgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGAGAGCCCTTTTTTTGAACTATTATTAAAAAATGCTGTCACTg                                                                                                                                                        >  2:97545‑M2/159‑251 (MQ=255)                                                             acaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttcgatggctccagatgacaaacatgatctcatatcagggacttgtgcgcgccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATGTTTTCATCGTGAGCCCTTTTGTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTc                                                                                                                                          <  1:161407‑M2/107‑1 (MQ=255)                                                                agaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCTATAAATATTCATc                                                                                                                                       <  1:115253‑M2/110‑1 (MQ=255)                                                                         cttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATc                                                                                                                                  <  2:77583‑M2/115‑1 (MQ=255)                                                                           tttctggtctgacggcgcttactgctgaattcactgtcggcgcaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAAtgtg                                                                                                                            >  1:135507‑M2/131‑251 (MQ=255)                                                                                ggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgaggactcattatgaaccctgtgctatggctccagatgacaaacatgatctcatatcaggggcttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAAGTTTTTCATCGTGAGCCCTGTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGatta                                                                                                                        <  1:342461‑M2/125‑1 (MQ=255)                                                                                     gacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTc                                                                                                                    <  1:54263‑M2/129‑1 (MQ=255)                                                                                              tactgctgaattcactgtcggcgaaggtacgttgatgactaattatgaaccctattctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGCACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGTGCCCTTTTTTTGTACTATTCTTCAAAAATGATGTCATTGCGCGATAATTATTCATCTAATCAATGTGATTATTTCAAGCTGTAAGt                                                                                                         >  2:115408‑M2/112‑251 (MQ=255)                                                                                               actgcttaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAAGTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATGATTTCAAGATGTAAGt                                                                                                         <  1:133517‑M2/140‑1 (MQ=255)                                                                                                  gctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCAAGATGTAAGttt                                                                                                       <  2:135507‑M2/142‑1 (MQ=255)                                                                                                          cactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTTAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCAAGATGTAAGTTTTGGTTTCTCg                                                                                             <  1:39351‑M2/152‑1 (MQ=255)                                                                                                            ctgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCAAGATGTAAGTTTTGGTTTCTCGtt                                                                                           >  1:46300‑M2/98‑251 (MQ=255)                                                                                                                                            gaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCAAGATGTAAGTTTTGGTTTCTCGTTTGATTTGTGAAGTAGATCTCTATTTTTATCt                                                            <  1:125978‑M2/185‑1 (MQ=255)                                                                                                                                                              ccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACGGCGCTATAAATATTCATCTAATCAATGTGATTATTTCAAGATGTAAGTTTTGGTTTCTCGTTTGATTTGTGAAGTAGATCTCTATTTTTATCTGAACTTTTTTCTATCGAAt                                         <  1:222765‑M2/204‑1 (MQ=255)                                                                                                                                                                    gacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCAAGATGTAAGTTTTGGTTTCTCGTTTGATTTGTGAAGTAGATCTCTATTTTTATCTGAACTTTTTTCTATCGAATCCTAtt                                   <  2:127859‑M2/210‑1 (MQ=255)                                                                                                                                                                                                    gcaccttccTTAGTACTGATGCAGCTACAGCCGCCATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCAAGATGTAAGTTTTGGTTTCTCGTTTGATTTGTGAAGTAGATCTCTATTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTAt   >  2:234996‑M2/10‑251 (MQ=255)                                                                                                                                                                                                     caccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCAAGATGTAAGTTTTGGTTTCTCGTTTGATTTGTGAAGTAGATCTCTATTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTAtt  >  2:26201‑M2/9‑251 (MQ=255)                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   CTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCAAGATGTAAGTTTTGGTTTCTCGTTTGATTTGTGAAGTAGATCTCTATTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATT  >  NC_000913/576573‑577019

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