Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	REL606	2,009,951	(GCA)3→2	intergenic (‑81/+23)	yeeX ← / ← yeeA	hypothetical protein/conserved inner membrane protein

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	REL606	2009944	2009946	3	11 [0]	[0] 12	yeeX/yeeA	hypothetical protein/conserved inner membrane protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	2009947 =	0 (0.000)	10 (0.510)	10/136	0.9	100%	intergenic (‑77/+29)	yeeX/yeeA	hypothetical protein/conserved inner membrane protein
	?	REL606	= 2009943	0 (0.000)					intergenic (‑73/+33)	yeeX/yeeA	hypothetical protein/conserved inner membrane protein

AGCTGCCGCAGATGACAATGCTTTTATCCCTTTTCACTATCATACCCTTTATCCTTGCTGAATCGAAGCAGCAGCAAGATGATTCTGAAGTTCAGGAATTATTTGCGCAAGGCCCGGCAAATCAGATTCGATAGCAACTCAAGCT  >  REL606/2009878‑2010022                                                                      |                                                                            catagctgccgcagatgacaatgcttttatcccttttcactatcataccctttatccttgctgaatcgaAgcagca                                                                       <  1:331373‑M1/7‑1 (MQ=255)   tagctgccgcagatgacaatgcttttatcccttttcactatcataccctttatccttgctgaatcgaAGCAGCAAg                                                                     >  1:193092‑M1/68‑76 (MQ=255)                     atgcttttatcccttttcactatcataccctttatccttgctgaatcgaAGCAGCAAGATGATTCTGAAGTTCAgg                                                   >  1:592929‑M1/50‑76 (MQ=255)                                   tttcactatcataccctttatccttgctgaatcgaAGCAGCAAGATGATTCTGAAGTTCAGGAATTATTTGCGCaa                                     >  2:29625‑M1/36‑76 (MQ=255)                                            cataccctttatccttgctgaatcgaAGCAGCAAGATGATTCTGAAGTTCAGGAATTATTTGCGCAAGGCCCGGCa                            >  1:664224‑M1/27‑76 (MQ=255)                                              taccctttatccttgctgaatcgaAGCAGCAAGATGATTCTGAAGTTCAGGAATTATTTGCGCAAGGCCCGGCaaa                          <  2:566191‑M1/52‑1 (MQ=255)                                                     tatccttgctgaatcgaAGCAGCAAGATGATTCTGAAGTTCAGGAATTATTTGCGCAAGGCCCGGCAAATCAGAtt                   <  2:239502‑M1/59‑1 (MQ=255)                                                       tccttgctgaatcgaAGCAGCAAGATGATTCTGAAGTTCAGGAATTATTTGCGCAAGGCCCGGCAAATCAGATTCg                 >  1:4803‑M1/16‑76 (MQ=255)                                                          ttgctgaatcgaAGCAGCAAGATGATTCTGAAGTTCAGGAATTATTTGCGCAAGGCCCGGCAAATCAGATTCGATa              <  1:154598‑M1/64‑1 (MQ=255)                                                                  tcgaAGCAGCAAGATGATTCTGAAGTTCAGGAATTATTTGCGCAAGGCCCGGCAAATCAGATTCGATAGCAACTCa      <  1:642182‑M1/72‑1 (MQ=255)                                                                     aAGCAGCAAGATGATTCTGAAGTTCAGGAATTATTTGCGCAAGGCCCGGCAAATCAGATTCGATAGCAACTCAAGc   <  1:666212‑M1/75‑1 (MQ=255)                                                                      AGCAGCAAGATGATTCTGAAGTTCAGGAATTATTTGCGCAAGGCCCGGCAAATCAGATTCGATAGCAACTCAAGCt  >  2:136149‑M1/1‑76 (MQ=255)                                                                      |                                                                            AGCTGCCGCAGATGACAATGCTTTTATCCCTTTTCACTATCATACCCTTTATCCTTGCTGAATCGAAGCAGCAGCAAGATGATTCTGAAGTTCAGGAATTATTTGCGCAAGGCCCGGCAAATCAGATTCGATAGCAACTCAAGCT  >  REL606/2009878‑2010022

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