Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,586,443	Δ555 bp		[ydeP]–[ydeQ]	[ydeP], [ydeQ]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1586443	1586997	555	28 [0]	[0] 28	[ydeP]–[ydeQ]	[ydeP],[ydeQ]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1586442	0 (0.000)	28 (0.620)	28/466	0.5	100%	coding (45/2280 nt)	ydeP	putative oxidoreductase YdeP
	?	NC_000913	1586998 =	0 (0.000)					coding (737/915 nt)	ydeQ	putative fimbrial adhesin protein YdeQ

CTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCA  >  NC_000913/1586754‑1587227                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                      aagtaatgattgaaccgtattctcagcaaagaaagaggcgtttacctgcttatccgtgacttcccaggcgattgcttttgcgccgttttcacaaatgtcgaatgacgcactgtgcttaggatctggccatgcacaacccggacagtcaaagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCg                                                                                                                                                                                                                                 <  2:348396‑M2/7‑1 (MQ=255)   gtaatgattgaaccgtattctcagcaaagaaagaggcgtttacctgcttatccgtgacttcccaggcgattgcttttgcgccgttttcacaaatgtcgaatgacgcactgtgcttaggatctggccatgcacaacccggacagtcaaagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGAt                                                                                                                                                                                                                               >  2:392323‑M2/243‑251 (MQ=255)                cgtattctcagcaaagaaagaggcgtttacctgcttatccgtgacttcccaggcgattgcttttgcgccgttttcacaaatgtcgaatgacgcactgtgcttaggatctggccatgcacaacccggacagtcaaagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCAt                                                                                                                                                                                                                  >  1:122369‑M2/230‑251 (MQ=255)                          gcaaagaaagaggcgttttcctgcgtatccgtgacttcccaggcgattgcttttgcgccgttttcacaaatgtcgaatgacgcactgtgcttaggatctggccatgcacaacccggacagtcaaagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccTCTCCCGATGCTTTTGTTGCATCAGGTGCtgt                                                                                                                                                                                                        <  2:293488‑M2/32‑1 (MQ=255)                                   gaggcgtttacctgcttatccgtgacttcccaggcgattgcttttgcgccgttttcacaaatgtcgaatgacgcactgtgcttaggatctggccatgcacaacccggacagtcaaagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCaaa                                                                                                                                                                                               <  1:221928‑M2/41‑1 (MQ=255)                                        gtttacctgcttatccgtgacttcccaggcgagtgcttttgcgccgttttcacaaatgtcgaatgactcactgtgcttaggatctggccatgcacaacccggacagtcaaagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgttctgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACtt                                                                                                                                                                                          <  1:121943‑M2/46‑1 (MQ=255)                                             cctgcttatccgtgacttcccaggcgattgcttttgcgccgttttcacaaatgtcgaatgacgcactgtgcttaggatctggccatgcacaacccggacagtcaaagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCgg                                                                                                                                                                                     >  1:314259‑M2/201‑251 (MQ=255)                                                  ttatccgtgacttcccaggcgattgcttttgcgccgttttcacaaatgtcgaatgacgcactgtgcttaggatctggccatgcacaacccggacagtcaaagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACt                                                                                                                                                                                >  2:146827‑M2/196‑251 (MQ=255)                                                            cttcccaggcgattgcttttgcgccgttttcacaaatgtcgaatgacgcactgtgcttaggatctggccatgcacaacccggaaagtcaaagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAgtggtg                                                                                                                                                                      >  2:278912‑M2/186‑251 (MQ=255)                                                                 caggcgattgcttttgcgccgttttcacaaatgtcgaatgacgcactgtgcttaggatctggccatgcacaacccggacagtcaaagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCAcc                                                                                                                                                                 <  2:188255‑M2/71‑1 (MQ=255)                                                                           cttttgcgccgttttcacaaatgtcgaatgacgcactgtgcttaggatctggccatgcacaacccggacagtcaaagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGtt                                                                                                                                                                                                     <  2:370035‑M2/35‑1 (MQ=255)                                                                           cttttgcgccgttttcacaaatgtcgaatgacgcactgtgcttaggatctggccatgcacaacccggacagtcaaagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGtt                                                                                                                                                                                                     >  1:370035‑M2/171‑205 (MQ=255)                                                                                                    gaatgacgcactgtgcttaggatctggccatgcacaacccggacagtcacagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGc                                                                                                                              >  1:158595‑M2/146‑251 (MQ=255)                                                                                                             actgtgcttaggatctggccatgcacaacccggacagtcaaagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgaactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGAGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGTCACCCGAAAGATAAAATGAGAGCTATTGTTCGCTTGAGCAATATACGc                                                                                                                     >  1:294183‑M2/137‑251 (MQ=255)                                                                                                                             ggccatgcacaacccggacagtcaaagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCGGTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTg                                                                                                      <  1:424019‑M2/130‑1 (MQ=255)                                                                                                                               ccatgcacaacccggacagtcaaagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTgcgc                                                                                                   >  2:347222‑M2/119‑251 (MQ=255)                                                                                                                                          ccggaaagtcaaagccctctggcttattaatgtcaaacatggcaataacatcctggcgtctatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGaa                                                                                          >  2:400757‑M2/108‑249 (MQ=255)                                                                                                                                                                  tattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGagag                                                                                                                                               <  2:451503‑M2/89‑1 (MQ=255)                                                                                                                                                                  tattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGagag                                                                                                                                               >  1:451503‑M2/84‑172 (MQ=255)                                                                                                                                                                                      ataacatcctggcgtatatccacctgcttacgtactgcagtcgctacggatttaacagcacccACTCCCTATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTGGTCGGGTACTGTCAGGGGTGGCTCCCGAAAGATAAAATGGGAGCTTTGGTGCGCTTGAGCAATATACGCCAAGGGGTATGTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACAGTTCGACCATCAACTCTGCAgcc                                            <  2:40657‑M2/188‑1 (MQ=255)                                                                                                                                                                                                ggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTg                                    <  1:53460‑M2/196‑1 (MQ=255)                                                                                                                                                                                                  cgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTGGGCAt                                <  1:392323‑M2/200‑1 (MQ=255)                                                                                                                                                                                                   gtatattcatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCAt                                >  1:136135‑M2/51‑250 (MQ=255)                                                                                                                                                                                                          catctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCGTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGTAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACAt                        <  1:374006‑M2/208‑1 (MQ=255)                                                                                                                                                                                                               tcttacgtactgcattcgctacgaatttaacagcacccACACCCGATGCTTTTGTTTCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGATCTTTTGTTCGCTTGAGCAATTTACGCCAAGGTGTATTTCTGCGCTTCCATGAAAGTCCGGGAGATCCTCTGATACATTTCGTCTTTCAACGGTGCAGCCGCCTGTTGGCATATCAACAtt                       <  1:278912‑M2/209‑1 (MQ=255)                                                                                                                                                                                                               gcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGAATCAGGTGCTGTGCAAGCAAACACTTGTCAGCTACAGTCAGTGGTGGCACCCGAAAGCTAAAATGAGCGATTTTGTTCGCTTGAGCAATATACGCCAAGGGGAATTTCTGCGCTTCCAGGAAAGTCCGGGAAATCAACGGTGACATTTCGACTCTCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGtt                    >  1:27874‑M2/39‑250 (MQ=255)                                                                                                                                                                                                                           cattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAataata       <  2:240389‑M2/225‑1 (MQ=255)                                                                                                                                                                                                                                gctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGACAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTTAGCAATATGCGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCTGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCa  <  1:76588‑M2/230‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                      CTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCA  >  NC_000913/1586754‑1587227

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: ‑