Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	W3110S.gb	4583689	4583850	162	153 [0]	[0] 49	mcrB	5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrB

GACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACT  >  W3110S.gb/4583851‑4583918 |                                                                    gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTCAAATTACCATAACCGAAACTCAATTTTACt  <  1:387062/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:608143/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:169241/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:175001/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:280986/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:301647/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:308250/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:335662/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:340899/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:399748/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:529647/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:547033/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:581935/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:169235/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:627277/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:633299/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:730962/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:741314/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:757328/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:796808/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:851324/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:86224/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:906841/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:940610/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:954915/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:1370939/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:102883/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:1077584/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:113641/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:1138564/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:1202045/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:1203908/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:1208467/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:1209471/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:1223206/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:1238492/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:128398/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:1024613/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:1410404/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:1423241/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:1471589/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:1484366/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:1488793/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:1490800/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:1518592/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:1525759/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:1557343/68‑1 (MQ=255) gACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:1567208/68‑1 (MQ=255) gACCTTCTCCAAGAAATCCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACt  <  1:1580101/68‑1 (MQ=255) |                                                                    GACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACT  >  W3110S.gb/4583851‑4583918

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