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	25 [0]	[0] 25	[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)	25 (0.440)	24/472	0.9	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

TAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAA  >  NC_000913/1586750‑1587233                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                            aggtaagtaatgattgaaccgtattctcagcaaagaaagaggcgtttacctgcttatccgtgacttcccaggcgattgcttttgcgccgttttcacaaatgtcgaatgacgcactgtgcttaggatctggccatgcacaacccggacagtcaaagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACt                                                                                                                                                                                                                                           <  2:154900‑M2/3‑1 (MQ=255)  ggtaagtaatgattgaaccgtattctcagcaaagaaagaggcgtttacctgcttatccgtgacttcccaggcgattgcttttgcgccgttttcacaaatgtcgaatgacgcactgtgcttaggatctggccatgcacaacccggacagtcaaagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTc                                                                                                                                                                                                                                          >  1:536159‑M2/248‑251 (MQ=255)       gtaatgattgaaccgtattctcagcaaagaaagaggcgtttacctgcttatccgtgacttcccaggcgattgcttttgcgccgttttcacaaatgtcgaatgacgcactgtgcttaggatctggccatgcacaacccggacagtcaaagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGAt                                                                                                                                                                                                                                     <  2:407618‑M2/9‑1 (MQ=255)       gtaatgattgaaccgtatcctcagcaaagaaagaggcgtttacctgcttatccgtgacttcccaggcgattgcttttgcgccgttttcacaaatgtcgaatgacgcactgtgcttaggatctggccatgcacaacccggacagtcaaagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGAt                                                                                                                                                                                                                                     <  2:220075‑M2/9‑1 (MQ=255)                    cgtattcgcagcaaagaaagaggcgtttacctgcttatccgtgactgcccaggcgattgcttttgcgccgttttcacaaatgtcgaatgacgcactgtgcttaggatctggccatgcacaacccggacagtcaaagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCAt                                                                                                                                                                                                                        <  2:365349‑M2/22‑1 (MQ=255)                       attctcagcaaagaaagaggcgtttacctgcttatccgtgacttcccaggcgattgcttttgcgccgttttcacaaatgtcgaatgacgcactgtgcttaggatctggccatgcacaacccggacagtcaaagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCa                                                                                                                                                                                                                      <  2:453271‑M2/24‑1 (MQ=255)                          ctcagcaaagaaagaggcgtttacctgcttatcagtgacttcccaggcgattgcttttgcgccgttttcacaaatgtcgaatgacgcactgtgcttaggatctggccatgcacaacccggacagtcaaagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTg                                                                                                                                                                                                                  >  2:400385‑M2/224‑251 (MQ=255)                             agcaaagaaagaggcgtttacctgcttatccgtgactgcccaggcgattgcttttgcgccgttttcacaaatgtcgaatgacgcactgtgcttaggatctggccatgcacaacccggacaggcaaagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTGTTGTTGCATCAGGTGCtg                                                                                                                                                                                                               <  2:575000‑M2/31‑1 (MQ=255)                                                           cgtgacttcccagtctattgcttgtgcgccgttttcacaaatggcgaatgtcgcactgtgcttaggatctggccatgcacaacccggacagtcaaagccctctggcttattcatggcaaacatggcaataatatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCGTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAg                                                                                                                                                                                 <  1:590586‑M2/61‑1 (MQ=255)                                                                     caggcgattgcttttgcgccgttttcacaaatgtcgaatgacgcactgtgcttaggatctggccatgcacaacccggacagtcaaagccctctggcttattcatgtcaaacatggcaagaacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGGTGCATCAGGTGCTGTTTTAGCAATCACTTGTCGGCTACTGTCAGTGGTGGCAcc                                                                                                                                                                       <  1:486241‑M2/71‑1 (MQ=255)                                                                              gcttttgcgccgttttcacaaatgtcgaatgacgcactgtgcttaggatctggccatgcacaacccggacagtcaaagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGAt                                                                                                                                                               >  1:170360‑M2/172‑250 (MQ=255)                                                                                     cgccgttttcacaaatgtcgaatgacgcactgtgcttaggatctggccatgcacaacccggacagtcaaagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATgag                                                                                                                                                       <  1:597265‑M2/87‑1 (MQ=255)                                                                                         gttttcacaaatgtcgaatgacgcactgtgcttaggatctggccatgcacaacccggacagtcaaagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCt                                                                                                                                                   <  1:251854‑M2/91‑1 (MQ=255)                                                                                                     gtcgaaagacgcactgtgcttaggattgggccatgcacaacccggacagtcaaagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatcctgcggcttacgtgctgcatgcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACATGTGTCGGCTACTGTCAGTGGTTGCACCCGAAAGATAAAATGAGAGCTTTTGGTCGCtt                                                                                                                                        <  2:519690‑M2/102‑1 (MQ=255)                                                                                                        gaatgacgcactgtgcttaggatctggccatgcacaacccggacagtcaaagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGCAAGATAAAATGAGAGCTTTTGTTCGCTTGAGc                                                                                                                                    >  2:22031‑M2/146‑251 (MQ=255)                                                                                                                    gtgcttaggatctggccatgcacaacccggacagtcaaagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGTTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCa                                                                                                                         <  2:508330‑M2/117‑1 (MQ=255)                                                                                                                                 ggccatgcacaacccggacagtcaaagccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTgc                                                                                                           <  1:424808‑M2/131‑1 (MQ=255)                                                                                                                                                          agccctctggcttattcatgtcaaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATcc                                                                                  >  1:338087‑M2/96‑251 (MQ=255)                                                                                                                                                                               caaacatggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCa                                                             <  1:376428‑M2/177‑1 (MQ=255)                                                                                                                                                                                      ggcaataacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGc                                                      <  1:400385‑M2/184‑1 (MQ=255)                                                                                                                                                                                           taacatcctggcgtatatccatctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAgccg                                                 >  1:261135‑M2/63‑251 (MQ=255)                                                                                                                                                                                                              catctgcttacgtactgcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGGTGGCATAACAACAt                              <  2:536159‑M2/208‑1 (MQ=255)                                                                                                                                                                                                                              gcattcgctacggatttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAaata                <  1:22031‑M2/222‑1 (MQ=255)                                                                                                                                                                                                                               cattcgctccggatttaacagcccccACTCCCGATGCTTTTGTTGCATCAGGTGCTGTTTTAGCAAACACTTGTCGGCTACTGTCAGTGGGGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTTAGCAATATACGCCAAGGTGTATGTCTGTGCTTCCAGGAAAGTCCTGGAGAGCCATGGTGACATGTCGACTATCAACGGTGCAGTCGCCTGTTGGCATAACAACATTGTTATGTGAAataata             <  2:350517‑M2/225‑1 (MQ=255)                                                                                                                                                                                                                                           atttaacagcacccACTCCCGATGCTTTTGTTGCATCAGGAGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATATTATTCCACGTaaa  >  1:322922‑M2/15‑250 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                            TAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAA  >  NC_000913/1586750‑1587233

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