Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
MC JC	NC_000913	3,366,751	Δ1,271 bp	100%	IS5‑mediated	[yhcE]–[yhcG]	[yhcE], yhcF, [yhcG]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3365662–3366754	3368021	1268–2360	29 [26]	[0] 118	insH1–[yhcG]	insH1,yhcE,yhcF,[yhcG]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3366750	NA (NA)	114 (1.320)	85/264	NT	100%	noncoding (1/1195 nt)	IS5	repeat region
	?	NC_000913	3368022 =	0 (0.000)					coding (196/1128 nt)	yhcG	DUF1016 family protein in the PD‑(D/E)XK nuclease superfamily

GCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCA  >  NC_000913/3365513‑3365697                                                                                                                                                     |                                     gCAAAAATGCCAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTAtt                                      <  8:198448/149‑1 (MQ=255) gCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTAtt                                      <  4:87851/149‑1 (MQ=255) gCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTAtt                                      <  3:329733/149‑1 (MQ=255)    aaaaTGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAg                                   <  5:378484/149‑1 (MQ=255)     aaaTGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTGCCGTAAATGCCTTGAATCAGCCTATTTAGa                                  <  2:161318/149‑1 (MQ=255)      aaTGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGAc                                 >  8:208038/1‑149 (MQ=255)       aTGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGATTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGAcc                                >  6:84816/1‑149 (MQ=255)         gACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGt                              >  7:36169/1‑149 (MQ=255)           cAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGATTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGttt                            <  5:84816/149‑1 (MQ=255)            aaTCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttc                           >  2:258539/1‑149 (MQ=255)             aTCCGGGCAGCACAGGAACAGCAGCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCATGAATCAGCCTATTTAGACAGTtact                          >  2:37457/1‑149 (MQ=255)                cGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttcttc                        <  3:108658/148‑1 (MQ=255)                  ggCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGcc                     <  1:368758/149‑1 (MQ=255)                   gcagcaCAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCa                    >  8:252363/1‑149 (MQ=255)                    cagcaCAGGAACAGCATCAATTATGGAAGGTGCGAATAAGCGGGGTAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCAt                   <  5:231003/149‑1 (MQ=255)                      gcaCAGGAACAGCATCAATTAGGGTAGGTGCGAATAAGCGGTGATATTCTTCTCGGCTGAATCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCAttt                 <  8:176783/149‑1 (MQ=255)                      gcaCAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCAttt                 <  8:30125/149‑1 (MQ=255)                      gcaCAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCAttt                 <  2:291584/149‑1 (MQ=255)                        acaGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCAtta                 >  8:141674/1‑146 (MQ=255)                         caGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAg              <  8:124082/149‑1 (MQ=255)                          ccgatccgCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAg              <  5:74877/141‑1 (MQ=25)                          aGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAgg             >  6:196468/1‑149 (MQ=40)                          aGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAgg             >  7:9059/1‑149 (MQ=40)                           ggAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTACTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGc            <  3:54953/149‑1 (MQ=255)                                  cATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGCCTCAGTCATTTCATTTCTTCATGTTTGAGCCGAGTTTTTCTCCCGTAAATGGCTTGAATAAGCCGATTTAGACACTTGCTTCGACATTTAAGGCgatcgg      >  6:74877/1‑143 (MQ=11)                                  cATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATcc     >  6:340476/1‑149 (MQ=34)                                   aTCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATccc    >  4:334820/1‑149 (MQ=33)                                    tCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTAATTTCTTCTTGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATcca    >  3:39616/1‑147 (MQ=17)                                      aaTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCa  <  3:334820/148‑1 (MQ=31)                                                                                                                                                     |                                     GCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCA  >  NC_000913/3365513‑3365697

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