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	68 [0]	[0] 64	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)	61 (1.570)	26/118	0.1	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

GGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAG  >  NC_000913/1978435‑1978568                                                                    |                                                                   ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                    >  1:885291‑M1/1‑68 (MQ=255) ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                    >  1:820125‑M1/1‑68 (MQ=255) ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                    >  1:1715346‑M1/1‑68 (MQ=255) ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                    >  1:39027‑M1/1‑68 (MQ=255)   ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGta                                                                  >  1:2589402‑M1/1‑66 (MQ=255)   ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGta                                                                  >  1:142008‑M1/1‑66 (MQ=255)      gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                               >  1:1038441‑M1/1‑63 (MQ=255)      gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                               >  1:267443‑M1/1‑63 (MQ=255)      gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                               >  1:3174235‑M1/1‑63 (MQ=255)           gATTTTCAATAGTGCGTGNTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcat                                                          >  1:2566819‑M1/1‑58 (MQ=255)           gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcat                                                          >  1:2798814‑M1/1‑58 (MQ=255)              tttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatact                                                       <  1:1790029‑M1/68‑14 (MQ=255)                  aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACAGAAATGtataagtcatacttttg                                                   <  1:2254378‑M1/68‑18 (MQ=255)                   ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgt                                                  <  1:3151520‑M1/68‑19 (MQ=255)                   ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgt                                                  <  1:2936855‑M1/68‑19 (MQ=255)                     aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttt                                                >  1:398133‑M1/1‑48 (MQ=255)                     aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttt                                                >  1:2276464‑M1/1‑48 (MQ=255)                     aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttt                                                >  1:3464103‑M1/1‑48 (MQ=255)                       tGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                              >  1:2259697‑M1/1‑46 (MQ=255)                          gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                           <  1:2932076‑M1/68‑26 (MQ=255)                          gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                           <  1:2538909‑M1/68‑26 (MQ=255)                          gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                           <  1:251093‑M1/68‑26 (MQ=255)                          gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                           >  1:1266994‑M1/1‑43 (MQ=255)                            gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgt                                         <  1:1576213‑M1/68‑28 (MQ=255)                              tGCAGATCACACAAAAGACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                       <  1:296241‑M1/68‑30 (MQ=255)                              tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                       <  1:2027909‑M1/68‑30 (MQ=255)                              tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                       <  1:1502702‑M1/68‑30 (MQ=255)                              tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                       <  1:1546575‑M1/68‑30 (MQ=255)                              tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                       <  1:3181117‑M1/68‑30 (MQ=255)                              tGCAGATCACACAAAACACGCAATTACTTAACATAAATGtataggtcatacttttgttttgggtgtat                                       <  1:879856‑M1/68‑30 (MQ=255)                               gCAGATCACACAAANCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                      >  1:1168417‑M1/1‑38 (MQ=255)                               gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                      >  1:177562‑M1/1‑38 (MQ=255)                               gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                      >  1:2360427‑M1/1‑38 (MQ=255)                               gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                      >  1:2273309‑M1/1‑38 (MQ=255)                                 aGATCACACAAAACACTCNATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                    >  1:2581671‑M1/1‑36 (MQ=255)                                    tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                 >  1:2523259‑M1/1‑33 (MQ=255)                                    tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                 >  1:3042987‑M1/1‑33 (MQ=255)                                    tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                 >  1:3390906‑M1/1‑33 (MQ=255)                                    tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                 >  1:158476‑M1/1‑33 (MQ=255)                                    tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                 >  1:1177573‑M1/1‑33 (MQ=255)                                      acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatct                               >  1:2302278‑M1/1‑31 (MQ=255)                                        acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                             >  1:1497849‑M1/1‑29 (MQ=255)                                        acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                             >  1:898534‑M1/1‑29 (MQ=255)                                        acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                             >  1:2216190‑M1/1‑29 (MQ=255)                                        acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                             <  1:1811736‑M1/68‑40 (MQ=255)                                         caAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                            >  1:2921542‑M1/1‑28 (MQ=255)                                           aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaa                          >  1:2748582‑M1/1‑26 (MQ=255)                                              cacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                       >  1:2679104‑M1/1‑23 (MQ=255)                                              cacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                       >  1:1615331‑M1/1‑23 (MQ=255)                                              cacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                       >  1:2565209‑M1/1‑23 (MQ=255)                                                    aTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                 >  1:919051‑M1/1‑17 (MQ=255)                                                    aTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                 >  1:1932614‑M1/1‑17 (MQ=255)                                                      tACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcg               >  1:2745550‑M1/1‑15 (MQ=255)                                                       aCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgc              >  1:2855673‑M1/1‑14 (MQ=255)                                                       aCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgc              >  1:126764‑M1/1‑14 (MQ=255)                                                       aCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgc              >  1:2934812‑M1/1‑14 (MQ=255)                                                       aCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgc              >  1:3109937‑M1/1‑14 (MQ=255)                                                         ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta            >  1:217253‑M1/1‑12 (MQ=255)                                                         ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta            >  1:3201541‑M1/1‑12 (MQ=255)                                                         ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta            >  1:1840654‑M1/1‑12 (MQ=255)                                                         ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta            >  1:133896‑M1/1‑12 (MQ=255)                                                         ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta            >  1:2527385‑M1/1‑12 (MQ=255)                                                         ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgataaaaagtttttcgcta            >  1:2999608‑M1/1‑12 (MQ=255)                                                          tAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctac           <  1:2901685‑M1/68‑58 (MQ=255)                                                          tAACATAAATGtataagtcatactttngttttgggtgtatttcaatctgttaaaaagtttttcgctac           <  1:2719595‑M1/68‑58 (MQ=255)                                                              aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta       >  1:2299931‑M1/1‑7 (MQ=255)                                                               tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctag      >  1:42458‑M1/1‑6 (MQ=255)                                                                   tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag  >  1:2777260‑M1/1‑2 (MQ=255)                                                                    |                                                                   GGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAG  >  NC_000913/1978435‑1978568

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