Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,978,503	Δ776 bp		insB‑5–insA‑5	insB‑5, insA‑5

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1978503–1979270	1979278	9–776	69 [0]	[0] 68	insB‑5–insA‑5	insB‑5,insA‑5

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1978502	0 (0.000)	57 (1.470)	30/112	0.0	100%	intergenic (‑305/+16)	flhD/insB‑5	DNA‑binding transcriptional dual regulator FlhD/IS1 protein InsB
	?	NC_000913	1979279 =	0 (0.000)					intergenic (‑64/‑474)	insA‑5/uspC	IS1 protein InsA/universal stress protein C

GGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATG  >  NC_000913/1978438‑1978565                                                                 |                                                                gggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                 <  1:836792‑M1/65‑1 (MQ=255)   gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGta                                                               >  1:3188787‑M1/1‑63 (MQ=255)   gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGta                                                               >  1:2162448‑M1/1‑63 (MQ=255)   gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGta                                                               >  1:2149713‑M1/1‑63 (MQ=255)     tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtata                                                             >  1:3236551‑M1/1‑61 (MQ=255)       cGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataag                                                           <  1:1248302‑M1/65‑7 (MQ=255)       cGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataag                                                           <  1:544259‑M1/65‑7 (MQ=255)        gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagt                                                          >  1:999211‑M1/1‑58 (MQ=255)        gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagt                                                          >  1:3291521‑M1/1‑58 (MQ=255)        gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagt                                                          >  1:2978940‑M1/1‑58 (MQ=255)                ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttt                                                  >  1:3446256‑M1/1‑50 (MQ=255)                ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttt                                                  >  1:3110071‑M1/1‑50 (MQ=255)                  aatGCGTGATGCAGATCACACAAAGCACTCAATTACTTAACATAAATGtataagtcatacttttg                                                >  1:2533456‑M1/1‑48 (MQ=255)                  aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttg                                                >  1:2738529‑M1/1‑48 (MQ=255)                     gCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttt                                             >  1:2407961‑M1/1‑45 (MQ=255)                       gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                           >  1:386202‑M1/1‑43 (MQ=255)                        tGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                          <  1:1890970‑M1/65‑24 (MQ=255)                         gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                         >  1:1995668‑M1/1‑41 (MQ=255)                          aTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                        <  1:1186801‑M1/65‑26 (MQ=255)                          aTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                        <  1:57255‑M1/65‑26 (MQ=255)                          aTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                        <  1:2746604‑M1/65‑26 (MQ=255)                           tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtg                                       >  1:3627996‑M1/1‑39 (MQ=255)                            gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgt                                      >  1:732035‑M1/1‑38 (MQ=255)                            gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgt                                      <  1:748563‑M1/65‑28 (MQ=255)                             cAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgta                                     <  1:3124748‑M1/65‑29 (MQ=255)                              aGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatccttttgttttgggtgtat                                    <  1:3329509‑M1/65‑30 (MQ=255)                              aGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                    <  1:3060024‑M1/65‑30 (MQ=255)                              aGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                    <  1:379447‑M1/65‑30 (MQ=255)                              aGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                    >  1:1258001‑M1/1‑36 (MQ=255)                                 tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                 >  1:311333‑M1/1‑33 (MQ=255)                                 tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                 >  1:2926090‑M1/1‑33 (MQ=255)                                 tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                 >  1:3581688‑M1/1‑33 (MQ=255)                                 tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                 >  1:2764228‑M1/1‑33 (MQ=255)                                 tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                 >  1:1477007‑M1/1‑33 (MQ=255)                                 tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                 >  1:2085845‑M1/1‑33 (MQ=255)                                 tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                 >  1:1610230‑M1/1‑33 (MQ=255)                                 tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                 >  1:1848823‑M1/1‑33 (MQ=255)                                      caAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatct                            >  1:1645500‑M1/1‑28 (MQ=255)                                      caAAACACTCAATTACTTAACATAAATGtataagtcatactcttgttttgggtgtatttcaatcc                            >  1:2409676‑M1/1‑28 (MQ=255)                                        aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                          >  1:2461372‑M1/1‑26 (MQ=255)                                        aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                          <  1:1540441‑M1/65‑40 (MQ=255)                                          acacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtta                        >  1:3341656‑M1/1‑24 (MQ=255)                                          acacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtta                        >  1:1891325‑M1/1‑24 (MQ=255)                                             cTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaa                     >  1:774396‑M1/1‑21 (MQ=255)                                               cAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaag                   >  1:283649‑M1/1‑19 (MQ=255)                                                aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagt                  <  1:655692‑M1/65‑48 (MQ=255)                                                 aTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtt                 >  1:320542‑M1/1‑17 (MQ=255)                                                 aTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtt                 >  1:197575‑M1/1‑17 (MQ=255)                                                 aTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtt                 >  1:2325339‑M1/1‑17 (MQ=255)                                                    aCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt              >  1:615796‑M1/1‑14 (MQ=255)                                                    aCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt              >  1:2900946‑M1/1‑14 (MQ=255)                                                    aCTTAACATAAATGtataagtcatacttttgttttgggcgtatttcaatctgttaaaaagttttt              >  1:1204374‑M1/1‑14 (MQ=255)                                                      ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcg            >  1:3132042‑M1/1‑12 (MQ=255)                                                      ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcg            >  1:1901849‑M1/1‑12 (MQ=255)                                                      ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcg            >  1:554679‑M1/1‑12 (MQ=255)                                                        gacATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgct          <  1:3520905‑M1/64‑56 (MQ=255)                                                        aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgct          <  1:1591232‑M1/65‑56 (MQ=255)                                                        aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgct          >  1:2662928‑M1/1‑10 (MQ=255)                                                          cATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctac        <  1:44895‑M1/65‑58 (MQ=255)                                                           aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacg       >  1:3472685‑M1/1‑7 (MQ=255)                                                           aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacg       >  1:3230667‑M1/1‑7 (MQ=255)                                                            tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgc      <  1:1638268‑M1/65‑60 (MQ=255)                                                             aaaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgct     <  1:1959152‑M1/65‑61 (MQ=255)                                                               aTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctag   <  1:2154984‑M1/65‑63 (MQ=255)                                                                tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagc  >  1:2511922‑M1/1‑2 (MQ=255)                                                                tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagc  >  1:2080159‑M1/1‑2 (MQ=255)                                                                tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagc  >  1:515143‑M1/1‑2 (MQ=255)                                                                tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagc  >  1:2836067‑M1/1‑2 (MQ=255)                                                                tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagc  >  1:92368‑M1/1‑2 (MQ=255)                                                                 |                                                                GGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATG  >  NC_000913/1978438‑1978565

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