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	574781–575790	576776	987–1996	7 [6]	[0] 16	[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)	16 (0.350)	15/488	1.6	100%	noncoding (1/1195 nt)	IS5	repeat region
	?	NC_000913	576777 =	0 (0.000)					pseudogene (49/1040 nt)	nmpC	putative outer membrane porin NmpC

GGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCAAGATGTAAGTTTTGGTTTCTCGTTTGATTTGTGAAGTAGATCTCTATTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCA  >  NC_000913/576567‑577024                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        ccttggggtaaaacggcccgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttccggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAAt                                                                                                                                                                                                                 >  2:437460‑M2/211‑251 (MQ=255)                     atggccttctccatgttttgccatggcagaatgttgtccatgggggacaagaaaatgtctgtacaggtctgacggcgcttactgctgaattcactgtcggcgaaggcaagttgatggctcatgatgacccctgttctatggctccagatgacaaacaagatctcatatcagggacttgttcgcaccttccATAGTACTGATGCAGCTACAGCAGCAATTGCCACTGTTAATTTTTTCATCGTGAGCcctaa                                                                                                                                                                                             <  2:14028‑M2/61‑3 (MQ=255)                       gacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgacttcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCttttt                                                                                                                                                                                           >  2:437416‑M2/189‑251 (MQ=255)                           tccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCttttttt                                                                                                                                                                                         <  1:60445‑M2/65‑1 (MQ=255)                              accatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcagagcagggacttgttcgcaccttccTGAGTACTGATGCAGCTACAGCTGAAAGTGCCACTGTTAATTGTTGCATCTTGGGCCCTTTTTTTGAACt                                                                                                                                                                                    <  1:178743‑M2/70‑1 (MQ=255)                                  tgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtacgttgatgactcatgatgaaccatgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGCACTGATGCCGCTACAGCAGAACTTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAActattt                                                                                                                                                                                >  2:309678‑M2/178‑250 (MQ=255)                                                               cgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAat                                                                                                                                                   >  2:96282‑M2/149‑251 (MQ=255)                                                                                             cggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCAAGa                                                                                                                     <  2:199189‑M2/133‑1 (MQ=255)                                                                                                                cactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTGTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCAAGATGTAAGTTTTGGTTtctc                                                                                                   <  1:43149‑M2/151‑1 (MQ=255)                                                                                                                      cggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGGGATTATTTCAAGATGTAAGTTTTGGTTTCTCGTTTGa                                                                                             <  1:34816‑M2/157‑1 (MQ=255)                                                                                                                                 agttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCAAGATGTAAGTTTTGGTTTCTCGTTTGATTTGTGAAGt                                                                                   >  1:302230‑M2/83‑249 (MQ=255)                                                                                                                                            catgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCAAGATGTAAGTTTTGGTTTCTCGTTTGATTTGTGAAGTAGATCTCTAttt                                                                       <  2:192547‑M2/179‑1 (MQ=255)                                                                                                                                                                                               gggacttgttcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCAAGATGTAAGTTTTGGTTTCTCGTTTGATTTGTGAAGTAGATCTCTATTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGaataaa                   <  1:437416‑M2/231‑1 (MQ=255)                                                                                                                                                                                                        tcgcaccttccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCAAGATGTAAGTTTTGGTTTCTCGTTTGATTTGTGAAGTAGATCTCTATTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATc           <  1:437460‑M2/239‑1 (MQ=255)                                                                                                                                                                                                                tccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAAGTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAAGCAATGTGATTATTTCAAGATTTAAGTTTTGGTTTCTCGTTGGAGTTGTGAAGTAGATCTCTATTTTTATCTGAACTTTTTTCTTGCGAATCTTATTCATGGCGCTTGGCTGAATAAAAATGAAGTTATTAGcc   <  2:441816‑M2/247‑1 (MQ=255)                                                                                                                                                                                                                 ccTTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCAAGATGTAAGTTTTGGTTTCTCGTTTGATTTGTGAAGTAGATCTCTATTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAATTCTATTAGCCa  >  2:190355‑M2/3‑250 (MQ=255)                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        GGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTTCATCGTGAGCCCTTTTTTTGAACTATTATTAAAAAATGATGTCACTGCGCGATAAATATTCATCTAATCAATGTGATTATTTCAAGATGTAAGTTTTGGTTTCTCGTTTGATTTGTGAAGTAGATCTCTATTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCA  >  NC_000913/576567‑577024

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