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	31 [0]	[0] 31	[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)	31 (0.540)	29/490	0.8	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

TCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATTTTAATAACGCTATCTGCGGATAAAGCAGAATAGGTGGTTAACCCCAGACATAAACCGAGGAAAATAATGTTATTGTATTTCATAATCTATTGTTCCTTAGCGACAGATTGCTGTCTGCTGGTTCAGTAAGGTACCAGGAGAAACTTCAGGAAGCTTGTACTCGACAATACAGTTTGAGTTTTTATCTTTGCCCCATGAAACCTGTAATTGCCCTGACTGTGGAAGTCCAGTCAGATAAA  >  NC_000913/1588522‑1588995                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                 tccttcATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatatta                                                                                                                                                                                                                                  <  1:546035‑M1/250‑17 (MQ=255)    ttcATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattgg                                                                                                                                                                                                                              >  1:93708‑M1/1‑231 (MQ=255)          cAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagc                                                                                                                                                                                                                        >  1:281630‑M1/1‑225 (MQ=255)           aGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagct                                                                                                                                                                                                                       >  1:191797‑M1/1‑224 (MQ=255)              ttCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgt                                                                                                                                                                                                                     >  1:543983‑M1/1‑221 (MQ=255)                caacaaAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaat                                                                                                                                                                                                                  <  1:299721‑M1/251‑33 (MQ=255)                caacaaAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaat                                                                                                                                                                                                                  <  1:457712‑M1/251‑33 (MQ=255)                  acaaAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTGGTCCCTTTTCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattc                                                                                                                                                                                                                <  2:316162‑M1/251‑35 (MQ=255)                        tattatTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatcttattggaatagctgtaattcagcga                                                                                                                                                                                                           >  2:400353‑M1/1‑211 (MQ=255)                            attTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACACCCGTTGTCCCTTTGATGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATggtttgccatatattattggaatagctgtaattcagcgaagtcg                                                                                                                                                                                                      >  1:206408‑M1/1‑207 (MQ=255)                                    cATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCGTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatca                                                                                                                                                                                              <  1:579575‑M1/251‑53 (MQ=255)                                        cTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatcacca                                                                                                                                                                                           >  2:262788‑M1/1‑195 (MQ=255)                                          cAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatcaccaa                                                                                                                                                                                          <  2:191875‑M1/249‑57 (MQ=255)                                          cAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatcacca                                                                                                                                                                                           <  1:262788‑M1/248‑56 (MQ=255)                                                cGTTAAATGCAACCCGAACCCCCGTTGTCCATTTGCTGCATTCACTTAACGAAAGCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATAGGGGACGCGCCCGATAATgttttgccatatattattggaatagctgtaattcagcgaagtcgaaatatcaccaaccggaac                                                                                                                                                                                   >  1:21438‑M1/1‑187 (MQ=255)                                                    aaaTGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGCTCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagttgtaattcagcgaagtcgtaatatcaccaaacggcacattt                                                                                                                                                                               >  2:370537‑M1/1‑183 (MQ=255)                                                    aaaTGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatcaccaaacggcacattt                                                                                                                                                                               <  1:400353‑M1/250‑68 (MQ=255)                                                                  cccccGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatcaccaaacggcacatttaatccaaatgatatt                                                                                                                                                                >  2:407610‑M1/1‑169 (MQ=255)                                                                             ccTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatcaccaaacggcacatttaatccaaatgatatttgctggtcgct                                                                                                                                                     <  1:532087‑M1/251‑94 (MQ=255)                                                                                         tCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatcaccaaacggcacatttaatccaaatgatatttgctggtcgctgcgtgacgtgt                                                                                                                                          <  1:507400‑M1/250‑105 (MQ=255)                                                                                              taaagtaatCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCAATGTTTTTTTGGAGATCTACGGAAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatcaccaaacggcacatttaatccaaatgatatttgctggtcgctgcgtgacgtgttcca                                                                                                                                      <  2:281630‑M1/245‑109 (MQ=255)                                                                                                  gTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatcaccaaacggcacatttaatccaaatgatatttgctggtcgctgcgtgacgtgttccagtaact                                                                                                                                <  1:250643‑M1/251‑115 (MQ=255)                                                                                                         gAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatcaccaaacggcacatttaatccaaatgatatttgctggtcgctgcgtgacgtgttccagtaactttggc                                                                                                                           <  1:190620‑M1/249‑120 (MQ=255)                                                                                                            aaGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatcaccaaacggcacatttaatccaaatgatatttgctggtcgctgcgtgacgtgttccagtaactttggcgactg                                                                                                                      <  1:82381‑M1/251‑125 (MQ=255)                                                                                                                                                                  tttGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatcaccaaacggcacatttaatccaaatgatatttgctggtcgctgcgtgacgtgttccagtaactttggcgactggcactgaaaaatgtcgtaccgtaatttccaagctgctgagagatgcttatttgt                                                                >  1:326151‑M1/1‑73 (MQ=255)                                                                                                                                                                                   aaaTTAAGCGAATCCGATGAGACTGGGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatcaccaaacggcacatttaatccaaatgatatttgctggtcgctgcgtgacgtgttccagtaactttggcgactggcactgaaaaatgtcgtaccgtaatttccaagctgctgagagatgcttatttgttcctgaccacgcttact                                               <  2:191797‑M1/251‑196 (MQ=255)                                                                                                                                                                                    aaTTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatcaccaaacggcacatttaatccaaatgatatttgctggtcgctgcgtgacgtgttccagtaactttggcgactggcactgaaaaatgtcgtaccgtaatttccaagctgctgagagatgcttatttgttcctgaccacgcttactg                                              >  1:118306‑M1/1‑55 (MQ=255)                                                                                                                                                                                          gCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatcaccaaacggcacatttaatccaaatgatatttgctggtcgctgcgtgacgtgttccagtaactttggcgactggcactgaaaaatgtcgtaccgtaatttccaagctgctgagagatgcttctttgttcctgaccacgcttactgtagaa                                         >  2:313785‑M1/1‑49 (MQ=255)                                                                                                                                                                                                                         gAGGACGCGCCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatcaccaaacggcacatttaatccaaatgatatttgctggtcgctgcgtgacgtgttccagtaactttggcgactggcactgaaaaatgtcgtaccgtaatttccaagctgctgagagatgcttatttgttcctgaccacgcttactgtagaacagattaaaataatcaataaattgtgtctgct         >  1:510023‑M1/1‑18 (MQ=255)                                                                                                                                                                                                                                cgcCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatcaccaaacggcacatttaatccaaatgatatttgctggtcgctgcgtgacgtgttccagtaactttggcgactggcactgaaaaatgtcgtaccgtaatttccaagctgctgagagatgcttatttgttcctgaccacgcttactgtagaacagattaaaataatcaataaattgtgtctgctcattggt  <  2:198878‑M1/251‑241 (MQ=255)                                                                                                                                                                                                                                cgcCCGCTAATgttttgccatatattattggaatagctgtaattcagcgaagtcgtaatatcaccaaacggcacatttaatccaaatgatatttgctggtcgctgcgtgacgtgttccagtaactttggcgactggcactgaaaaatgtcgtaccgtaatttccaagctgctgagagatgcttatttgttcctgaccacgcttactgtagaacagattaaaataatcaataaattgtgtctgctcattggt  <  2:93708‑M1/251‑241 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                 TCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTGCTACCGGTCGTTGGAAATTGTCTGGCACTGTTTTTTTGGAGATCTACGGTAAAATTAAGCGAATCCGATGAGACTGTGCAGCCATAATCGAGGACGCGCCCGCTAATTTTAATAACGCTATCTGCGGATAAAGCAGAATAGGTGGTTAACCCCAGACATAAACCGAGGAAAATAATGTTATTGTATTTCATAATCTATTGTTCCTTAGCGACAGATTGCTGTCTGCTGGTTCAGTAAGGTACCAGGAGAAACTTCAGGAAGCTTGTACTCGACAATACAGTTTGAGTTTTTATCTTTGCCCCATGAAACCTGTAATTGCCCTGACTGTGGAAGTCCAGTCAGATAAA  >  NC_000913/1588522‑1588995

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