Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,588,756	Δ970 bp		[ydeS]–[ydeT]	[ydeS], [ydeT]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1588756	1589725	970	18 [0]	[0] 18	[ydeS]–[ydeT]	[ydeS],[ydeT]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1588755	0 (0.000)	18 (0.520)	16/490	1.0	100%	coding (85/531 nt)	ydeS	putative fimbrial protein YdeS
	?	NC_000913	1589726 =	0 (0.000)					coding (276/1149 nt)	ydeT	fimbrial usher domain‑containing protein YdeT

ATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATTTTAATAACGCTATCTGCGGATAAAGCAGAATAGGTGGTTAACCCCAGACATAAACCGAGGAAAATAATGTTATTGTATTTCATAATCTATTGTTCCTTAGCGACAGATTGCTGTCTGCTGGTTCAGTAAGGTACCAGGAGAAACTTCAGGAAGCTTGTACTCGACAATACAGTTTGAGTTTTTATCTTTGCCCCATGAAACCTGTAATTGCCCTGACTGTGGAAGTCCAGTCA  >  NC_000913/1588516‑1588989                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                           aTTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccat                                                                                                                                                                                                                                  >  2:6421‑M1/1‑240 (MQ=255)                      caacaaAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaat                                                                                                                                                                                                            >  2:356583‑M1/1‑219 (MQ=255)                       aacaaAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaatt                                                                                                                                                                                                           <  1:149754‑M1/251‑34 (MQ=255)                              tattatTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGAAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcga                                                                                                                                                                                                     >  1:142689‑M1/1‑211 (MQ=255)                                             ccTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatcacc                                                                                                                                                                                      >  2:311668‑M1/1‑196 (MQ=255)                                                                                                      aCGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatcaccaaacggcacatttaatccaaatgatatttgctggtcgctgcgtgacgtgttccagta                                                                                                                             >  1:83801‑M1/1‑139 (MQ=255)                                                                                                           aTCTGAAAAGGGCCGGCTGGACTTGTGCTACCGGTCGTTGGACATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATACGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatcaccaaacggcacatttaatccaaatgatatttgctggtcgctgcgtgacgtgttccagtaactttg                                                                                                                       >  2:223581‑M1/1‑134 (MQ=255)                                                                                                                    gggACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatcaccaaacggcacatttaatccaaatgatatttgctggtcgctgcgtgacgtgttccagtaactttggcgactggc                                                                                                              <  1:311668‑M1/251‑127 (MQ=255)                                                                                                                                   gtgCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatcaccaaacggcacatttaatccaaatgatatttgctggtcgctgcgtgacgtgttccagtaactttggcgactggcactgaaaaatgtcgt                                                                                               <  2:38621‑M1/251‑142 (MQ=255)                                                                                                                                      cTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatcaccaaacggcacatttaatccaaatgatatttgctggtcgctgcgtgacgtgttccagtaactttggcgactggcactgaaaaatgtcgtacc                                                                                            >  2:142971‑M1/1‑107 (MQ=255)                                                                                                                                      cTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatcaccaaacggcacatttaatccaaatgatatttgctggtcgctgcgtgacgtgttccagcaactttggcgactggcactgaaaaatgtcgtacc                                                                                            >  2:142983‑M1/1‑107 (MQ=255)                                                                                                                                                 tGGAAATTGCATGGCACTTTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatatgattggaatagctgtgattcagcgaagtcgtaatatcaccaaacggcacatttaatccaaatgatatttgctggtcgctgcgtgacgtgttccagtacctttggcgactggcactgagaaatgtcgtacagtaagttcca                                                                                  <  2:248192‑M1/250‑155 (MQ=255)                                                                                                                                                             ggCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatcaccaaacggcacatttaatccaaattatatttgctggtcgctgcgtgacgtgttccagtaactttggcgactggcactgaaaaatgtcgtaccgtaatttccaagctgctgag                                                                        <  1:332248‑M1/248‑165 (MQ=255)                                                                                                                                                                                           aTTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagacgtaatatcaccaaacggcacatttaatccaaatgatatttgctggtcgctgcgtgacgtgttccagtaactttggcgactggcactgaaaaatgtcgtaccgtaatttccaagctgctgagagatgcttatttgttcctgaccacgcttactgt                                       >  1:43723‑M1/1‑54 (MQ=255)                                                                                                                                                                                             tAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatcaccaaacggcacatttaatccaaatgatatttgctggtcgctgcgtgacgtgttccagtaactttggcgactggcactgaaaaatgtcgtaccgtaatttccaagctgctgagagatgcttatttgttcctgaccacgcttactgta                                      <  1:66935‑M1/250‑199 (MQ=255)                                                                                                                                                                                              aaGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatcaccaaacggcacatttaatccaaatgatatttgctggtcgctgcgtgacgtgttccagtaactttggcgactggcactgaaaaatgtcgtaccgtaatttccaagctgctgagagatgcttatttgttcctgaccccgcttactgtaga                                    <  2:310631‑M1/251‑201 (MQ=255)                                                                                                                                                                                              aaGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatcaccaaacggcacatttaatccaaatgatatttgctggtcgctgcgtgacgtgttccagtaactttggcgactggcactgaaaaatgtcgtaccgtaatttccaagctgctgagagatgcttatttgttcctgaccacgcttactgtaga                                    <  1:210288‑M1/251‑201 (MQ=255)                                                                                                                                                                                                                                aGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatcaccaaacggcacatttaatccaaatgatatttgctggtcgctgcgtgacgtgttccagtaactttggcgactggcactgaaaaatgtcgtaccgtaatttccaagctgctgagagatgcttatttgttcctgaccacgcttactgtagaacagattaaaataatcaataaattgtgtctgctc  <  1:322807‑M1/251‑235 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                           ATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATTTTAATAACGCTATCTGCGGATAAAGCAGAATAGGTGGTTAACCCCAGACATAAACCGAGGAAAATAATGTTATTGTATTTCATAATCTATTGTTCCTTAGCGACAGATTGCTGTCTGCTGGTTCAGTAAGGTACCAGGAGAAACTTCAGGAAGCTTGTACTCGACAATACAGTTTGAGTTTTTATCTTTGCCCCATGAAACCTGTAATTGCCCTGACTGTGGAAGTCCAGTCA  >  NC_000913/1588516‑1588989

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