Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	W3110S.gb	2888581	2888618	38	42 [0]	[0] 15	cysI	sulfite reductase, beta subunit, NAD(P)‑binding

CATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCT  >  W3110S.gb/2888516‑2888580                                                                 | caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:435866/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:1705332/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:1065367/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:1869610/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:208910/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:216436/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:275853/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:292838/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:35845/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:396316/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:1829945/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:551796/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:561578/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:637014/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:771720/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:7862/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:899456/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:90881/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:959567/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:978508/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:1694343/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:1187538/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:120729/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:1249513/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:1358091/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:1406017/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:1435988/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:1172302/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:1464474/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:1538748/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:1540951/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:1568146/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:1594023/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:1629188/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:16297/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:164379/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:1664408/65‑1 (MQ=255) caTCGCGTGGCGCGGCTCCAGCTTCTGTTCAACACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:562799/65‑1 (MQ=255) caTCGCGGGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:1353396/65‑1 (MQ=255)  aTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:1067619/64‑1 (MQ=255)  aTCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:486911/64‑1 (MQ=255)                              cAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct  <  1:493577/36‑1 (MQ=255)                                                                 | CATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCT  >  W3110S.gb/2888516‑2888580

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