Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	687,860	Δ1,185 bp	IS5‑mediated	insH‑3	insH‑3

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	687860	689044	1185	25 [0]	[0] 27	insH‑3	insH‑3

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 687859	0 (0.000)	24 (0.890)	22/94	0.1	100%	noncoding (1187/1195 nt)	IS5	repeat region
	?	NC_000913	689045 =	0 (0.000)					noncoding (1/1195 nt)	IS5	repeat region

TATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCAT  >  NC_000913/687811‑687902                                                 |                                            tATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGC                                             >  1:724012‑M1/1‑49 (MQ=255)  aTCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCc                                            >  1:1612562‑M1/1‑48 (MQ=255)    cGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCct                                           >  1:1600878‑M1/1‑46 (MQ=255)     gTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCctta                                         <  1:1095482‑M1/49‑5 (MQ=255)     gTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCctt                                          <  1:1276929‑M1/48‑4 (MQ=255)     gTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCctt                                          <  1:573909‑M1/48‑4 (MQ=255)     gTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCct                                           >  1:489746‑M1/1‑45 (MQ=255)      tCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCcttag                                        <  1:1177861‑M1/49‑6 (MQ=255)       cTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCcttagg                                       >  1:585098‑M1/1‑43 (MQ=255)       cTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCcttagg                                       >  1:1992725‑M1/1‑43 (MQ=255)        tGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCcttagg                                       >  1:213373‑M1/1‑42 (MQ=255)                tATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCcttaggtaacattt                               <  1:687488‑M1/48‑15 (MQ=255)                     tgCAGTGTTGTGCCTGTTAGGGAAGGTGCcttaggtaacatttagtttg                         <  1:1026648‑M1/49‑21 (MQ=255)                       cAGTGTTGTGCCTGTTAGGGAAGGTGCcttaggtaacatttagtttggc                       >  1:624456‑M1/1‑27 (MQ=255)                        aGTGTTGTGCCTGTTAGGGAAGGTGCcttaggtaacatttagtttggc                       >  1:418699‑M1/1‑26 (MQ=255)                          tgttgtGCCTGTTAGGGAAGGTGCcttaggtaacatttagtttggcta                     >  1:1132030‑M1/1‑24 (MQ=255)                               tgCCTGTTAGGGAAGGTGCcttaggtaacatttagtttggctaaatgta               <  1:866641‑M1/49‑31 (MQ=255)                                gCCTGTTAGGGAAGGTGCcttaggtaacatttagtttggctaaatgta               >  1:1871609‑M1/1‑18 (MQ=255)                                  cTGTTAGGGAAGGTGCcttaggtaacatttagtttggctaaatgtaaag            <  1:1679333‑M1/49‑34 (MQ=255)                                   tGTTAGGGAAGGTGCcttaggtaacatttagtttggctaaatgta               >  1:1071217‑M1/1‑15 (MQ=255)                                    gTTAGGGAAGGTGCcttaggtaacatttagtttggctaaatgtaaaga           >  1:836117‑M1/1‑14 (MQ=255)                                     ttAGGGAAGGTGCcttaggtaacatttagtttggctaaatgtaaaga           <  1:1907091‑M1/47‑35 (MQ=255)                                          gAAGGTGCcttaggtaacatttagtttggctaaatgtaaagatattgct    <  1:87629‑M1/49‑42 (MQ=255)                                           aaGGTGCcttaggtaacatttagtttggctaaatgtaaagatattgctg   <  1:1097195‑M1/49‑43 (MQ=255)                                            aGGTGCcttaggtaacatttagtttggctaaatgtaaagatattgctgt  >  1:1535515‑M1/1‑6 (MQ=255)                                                 |                                            TATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCAT  >  NC_000913/687811‑687902

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