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	71 [0]	[0] 65	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)	56 (1.270)	30/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:295518‑M1/1‑68 (MQ=255) ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                    >  1:1976008‑M1/1‑68 (MQ=255) ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                    >  1:3821769‑M1/1‑68 (MQ=255) ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                    >  1:3257759‑M1/1‑68 (MQ=255) ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                    >  1:384430‑M1/1‑68 (MQ=255) ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                    >  1:401334‑M1/1‑68 (MQ=255) ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                    >  1:1508504‑M1/1‑68 (MQ=255)  gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGt                                                                   <  1:291053‑M1/68‑2 (MQ=255)  gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGt                                                                   <  1:3501877‑M1/68‑2 (MQ=255)  gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGt                                                                   <  1:4068829‑M1/68‑2 (MQ=255)   ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGta                                                                  >  1:3628941‑M1/1‑66 (MQ=255)      gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                               >  1:2686534‑M1/1‑63 (MQ=255)      gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                               >  1:3922231‑M1/1‑63 (MQ=255)      gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                               >  1:3349144‑M1/1‑63 (MQ=255)      gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                               >  1:3434601‑M1/1‑63 (MQ=255)       cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataag                                                              <  1:3039973‑M1/68‑7 (MQ=255)       cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataag                                                              <  1:1220887‑M1/68‑7 (MQ=255)        tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagt                                                             <  1:1021292‑M1/68‑8 (MQ=255)        tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagt                                                             <  1:3785465‑M1/68‑8 (MQ=255)           gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcat                                                          >  1:3377524‑M1/1‑58 (MQ=255)           gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcat                                                          >  1:1382517‑M1/1‑58 (MQ=255)              tttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatact                                                       <  1:1580218‑M1/68‑14 (MQ=255)              tttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatact                                                       <  1:505932‑M1/68‑14 (MQ=255)                   ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgt                                                  <  1:3792705‑M1/68‑19 (MQ=255)                   ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgt                                                  >  1:3258802‑M1/1‑50 (MQ=255)                   ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttcgt                                                  >  1:89473‑M1/1‑50 (MQ=255)                    taatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtt                                                 <  1:195396‑M1/68‑20 (MQ=255)                        gCGTAATGCAGATCACACAAAACACTCAATTACTTAACAGAAATGtataagtcatacttttgttttgg                                             <  1:2411243‑M1/68‑24 (MQ=255)                         cGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                            <  1:1762946‑M1/68‑25 (MQ=255)                          gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                           >  1:3025318‑M1/1‑43 (MQ=255)                          gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                           >  1:662476‑M1/1‑43 (MQ=255)                          gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                           <  1:2819573‑M1/68‑26 (MQ=255)                          gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                           <  1:1684451‑M1/68‑26 (MQ=255)                            gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgt                                         <  1:1970043‑M1/68‑28 (MQ=255)                              tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                       <  1:340972‑M1/68‑30 (MQ=255)                              tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                       <  1:1819304‑M1/68‑30 (MQ=255)                              tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                       <  1:538634‑M1/68‑30 (MQ=255)                               gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                      >  1:3667219‑M1/1‑38 (MQ=255)                               gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                      >  1:3647983‑M1/1‑38 (MQ=255)                               gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                      >  1:2134031‑M1/1‑38 (MQ=255)                               gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                      >  1:1260110‑M1/1‑38 (MQ=255)                                 aGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                    >  1:2776919‑M1/1‑36 (MQ=255)                                 aGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                    >  1:297762‑M1/1‑36 (MQ=255)                                   aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaa                                  >  1:1959463‑M1/1‑34 (MQ=255)                                    tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcctg                                 >  1:2082084‑M1/1‑33 (MQ=255)                                    tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                 >  1:2808441‑M1/1‑33 (MQ=255)                                    tCACACAAAACACTCAATTACTTAACATAAATGtataagtcacacttttgttttgggtgtatttcaat                                 >  1:2742599‑M1/1‑33 (MQ=255)                                        acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                             <  1:1452094‑M1/68‑40 (MQ=255)                                        acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                             <  1:1168441‑M1/68‑40 (MQ=255)                                        acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggctgtatttcaatctgt                             >  1:43615‑M1/1‑29 (MQ=255)                                         caAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                            >  1:2027659‑M1/1‑28 (MQ=255)                                             acacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaa                        >  1:293794‑M1/1‑24 (MQ=255)                                              cacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                       >  1:448028‑M1/1‑23 (MQ=255)                                                cTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagt                     <  1:3254466‑M1/68‑48 (MQ=255)                                                cTCAATTACTTAACATAAATGtataagtcatacttttgttgtgggtgtatttcaatctgttaaaaagt                     >  1:3398016‑M1/1‑21 (MQ=255)                                                    aTTACTTAACATAAATNtataagtcatacttttgttttgggtgtatntcaatatgttaaaaagctttt                 >  1:3681185‑M1/1‑17 (MQ=255)                                                    aTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                 >  1:1699010‑M1/1‑17 (MQ=255)                                                       aCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgc              >  1:1829578‑M1/1‑14 (MQ=255)                                                       aCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgc              >  1:72046‑M1/1‑14 (MQ=255)                                                       aCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgc              >  1:179335‑M1/1‑14 (MQ=255)                                                         ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta            >  1:520572‑M1/1‑12 (MQ=255)                                                         ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta            >  1:3042341‑M1/1‑12 (MQ=255)                                                         ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta            >  1:1069190‑M1/1‑12 (MQ=255)                                                               tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctcg      >  1:236699‑M1/1‑6 (MQ=255)                                                                   tGtataagtcatacttttgttttgggtgtatttcaatctgttaaacagtttttcgctacgctagcaag  >  1:3952219‑M1/1‑2 (MQ=255)                                                                   tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag  >  1:3978944‑M1/1‑2 (MQ=255)                                                                   tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag  <  1:2210541‑M1/68‑67 (MQ=255)                                                                   tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag  >  1:3340116‑M1/1‑2 (MQ=255)                                                                   tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag  >  1:628492‑M1/1‑2 (MQ=255)                                                                   tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag  >  1:3024571‑M1/1‑2 (MQ=255)                                                                   tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag  >  1:824647‑M1/1‑2 (MQ=255)                                                                   tGtataagtcacacttttgttttgggtgcatttcaatctgttaaaaagtttttcgctacgctagcaag  >  1:1194639‑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 < 16 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑