Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,299,499	Δ1,199 bp		insH21	insH21

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1299499–1300693	1300697	5–1199	54 [0]	[0] 54	insH21	insH21

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1299498	0 (0.000)	51 (1.220)	38/268	0.1	100%	intergenic (+253/‑33)	ychE/insH21	putative inner membrane protein/IS5 transposase and trans‑activator
	?	NC_000913	1300698 =	0 (0.000)					intergenic (+151/‑484)	insH21/oppA	IS5 transposase and trans‑activator/oligopeptide ABC transporter periplasmic binding protein

AGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTT  >  NC_000913/1299363‑1299637                                                                                                                                        |                                                                                                                                            agttagttATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataa                                                                                                                                         <  1:542109‑M1/140‑5 (MQ=255)               aTTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaa                                                                                                                           >  1:60100‑M1/1‑122 (MQ=255)                tttCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaa                                                                                                                          >  1:270528‑M1/1‑121 (MQ=255)                 ttCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaat                                                                                                                         >  2:67944‑M1/1‑120 (MQ=255)                 ttCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaat                                                                                                                         >  1:492498‑M1/1‑120 (MQ=255)                 ttCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaat                                                                                                                                 <  2:449480‑M1/132‑13 (MQ=255)                 ttCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaat                                                                                                                                 >  1:449480‑M1/1‑120 (MQ=255)                           gCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcat                                                                                                               >  2:196541‑M1/1‑110 (MQ=255)                           gCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcat                                                                                                               >  1:459508‑M1/1‑110 (MQ=255)                                      ccATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactgg                                                                                                    <  2:459508‑M1/140‑42 (MQ=255)                                       catcatGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactgga                                                                                                   >  1:192060‑M1/1‑98 (MQ=255)                                          catGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaact                                                                                                <  2:60100‑M1/140‑46 (MQ=255)                                                 ttttCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaaatgttaaattattgtgcatttcactactggaactgtaatca                                                                                         >  2:683737‑M1/1‑88 (MQ=255)                                                      gCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaa                                                                                    >  2:465873‑M1/1‑83 (MQ=255)                                                       ccACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaag                                                                                   >  1:647730‑M1/1‑82 (MQ=255)                                                            tcatAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagataga                                                                              <  1:465873‑M1/140‑64 (MQ=255)                                                                  tCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgct                                                                        <  2:647730‑M1/140‑70 (MQ=255)                                                                     tGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttag                                                                     <  2:140524‑M1/140‑73 (MQ=255)                                                                       ggCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagcc                                                                   >  2:48220‑M1/1‑66 (MQ=255)                                                                       ggCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagcc                                                                   >  2:141680‑M1/1‑66 (MQ=255)                                                                       ggCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagcc                                                                   >  2:259907‑M1/1‑66 (MQ=255)                                                                          ttGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaat                                                                >  2:730892‑M1/1‑63 (MQ=255)                                                                            gCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatct                                                              >  2:599139‑M1/1‑61 (MQ=255)                                                                            gCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatct                                                              >  2:497018‑M1/1‑61 (MQ=255)                                                                                  gaaTAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctattt                                                        >  1:23995‑M1/1‑55 (MQ=255)                                                                                      aatTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgatt                                                    <  2:492498‑M1/140‑90 (MQ=255)                                                                                      aatTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgatt                                                    <  2:23995‑M1/140‑90 (MQ=255)                                                                                         tGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaa                                                 <  1:48220‑M1/140‑93 (MQ=255)                                                                                         tGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaa                                                 <  2:192060‑M1/140‑93 (MQ=255)                                                                                          gAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaat                                                >  2:674481‑M1/1‑47 (MQ=255)                                                                                                 tattatTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaaga                                         >  2:607656‑M1/1‑40 (MQ=255)                                                                                                  attattAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagat                                        <  2:95078‑M1/140‑102 (MQ=255)                                                                                                  attattAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagat                                        >  1:249864‑M1/1‑39 (MQ=255)                                                                                                  attattAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagat                                        >  1:560442‑M1/1‑39 (MQ=255)                                                                                                   ttattaATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatg                                       <  1:196541‑M1/140‑103 (MQ=255)                                                                                                          ttCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgtta                                <  2:560442‑M1/140‑110 (MQ=255)                                                                                                           tgcacTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgcgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaa                               >  2:424709‑M1/3‑30 (MQ=255)                                                                                                           tCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaa                               >  2:140059‑M1/1‑30 (MQ=255)                                                                                                           tCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaa                               >  2:215442‑M1/1‑30 (MQ=255)                                                                                                           tCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaa                               >  2:93658‑M1/1‑30 (MQ=255)                                                                                                               cTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggca                           >  1:62979‑M1/1‑26 (MQ=255)                                                                                                               cTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggca                           >  1:285851‑M1/1‑26 (MQ=255)                                                                                                                tGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcat                          >  1:369097‑M1/1‑25 (MQ=255)                                                                                                                     ttGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatg                     <  2:678138‑M1/140‑121 (MQ=255)                                                                                                                      tGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgc                    <  1:607656‑M1/140‑122 (MQ=255)                                                                                                                        gTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaa                  <  1:215442‑M1/140‑124 (MQ=255)                                                                                                                           gagGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagct               <  1:424709‑M1/140‑127 (MQ=255)                                                                                                                              gAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctata            <  1:461508‑M1/140‑130 (MQ=255)                                                                                                                                aaGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctataga          <  1:599139‑M1/140‑132 (MQ=255)                                                                                                                                  gTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagatt        >  2:671643‑M1/1‑7 (MQ=255)                                                                                                                                   tGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattc       >  2:280302‑M1/1‑6 (MQ=255)                                                                                                                                     cTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaacctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctg     >  2:727022‑M1/1‑4 (MQ=255)                                                                                                                                       aaataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgat   <  2:177231‑M1/140‑139 (MQ=255)                                                                                                                                        aataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgata  <  2:62979‑M1/140‑140 (MQ=255)                                                                                                                                        |                                                                                                                                            AGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTT  >  NC_000913/1299363‑1299637

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