Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	REL606	475,285	Δ2,907 bp		[ybaL]–[ushA]	[ybaL], fsr, [ushA]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	REL606	475285	478191	2907	31 [0]	[0] 32	[ybaL]–[ushA]	[ybaL],fsr,[ushA]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 475284	0 (0.000)	27 (0.800)	23/138	0.5	100%	coding (22/1677 nt)	ybaL	predicted transporter with NAD(P)‑binding Rossmann‑fold domain
	?	REL606	478192 =	0 (0.000)					coding (1212/1653 nt)	ushA	UDP‑sugar hydrolase

TTGTACAGACCAATATGGGGCGGTTGATTCTGGCAGCCCAAATGGATCGCACTGGTGCCGACTTTGCGGTGATGAGCGGAGGCGGAATTCGTGATTCTATCGAAGCAGGCGATATCAGCTATAAAAACGTGCTGAAAGTGCAGCCATTC  >  REL606/478119‑478267                                                                          |                                                                            gagaaatacgtagtttattggccagcatgccgaggataaaggcgagcacaaggccgccaacaatggtggtgatGAg                                                                           <  1:635967‑M2/3‑1 (MQ=255)  agaaatacgtagtttattggccagcatgccgaggataaaggcgagcacaaggccgccaacaatggtggtgatGAGc                                                                          >  2:1079074‑M2/73‑76 (MQ=255)    aaatacgtagtttattggccagcatgccgaggataaaggcgagcacaaggccgccaacaatggtggtgatGAGCgg                                                                        <  2:464999‑M2/6‑1 (MQ=255)          gtagtttattggccagcatgccgaggataaaggcgagcacaaggccgccaacaatggtggtgatGAGCGGAGGCgg                                                                  <  2:194118‑M2/12‑1 (MQ=255)          gtagtttattggccagcatgccgaggataaaggcgagcacaaggccgccaacaatggtggtgatGAGCGGAGGCgg                                                                  <  2:729253‑M2/12‑1 (MQ=255)           tagtttattggccagcatgccgaggataaaggcgagcacaaggccgccaacaatggtggtgatGAGCGGAGGCGGa                                                                 >  2:940058‑M2/64‑76 (MQ=255)               ttattggccagcatgccgaggataaaggcgagcacaaggccgccaacaatggtggtgatGAGCGGAGGCGGAATTc                                                             <  1:191341‑M2/17‑1 (MQ=255)               ttattggccagcatgccgaggataaaggcgagcacaaggccgccaacaatggtggtgatGAGCGGAGGCGGAATTc                                                             >  2:894125‑M2/60‑76 (MQ=255)                  ttggccagcatgccgaggataaaggcgagcacaaggccgccaacaatggtggtgatGAGCGGAGGCGGAATTCGTg                                                          >  2:769197‑M2/57‑76 (MQ=255)                   tggccagcatgccgaggataaaggcgagcacaaggccgccaacaatggtggtgatGAGCGGAGGCGGAATTCGTGa                                                         >  1:426190‑M2/56‑76 (MQ=255)                   tggccagcatgccgaggataaaggcgagcacaaggccgccaacaatggtggtgatGAGCGGAGGCGGAATTCGTGa                                                         >  1:51829‑M2/56‑76 (MQ=255)                   tggccagcatgccgaggataaaggcgagcacaaggccgccaacaatggtggtaatGAGCGGATGCGGAATTCGTGa                                                         >  1:996388‑M2/56‑76 (MQ=255)                    ggccagcatgccgaggataaaggcgagcacaaggccgccaacaatggtggtgatGAGCGGAGGCGGAATTCGTGAt                                                        <  1:841773‑M2/22‑1 (MQ=255)                           atgccgaggataaaggcgagcacaaggccgccaacaatggtggtgatGAGCGGAGGCGGAATTCGTGATTCTATCg                                                 <  1:482486‑M2/29‑1 (MQ=255)                            tgccgaggataaaggcgagcacaaggccgccaacaatggtggtgatGAGCGGAGGCGGAATTCGTGATTCTATCGa                                                <  2:91609‑M2/30‑1 (MQ=255)                              ccgaggataaaggcgagcacaaggcctccaacaatggtggtgatGAGCGGAGGCGGAATTCGTGATTCTATCGAAg                                              <  2:866777‑M2/32‑1 (MQ=255)                              ccgaggataaaggcgagcacaaggccgccaacaatggtggtgatGAGCGGAGGCGGAATTCGTGATTCTATCGAAg                                              >  1:690184‑M2/45‑76 (MQ=255)                              ccgaggataaaggcgagcacaaggccgccaacaatggtggtgatGAGCGGAGGCGGAATTCGTGATTCTATCGAAg                                              <  2:392442‑M2/32‑1 (MQ=255)                                       aaggcgagcacaaggccgccaacaatggtggtgatGAGCGGAGGCGGAATTCGTGATTCTATCGAAGCAGGCGata                                     >  2:120093‑M2/36‑76 (MQ=255)                                       aaggcgagcacaaggccaccaacaatggtggtgatGAGCGGAGGCGGAATTCGTGANTCTATCGAAGCAGGCGata                                     <  1:1089072‑M2/41‑1 (MQ=255)                                        aggcgagcacaaggccgccaacaatggtggtgatGAGCGGAGGCGGAATTCGTGATTCTATCGAAGCAGGCGatat                                    <  1:646048‑M2/42‑1 (MQ=255)                                             agcacaaggccgccaacaatggtggtgatGAGCGGAGGCGGAATTCGTGATTCTATCGAAGCAGGCGATATCAGCt                               <  2:285424‑M2/47‑1 (MQ=255)                                               canaagtccgcgaacaatggtggtggtGAGCGGAGGCGGAATTCGTGATTCTATCGAAGCAGGCGATATCAGCtat                             <  2:1102835‑M2/49‑1 (MQ=255)                                                acaaggccgccaacaatggtggtgatGAGCGGAGGCGGAATTCGTGATTCTATCGAAGCAGGCGATATCAGCtata                            <  1:247479‑M2/50‑1 (MQ=255)                                                      ccgccaacaatggtggtgatGAGCGGAGGCGGAATTCGTGATTCTATCGAAGCAGGCGATATCAGCTATAAAAACg                      <  2:896895‑M2/56‑1 (MQ=255)                                                           aacaatggtggtgatGAGCGGAGGCGGAATTCGTGATTCTATCGAAGCAGGCGATATCAGCTATAAAAACGTGCTg                 <  2:786409‑M2/61‑1 (MQ=255)                                                            acaatggtggtgatGAGCGGAGGCGGAACTCGTGATCATATCGAAGCNGGAGATATCAGCTATAAAAACGTGCTGa                >  1:96860‑M2/15‑76 (MQ=255)                                                                    ggtgatGAGCGGAGGCGGAATTCGTGATTCTATCGAAGCAGGCGATATCAGCTATAAAAACGTGCTGAAAGTGCAg        <  2:588379‑M2/70‑1 (MQ=255)                                                                    ggtgatGAGCGGAGGCGGAATTCGTGATTCTATCGAAGCAGGCGATATCAGCTATAAAAACGTGCTGAAAGTGCAg        <  1:376645‑M2/70‑1 (MQ=255)                                                                         tGAGCGGANGCGGAGTTCGNGATTCTATCGAAGCAGGCGATANCAGCTATAAAAACGTGCTGAAAGTGCAGCCAtt   <  2:971857‑M2/75‑1 (MQ=255)                                                                         tGAGCGGAGGCGGAATTCGTGATTCTATCGAAGCAGGCGATATCAGCTATAAAAACGTGCTGAAAGTGCAGCCAtt   <  1:529705‑M2/75‑1 (MQ=255)                                                                          gaGCGGAGGCGGAATTCGTGATTCTATCGAAGCAGGCGATATCAGCTATAAAAACGTGCTGAAAGTGCAGCCATTc  <  2:552211/76‑1 (MQ=255)                                                                          |                                                                            TTGTACAGACCAATATGGGGCGGTTGATTCTGGCAGCCCAAATGGATCGCACTGGTGCCGACTTTGCGGTGATGAGCGGAGGCGGAATTCGTGATTCTATCGAAGCAGGCGATATCAGCTATAAAAACGTGCTGAAAGTGCAGCCATTC  >  REL606/478119‑478267

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