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

AAGTAATGATTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGCACAGTACATACAAAGGGATTATCGCCAACATGGATATATAAGAGTAGGGGACTACTCTCTTTCTATTTTTAATCTATTATCGTTAAACAATAGATTAAATACGATAAGTTAAGCATCATCTCTCATAGTGTTTA  >  NC_000913/1586199‑1586672                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                       aagtaaTGATTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccg                                                                                                                                                                                                                                 <  2:348396‑M1/251‑8 (MQ=255)   gtaaTGATTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgat                                                                                                                                                                                                                               >  2:392323‑M1/1‑242 (MQ=255)                cGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcat                                                                                                                                                                                                                  >  1:122369‑M1/1‑229 (MQ=255)                          gCAAAGAAAGAGGCGTTTTCCTGCGTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCtctcccgatgcttttgttgcatcaggtgctgt                                                                                                                                                                                                        <  2:293488‑M1/251‑33 (MQ=255)                                   gagGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaa                                                                                                                                                                                               <  1:221928‑M1/251‑42 (MQ=255)                                        gTTTACCTGCTTATCCGTGACTTCCCAGGCGAGTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACTCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTTCTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacactt                                                                                                                                                                                          <  1:121943‑M1/251‑47 (MQ=255)                                             ccTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcgg                                                                                                                                                                                     >  1:314259‑M1/1‑200 (MQ=255)                                                  ttATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctact                                                                                                                                                                                >  2:146827‑M1/1‑195 (MQ=255)                                                            cTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGAAAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtg                                                                                                                                                                      >  2:278912‑M1/1‑185 (MQ=255)                                                                 cAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacc                                                                                                                                                                 <  2:188255‑M1/251‑72 (MQ=255)                                                                           cTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgtt                                                                                                                                                                                                     <  2:370035‑M1/205‑36 (MQ=255)                                                                           cTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgtt                                                                                                                                                                                                     >  1:370035‑M1/1‑170 (MQ=255)                                                                                                    gAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCACAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagc                                                                                                                              >  1:158595‑M1/1‑145 (MQ=255)                                                                                                             aCTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGAACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggagctgtgttagcaaacacttgtcggctactgtcagtggtgtcacccgaaagataaaatgagagctattgttcgcttgagcaatatacgc                                                                                                                     >  1:294183‑M1/1‑136 (MQ=255)                                                                                                                             ggCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcggttgttcgcttgagcaatatacgccaaggggtatttctg                                                                                                      <  1:424019‑M1/250‑131 (MQ=255)                                                                                                                               ccATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgccaaggggtatttctgcgc                                                                                                   >  2:347222‑M1/1‑118 (MQ=255)                                                                                                                                          ccGGAAAGTCAAAGCCCTCTGGCTTATTAATGTCAAACATGGCAATAACATCCTGGCGTCTATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgccaaggggtatttctgcgcttccaggaa                                                                                          >  2:400757‑M1/1‑107 (MQ=255)                                                                                                                                                                  tATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagag                                                                                                                                               <  2:451503‑M1/172‑90 (MQ=255)                                                                                                                                                                  tATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagag                                                                                                                                               >  1:451503‑M1/1‑83 (MQ=255)                                                                                                                                                                                      aTAACATCCTGGCGTATATCCACCTGCTTACGTACTGCAGTCGCTACGGATTTAACAGCACCCactccctatgcttttgttgcatcaggtgctgtgttagcaaacactggtcgggtactgtcaggggtggctcccgaaagataaaatgggagctttggtgcgcttgagcaatatacgccaaggggtatgtctgcgcttccaggaaagtccgggagatccacggtgacagttcgaccatcaactctgcagcc                                            <  2:40657‑M1/251‑189 (MQ=255)                                                                                                                                                                                                ggCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgccaaggggtatttctgcgcttccaggaaagtccgggagatccacggtgacatttcgactatcaacggtgcagccgcctgttg                                    <  1:53460‑M1/249‑197 (MQ=255)                                                                                                                                                                                                  cGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgccaaggggtatttctgcgcttccaggaaagtccgggagatccacggtgacatttcgactatcaacggtgcagccgcctgtgggcat                                <  1:392323‑M1/251‑201 (MQ=255)                                                                                                                                                                                                   gTATATTCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgccaaggggtatttctgcgcttccaggaaagtccgggagatccacggtgacatttcgactatcaacggtgcagccgcctgttggcat                                >  1:136135‑M1/1‑50 (MQ=255)                                                                                                                                                                                                          cATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcgtttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgccaaggggtatttctgcgcttccagtaaagtccgggagatccacggtgacatttcgactatcaacggtgcagccgcctgttggcataacaacat                        <  1:374006‑M1/251‑209 (MQ=255)                                                                                                                                                                                                               tctTACGTACTGCATTCGCTACGAATTTAACAGCACCCacacccgatgcttttgtttcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagatcttttgttcgcttgagcaatttacgccaaggtgtatttctgcgcttccatgaaagtccgggagatcctctgatacatttcgtctttcaacggtgcagccgcctgttggcatatcaacatt                       <  1:278912‑M1/246‑210 (MQ=255)                                                                                                                                                                                                               gCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgaatcaggtgctgtgcaagcaaacacttgtcagctacagtcagtggtggcacccgaaagctaaaatgagcgattttgttcgcttgagcaatatacgccaaggggaatttctgcgcttccaggaaagtccgggaaatcaacggtgacatttcgactctcaacggtgcagccgcctgttggcataacaacattgtt                    >  1:27874‑M1/1‑38 (MQ=255)                                                                                                                                                                                                                           cATTCGCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgtcagtggtggcacccgaaagataaaatgagagcttttgttcgcttgagcaatatacgccaaggggtatttctgcgcttccaggaaagtccgggagatccacggtgacatttcgactatcaacggtgcagccgcctgttggcataacaacattgttatttgaaataata       <  2:240389‑M1/251‑226 (MQ=255)                                                                                                                                                                                                                                gCTACGGATTTAACAGCACCCactcccgatgcttttgttgcatcaggtgctgtgttagcaaacacttgtcggctactgacagtggtggcacccgaaagataaaatgagagcttttgttcgctttagcaatatgcgccaaggggtatttctgcgcttccaggaaagtccgggagatccacggtgacatttcgactatcaacggtgctgccgcctgttggcataacaacattgttatttgaaataatattcca  <  1:76588‑M1/251‑231 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                       AAGTAATGATTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTACGTACTGCATTCGCTACGGATTTAACAGCACCCCAACCACCTGCAGCACCCTGGTAGGATTCAATTTTTTTCTTCATCTACTTATCCTGTGTGGAAATGAATTATTAGGATTTGTAATAGCGATGATGCACAGTACATACAAAGGGATTATCGCCAACATGGATATATAAGAGTAGGGGACTACTCTCTTTCTATTTTTAATCTATTATCGTTAAACAATAGATTAAATACGATAAGTTAAGCATCATCTCTCATAGTGTTTA  >  NC_000913/1586199‑1586672

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