Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,335,153	Δ593 bp		[ymiA]–yciX	[ymiA], yciX

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1335153	1335745	593	22 [0]	[2] 26	[ymiA]–yciX	[ymiA],yciX

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1335152	0 (0.000)	22 (0.590)	20/236	1.0	95.7%	coding (29/165 nt)	ymiA	uncharacterized protein
	?	NC_000913	1335746 =	2 (0.050)					intergenic (+288/‑85)	yciX/acnA	uncharacterized protein/aconitate hydratase 1; aconitase A

AAAAGATTGAATTATTCACTGGAGACGATTCGTAAATTCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGGTTAGCAATGCCTTCTGGAAATCAGGAACCCCGCCGCGATCCTGAATTAAAACGTAAAGCCTGGCTGGCGGTTTTTCTTGGTTCTGCACTTTTCTGGGTGGTTGTCGCACT  >  NC_000913/1335013‑1335263                                                                                                                                            |                                                                                                                aaaaGATTGAATTATTCACTGGAGACGATTCGTAAATTCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGaaatcgtta                                                                                                        >  5:15855‑M1/1‑140 (MQ=255)                  aCTGGAGACGATTCGTAAATTCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGaaatcgttac                                                                                                       >  3:176696‑M1/1‑123 (MQ=255)                  aCTGGAGACGATTCGTAAATTCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGaaatcgttac                                                                                                       <  4:176696‑M1/133‑11 (MQ=255)                                        gCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGaaatcgttacgcgatgtttgtgttatctttaatattcaccctgaa                                                                    >  6:137313‑M1/1‑101 (MQ=255)                                                     cTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGaaatcgttacgcgatgtttgtgttatctttaatattcaccctgaagagaatcagggcttcg                                                    >  7:116474‑M1/1‑88 (MQ=255)                                                          ggttctTTAGGATTTGTCTCATCGTATGATTAATTACACAAATCGTTCAACGTCTTGCTGATAAGTGATTGCGATTATATGAGaaatcgttacgcgatgtttgtgatatctttaatattcaccctgaagagaatcagatcttcgcaacc                                               >  7:38841‑M1/3‑83 (MQ=255)                                                                 tAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGaaatcgttacgcgatgtttgtgttatctttaatattcaccctgaagagaatcagggcttcgcaaccctgtca                                         <  6:15855‑M1/148‑73 (MQ=255)                                                                    gATTTGTGTCATCGTATGATTAATTACTCCAATCGTTCAACGTGTTTCTGATATGTGATGGCGATTATATGAGtaatcgttacgctatgtttgtgttatctttaatattcaccctgaagagaatctgggcttcgcaacgctgtcattaa                                     <  8:38841‑M1/149‑77 (MQ=255)                                                                    gATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGaaatcgttacgcgatgtttgtgttatctttaatattcaccctgaagagaatc                                                             <  5:162828‑M1/125‑53 (MQ=255)                                                                    gATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGaaatcgttacgcgatgtttgtgttatctttaatattcaccctgaagagaatc                                                             >  6:162828‑M1/1‑73 (MQ=255)                                                                          tctcATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGaaatcgttacgcgatgtttgtgttatctttaatattcaccctgaagagaatcagggcttcgcaaccctgtcattaa                                     >  8:19807‑M1/1‑67 (MQ=255)                                                                          tctcATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGaaatcgttacgcgatgtttgtgttatctttaatattcaccctgaagagaatcagggcttcgcaaccctgtcattaa                                     <  7:19807‑M1/143‑77 (MQ=255)                                                                              aTCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGaaatcgttacgcgatgtttgtgttatctttaatattcaccctgaagagaatcagggcttcgcaaccctgtcattaaggaggagct                            <  8:116474‑M1/148‑86 (MQ=255)                                                                                      attaattaCACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGaaatcgttacgcgatgtttgtgttatctttaatattcaccctgaagagaatcagggcttcgcaaccctgtcattaaggaggagctatgtcgtc                    <  2:86083‑M1/148‑94 (MQ=255)                                                                                      attaattaCACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGaaatcgttacgcgatgtttgtgttatctttaatattcaccctgaagagaatcagggcttcgcaaccctgtcattaaggaggagctatgtcgtc                    >  1:86083‑M1/1‑55 (MQ=255)                                                                                                       ttCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGaaatcgttacgcgatgtttgtgttatctttaatattcaccctgaagagaatcagggcttcgcaaccctgtcattaaggaggagctatgtcgtcaaccctacgagaagccag  <  5:137313‑M1/149‑112 (MQ=255)                                                                                                                           tGATGGCGATTATATGAGaaatcgttacgcgatgtttgtgttatctttaatattcaccctgaagagaatcagggcttcgcaaccctgtcattaaggaggagctatgtcgtcaaccctacgagaag      <  7:27876‑M1/125‑108 (MQ=255)                                                                                                                           tGATGGCGATTATATGAGaaatcgttacgcgatgtttgtgttatctttaatattcaccctgaagagaatcagggcttcgcaaccctgtcattaaggaggagctatgtcgtcaaccctacgagaag      >  8:27876‑M1/1‑18 (MQ=255)                                                                                                                               ggCGATTATATGAGaaatcgttacgcgatgtttgtgttatctttaatattcaccctgaagagaatcagggcttcgcaaccctgtcattaaggaggagctatgtcgtcaa                  <  6:69361‑M1/109‑96 (MQ=255)                                                                                                                               ggCGATTATATGAGaaatcgttacgcgatgtttgtgttatctttaatattcaccctgaagagaatcagggcttcgcaaccctgtcattaaggaggagctatgtcgtcaa                  >  5:69361‑M1/1‑14 (MQ=255)                                                                                                                                          gagaaatcgttacgcgatgtttgtgttatctttaatattcaccctgaagagaatcagggcttcgcaaccctgtcattaaggagga                               >  7:115822‑M1/1‑3 (MQ=255)                                                                                                                                          gagaaatcgttacgcgatgtttgtgttatctttaatattcaccctgaagagaatcagggcttcgcaaccctgtcattaaggagga                               <  8:115822‑M1/85‑83 (MQ=255)                                                                                                                                            |                                                                                                                AAAAGATTGAATTATTCACTGGAGACGATTCGTAAATTCGCTGGATTTTTGACTAAAGTTTCTTTAGGATTTGTCTCATCGTATGATTAATTACACCAATCGTTCAACGTCTTGCTGATAAGTGATGGCGATTATATGAGGTTAGCAATGCCTTCTGGAAATCAGGAACCCCGCCGCGATCCTGAATTAAAACGTAAAGCCTGGCTGGCGGTTTTTCTTGGTTCTGCACTTTTCTGGGTGGTTGTCGCACT  >  NC_000913/1335013‑1335263

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