Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC	REL606	1,609,175	Δ7,554 bp	between IS3	ybcW–[ECB_01526]	ybcW, gnsB, ynfN, ECB_01516, cspI, ydfP, ydfQ, ydfR, essQ, ECB_01522, ECB_01523, insE‑4, insF‑4, [ECB_01526]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	REL606	1607956–1609176	1616669–1615469	6294–8714	22 [21]	[20] 22	insE‑3–[insF‑4]	insE‑3,insF‑3,ybcW,gnsB,ynfN,ECB_01516,cspI,ydfP,ydfQ,ydfR,essQ,ECB_01522,ECB_01523,insE‑4,[insF‑4]

GAACTCAAACTGTACCGTCGGATCATAAACCGCAAGATGCGGCGTGGCGGTTATCTGAAAATTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACC  >  REL606/1607855‑1608013                                                                                                     |                                                           gAACTCAAACTGTACCGTCGGATCATAAACCGCAAGATGCGGCGTGGCGGTTATCTGAAAATTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAg                                                            >  2:893143/1‑101 (MQ=255)  aaCTCAAACTGTACCGTCGGATCATAAACCGCAAGATGCGGCGTGGCGGTTATCTGAAAATTGATCTTACCCAGCAATAGCGGACACGCGGCTAAGTGAGt                                                           >  1:1158364/1‑101 (MQ=255)   aCTCAAACTGTACCGTCGGATCATAAACCGCAAGATGCGGCGTGGCGGTTATCTGAAAATTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTa                                                          <  1:894162/101‑1 (MQ=255)      cAAACTGTACCGTCGGATCATAAACCGCAAGATGCGGCGTGGCGGTTATCTGAAAATTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAAc                                                       <  1:315349/101‑1 (MQ=255)      cAAACTGTACCGTCGGATCATAAACCGCAAGATGCGGCGTGGCGGTTATCTGAAAATTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAAc                                                       <  2:533358/101‑1 (MQ=255)                 gTCGGATCATAAACCGCAAGATGCGGCGTGGCGGTTATCTGAAAATTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCaga                                            >  1:763398/1‑101 (MQ=255)                  tCGGATCATAATCCGCAAGATGCGGCGTGGCGGTTATCTGAAAATTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCagag                                           <  2:751498/101‑1 (MQ=255)                     gATCATAAACCGCAAGATGCGGCGTGGCGGTTATCTGAAAATTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTg                                        >  2:1630176/1‑101 (MQ=255)                     gATCATAAACCGCAAGATGCGGCGTGGCGGTTATCTGAAAATTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTg                                        <  2:493092/101‑1 (MQ=255)                          tAAACCGCAAGATGCGGCGTGGCGGTTATCTGAAAATTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTca                                   <  2:1588453/101‑1 (MQ=255)                          tAAACCGCAAGATGCGGCGTGGCGGTTATCTGAAAATTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTca                                   <  1:1474692/101‑1 (MQ=255)                                     aTGCGGCGTGGCGGTTATCTGAAAATTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaaaa                        <  2:888639/101‑1 (MQ=255)                                      tGCGGCGTGGCGGTTATCTGAAAATTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAAc                       <  2:1146559/101‑1 (MQ=255)                                         ggcgTGGCGGTTATCTGAAAATTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGt                    <  1:1718760/101‑1 (MQ=255)                                               gCGGTTATCTGAAAATTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAAc              <  2:1417452/101‑1 (MQ=255)                                                 ggTTATCTGAAAATTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGGGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCa            >  1:1676828/1‑101 (MQ=255)                                                  gTTATCTGAAAATTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAg           <  1:893148/101‑1 (MQ=40)                                                   ttATCTGAAAATTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGt          <  1:142119/101‑1 (MQ=40)                                                    tATCTGAAAATTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTa         >  2:187718/1‑101 (MQ=38)                                                     aTCTGAAAATTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTaa        <  1:120451/101‑1 (MQ=38)                                                          aaaaTTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAAc   <  1:109547/101‑1 (MQ=32)                                                           aaaTTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAAcc  <  1:893874/101‑1 (MQ=31)                                                                                                     |                                                           GAACTCAAACTGTACCGTCGGATCATAAACCGCAAGATGCGGCGTGGCGGTTATCTGAAAATTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACC  >  REL606/1607855‑1608013

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑