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

AGGTAAGTAATGATTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGCACAGTACATACAAAGGGATTATCGCCAACATGGATATATAAGAGTAGGGGACTACTCTCTTTCTATTTTTAATCTATTATCGTTAAACAATAGATTAAATACGATAAGTTAAGCATCATCTCTCATAGTGTTTAAATTAA  >  NC_000913/1586195‑1586678                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             aGGTAAGTAATGATTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCact                                                                                                                                                                                                                                           <  2:154900‑M1/251‑4 (MQ=255)  ggTAAGTAATGATTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactc                                                                                                                                                                                                                                          >  1:536159‑M1/1‑247 (MQ=255)       gtaaTGATTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgat                                                                                                                                                                                                                                     <  2:407618‑M1/251‑10 (MQ=255)       gtaaTGATTGAACCGTATCCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgat                                                                                                                                                                                                                                     <  2:220075‑M1/251‑10 (MQ=255)                    cGTATTCGCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTGCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcat                                                                                                                                                                                                                        <  2:365349‑M1/251‑23 (MQ=255)                       aTTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatca                                                                                                                                                                                                                      <  2:453271‑M1/250‑25 (MQ=255)                          ctcAGCAAAGAAAGAGGCGTTTACCTGCTTATCAGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtg                                                                                                                                                                                                                  >  2:400385‑M1/1‑223 (MQ=255)                             aGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTGCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGGCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgctgttgttgcatcaggtgctg                                                                                                                                                                                                               <  2:575000‑M1/251‑32 (MQ=255)                                                           cGTGACTTCCCAGTCTATTGCTTGTGCGCCGTTTTCACAAATGGCGAATGTCGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGGCAAACATGGCAATAATATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcgtttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcag                                                                                                                                                                                 <  1:590586‑M1/251‑62 (MQ=255)                                                                     cAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAAGAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttggtgcatcaggtgctgttttagcaatcacttgtcggctactgtcagtggtggcacc                                                                                                                                                                       <  1:486241‑M1/251‑72 (MQ=255)                                                                              gCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagat                                                                                                                                                               >  1:170360‑M1/1‑171 (MQ=255)                                                                                     cgcCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgag                                                                                                                                                       <  1:597265‑M1/251‑88 (MQ=255)                                                                                         gTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagct                                                                                                                                                   <  1:251854‑M1/251‑92 (MQ=255)                                                                                                     gTCGAAAGACGCACTGTGCTTAGGATTGGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCTGCGGCTTACGTGCTGCATGCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacatgtgtcggctactgtcagtggttgcacccgaaagataaaatgagagcttttggtcgctt                                                                                                                                        <  2:519690‑M1/250‑103 (MQ=255)                                                                                                        gAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgcaagataaaatgagagcttttgttcgcttgagc                                                                                                                                    >  2:22031‑M1/1‑145 (MQ=255)                                                                                                                    gtgCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtgttggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgcca                                                                                                                         <  2:508330‑M1/250‑118 (MQ=255)                                                                                                                                 ggCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgccaaggggtatttctgc                                                                                                           <  1:424808‑M1/251‑132 (MQ=255)                                                                                                                                                          aGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgccaaggggtatttctgcgcttccaggaaagtccgggagatcc                                                                                  >  1:338087‑M1/1‑95 (MQ=255)                                                                                                                                                                               cAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgccaaggggtatttctgcgcttccaggaaagtccgggagatccacggtgacatttcgactatca                                                             <  1:376428‑M1/251‑178 (MQ=255)                                                                                                                                                                                      ggCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgccaaggggtatttctgcgcttccaggaaagtccgggagatccacggtgacatttcgactatcaacggtgc                                                      <  1:400385‑M1/251‑185 (MQ=255)                                                                                                                                                                                           tAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgccaaggggtatttctgcgcttccaggaaagtccgggagatccacggtgacatttcgactatcaacggtgcagccg                                                 >  1:261135‑M1/1‑62 (MQ=255)                                                                                                                                                                                                              cATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgccaaggggtatttctgcgcttccaggaaagtccgggagatccacggtgacatttcgactatcaacggtgcagccgcctggtggcataacaacat                              <  2:536159‑M1/251‑209 (MQ=255)                                                                                                                                                                                                                              gCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgccaaggggtatttctgcgcttccaggaaagtccgggagatccacggtgacatttcgactatcaacggtgcagccgcctgttggcataacaacattgttatttgaaata                <  1:22031‑M1/249‑223 (MQ=255)                                                                                                                                                                                                                               cATTCGCTCCGGATTTAACAGCCCCCactcccgatgcttttgttgcatcaggtgctgttttagcaaacacttgtcggctactgtcagtgggggcacccgaaagataaaatgagagcttttgttcgctttagcaatatacgccaaggtgtatgtctgtgcttccaggaaagtcctggagagccatggtgacatgtcgactatcaacggtgcagtcgcctgttggcataacaacattgttatgtgaaataata             <  2:350517‑M1/251‑226 (MQ=255)                                                                                                                                                                                                                                           aTTTAACAGCACCCactcccgatgcttttgttgcatcaggagctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgccaaggggtatttctgcgcttccaggaaagtccgggagatccacggtgacatttcgactatcaacggtgcagccgcctgttggcataacaacattgttatttgaaatattattccacgtaaa  >  1:322922‑M1/1‑14 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             AGGTAAGTAATGATTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGCACAGTACATACAAAGGGATTATCGCCAACATGGATATATAAGAGTAGGGGACTACTCTCTTTCTATTTTTAATCTATTATCGTTAAACAATAGATTAAATACGATAAGTTAAGCATCATCTCTCATAGTGTTTAAATTAA  >  NC_000913/1586195‑1586678

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