Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,366,751	Δ295 bp	IS5‑mediated	yhcE–[yhcF]	yhcE, [yhcF]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3365727–3366749	3367045	297–1319	18 [15]	[0] 95	[insH‑10]–[yhcF]	[insH‑10],yhcE,[yhcF]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3366750	NA (NA)	95 (1.440)	71/426	0.0	100%	noncoding (1/1195 nt)	IS5	repeat region
	?	NC_000913	3367046 =	0 (0.000)					coding (121/717 nt)	yhcF	DUF1120 domain‑containing protein YhcF

AATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTT  >  NC_000913/3365488‑3365789                                                                                                                                                                                                                                               |                                                                aaTCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGt                                                                 <  1:197234/239‑1 (MQ=255)       gACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGGCTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCAttt                                                          >  2:222096/1‑240 (MQ=255)       gACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGt                                                                 <  2:587933/233‑1 (MQ=255)       gACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGt                                                                 >  1:587941/1‑233 (MQ=255)              cGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCg                                                    <  1:276234/239‑1 (MQ=255)                       gTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAGCAg                                          >  1:550899/1‑240 (MQ=255)                        tGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAgg                                         >  1:682009/1‑240 (MQ=255)                         ggCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGt                                        <  2:589153/240‑1 (MQ=255)                             aaaaTGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGc                                    <  1:445200/240‑1 (MQ=255)                             aaaaTGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGc                                    <  2:523215/240‑1 (MQ=255)                                  gACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGc                                <  1:536959/239‑1 (MQ=255)                                                     gAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGtt             >  2:581716/1‑239 (MQ=35)                                                          gCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGaa                                              >  2:27155/1‑201 (MQ=33)                                                          gCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTAt       >  1:159181/1‑240 (MQ=33)                                                          gCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTAt       >  2:31482/1‑240 (MQ=33)                                                          gCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTAt       >  1:210495/1‑240 (MQ=33)                                                           cATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATc      >  1:600304/1‑240 (MQ=32)                                                               aaTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGttt  >  1:754744/1‑240 (MQ=11)                                                                                                                                                                                                                                               |                                                                AATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTT  >  NC_000913/3365488‑3365789

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑