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

ACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATGCCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCAATGACTGAAGGGAACATTAAGCGTAAAAGCGAGTAAATGATCCCGATCGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCACGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGT  >  NC_000913/1589492‑1589965                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                tccttcatccagtttcaacaaagtattattttctgcatcctcaataccgttaaatgcaacccgaacccccgttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATAtatta                                                                                                                                                                                                                                  <  1:546035‑M2/16‑1 (MQ=255)    ttcatccagtttcaacaaagtattattttctgcatcctcaataccgttaaatgcaacccgaacccccgttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTgg                                                                                                                                                                                                                              >  1:93708‑M2/232‑251 (MQ=255)          cagtttcaacaaagtattattttctgcatcctcaataccgttaaatgcaacccgaacccccgttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGc                                                                                                                                                                                                                        >  1:281630‑M2/226‑251 (MQ=255)           agtttcaacaaagtattattttctgcatcctcaataccgttaaatgcaacccgaacccccgttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCt                                                                                                                                                                                                                       >  1:191797‑M2/225‑251 (MQ=255)              ttcaacaaagtattattttctgcatcctcaataccgttaaatgcaacccgaacccccgttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGt                                                                                                                                                                                                                     >  1:543983‑M2/222‑250 (MQ=255)                caacaaagtattattttctgcatcctcaataccgttaaatgcaacccgaacccccgttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAAt                                                                                                                                                                                                                  <  1:299721‑M2/32‑1 (MQ=255)                caacaaagtattattttctgcatcctcaataccgttaaatgcaacccgaacccccgttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAAt                                                                                                                                                                                                                  <  1:457712‑M2/32‑1 (MQ=255)                  acaaagtattattttctgcatcctcaataccgttaaatgcaacccgaacccccgtggtcccttttctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTc                                                                                                                                                                                                                <  2:316162‑M2/34‑1 (MQ=255)                        tattattttctgcatcctcaataccgttaaatgcaacccgaacccccgttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATCTTATTGGAATAGCTGTAATTCAGCGa                                                                                                                                                                                                           >  2:400353‑M2/212‑250 (MQ=255)                            attttctgcatcctcaataccgttaaatgcaacccgaacacccgttgtccctttgatgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGGTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCg                                                                                                                                                                                                      >  1:206408‑M2/208‑251 (MQ=255)                                    catcctcaataccgttaaatgcaacccgaacccccgttgtcccgttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCa                                                                                                                                                                                              <  1:579575‑M2/52‑1 (MQ=255)                                        ctcaataccgttaaatgcaacccgaacccccgttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCa                                                                                                                                                                                           >  2:262788‑M2/196‑250 (MQ=255)                                          caataccgttaaatgcaacccgaacccccgttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCaa                                                                                                                                                                                          <  2:191875‑M2/56‑1 (MQ=255)                                          caataccgttaaatgcaacccgaacccccgttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCa                                                                                                                                                                                           <  1:262788‑M2/55‑1 (MQ=255)                                                cgttaaatgcaacccgaacccccgttgtccatttgctgcattcacttaacgaaagctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataataggggacgcgcccgataatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGAAATATCACCAaccggaac                                                                                                                                                                                   >  1:21438‑M2/188‑247 (MQ=255)                                                    aaatgcaacccgaacccccgttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagctctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGTTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACAttt                                                                                                                                                                               >  2:370537‑M2/184‑250 (MQ=255)                                                    aaatgcaacccgaacccccgttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACAttt                                                                                                                                                                               <  1:400353‑M2/67‑1 (MQ=255)                                                                  cccccgttgtccctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATAtt                                                                                                                                                                >  2:407610‑M2/170‑251 (MQ=255)                                                                             cctttgctgcattcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCt                                                                                                                                                     <  1:532087‑M2/93‑1 (MQ=255)                                                                                         tcacttaacgtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACgtgt                                                                                                                                          <  1:507400‑M2/104‑1 (MQ=255)                                                                                              taaagtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcaatgtttttttggagatctacggaaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCa                                                                                                                                      <  2:281630‑M2/108‑1 (MQ=255)                                                                                                  gtaatctgaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACt                                                                                                                                <  1:250643‑M2/114‑1 (MQ=255)                                                                                                         gaaaagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGc                                                                                                                           <  1:190620‑M2/119‑1 (MQ=255)                                                                                                            aagggacggctggacttgtgctaccggtcgttggaaattgtctggcactgtttttttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTg                                                                                                                      <  1:82381‑M2/124‑1 (MQ=255)                                                                                                                                                                  tttggagatctacggtaaaattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGt                                                                >  1:326151‑M2/74‑251 (MQ=255)                                                                                                                                                                                   aaattaagcgaatccgatgagactgggcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCACGCTTACt                                               <  2:191797‑M2/195‑1 (MQ=255)                                                                                                                                                                                    aattaagcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCACGCTTACTg                                              >  1:118306‑M2/56‑251 (MQ=255)                                                                                                                                                                                          gcgaatccgatgagactgtgcagccataatcgaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTCTTTGTTCCTGACCACGCTTACTGTAGaa                                         >  2:313785‑M2/50‑250 (MQ=255)                                                                                                                                                                                                                         gaggacgcgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCACGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCt         >  1:510023‑M2/19‑251 (MQ=255)                                                                                                                                                                                                                                cgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCACGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGgt  <  2:198878‑M2/240‑1 (MQ=255)                                                                                                                                                                                                                                cgcccgctaatGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCACGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGgt  <  2:93708‑M2/240‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                ACTGGTGCCAGACGATGTATTACCTCCGTGGGTGTTACCGACCTGAACGCTATAATTCAGGTTATTATCCGGCAGCAGAGTGCCATAAACCCCCGATAGATTGGTCATGCCGCCTTTCAAATCGTTTGACATACTGTAACTGGCGTTTGAATTACGAAATGCCGACTGACTGTCTGTACGCATCCAATGACTGAAGGGAACATTAAGCGTAAAAGCGAGTAAATGATCCCGATCGTTTTGCCATATATTATTGGAATAGCTGTAATTCAGCGAAGTCGTAATATCACCAAACGGCACATTTAATCCAAATGATATTTGCTGGTCGCTGCGTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCACGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGT  >  NC_000913/1589492‑1589965

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